forked from nthiros/ER_PLM_ParFlow
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathRunLog.txt
More file actions
144 lines (112 loc) · 5.83 KB
/
Copy pathRunLog.txt
File metadata and controls
144 lines (112 loc) · 5.83 KB
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
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
Directories:
NOTE: RunB.0.mint is the mint directly as of 04/19/2022
RunA.0 -- 'Base' case model from Erica work -- wshale=5.605e-6, fshale=5.605e-08
RunA.0.clmup.sm -- 'Base' model including smoothed DEM and CLM vegetation types
RunA.1 -- High K -- wshale=5.605e-6, fshale=5.605e-07
RunA.2 -- Low K -- wshale=5.605e-6, fshale=5.605e-09
RunA.3 -- Exponential K with low decay rate -- wshale=5.605e-6, fshale_top=5.605e-07, fshale_bottom=5.605e-9
RunA.4 -- Exponential K with high decay rate -- wshale=5.605e-6, fshale_top=5.605e-07, fshale_bottom=5.605e-12
RunB.0 -- 'Base' model including smoothed DEM and CLM vegetation types -- same as RunA.0.clmup.sm (I think), Use this dir as base case
RunB.0.mint -- same as above (Normal RunB.0 was corrupted somehow)
RunB.1 -- Exponential K with medium decay rate -- wshale=5.605e-06, fshale_top=5.605e-07, fshale_bottom=5.605e-10
RunB.2 -- Exponential K with low weathered shale K -- wshale=5.605e-07, fshale_top=5.605e-07, fshale_bottom=5.605e-10
RunB.3 -- Exponential K with high weathered shale K -- wshale=5.605e-06, fshale_top=5.605e-06, fshale_bottom=5.605e-10
RunB.4 -- Exponential K with high weathered shale K and high exponential decay rate
-- wshale=5.605e-6, fshale_top=5.605e-6, fshale_bottom=5.605e-12
RunC.0 -- 'Base' model with shallower (3m) soil layer, bedrock at (9m) depth
RunC.1 -- Exponential K with med decay rate -- wshale=5.605e-06, fshale_top=5.605e-07, fshale_bottom=5.605e-10
RunC.2 -- Exponential K with low weathered shale K -- wshale=5.605e-07, fshale_top=5.605e-07, fshale_bottom=5.605e-10
RunD.0 -- (2 layer model) Thin soil layer (2mbls) and bedrock at 4mbls
-- med K bedrock -- wshale=5.605e-7, fshale=5.605e-7
RunD.1 -- low K bedrock -- wshale=5.605e-8, fshale=5.605e-8
RunD.2 -- high K bedrock -- wshale=5.605e-6, fhsale=5.605e-6
Completion Notes:
Done: A.1, A.2, A.3, A.4 -- Note that A.0 and B.0 are the same, but B.0 is more up-to-date
B.0, B.1, B.2, B.3, B.4
C.0, C.1, C.2
Not Done: D.0, D.1, D.2
Workflow Outline:
ParFlow-CLM:
(1): Initial Conditions Pressure Field
(a): 1000 year initial conditions run with constant recharge (transect_ic.tcl)
(b): Restart for 50 year simulation with CLM turned on (transect_spinupA.tcl)
NOTE: only doing these runs once, then will use this as IC pressure throughout
(2): Spinup:
(a): 10 years using average MET data (transect_spinup.tcl)
(b): Dump at 24 hours
(3): Transient 2000-2016
(a): Dump at 24 hours
(4): Transient 2017-2021
(a): Dump at 3 hours
Ecoslim:
(1): Spinup1
(a): ParFlow spinup at 240 hour velocity fields
(b): Run for 1000 years
(2): Spinup2
(a): Restart ecoslim spinup1 using ParFlow spinup at 24 hour velocity fields
(b): Run for 50 years
(3): 2000-2016
(a): 24 hour velocity fields
(4): 2016-2021
(a): 24 hour velocity fields
Practical Workflow with Scripts:
Simulation Setup:
Running Simulation:
ParFlow
run_both.sc -- runs both parflow and ecoslim
run_tcl.sc -- runs parflow
EcoSLIM:
run_both.sc -- runs both parflow and ecoslim
Or
run_eco.sc -- runs just ecoslim
Or
loop_eco.sc -- perfrom run_eco.short.sc
run_eco.short.sc -- only re-runs slimin.txt.2017_2021
Local Analysis -- after the full parflow/ecoslim is completed
Run: python pull_parflow_wl_v5.py
Run: python read_vtk_v4.py
Run: tclsh vis_vtk.tcl
Transferring Files -- only select files because .pfb are huge
(1) rsync full directory minus big files from server to local
- use 'utils/rsync_helper.txt'
(2) manually copy files listed below
-----------
Plotting
-----------
Files to Copy:
Use utils 'rsync_helper.txt'
ParFlow Fields in 'wy_2017_2021' dir
wy_2017_2021.out.velx.01683.pfb
wy_2017_2021.out.vely.01683.pfb
wy_2017_2021.out.velz.01683.pfb
wy_2017_2021.out.satur.01683.pfb
wy_2017_2021.out.evaptrans.01683.pfb
wy_2017_2021.out.porosity.pfb
wy_2017_2021.out.press.01683.pfb
wy_2017_2021.out.specific_storage.pfb
Use: rsync -auvz --include='*01683.pfb*' --exclude='*.*' nt126396@skypilot.geol.umt.edu:projects/PLM_transect.v5/RunA.1/wy_2017_2021/ ./
EcoSLIM from 'ecoslim_2017_2021' dir
Note: needs to be updated EcoSLIM)
wy2017_2021_eco_pnts.00001683.vtk
Files that should be handled with 'rsync_helper.txt'
Permeability and Saturation Fields (from vis_vtk.tcl)
tfg.out.Perm.vtk
wy_2017_2021.out.satur.01683.pfb
From read_vtk_v4.py
ecoslim_rtd.pk
parflow_out/wy_2017_2021_wt_bls.csv
Up-to-date Scripts:
pull_parflow_wl_v5.py - Reads water table simulations
read_vtk_v4.py - Reads Ecoslim output
plot_results_v5.py - Post run plotting and analysis
sparse_vel_files.py - Takes every 10th parflow timestep file and writes to new dir. Using this as way to run Ecoslim and 10 day resolution for spinup.
mangle_clm.py - utilities script to update the 1D clm vegetation file
mangle_vel.py - post analysis script to write velocity vectors to vtk file for Paraview vis
send_pnts.py - utilities script to put select ecoslim pnts files into a new directory
vel_decomp.py - # Script to process the Parflow velocity field at a single timestep.
# Plots velocity field angle and magnitude variations along the hillslope at the saprolite-bedrock interface.
# Plots the Ecoslim particle infiltration locations along the hillslope.
vel_decomp_2_rech.py - # Script to read in and process the Parflow .pfb files for all timesteps
# Saves output to 'parflow_out/pf_out_dict.pk'
# Outputs in pf_out_dict.pk: WTD, specific storage, gw velocities, ET, and pressures
# Notes: Need to manually specify ParFlow time and domain details below