Code associated with 'Aquifer conditions, not irradiance determine the potential of photovoltaic energy for groundwater pumping across Africa'

Matlab code associated with the journal article 'Aquifer conditions, not irradiance determine the potential of photovoltaic energy for groundwater pumping across Africa' (https://www.nature.com/articles/s43247-023-00695-8)

Reference for the code: Meunier, S., Kitanidis, P.K., Cordier, A., MacDonald, A.M. Aquifer conditions, not irradiance determine the potential of photovoltaic energy for groundwater pumping across Africa. Communications Earth & Environment 4, 52 (2023).

For any questions about the code or the journal article, please contact Simon Meunier (simon.meunier@centralesupelec.fr)

Note : the matlab relation to obtain the pixel coordinates (row and column, which are integers) from the geographic coordinates (latitude and longitude in degrees) is :
row = round((latitude+37.1)/0.2 + 0.5); & column = round((longitude+20.1)/0.2 + 0.5);
For instance, the point {latitude = 17.5°, longitude = 9.00°}, which is in Niger, is at the coordinates {row (Y) = 274, column (X) = 146} on the maps. So, in the 'main.m' function, if you want to obtain results for the location {latitude = 17.5°, longitude = 9.00°}, input 'II = 274;' and 'JJ = 146;'.
Note that in the plotted maps, the point {row (Y) = 1, column (X) = 1} is at the bottom left of the map, the point {row (Y) = 386, column (X) = 1} is at the top left, the point {row (Y) = 1, column (X) = 371} is at the bottom right.