Property:Describe input parameters model
From CSDMS
This is a property of type Text.
A
The area to be simulated has to be described (DEM, landuse). The meteorological input data (air temperature, relative humidity, precipitations...) have to be described (units, interpolations types). Some parameters about the model itself must be given (precision of the radiation ray tracing algorithms, characteristic lengths, parameters for a bucket model of runoff...) +
S
The bathymetry, current, water level, bottom friction and wind (if spatially variable) need to be provided to SWAN on so-called input grids. It is best to make an input grid so large that it completely covers the computational grid. +
T
The behavior of this component is controlled with a configuration (CFG) file, which may point to other files that contain input data. Here is a sample configuration (CFG) file for this component:
Method code: 1
Method name: Read_from_binary_file
Time step: Scalar 10800.00000000 (sec)
ET rate: Grid_Sequence Space-time_Rain_Test/Rain_TEST.rts (mm/hr) +
O
The default climate dataset used by OGGM is the Climatic Research Unit (CRU) TS v4.01 dataset +
W
The input data for WOFOST consists of three categories:
1. Daily weather variables (temperature, radiation, precipitation, humidity, windspeed)
2. Parameters for the crop, soil and site
3. Agromanagement information related to the cropping practices: sowing, harvesting, irrigation, nutrient application, etc.
How these inputs are provided to the model depends on the implementation. +
S
The input file is a DEM in .flt format. A driver text file is also required which contains the parameters used for the extraction. +
D
The input file is a DEM in .flt format. A driver text file is also required which contains the parameters used for the extraction. Information on the parameters needed in the driver file is available in the documentation (http://www.geos.ed.ac.uk/~smudd/LSDTT_docs/html/channel_heads.html). +
G
The input file is a text file and users are required to input:
Time (in model years)
dT (the time step in fractions of a year)
tauc (Critical Shear stress for portions of the channel that are vegetated in Pascals)
taucWepp (Critical Shear stress for portions of the channel that are soil in Pascals)
lenzone (the length of channel that is bare soil in Meters)
Pmmphr (Rainfall to be used for erosion in Millimeters per Hour)
tval (this is the number of loop iterations before a profile is saved as output)
Immphr (Infiltration to be used for erosion in Millimeters per Hour)
One additional input:
One must supply the input length of the channel as a matlab data array called xcell.mat +
A
The input file is required for each run and provides the basic simulation parameters. It must conform to the ACADIA.inp standard format as described below
Line1: should read the character string 'Comment�:' which is the label for the next line and is ignored during data input
Line2: inputs a comment of maximum 72 characters about the current simulation
Line3: should read 'Case name', which is the label for the next line and is ignored during data input
Line4: specifes the case name (including the directory location if it is not soft-linked to the current directory). The code will append the suffixes '.nod',
'.ele',and '.bat' to this string in order to know where to find the NML standard files
Line5: should read 'Simulation Parameters', which is the label for the next line and is ignored during data input
Line6: inputs the number of transport variables
Line7: inputs the number of �uid layers
Line8: inputs the degrees latitude of the center of the mesh
Line9: inputs the time step in seconds
Line10: inputs the starting date (day, month, year), time (in seconds after midnight). The value of time can be greater than 8.64E4 (equivalent of 1 day) as the code will automatically adjust the day and time accordingly
Last Line: inputs the quitting criteria for the code as day, month, year and time in seconds after midnight. Again the value of time can be greater than 8.64E4. This allows the user to specify an overall length of the simulation instead of the exact date and time of the end of the simulation +
H
The input files are a DEM in .flt format, a channel heads file generated using the DrEICH algorithm (https://csdms.colorado.edu/wiki/Model:DrEICH_algorithm) and an optional floodplain mask in .flt format. Input parameters are also supplied at the command line. Information on the parameters is available in the documentation (http://www.geos.ed.ac.uk/~smudd/LSDTT_docs/html/channel_heads.html). +
R
The input includes file names for output, subgrid information, physical parameters and wave conditions. Water depth is obtained from the master program though a file name for water depth input is still kept in "indat.dat". +
C
The input is a 'channel file' and a 'driver file'. The channel file contains data on channel profiles within a channel network composed of a main stem and tributaries flowing into that main stem (that is, there are no tributaries of tributaries). The driver file contains parameters for the model run. The format of these files is described in the documentation that accompanies the model source code. +
The input parameters to a model run consist of an initial shoreline, a wave file, and a set of configuration parameters. The initial shoreline is stored within a custom binary formatted-file. Since CEM has been used for abstract simulations of coastline evolution, the initial model condition consists either of a mostly-smooth shoreline with initial perturbations to the shoreline position (generated by a tool provided with the model), or using a shoreline that resulted from a previous model run. The wave file consists of a set of wave approach angles and wave heights that are used during the model run. This wave file is also generated by a tool provided with the model, and takes as input the statistical distribution of wave-approach angles. Finally, basic model parameters (e.g. number of time steps to simulate, etc.) are specified within an XML-formatted text file. An example is provided with the model. +
T
The input variables for modeling the net flux of shortwave radiation are defined as follows:
Tair = air temperature (deg C)
RH = relative humidity (unitless) (in (0,1))
albedo = surface albedo (unitless) (in (0,1))
dust att. = dust attenuation factor (unitless) (in (0,1))
factor = cloud or canopy cover factor (unitless) (in (0,1))
slope = topographic slope (unitless, m/m) (in (0,Infinity))
aspect = aspect angle (radians) (in (0,1))
The behavior of this component is controlled with a configuration (CFG) file, which may point to other files that contain input data. Here is a sample configuration (CFG) file for this component:
Method code: 1
Method name: Standard
Time step: Scalar 3600.0 (sec)
rho_H2O: Scalar 1000.00000000 (kg/m^3)
Cp_air: Scalar 1005.70001221 (J/kg/K)
rho_air: Scalar 1.26139998 (kg/m^3)
P: Time_series Case5_rain_rates.txt (mm/hr)
T_air: Scalar 20.00000000 (deg C)
T_surf: Scalar -5.00000000 (deg C)
RH: Scalar 0.50000000 (none)
p0: Scalar 1000.00000000 (mbar)
uz: Scalar 3.00000000 (m/s)
z: Scalar 10.00000000 (m)
z0_air: Scalar 0.02000000 (m)
Qn_SW: Scalar 100.00000000 (W/m^2)
Qn_LW: Scalar 10.00000000 (W/m^2)
Save grid timestep: Scalar 60.00000000 (sec)
Save ea grids: 0 Case5_2D-ea.rts (mbar)
Save es grids: 0 Case5_2D-es.rts (mbar)
Save pixels timestep: Scalar 60.00000000 (sec)
Save ea pixels: 0 Case5_0D-ea.txt (mbar)
Save es pixels: 0 Case5_0D-es.txt (mbar)
The input variables for the Energy Balance method of estimating losses due to evaporation are defined as follows:
Q_SW = net shortwave radiation (W / m^2)
Q_LW = net longwave radiation (W / m^2)
T_air = air temperature (deg C)
T_surf = surface (snow) temperature (deg C)
T_soil_x = soil temperature at depth x (deg C)
x = reference depth in soil (m)
K_soil = thermal conductivity of soil (W / (m deg_C))
u_z = wind velocity at height z (m / s)
z = reference height for wind (m) (above land surface)
z_0 = surface roughness height (m) (with no snow)
h0_snow = initial snow depth (m)
ρ_air = density of the air (kg / m^3)
c_air = specific heat capacity of air (J / (kg deg_C))
L_v = latent heat of vaporization, water (J / kg) (2500000)
g = gravitational constant, Earth = 9.81 (m / s^2)
κ = von Karman's constant = 0.41 (unitless)
The behavior of this component is controlled with a configuration (CFG) file, which may point to other files that contain input data. Here is a sample configuration (CFG) file for this component:
Method code: 2
Method name: Energy_Balance
Time step: Scalar 3600.00000000 (sec)
alpha: Scalar 1.20000000 (none)
K_soil: Scalar 0.44999999 (W/m/deg_C)
soil_x: Scalar 0.05000000 (m)
T_soil_x: Scalar 0.00000000 (deg C)
Save grid timestep: Scalar 60.00000000 (sec)
Save er grids: 0 Case5_2D-ETrate.rts (m/s)
Save pixels timestep: Scalar 60.00000000 (sec)
Save er pixels: 0 Case5_0D-ETrate.txt (m/s) +
The input variables for the Energy Balance method of estimating runoff due to snowmelt are defined as follows:
Q_SW = net shortwave radiation (W / m^2)
Q_LW = net longwave radiation (W / m^2)
T_air = air temperature (deg C)
T_surf = surface (snow) temperature (deg C)
RH = relative humidity (none) (in (0,1))
p_0 = atmospheric pressure (mbar)
u_z = wind velocity at height z (m / s)
z = reference height for wind (m)
z0_air = surface roughness height (m)
h0_snow = initial snow depth (m)
h0_swe = initial depth, snow water equivalent (m)
ρ_snow = density of the snow (kg / m^3)
c_snow = specific heat capacity of snow (J / (kg deg_C))
ρ_air = density of the air (kg / m^3)
c_air = specific heat capacity of air (J / (kg deg_C))
L_f = latent heat of fusion, water (J / kg) (334000)
L_v = latent heat of vaporization, water (J / kg) (2500000)
e_air = air vapor pressure at height z (mbar)
e_surf = vapor pressure at the surface (mbar)
g = gravitational constant = 9.81 (m / s^2)
κ = von Karman's constant = 0.41 (unitless)
The behavior of this component is controlled with a configuration (CFG) file, which may point to other files that contain input data. Here is a sample configuration (CFG) file for this component:
Method code: 2
Method name: Energy_Balance
Time step: Scalar 3600.00000000 (sec)
Cp_snow: Scalar 2090.00000000 (J/kg/K)
rho_snow: Scalar 300.00000000 (kg/m^3)
c0: Scalar 2.70000005 (mm/day/deg C)
T0: Scalar -0.20000000 (deg C)
h0_snow: Scalar 0.50000000 (m)
h0_swe: Scalar 0.15000000 (m)
Save grid timestep: Scalar 60.00000000 (sec)
Save mr grids: 0 Case5_2D-SMrate.rts (m/s)
Save hs grids: 0 Case5_2D-hsnow.rts (m)
Save sw grids: 0 Case5_2D-hswe.rts (m)
Save cc grids: 0 Case5_2D-Ecc.rts (J/m^2)
Save pixels timestep: Scalar 60.00000000 (sec)
Save mr pixels: 0 Case5_0D-SMrate.txt (m/s)
Save hs pixels: 0 Case5_0D-hsnow.txt (m)
Save sw pixels: 0 Case5_0D-hswe.txt (m)
Save cc pixels: 0 Case5_0D-Ecc.txt (J/m^2)
The input variables for the Priestley-Taylor method of estimating losses due to evaporation are defined as follows:
Q_SW = net shortwave radiation (W / m^2)
Q_LW = net longwave radiation (W / m^2)
T_air = air temperature (deg C)
T_surf = surface (snow) temperature (deg C)
T_soil_x = soil temperature at depth x (deg C)
x = reference depth in soil (m)
K_soil = thermal conductivity of soil (W / (m deg_C))
α = coefficient (unitless)
L_v = latent heat of vaporization, water (J / kg) (2500000)
The behavior of this component is controlled with a configuration (CFG) file, which may point to other files that contain input data. Here is a sample configuration (CFG) file for this component:
Method code: 0
Method name: Priestley-Taylor
Time step: Scalar 3600.00000000 (sec)
alpha: Scalar 1.20000000 (none)
K_soil: Scalar 0.44999999 (W/m/deg_C)
soil_x: Scalar 0.05000000 (m)
T_soil_x: Scalar 0.00000000 (deg C)
Save grid timestep: Scalar 60.00000000 (sec)
Save er grids: 0 Case5_2D-ETrate.rts (m/s)
Save pixels timestep: Scalar 60.00000000 (sec)
Save er pixels: 0 Case5_0D-ETrate.txt (m/s) +
The input variables used by the Smith-Parlange 3-parameter method for modeling infiltration are defined as follows:
K_s = saturated hydraulic conductivity (m / s)
K_i = initial hydraulic conductivity (m / s) (typically much less than K_s)
θ_s = soil water content at ψ=0 (unitless) (typically set to the porosity, φ)
θ_i = initial soil water content (unitless)
G = capillary length scale (meters)
= integral over all ψ of K(ψ) / K_s
= (almost always between ψ_B and 2*ψ_B)
P = precipitation rate (mm / sec)
M = snowmelt rate (mm / sec)
γ = Smith-Parlange method parameter (between 0 and 1, near 0.8)
The behavior of this component is controlled with a configuration (CFG) file, which may point to other files that contain input data. Here is a sample configuration (CFG) file for this component:
Method code: 2
Method name: Smith_Parlange
Number of layers: 1
Time step: Scalar 60.0000000 (sec)
Ks: Scalar 0.00000720 (m/s)
Ki: Scalar 0.00000010 (m/s)
qs: Scalar 0.48500001 (none)
qi: Scalar 0.37580763 (none)
G: Scalar 0.72400000 (m)
gamma: Scalar 0.82000000 (none)
Closest soil_type: silt_loam
Save grid timestep: Scalar 60.00000000 (sec)
Save v0 grids: 0 Case5_2D-v0.rts (m/s)
Save I grids: 0 Case5_2D-I.rts (m)
Save pixels timestep: Scalar 60.00000000 (sec)
Save v0 pixels: 0 Case5_0D-v0.txt (m/s)
Save I pixels: 0 Case5_0D-I.txt (m) +
The input variables used by the simple Green-Ampt method for modeling infiltration are defined as follows:
K_s = saturated hydraulic conductivity (m / s)
K_i = initial hydraulic conductivity (m / s) (typically much less than Ks)
θ_s = soil water content at ψ=0 (unitless) (typically set to the porosity, φ)
θ_i = initial soil water content (unitless)
G = capillary length scale (meters)
= integral over all ψ of K(ψ) / K_s
= (almost always between ψ_B and 2*ψ_B)
P = precipitation rate (mm / sec)
M = snowmelt rate (mm / sec)
The behavior of this component is controlled with a configuration (CFG) file, which may point to other files that contain input data. Here is a sample configuration (CFG) file for this component:
Method code: 0
Method name: Green_Ampt
Number of layers: 1
Time step: Scalar 60.0000000 (sec)
Ks: Scalar 0.00000720 (m/s)
Ki: Scalar 0.00000010 (m/s)
qs: Scalar 0.48500001 (none)
qi: Scalar 0.37580763 (none)
G: Scalar 0.72400000 (m)
Closest soil_type: silt_loam
Save grid timestep: Scalar 60.00000000 (sec)
Save v0 grids: 0 Case5_2D-v0.rts (m/s)
Save I grids: 0 Case5_2D-I.rts (m)
Save pixels timestep: Scalar 60.00000000 (sec)
Save v0 pixels: 0 Case5_0D-v0.txt (m/s)
Save I pixels: 0 Case5_0D-I.txt (m) +
The input variables used for modeling horizontal subsurface flow in the saturated zone via Darcy's Law are defined as follows:
K_s = saturated hydraulic conductivity (m / s)
S_y = specific yield or drainable porosity (unitless)
(less than or equal to the porosity, φ, see Notes)
thickness = soil layer thickness (meters) +
