-
Notifications
You must be signed in to change notification settings - Fork 4
/
settings.txt
99 lines (90 loc) · 2.46 KB
/
settings.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
#-------------------------------------------------------------------------------
#grid settings
#depth of domain (m)
depth = 1
#cell depth at surface (m)
delz0 = 0.02
#fraction increase in cell depth
delzfrac = 1.025
#maximum cell depth
delzmax = 0.1
#output grid variables to output directory?
save_grid = True
#-------------------------------------------------------------------------------
#model setup and integration settings
#integration time (a multiple of tunit seconds)
tint = 10
#number of seconds to use for tint (seconds) an Earth year is 31557600
tunit = 3600
#number of snapshots to write to output directory
nsnap = 5
#maximum length of tracker variables which capture every step
nmaxout = 1e8
#safety factor applied to maximum stable time step
dtfac = 0.3
#-------------------------------------------------------------------------------
#physical parameters
#porosity
poro = 0.4
#permeability (m^2)
perm = 3e-13
#gravity (m/s^2)
g = 9.8
#water viscosity (Pa*s)
mu = 9e-4
#water density (kg/m^3)
rho = 1e3
#Brooks-Corey parameter
b = 4
#wilting saturation or minimum saturation fraction (as fraction of porosity)
wilt = 0.1
#-------------------------------------------------------------------------------
#periodic forcing
#surface evaporation time scale for exponential (seconds)
tauevap = 1800
#surface evaporation length scale (m)
Levap = 0.01
#infiltration period (seconds)
infper = 172800
#infiltration duration (seconds)
infdur = 3600
#-------------------------------------------------------------------------------
#tracker and output settings
#porosity at cell centers
poroc = True
#porosity at cell edges
poroe = True
#saturated hydraulic conductivity at cell edges (m/s)
Ksat = True
#saturated matric head at cell edges (m)
psisat = True
#derivative of matric head w/r/t water fraction at cell edges
dpsidw = True
#derivative of water fraction w/r/t depth
dwdz = True
#unsaturated hydraulic conductivity at cell edges (m/s)
K = True
#unsaturated "diffusivity" at cell edges (m^2)
D = True
#water fraction at cell centers
w = True
#water fraction at cell edges
we = True
#all water fractions
wall = True
#water flux at cell edges at snapshots
q = True
#water flux at surface at all time steps
qtop = True
#water flux at approximate middle of domain at all time steps
qmid = True
#water flux at bottom of domain at all time steps
qbot = True
#water flux at all edges at all time steps
qall = True
#whether to track infiltration flag
infil = True
#time
t = True
#snapshot times
tsnap = True