Skip to content

Commit 1b79a86

Browse files
authored
document sipnet.in, sipnet.param, add context to use of events.out (#73)
1 parent db542ed commit 1b79a86

File tree

1 file changed

+69
-6
lines changed

1 file changed

+69
-6
lines changed

docs/parameters.md

Lines changed: 69 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -327,7 +327,65 @@ Run-time parameters can change from one run to the next, or when the model is st
327327
| E\_STAR\_SNOW | 0.6 | approximate saturation vapor pressure at 0°C (kPa) |
328328
| | | |
329329

330-
## Drivers
330+
331+
## Input Files
332+
333+
### Run Settings
334+
335+
The `sipnet.in` file specifies run settings for SIPNET, including the run type, input file names, and output options. The file is self-documenting, with comments describing each option. Key features of interest include
336+
337+
- `FILENAME` Set output filenames.
338+
- `PRINT_HEADER` Print header information to output files
339+
- `DO_SINGLE_OUTPUTS` Write each variable to a separate file
340+
341+
Multi-site runs, sensitivity tests, and Monte Carlo runs are no longer supported. Typically these analyses are handled
342+
using the [PEcAn Framework](https://pecanproject.org/).
343+
344+
### Parameters and Initial Conditions
345+
346+
347+
Both initial conditions and parameters are specified in a file named `sipnet.param`.
348+
349+
_Note: A `sipnet.param-spatial` can be used for multi-site runs, but this is not currently used
350+
or documented._
351+
352+
353+
The SIPNET parameter file (`sipnet.param`) specifies model parameters and their properties for each simulation.
354+
Each line in the file corresponds to a single parameter and contains five or six space-separated values.
355+
356+
| Column | Description |
357+
| -------------- | ----------------------------------------------------------------------------------------- |
358+
| Parameter Name | Name of the parameter |
359+
| Value | Value of the parameter to use in the model |
360+
| Flag | Parameter estimation flag: <br> `0` = fixed (not estimated), <br> `1` = estimated (free) |
361+
| Min | Minimum value for optimization |
362+
| Max | Maximum value for optimization |
363+
| Stddev | (Optional) Standard deviation for prior distribution |
364+
365+
#### Example `sipnet.param` file
366+
367+
Column names are not used, but are:
368+
```
369+
param_name value flag min max stddev
370+
```
371+
372+
373+
```
374+
plantWoodInit 110 0 6600 14000 200
375+
laiInit 0 0 0 5.2 0.2
376+
litterInit 200 0 130 1200 25
377+
soilInit 7000 1 3300 19000 3000
378+
litterWFracInit 0.5 0 0 1 0.1
379+
soilWFracInit 0.6 1 0 1 0.1
380+
snowInit 1 0 0 0 1
381+
microbeInit 0.5 0 0.02 1 0.001
382+
fineRootFrac 0.2 0 0 1 0.001
383+
coarseRootFrac 0.2 0 0 1 0.001
384+
aMax 95 1 0 200 0.2
385+
aMaxFrac 0.85 0 0.66 0.86 0.005
386+
...
387+
```
388+
331389

332390
### Climate
333391

@@ -350,7 +408,7 @@ For each step of the model, for each location, the following inputs are needed.
350408
| 13 | wspd | avg. wind speed | m/s | |
351409
| 14 | soilWetness | fractional soil wetness | unitless (0-1) | $f_\text{WHC}$; Used if `MODEL_WATER=0`; if `MODEL_WATER=1`, soil wetness is simulated|
352410

353-
#### Examples of climate files:
411+
#### Example `sipnet.clim` file:
354412

355413
Column names are not used, but are:
356414

@@ -376,7 +434,6 @@ loc year day time length tair tsoil par precip vpd vpdSoil vPress wspd s
376434
0 1998 306 7.00 0.417 2.2 1.4 2.7104 1.0000 114.1 71.6 741.8 0.9690 0.0000
377435
```
378436

379-
380437
### Agronomic Events
381438

382439
For managed ecosystems, the following inputs are provided in a file named `events.in` with the following columns:
@@ -469,7 +526,7 @@ Notes:
469526
- for perennials, some biomass may remain (col 5 + col 7 <= 1 and col 6 + col 8 <= 1; remainder is living).
470527
- root biomass is only removed for root crops
471528

472-
#### Example of events.in file:
529+
#### Example of `events.in` file:
473530

474531
```
475532
1 2022 40 irrig 5 1 # 5cm canopy irrigation on day 40
@@ -481,11 +538,14 @@ Notes:
481538

482539
#### Events output
483540

484-
SIPNET will create a file named `events.out` when event handling is enabled. This file
485-
will have one row for each agronomic event that is processed. Each row lists location,
541+
SIPNET will create a file named `events.out` when event handling is enabled.
542+
543+
544+
This file will have one row for each agronomic event that is processed. Each row lists location,
486545
time, event type, and parameter name/value pairs for all state variables that the event
487546
affects.
488547

548+
489549
Example events.out file below, with header enabled for clarity. Note the delimiters: spaces
490550
up to the param-values pairs, commas separating PV pairs, and `=` between param and value.
491551
```
@@ -497,6 +557,9 @@ loc year day type param_name=delta[,param_name=delta,...]
497557
1 2024 70 irrig fluxes.immedEvap=2.50,envi.soilWater=2.50
498558
1 2024 200 harv env.litter=4.25,envi.plantLeafC=-1.39,envi.plantWoodC=-1.63,envi.fineRootC=-2.52,envi.coarseRootC=-2.97
499559
```
560+
561+
_Note: `events.out` logs all parameters changed by an event for debugging and testing purposes;
562+
For downstream analyses that only need the date and event type, `events.in` is equivalent and easier to parse._
500563
## Outputs
501564

502565
| | Symbol | Parameter Name | Definition | Units |

0 commit comments

Comments
 (0)