ZZ_M2GD Inputs, Outputs, and Parameters
We are in the process of re-working the I/O format for dvmdostem. The reasons for the redesigning the format are:
- reduce pre- and post-processing effort
- use standard tools for viewing NetCDF
- use standard tools for pre- and post-processing
Note: We think the reason for the old structure was to reduce the space needed to store the data. This was done by "indexing" all the pixels in the IEM domain that have data. December 2014, we did a test that showed that the compression available with the NetCDF format would allow us to get acceptable sizes without having to change formats.
We currently have a scripted process for building dvmdostem input data
from tifs provided by SNAP (for the IEM project).
- Start by collecting all the data you need from SNAP and or your colleagues.
As of 09/18/2017 this includes:
- tif files for the 4 climate variables for the historic period
- tif files for the 4 climate variables for the projected period
- tif files for vegetation cover, soil type, drainage, and topography All the files you use as source material must have the projection information with them.
- Subset the original files (rectangular subset, pixel coordinates). Use the
create_region_input.pyscript to carry out the subsetting. - All files must have the same spatial extent and the same projection. The model is currently "dumb" in that it does not figure out or guess anything about the inputs.
- Run the
gapfill.pyscript to interpolate over single missing data points along the time axis - Run the
--conform-mask-to-inputsoption
Some details are described in the thread for Issue 135
NetCDF outputs are specified in a csv file named in config/config.js. The csv file specifies a variable name (for identification only - it does not correspond to the variable name in the code), a short description, units, and what level of detail to output on (timestep and variable part).
[Link to default file after PR merge]
Variable name, Description, Units, Yearly, Monthly, Daily, PFT, Compartment, Layer,
Example entry: VEGC,Total veg. biomass C,gC/m2,y,m,,p,c,,
This will output VegC every month, and provide both PFT and PFT compartment values.
The file is more user-friendly when viewed in a spreadsheet.
[example]
A complete list of output combinations is below
The initial list of outputs can be found at Issue #252
LAYERDEPTH, LAYERDZ, and LAYERTYPE should be automatically output if the user specifies any by-layer output. They are not currently, so ensure that they are specified on the same timestep as the desired output.
HKLAYER, LAYERDEPTH, LAYERDZ, LAYERTYPE, TCLAYER, TLAYER, and VWCLAYER must have the layer option specified or they will generate NetCDF dimension bound errors.
A single output file will be produced for each entry in the specifying file, based on variable name, timestep, and run stage.
VEGC_monthly_eq.nc
At the beginning of the model run, an output file will be constructed for each variable specified, for each run stage where NetCDF output is indicated and that has more than 0 years of run time.
Currently the model tracks the variables specified for each timestep as separate sets (i.e. monthly separate from yearly, etc). This reduces the number of map lookups every time the output function is called, but increases the number of monthly vs. yearly string comparisons.
'-' indicates that the combination is not an option 'x' indicates that the combination has been implemented in the code '?' indicates that it is undecided if the combination should be made available, or that structure in the code needs to be modified to make data available for output.
Three variables should be automatically written out if any by-layer variable is specified: Layer type Layer depth Layer thickness Currently they are written out like standard variables. Automation will need to be added in the future.
| Var | Yearly | Monthly | Daily | |
|---|---|---|---|---|
| ALD | Max | x | - | - |
| AVLN | Layer | x | x | - |
| Total | x | x | - | |
| BURNAIR2SOIN | x | x | - | |
| BURNSOIC | Layer | ? | ? | - |
| Total | x | x | - | |
| BURNSOILN | Layer | ? | ? | - |
| Total | x | x | - | |
| BURNTHICK | Total | x | x | - |
| BURNVEG2AIRC | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| BURNVEG2AIRN | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| BURNVEG2DEADC | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| BURNVEG2DEADN | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| BURNVEG2SOIABVC | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| BURNVEG2SOIABVN | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| BURNVEG2SOIBLWC | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| BURNVEG2SOIBLWN | Total | x | x | - |
| PFT | ? | ? | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| DEADC | x | x | - | |
| DEADN | x | x | - | |
| DEEPC | Total | x | x | - |
| DEEPDZ | Total | x | - | - |
| DWDC | x | x | - | |
| DWDN | x | x | - | |
| EET | PFT | x | x | x |
| Total | x | x | x | |
| GPP | Total | x | x | - |
| PFT | x | x | - | |
| Comp | x | x | - | |
| Both | x | x | - | |
| GROWEND | x | - | - | |
| GROWSTART | x | - | - | |
| HKDEEP | x | x | x | |
| HKLAYER | Layer | x | x | - |
| HKMINEA | x | x | x | |
| HKMINEB | x | x | x | |
| HKMINEC | x | x | x | |
| HKSHLW | x | x | x | |
| LAI | Total | x | x | - |
| PFT | x | x | - | |
| LAYERDEPTH | Layer | x | x | - |
| LAYERDZ | Layer | x | x | - |
| LAYERTYPE | Layer | x | x | - |
| LTRFALC | Total | x | x | - |
| PFT | x | x | - | |
| Comp | x | x | - | |
| Both | x | x | - | |
| LTRFALN | Total | x | x | - |
| PFT | x | x | - | |
| Comp | x | x | - | |
| Both | x | x | - | |
| MINEC | Total | x | x | - |
| MOSSDZ | Total | x | - | - |
| NDRAIN | Layer | ? | ? | - |
| Total | ? | ? | - | |
| NETNMIN | Layer | x | x | - |
| Total | x | x | - | |
| NIMMOB | Layer | x | x | - |
| Total | x | x | - | |
| NINPUT | Layer | ? | ? | - |
| Total | x | x | - | |
| NLOST | Layer | ? | ? | - |
| Total | x | x | - | |
| NPP | Total | x | x | - |
| PFT | x | x | - | |
| Comp | x | x | - | |
| Both | x | x | - | |
| NUPTAKEIN | Total | x | x | - |
| PFT | x | x | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| NUPTAKELAB | Total | x | x | - |
| PFT | x | x | - | |
| Comp | ? | ? | - | |
| Both | ? | ? | - | |
| NUPTAKEST | Total | x | x | - |
| PFT | x | x | - | |
| Comp | x | x | - | |
| Both | x | x | - | |
| ORGN | Layer | x | x | - |
| Total | x | x | - | |
| PERMAFROST | x | - | - | |
| PET | PFT | x | x | x |
| Total | x | x | x | |
| QDRAINAGE | Ratio | x | x | - |
| QINFILTRATION | Ratio | x | x | - |
| QRUNOFF | Ratio | x | x | - |
| RAINFALL | Total | - | - | - |
| RH | Layer | x | x | - |
| Total | x | x | - | |
| ROLB | x | x | - | |
| SHLWC | Total | x | x | - |
| SHLWDZ | Total | x | - | - |
| SNOWEND | x | - | - | |
| SNOWFALL | Total | - | - | - |
| SNOWSTART | x | - | - | |
| SNOWTHICK | Total | x | x | x |
| SOC | Layer | x | x | - |
| Total | x | x | - | |
| SWE | Total | x | x | x |
| TCDEEP | x | x | x | |
| TCLAYER | Layer | x | x | - |
| TCMINEA | x | x | x | |
| TCMINEB | x | x | x | |
| TCMINEC | x | x | x | |
| TCSHLW | x | x | x | |
| TDEEP | Total | x | x | x |
| TLAYER | Layer | x | x | - |
| TMINEA | Total | x | x | x |
| TMINEB | Total | x | x | x |
| TMINEC | Total | x | x | x |
| TSHLW | Total | x | x | x |
| VEGC | Total | x | x | - |
| PFT | x | x | - | |
| Comp | x | x | - | |
| Both | x | x | - | |
| VEGN | Total | x | x | - |
| PFT | x | x | - | |
| Comp | x | x | - | |
| Both | x | x | - | |
| VWCDEEP | Total | x | x | x |
| VWCLAYER | Layer | x | x | - |
| VWCMINEA | Total | x | x | x |
| VWCMINEB | Total | x | x | x |
| VWCMINEC | Total | x | x | x |
| VWCSHLW | Total | x | x | x |
| WATERTAB | x | x | x | |
| WDRH | x | x | - | |
| YSD | x | - | - |