Skip to content

Commit cc52152

Browse files
committed
Added 4 csub package example problems
1 parent 0876241 commit cc52152

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

63 files changed

+8758
-0
lines changed
Binary file not shown.
Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,2 @@
1+
This is the Jacob (1939) train problem simulates elastic compaction of aquifer materials in response to the loading of an aquifer by
2+
a passing train. The problem is described in the MODFLOW 6 CSUB Package Example Problems document released with MODFLOW 6 (version 6.1.0 or higher). This problem demonstrates elastic compaction of coarse-grained aquifer materials and time-varying loading at the surface to simulate the train passing an observation well in MODFLOW 6. The default effective-stress formulation is used to simulate coarse-grained material compaction.

mf6/test061_csub_jacob/fig4_base.csub

Lines changed: 34 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,34 @@
1+
BEGIN options
2+
PRINT_INPUT
3+
SAVE_FLOWS
4+
GAMMAW 9806.65000000
5+
BETA 4.65120000E-10
6+
UPDATE_MATERIAL_PROPERTIES
7+
TS6 FILEIN fig4_base.load.ts
8+
OBS6 FILEIN fig4_base.csub.obs
9+
END options
10+
11+
BEGIN dimensions
12+
NINTERBEDS 0
13+
MAXSIG0 1
14+
END dimensions
15+
16+
BEGIN griddata
17+
cg_ske_cr LAYERED
18+
CONSTANT 3.30000000E-05
19+
CONSTANT 6.60000000E-04
20+
CONSTANT 4.49913843E-07
21+
cg_theta LAYERED
22+
CONSTANT 0.25000000
23+
CONSTANT 0.50000000
24+
CONSTANT 0.30000000
25+
sgm
26+
CONSTANT 1.70000000
27+
sgs
28+
CONSTANT 2.00000000
29+
END griddata
30+
31+
BEGIN period 1
32+
1 1 1 load
33+
END period 1
34+
Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
BEGIN options
2+
DIGITS 10
3+
PRINT_INPUT
4+
END options
5+
6+
BEGIN continuous FILEOUT csub_obs.csv
7+
tc1 compaction-cell 1 1 1
8+
tc2 compaction-cell 2 1 1
9+
tc3 compaction-cell 3 1 1
10+
wc1 compaction-cell 1 1 12
11+
wc2 compaction-cell 2 1 12
12+
wc3 compaction-cell 3 1 12
13+
tes3 estress-cell 3 1 1
14+
tgs3 gstress-cell 3 1 1
15+
END continuous FILEOUT csub_obs.csv
16+

mf6/test061_csub_jacob/fig4_base.dis

Lines changed: 24 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,24 @@
1+
BEGIN options
2+
END options
3+
4+
BEGIN dimensions
5+
NLAY 3
6+
NROW 1
7+
NCOL 35
8+
END dimensions
9+
10+
BEGIN griddata
11+
delr
12+
INTERNAL FACTOR 1.0
13+
0.50000000 0.60000000 0.72000000 0.86400000 1.03680000 1.24416000 1.49299200 1.79159040 2.14990848 2.57989018 3.09586821 3.71504185 4.45805022 5.34966027 6.41959232 7.70351079 9.24421294 11.09305553 13.31166664 15.97399997
14+
19.16879996 23.00255995 27.60307195 33.12368633 39.74842360 47.69810832 57.23772999 68.68527598 82.42233118 98.90679742 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000
15+
delc
16+
CONSTANT 100.58491653
17+
top
18+
CONSTANT 0.00000000
19+
botm LAYERED
20+
CONSTANT -12.20000000
21+
CONSTANT -21.30000000
22+
CONSTANT -30.50000000
23+
END griddata
24+

mf6/test061_csub_jacob/fig4_base.ic

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,5 @@
1+
BEGIN griddata
2+
strt
3+
CONSTANT -10.66809721
4+
END griddata
5+

mf6/test061_csub_jacob/fig4_base.ims

Lines changed: 20 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
BEGIN options
2+
PRINT_OPTION summary
3+
END options
4+
5+
BEGIN nonlinear
6+
OUTER_HCLOSE 1.00000000E-09
7+
OUTER_MAXIMUM 500
8+
UNDER_RELAXATION none
9+
END nonlinear
10+
11+
BEGIN linear
12+
INNER_MAXIMUM 300
13+
INNER_HCLOSE 1.00000000E-09
14+
inner_rclose 1.00000000E-06
15+
LINEAR_ACCELERATION bicgstab
16+
RELAXATION_FACTOR 1.00000000
17+
SCALING_METHOD none
18+
REORDERING_METHOD none
19+
END linear
20+
Lines changed: 128 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,128 @@
1+
BEGIN attributes
2+
NAMES load
3+
METHODS linear
4+
SFACS 1.05000000
5+
END attributes
6+
7+
BEGIN timeseries
8+
0.00000000 0.80960210
9+
0.50000000 0.80960210
10+
1.00000000 0.80960211
11+
1.50000000 0.80960212
12+
2.00000000 0.80960210
13+
2.50000000 0.80960209
14+
3.00000000 0.80960216
15+
3.50000000 0.80960220
16+
4.00000000 0.80960191
17+
4.50000000 0.80960183
18+
5.00000000 0.80960287
19+
5.50000000 0.80960331
20+
6.00000000 0.80959984
21+
6.50000000 0.80959806
22+
7.00000000 0.80961312
23+
7.50000000 0.80961853
24+
8.00000000 0.80956902
25+
8.50000000 0.80954718
26+
9.00000000 0.80976097
27+
9.50000000 0.80982555
28+
10.00000000 0.80912043
29+
10.50000000 0.80885893
30+
11.00000000 0.81258710
31+
11.50000000 0.82110941
32+
12.00000000 0.83458525
33+
12.50000000 0.85302994
34+
13.00000000 0.87985583
35+
13.50000000 0.91277717
36+
14.00000000 0.94747637
37+
14.50000000 0.98720951
38+
15.00000000 1.03806423
39+
15.50000000 1.09509371
40+
16.00000000 1.15129288
41+
16.50000000 1.21259107
42+
17.00000000 1.28778579
43+
17.50000000 1.36923460
44+
18.00000000 1.44721035
45+
18.50000000 1.53035020
46+
19.00000000 1.61866946
47+
19.50000000 1.72436087
48+
20.00000000 1.83640435
49+
20.50000000 1.94166246
50+
21.00000000 2.05216507
51+
21.50000000 2.18289205
52+
22.00000000 2.31966290
53+
22.50000000 2.44340071
54+
23.00000000 2.56387296
55+
23.50000000 2.67594841
56+
24.00000000 2.78538404
57+
24.50000000 2.87502349
58+
25.00000000 2.93809402
59+
25.50000000 2.98723141
60+
26.00000000 3.02712313
61+
26.50000000 3.05294024
62+
27.00000000 3.06447585
63+
27.50000000 3.06648920
64+
28.00000000 3.05838992
65+
28.50000000 3.04137486
66+
29.00000000 3.02001031
67+
29.50000000 2.99453412
68+
30.00000000 2.96647691
69+
30.50000000 2.93288966
70+
31.00000000 2.89703581
71+
31.50000000 2.86273273
72+
32.00000000 2.82580058
73+
32.50000000 2.78059800
74+
33.00000000 2.73129875
75+
33.50000000 2.68344228
76+
34.00000000 2.63147749
77+
34.50000000 2.56682919
78+
35.00000000 2.49376546
79+
35.50000000 2.41921903
80+
36.00000000 2.33401911
81+
36.50000000 2.23878921
82+
37.00000000 2.12369112
83+
37.50000000 2.00522986
84+
38.00000000 1.90071814
85+
38.50000000 1.79813777
86+
39.00000000 1.68437695
87+
39.50000000 1.57144785
88+
40.00000000 1.47110748
89+
40.50000000 1.37066923
90+
41.00000000 1.27084806
91+
41.50000000 1.16039656
92+
42.00000000 1.05289011
93+
42.50000000 0.96104275
94+
43.00000000 0.87394114
95+
43.50000000 0.78345589
96+
44.00000000 0.70211315
97+
44.50000000 0.63826703
98+
45.00000000 0.58114009
99+
45.50000000 0.52349548
100+
46.00000000 0.47102555
101+
46.50000000 0.42831718
102+
47.00000000 0.38899887
103+
47.50000000 0.35286725
104+
48.00000000 0.31595704
105+
48.50000000 0.28288943
106+
49.00000000 0.25705666
107+
49.50000000 0.23464476
108+
50.00000000 0.21370591
109+
50.50000000 0.19623844
110+
51.00000000 0.18196131
111+
51.50000000 0.16698585
112+
52.00000000 0.14726014
113+
52.50000000 0.12500568
114+
53.00000000 0.10543738
115+
53.50000000 0.08853821
116+
54.00000000 0.07624396
117+
54.50000000 0.06685801
118+
55.00000000 0.05845986
119+
55.50000000 0.04878943
120+
56.00000000 0.03599087
121+
56.50000000 0.02000816
122+
57.00000000 0.00514449
123+
57.50000000 0.00000000
124+
58.00000000 0.00000000
125+
58.50000000 0.00000000
126+
59.00000000 0.00000000
127+
END timeseries
128+

mf6/test061_csub_jacob/fig4_base.nam

Lines changed: 14 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,14 @@
1+
BEGIN options
2+
NEWTON
3+
END options
4+
5+
BEGIN packages
6+
dis6 fig4_base.dis dis
7+
ic6 fig4_base.ic ic
8+
npf6 fig4_base.npf npf
9+
obs6 fig4_base.obs head_obs
10+
sto6 fig4_base.sto sto
11+
csub6 fig4_base.csub csub
12+
oc6 fig4_base.oc oc
13+
END packages
14+

mf6/test061_csub_jacob/fig4_base.npf

Lines changed: 20 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
BEGIN options
2+
SAVE_FLOWS
3+
SAVE_SPECIFIC_DISCHARGE
4+
END options
5+
6+
BEGIN griddata
7+
icelltype LAYERED
8+
CONSTANT 1
9+
CONSTANT 0
10+
CONSTANT 0
11+
k LAYERED
12+
CONSTANT 1.80000000E-05
13+
CONSTANT 3.50000000E-10
14+
CONSTANT 3.10070996E-05
15+
k33 LAYERED
16+
CONSTANT 1.80000000E-05
17+
CONSTANT 3.50000000E-10
18+
CONSTANT 3.10070996E-05
19+
END griddata
20+

mf6/test061_csub_jacob/fig4_base.obs

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,9 @@
1+
BEGIN options
2+
DIGITS 10
3+
PRINT_INPUT
4+
END options
5+
6+
BEGIN continuous FILEOUT gwf_calib_obs.csv
7+
w3_1_1 head 3 1 12
8+
END continuous FILEOUT gwf_calib_obs.csv
9+

mf6/test061_csub_jacob/fig4_base.oc

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
BEGIN options
2+
BUDGET FILEOUT fig4_base.cbc
3+
HEAD FILEOUT fig4_base.hds
4+
HEAD PRINT_FORMAT COLUMNS 10 WIDTH 15 DIGITS 6 general
5+
END options
6+
7+
BEGIN period 1
8+
SAVE head all
9+
SAVE budget all
10+
PRINT head all
11+
PRINT budget all
12+
END period 1
13+

mf6/test061_csub_jacob/fig4_base.sto

Lines changed: 22 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,22 @@
1+
BEGIN options
2+
END options
3+
4+
BEGIN griddata
5+
iconvert
6+
CONSTANT 1
7+
ss
8+
CONSTANT 0.00000000
9+
sy LAYERED
10+
CONSTANT 0.10000000
11+
CONSTANT 0.05000000
12+
CONSTANT 0.25000000
13+
END griddata
14+
15+
BEGIN period 1
16+
STEADY-STATE
17+
END period 1
18+
19+
BEGIN period 2
20+
TRANSIENT
21+
END period 2
22+

mf6/test061_csub_jacob/fig4_base.tdis

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
BEGIN options
2+
TIME_UNITS seconds
3+
END options
4+
5+
BEGIN dimensions
6+
NPER 2
7+
END dimensions
8+
9+
BEGIN perioddata
10+
0.50000000 1 1.00000000
11+
58.50000000 117 1.00000000
12+
END perioddata
13+

mf6/test061_csub_jacob/mfsim.nam

Lines changed: 18 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,18 @@
1+
BEGIN options
2+
END options
3+
4+
BEGIN timing
5+
TDIS6 fig4_base.tdis
6+
END timing
7+
8+
BEGIN models
9+
gwf6 fig4_base.nam fig4_base
10+
END models
11+
12+
BEGIN exchanges
13+
END exchanges
14+
15+
BEGIN solutiongroup 1
16+
ims6 fig4_base.ims fig4_base
17+
END solutiongroup 1
18+
7.42 KB
Binary file not shown.
Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
This is the Problem 1 from the Subsidence Package for MODFLOW-2000 documentation (Hoffman and others, 2003), which is described in the MODFLOW 6 CSUB Package Example Problems document released with MODFLOW 6 (version 6.1.0 or higher). This simulates the drainage of a thick interbed caused by a step decrease in hydraulic head in the aquifer in MODFLOW 6. The thick interbed is simulated using a delay-bed interbed and the default effective-stress formulation is used to simulate thick interbed compaction.

mf6/test062_csub_sub01/mfsim.nam

Lines changed: 18 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,18 @@
1+
BEGIN options
2+
END options
3+
4+
BEGIN timing
5+
TDIS6 sub01es.tdis
6+
END timing
7+
8+
BEGIN models
9+
gwf6 sub01es.nam sub01es
10+
END models
11+
12+
BEGIN exchanges
13+
END exchanges
14+
15+
BEGIN solutiongroup 1
16+
ims6 sub01es.ims sub01es
17+
END solutiongroup 1
18+

mf6/test062_csub_sub01/sub01es.chd

Lines changed: 12 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,12 @@
1+
BEGIN options
2+
END options
3+
4+
BEGIN dimensions
5+
MAXBOUND 2
6+
END dimensions
7+
8+
BEGIN period 1
9+
1 1 1 0.00000000
10+
1 1 3 0.00000000
11+
END period 1
12+

0 commit comments

Comments
 (0)