Outputs¶
There are two main output files generated by SIPNET:
1. sipnet.out: Model state variables and fluxes at each timestep.
2. events.out (name configurable): Contains a record of events processed during the simulation if event handling is enabled.
Model Outputs¶
The sipnet.out file contains a time series of state variables and fluxes from the simulation.
| # | Symbol | Output Name | Definition / Notes | Units |
|---|---|---|---|---|
| 1 | year | Year of start of timestep | (integer) | |
| 2 | day | Day-of-year of start of timestep | (integer) | |
| 3 | time | Hour-of-day (fractional) at start of timestep | hours | |
| 4 | plantWoodC | Woody plant carbon | g C m\(^{-2}\) | |
| 5 | plantLeafC | Leaf carbon | g C m\(^{-2}\) | |
| 6 | woodCreation | Carbon allocated to wood during timestep | g C m\(^{-2}\) | |
| 7 | soil | (Single) soil organic carbon (or mineral soil C pool) | g C m\(^{-2}\) | |
| 8 | coarseRootC | Coarse root carbon | g C m\(^{-2}\) | |
| 9 | fineRootC | Fine root carbon | g C m\(^{-2}\) | |
| 10 | litter | Litter carbon (0 if litter pool disabled) | g C m\(^{-2}\) | |
| 11 | soilWater | Soil water content | cm | |
| 12 | \(\overline{f}_{\text{WHC}}\) | soilWetnessFrac | Mean soil wetness as fraction of water holding capacity over timestep 1 | unitless |
| 13 | snow | Snow water equivalent | cm | |
| 14 | npp | Net primary production for timestep | g C m\(^{-2}\) | |
| 15 | nee | Net ecosystem exchange (sign convention per code: - (NPP - RH)) | g C m\(^{-2}\) | |
| 16 | cumNEE | Cumulative NEE since simulation start | g C m\(^{-2}\) | |
| 17 | \(GPP\) | gpp | Gross primary production | g C m\(^{-2}\) |
| 18 | rAboveground | Aboveground autotrophic respiration (leaves + wood) | g C m\(^{-2}\) | |
| 19 | \(R_H\) | rSoil | Heterotrophic respiration | g C m\(^{-2}\) |
| 20 | rRoot | Root (autotrophic) respiration | g C m\(^{-2}\) | |
| 21 | ra | Total autotrophic respiration (rAboveground + rRoot) | g C m\(^{-2}\) | |
| 22 | rh | Total heterotrophic respiration (litter + soil components) | g C m\(^{-2}\) | |
| 23 | \(R\) | rtot | Total ecosystem respiration (ra + rh) | g C m\(^{-2}\) |
| 24 | evapotranspiration | ET (transpiration + immedEvap + evaporation + sublimation) for timestep | cm | |
| 25 | fluxestranspiration | Transpiration component only (as a flux) | cm day\(^{-1}\) | |
| 27 | \(N_\textrm{min}\) | minN | Soil mineral nitrogen | g N m\(^{-2}\) |
| soilOrgN | Soil organic nitrogen | g N m\(^{-2}\) | ||
| litterN | Litter nitrogen | g N m\(^{-2}\) | ||
| n2o | Nitrous oxide production | g N m\(^{-2}\) | ||
| nLeaching | Nitrogen lost to leaching | g N m\(^{-2}\) | ||
| nFixation | Nitrogen demand met by fixation | g N m\(^{-2}\) | ||
| nUptake | Nitrogen demand met by uptake from soil | g N m\(^{-2}\) | ||
| ch4 | Methane production | g C m\(^{-2}\) |
An example output file can be found in tests/smoke/sipnet.out.
year day time plantWoodC plantLeafC woodCreation soil coarseRootC fineRootC litter soilWater soilWetnessFrac snow npp nee cumNEE gpp rAboveground rSoil rRoot ra rh rtot evapotranspiration fluxestranspiration minN soilOrgN litterN n2o nLeaching nFixation nUptake ch4
2016 1 0.00 1313.6134 0.0000 0.0000 2688.1376 437.8735 437.8621 279.9748 3.3419 0.2701 0.7981 -0.0258 0.0871 0.0871 0.0000 0.0226 0.0645 0.0032 0.0258 0.0612 0.0871 0.00157021 0.0000 1.0014 134.9991 13.9985 0.013744 0.0000 0.0000 0.0000 0.0000
2016 1 3.00 1313.5850 0.0000 0.0000 2688.1342 437.8651 437.8423 279.9475 3.4565 0.2833 0.6794 -0.0221 0.0872 0.1743 0.0000 0.0189 0.0683 0.0032 0.0221 0.0651 0.0872 0.00411191 0.0000 1.0029 134.9981 13.9969 0.013748 0.0000 0.0000 0.0000 0.0000
2016 1 6.00 1313.5586 0.0000 0.0000 2688.1296 437.8567 437.8225 279.9191 3.5189 0.2906 0.6125 -0.0201 0.0877 0.2619 0.0000 0.0169 0.0708 0.0033 0.0201 0.0675 0.0877 0.00442997 0.0000 1.0046 134.9970 13.9952 0.013810 0.0000 0.0000 0.0000 0.0000
Events output¶
When event handling is enabled, SIPNET will create events.out by default, or
<EVENTS_PREFIX>.out when a custom events prefix is configured.
This file is designed primarily for testing and debugging.
It contains one row for each agronomic event that is processed.
Each row lists the year, day, event type, and parameter name/value pairs.
The name/value pairs represent the state variables that are directly changed by an event, recording the change (delta) applied to each.
Information in the events output file can, in principle, be reconstructed or
inferred from the corresponding events input file and sipnet.out though this
may be confounded if simultaneous events affect the same variable.
Still, sipnet.out is the authoritative source for information about system state and evolution in time, including responses to events.
Below is an example events output file, with header enabled for clarity.
Note the delimiters: spaces separate columns, commas separate name/value pairs, and = map names with their values (deltas).
year day type param_name=delta[,param_name=delta,...]
2023 65 plant envi.plantLeafC=3.00,envi.plantWoodC=4.00,envi.fineRootC=5.00,envi.coarseRootC=6.00
2023 70 irrig envi.soilWater=5.00
2023 200 harv env.litter=5.46,envi.plantLeafC=-5.93,envi.plantWoodC=-4.75,envi.fineRootC=-3.73,envi.coarseRootC=-3.89
2024 65 plant envi.plantLeafC=3.00,envi.plantWoodC=5.00,envi.fineRootC=7.00,envi.coarseRootC=9.00
2024 70 irrig fluxes.immedEvap=2.50,envi.soilWater=2.50
2024 200 harv env.litter=4.25,envi.plantLeafC=-1.39,envi.plantWoodC=-1.63,envi.fineRootC=-2.52,envi.coarseRootC=-2.97
-
Mean soilWetnessFrac (ratio of soil water / water holding capacity) calculated as average between previous and current time step. Reported for diagnostics only. Internal moisture dependency functions use instantaneous \(W_{soil}/W_{WHC}\) (not this average), and clip that ratio to [0,1] where those dependency functions are defined. ↩