How to handle time dimension in export_gridded_model() #22
Description
The ISMIP6 data request requires:
"Two-dimensional state variables (ST) are requested as a yearly snapshot corresponding to the end of the year in a simulation for state variables (such as ice thickness) and as a yearly flux average for flux variables (FL, such as surface mass balance). Fields such as surface mass balance flux should be what was applied as a boundary conditions to the ice sheet model and may be different from the input forcing file"
If ISSM models output data sub-annually, or at a different timesteps, we need to deal with the time aspect here by, potentially by aggregating or disaggregating multiple timesteps. Should this be done within the export_gridded_modelI() function? We could provide an output_frequency or something to specify how the time dimension should be handled. If only select timesteps are needed, we could pass a list of times to explicitly subset the data? Perhaps we can trim/time average the model externally to the function?