#-------------------------------------------------------------------------------
#
# Mesh009.in: Test of Child mesh synthesis from xyz file
# Sky Coyote, February 2008
#
# --- Inlet and 'no' rainfall ---
#     See notes in climate and inlet sections below
#
#-------------------------------------------------------------------------------

#-------------------------------------------------------------------------------
# Run control parameters
#
# The following parameters control the name and duration of the run along
# with a couple of other general settings.
#-------------------------------------------------------------------------------

OptInitMeshAndExit: Initialize and write mesh files and then exit
0

OUTFILENAME: Base name of output files
Mesh009/Mesh009

RUNTIME: Duration of run (years)
1000

OPINTRVL: Output interval (years)
100

#-------------------------------------------------------------------------------
# Mesh setup parameters
#
# These parameters control the initial configuration of the mesh. Here you
# specify whether a new or existing mesh is to be used; the geometry and
# resolution of a new mesh (if applicable); the boundary settings; etc.
#
# Notes:
#
#    OPTREADINPUT - controls the source of the initial mesh setup:
#                      0 = create a new mesh in a rectangular domain
#                      1 = read in an existing triangulation (eg, earlier run)
#                      2 = create a new mesh by triangulating a given set
#                          of (x,y,z,b) points
#    INPUTDATAFILE - use this only if you want to read in an existing
#                    triangulation, either from an earlier run or from
#                    a dataset.
#    INPUTTIME - if reading in a mesh from an earlier run, this specifies
#                the time slice number
#-------------------------------------------------------------------------------

OPTREADINPUT: reads data from inputdatafile if = 1
2

POINTFILENAME: xyz text file
Mesh009/xyz009.txt

INPUTTIME: which time to read from inputdatafile
0

OPTINITMESHDENS
0

#-------------------------------------------------------------------------------
# Climate parameters
#-------------------------------------------------------------------------------

OPTVAR: Option for rainfall variation
0

SEED: Random seed
021508

ST_PMEAN: Mean rainfall intensity (m/yr) (16.4 m/yr = Atlanta, GA)(instantaneous)
1.0e-12
# must have rainfall for an inlet to produce water

ST_STDUR: Mean storm duration (yr)
1.0

ST_ISTDUR: Mean time between storms (yr)
0.0

CRITICAL_FLOW: minimum flow for which we calc. meandering (m3/yr)
1.0e-6
# also used to plot discharge of mesh

#-------------------------------------------------------------------------------
# Inlet parameters
#-------------------------------------------------------------------------------

OPTINLET: 1 = add an "inlet" discharge boundary condition (0=none)
1

INDRAREA: drainage area at an inlet (if greater than zero ) (m2)
1e18
# inlet discharge is drain area * rain rate

INSEDLOAD: sediment load at inlet (m3/yr)
1000.0

INLET_X: x-coordinate of inlet (m)
484.974226119285
# at an existing node

INLET_Y: y-coordinate of inlet (m)
960.0
# at an existing node

#-------------------------------------------------------------------------------
# Erosion and sediment transport parameters
# (Note: choice of sediment-transport law is dictated at compile-time;
# see tErosion.h)
#
# Important notes on parameters:
#
#    (1) kb, kt, mb, nb and pb are defined as follows:
#           E = kb * ( tau - taucrit ) ^ pb,
#           tau = kt * q ^ mb * S ^ nb,
#           q = Q / W,  W = Wb ( Q / Qb ) ^ ws,  Wb = kw Qb ^ wb
#        where W is width, Q total discharge, Qb bankfull discharge,
#        Wb bankfull width. Note that kb, mb and nb are NOT the same as the
#        "familiar" K, m, and n as sometimes used in the literature.
#
#    (2) For power-law sediment transport, parameters are defined as follows:
#           capacity (m3/yr) = kf * W * ( tau - taucrit ) ^ pf
#           tau = kt * q ^ mf * S ^ nf
#        q is as defined above
#
#    (3) KT and TAUC are given in SI units -- that is, time units of seconds
#        rather than years. The unit conversion to erosion rate or capacity
#        is made within the code.
#-------------------------------------------------------------------------------

DETACHMENT_LAW:
1

KR: regolith erodibility coefficient (dimensions same as KB)
1e-5

KT:  Shear stress (or stream power) coefficient (in SI units)
9800.0

KB: bedrock erodibility coefficient (dimensions in m, kg, yr)
1e-5

MB: bedrock erodibility discharge exponent
1.0

NB: bedrock erodibility slope exponent
1.0

PB: Exponent on excess erosion capacity (e.g., excess shear stress)
1.0

TRANSPORT_LAW:
2

KF: sediment transport efficiency factor (dimensions vary)
0

MF: sediment transport capacity discharge exponent
2

NF: sed transport capacity slope exponent (ND)
2

TAUCD: critical shear stress for detachment-limited-erosion (kg/m/s^2)
0

KD: diffusivity coef
0.01

OPTDIFFDEP: if = 1 then diffusion only erodes, never deposits
1

DIFFUSIONTHRESHOLD: contributing area above which diffusive flux out is set to zero
1e8

#-------------------------------------------------------------------------------
# Bedrock and regolith
#-------------------------------------------------------------------------------

MAXREGDEPTH: maximum depth of a single regolith layer (also "active layer")
1000000.0

BEDROCKDEPTH: initial depth of bedrock (make this arbitrarily large)
1000000.0

REGINIT: initial regolith thickness
0.0

#-------------------------------------------------------------------------------
# Tectonics / baselevel boundary conditions
#-------------------------------------------------------------------------------

UPTYPE: type of uplift (0=none, 1=uniform, 2=block, etc)
0

#-------------------------------------------------------------------------------
# Grain size parameters
#
# (Note: for Wilcock sand-gravel transport formula, NUMGRNSIZE must be 2;
# otherwise, NUMGRNSIZE must be 1. Grain diameter has no effect if the
# Wilcock model is not used.)
#-------------------------------------------------------------------------------

NUMGRNSIZE: number of grain size classes 
1

GRAINDIAM1: representative diameter of first grain size class [m]
0.0005

REGPROPORTION1: proportion of sediments of grain size diam1 in regolith [.]
1.0

BRPROPORTION1: proportion of sediments of grain size diam1 in bedrock [.]
1.0

#-------------------------------------------------------------------------------
# Hydraulic geometry parameters
#
# Width is the most critical parameter as it is used in erosion and
# transport capacity calculations. HYDR_WID_COEFF_DS is the "kw" parameter
# referred to above (equal to bankfull width in m at unit bankfull discharge
# in cms)
#
# CHAN_GEOM_MODEL options are:
#    1 = empirical "regime" model: Wb = Kw Qb ^ wb, W / Wb = ( Q / Qb ) ^ ws
#    2 = Parker width closure: tau / tauc = const
#-------------------------------------------------------------------------------

CHAN_GEOM_MODEL: option for channel width closure
1

HYDR_ROUGH_COEFF_DS: coeff. on downstrm hydraulic roughness reln. (manning n)
5e-2

BANKFULLEVENT: precipitation rate of a bankfull event, in m/year
31.6

HYDR_WID_EXP_STN: exp. on at-a-station hydraulic width relation
0.25

HYDR_ROUGH_EXP_DS: exp. on downstream hydraulic roughness
-0.01

HYDR_ROUGH_EXP_STN: exp on at-a-station hydr. rough.
-0.21

BANK_ROUGH_COEFF: coeff. on downstream bank roughness relation 
3.0

BANK_ROUGH_EXP: exp on discharge for downstream bank roughness relation 
0.50

HYDR_WID_COEFF_DS: coeff. on downstream hydraulic width relation (m/(m3/s)^exp)
3.0

HYDR_DEP_COEFF_DS: coeff. on downstream hydraulic depth relation (m/(m3/s)^exp)
0.3417

HYDR_WID_EXP_DS: exponent on downstream hydraulic width relation 
0.5

HYDR_DEP_EXP_DS: exponent on downstream hydraulic depth relation 
0.4

HYDR_DEP_EXP_STN: exp. on at-a-station hydraulic depth relation
0.4

#-------------------------------------------------------------------------------
# Other various options
#-------------------------------------------------------------------------------

OPTFLOODPLAIN: option for overbank deposition using modified Howard 1992 model
0

OPTLOESSDEP: space-time uniform surface accumulation of sediment (loess)
0

OPTVEG: option for dynamic vegetation growth and erosion
0

OPTMESHADAPTAREA: dynamic adaptive meshing based on drainage area
0

OPTMESHADAPTDZ: dynamic adaptive meshing based on erosion rates
0

OPTINTERPLAYER: for node moving, do we care about tracking the layers? yes=1
0

OPTMEANDER: Option for meandering
0

OPTDETACHLIM: for detachment limited erosion = 1
1

FLOWGEN: flow generation, 1 = saturated flow
0

INFILTRATION: infiltration capacity (for Hortonian option) (m/yr)
0

OPTSINVARINFILT: option for sinusoidal variations in infiltration capacity
0

LAKEFILL: fill lakes if = 1
1

OPTTSOUTPUT: option for writing mean erosion rates, etc, at each time step
0

OPTSTRATGRID
0

OPTREADLAYER: option to read layer information from file (only if reading mesh)
0