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.