trun.py: change hydrology model and shorten default simulation time - pism-exp-gsw - ice stream and sediment transport experiments
(HTM) git clone git://src.adamsgaard.dk/pism-exp-gsw
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
(DIR) LICENSE
---
(DIR) commit 01c7bbb7016088050c404dc4b8db89934fdc6142
(DIR) parent a667852e299cec83941bdecfbda11373b96a3fa9
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 26 Nov 2021 14:51:27 +0100
run.py: change hydrology model and shorten default simulation time
Diffstat:
M run.py | 21 ++++++---------------
1 file changed, 6 insertions(+), 15 deletions(-)
---
(DIR) diff --git a/run.py b/run.py
t@@ -144,19 +144,20 @@ class Experiment:
"-subgl", # PIK: https://pism-docs.org/sphinx/manual/modeling-choices/marine/pik.html#sub-grid-treatment-of-the-grounding-line-position
#"-ssafd_ksp_rtol 1e-7",
"-ys 0",
- "-ye 1e4",
+ "-ye 1e3",
+ #"-ye %d" % MISMIP.run_length(self.experiment, step),
"-options_left",
- #"-hydrology routing",
+ "-hydrology routing",
"-plastic_phi 24", # Tulaczyk et al 2000
"-till_cohesion 0", # Tulaczyk et al 2000
- #"-geothermal_flux 70e-3", # Feldmann and Levermann 2017
+ "-geothermal_flux 70e-3", # Feldmann and Levermann 2017
"-sia_e 4.5", # Martin et al 2010
"-ssa_e 0.512", # Martin et al 2010
"-bed_def lc", # Lingle and Clark 1985, Bueler et al 2007
"-stress_balance.sia.max_diffusivity 1e4",
+ "-sea_level constant,delta_sl -ocean_delta_sl_file sealvl.nc",
]
-
#options = [,
#"-yield_stress mohr_coulomb",
#"-Lz 1e4",
t@@ -169,15 +170,6 @@ class Experiment:
#"-options_left",
#]
-
- #if self.model == 1:
- #options.extend(["-stress_balance ssa"])
- #else:
- #options.extend(["-stress_balance ssa+sia",
- #"-sia_flow_law isothermal_glen", # isothermal setup
- #])
-
-
return options
def config(self, step):
t@@ -240,8 +232,7 @@ class Experiment:
"-extra_vars thk,topg,velbar_mag,flux_mag,mask,dHdt,usurf,hardav,velbase_mag,nuH,tauc,taud,taub,flux_divergence,cell_grounded_fraction",
"-ts_file %s" % ts_file,
"-ts_times 0:50:3e4",
- #"-backup_size big",
- "-backup_size big_2d",
+ "-backup_size big",
"-o %s" % output_file,
"-o_order zyx",
]