|
Chapter 5 ....continued
|
5.13 |
Simulation control inputs |
The final step in preparing a Topog_dynamic simulation is the building of a 'parameter file', usually referred to as a .par file. This contains:
- the names of all the input files to be used
- values for selected input parameters
- flags to invoke particular computational routines in the model
- a list of requested data output dumps
The .par file is built (and can be edited) using Topog_simgen.
To run Topog_simgen, type: $ _simgenA large template will appear, such as shown below:
Not all of these fields need to be filled in, as the input requirements for the model vary, depending on:
- which modules you choose to run
- whether you distribute input parameters spatially or not
- what types of data ouput dumps you desire
The fields that must be filled in under all circumstances include:
- basename
- run identifier
- start time
- end time
- climate files
- soil files
- nodes files
- either restart .fin file -or - initial theta values
Note that many of the fields are set to default values and may need to be modified by you.
There are four basic control buttons, located on the left hand side of the template. These have the following functions:
OPEN PARAMETER FILE loads an existing .par file
SAVE PARAMETER FILE saves the screen entries to a new .par file
RUN SIMULATION runs Topog_dynamic, using the screen data as input
CLOSE ends the current Topog_simgen session
5.14 |
Outputs from Topog_dynamic |
There are three main kinds of output generated by Topog_dynamic:
- time series data: such as daily recharge over time for a given element, viewed using Topog_chart
- spatial data: such as a map of soil moisture storage across a catchment on a given timestep, viewed using Topog_display
- profile data: such as the vertical distribution of soil moisture for every catchment element at a given timestep, viewed using Topog_profile
All output files generated by Topog_dynamic are given:
- a basename
- a run number
- an extension
eg: myrtle.run1.log
myrtle.run1.finTime series data files for a specific element are also allocated:
- a locator (contour and element number)
eg: myrtle.run1.tstm.001.003
myrtle.run3.tirm.015.046Spatial data output files are also given:
- a qualifier
- a day number
eg: myrtle_y.run1.scar.0123
myrtle_y.run7.sevt.2539Profile data files are also given:
- a day number
eg: myrtle.run1.pprf.0123
myrtle.run4.pprf.1745
5.14.1 Time series data output by Topog_dynamic
There are ten different time series dumps, each viewable using program Topog_chart:
File extension Contents .tbal time series of the following attributes for a single element (or catchment average/total):
- gross rainfall (mm)
- net rainfall (mm)
- overstorey interception (mm)
- understorey interception (mm)
- litter interception (mm)
- total interception (mm)
- overstorey transpiration (mm)
- understorey transpiration (mm)
- total transpiration (mm)
- soil/litter evaporation (mm)
- total evapotranspiration (mm)
- subsurface flow (mm)
- overland flow (mm)
- total flow (mm)
- soil moisture storage (mm)
- soil salt storage (kg/m2)
- watertable depth (m)
- percent saturated area (%)
- overstorey leaf area index (m/m)
- understorey leaf area index (m/m)
- litter leaf area index (m/m)
- total leaf area index (m/m)
.tstm time series of the following attributes for a single element (or catchment average/total):
- net rainfall (mm)
- watertable depth (m)
- percent saturated area (%)
- soil moisture storage (mm)
- overland flow height (mm)
- overland flow velocity (m/s)
- sediment export (kg/m2)
- sediment entrainment (kg/m2)
- sediment deposition (kg/m2)
- net sediment flux (+/- kg/m2)
.tvmc time series of volumetric soil moisture content for every node in a specific catchment element (cm3/cm3) .tpsi time series of soil water potential for every node in a specific catchment element (m) .tsat time series of percent saturation for every node in a specific catchment element (0-1) .tsol time series of soil salt content content for every node in a specific catchment element (kg/m2) .troc time series of root carbon content for every node in a specific catchment element (kg/m2) .tirm time series of the following parameter values in the plant growth equation for a specific catchment element:
- overstorey soil water availability (0-1)
- overstorey light availability (0-1)
- overstorey temperature favourability (0-1)
- overstorey growth rate (0-1)
- overstorey carbon assimilation rate (kg/m2)
- understorey soil water availability (0-1)
- understorey light availability (0-1)
- understorey temperature favourability (0-1)
- understorey growth rate (0-1)
- understorey carbon assimilation rate (kg/m2)
.tcan time series of the following canopy values for a specific catchment element:
- overstorey radiation (KJ/m2/d)
- understorey radiation (KJ/m2/d)
- soil/litter layer radiation (KJ/m2/d)
- overstorey VPD (mbars)
- understorey VPD (mbars)
- soil/litter layer VPD (mbars)
- overstorey canopy conductance (m/s)
- understorey canopy conductance (m/s)
- soil/litter layer conductance (m/s)
.tcar time series of the following carbon balance values for a specific catchment element:
- overstorey leaf carbon (kg/m2)
- overstorey stem carbon (kg/m2)
- overstorey root carbon (kg/m2)
- understorey leaf carbon (kg/m2)
- understorey stem carbon (kg/m2)
- understorey root carbon (kg/m2)
- litter pool (kg/m2)
5.14.2 Spatial data output by Topog_dynamic
There are four different spatial dumps, each viewable using program Topog_display:
File extension Contents .sbal map of the following attributes for each element for a specified timestep:
- total evapotranspiration (mm)
- recharge to groundwater (mm)
- soil moisture storage (mm)
- soil salt storage (kg/m2)
- water table depth (m)
.sevt map of the following attributes for each element for a specified timestep:
- overstorey transpiration (mm)
- understorey transpiration (mm)
- total transpiration (mm)
- overstorey interception (mm)
- understorey interception (mm)
- litter interception (mm)
- total interception (mm)
- soil/litter evaporation (mm)
- overstorey evapotranspiration (mm)
- understorey evapotranspiration (mm)
- total evapotranspiration (mm)
.sstm map of the following attributes for each element for a specified timestep:
- overland flow height (mm)
- overland flow velocity (m/s)
- sediment export (kg/m2)
- sediment entrainment (kg/m2)
- sediment deposition (kg/m2)
- net sediment flux (+/- kg/m2)
.scar map of the following attributes for each element for a specified timestep:
- overstorey leaf area index (m/m)
- understorey leaf area index (m/m)
- litter leaf area index (m/m)
- total leaf area index (m/m)
- overstorey leaf carbon (kg/m2)
- understorey leaf carbon (kg/m2)
- litter leaf carbon (kg/m2)
- total leaf carbon (kg/m2)
- overstorey stem carbon (kg/m2)
- understorey stem carbon (kg/m2)
- total stem carbon (kg/m2)
- overstorey root carbon (kg/m2)
- understorey root carbon (kg/m2)
- total root carbon (kg/m2)
5.14.3 Profile data output by Topog_dynamic
There is one type of profile dump, viewable using program Topog_profile
File extension Contents .pprf profile of the following attributes for each element for a specified timestep:
- soil moisture content (cm3/cm3)
- soil water potential (m)
- soil salt content (kg/m2)
- root carbon (kg/m2)
5.14.4 Other data output by Topog_dynamic
There are two other types of output files:
File extension Contents .fin dump of all sub-surface state parameters at end of a simulation, including the following parameters for every node in every element: soil moisture content soil moisture solute concentration root carbon .log log of all simulation inputs and basic outputs The .fin file can be used to start a new simulation by listing it in the .par file, built using Topog_simgen (see section 5.13).
5.15 |
Running Topog_dynamic |
To run Topog_dynamic, type: $ _dynamic < parameter filename
| Take me out of frames | on to Chapter 6 .... |
