Compute and save maps data separated from plots
Right now we do not save the data generated to plot the maps. A more flexible approach would be to first save the generated data and all the metadata necessary to understand the data of the maps and then load this data to plot the maps in whatever subplots structure. This would facilitate things like averaging maps from different models, computing the average of maps for certain regions of the space, and plotting maps in subplots arrangement different than the default one. Also, saving the metadata used to generate the maps it is useful to report the fixed values of the variables that do not change.