Additional "annual" (per growing season) crop outputs #1649
Labels
enhancement
new capability or improved behavior of existing capability
priority: low
Background task that doesn't need to be done right away.
Over email, @danicalombardozzi, @billsacks, and I talked about how it might be a good idea to use the "annual" (technically, per growing season) crop output infrastructure I'm developing for #1537 (PR #1616) for additional crop output variables. This would reduce the data storage/transfer requirements and postprocessing overhead needed (as well as probability of postprocessing mistakes) for people working with crop outputs.
For example, currently, the only way to get crop yield is to sum up the daily grain C addition over the course of a growing season. It would be much easier and less storage-intensive to simply output yield once per harvest, on my new
mxharvests
dimension.I thought we could use this issue to discuss the various variables that should have such output enabled, as well as implementation details such as in what data structure(s) the values should be stored.
As a starting point, I took the list of output variables requested for GGCMI phase 3. I've included the units specified by the GGCMI protocol, but these could be changed for consistency with standard CLM output units as desired.
SDATES
(taken fromcrop_inst%sdates_thisyr
).HDATES
(taken fromcrop_inst%hdates_thisyr
), which has units "day of year" likeSDATES
.1Soil moisture of upper 1m (monthly)kg m-3monthly average actual absolute available soil water content of the upper 1m soil column (above wilting point);All of these could be output on the
mxharvests
dimension. Currently, sowing date (SDATES
) is output on themxsowings
dimension, but because by definitionmxharvests
>mxsowings
, it could go onmxharvests
too.As far as implementation goes, right now I have
SDATES
andHDATES
being set up in theCropType
module:(Note the
avgflag='I'
, as these are instantaneous variables.) I could see an argument being made either way: to put all the new variables intoCropType
, or to put them into more thematically related structures (water, C, etc.).Footnotes
Instead, we could think about a "growing season length" variable for this, with
HDATES
remaining as the actual "harvest date." ↩The text was updated successfully, but these errors were encountered: