Changeset 24285

Timestamp:

10/28/19 19:55:36 (5 years ago)

Author:

jdquinn

Message:

BUG: Missing slr.hydro_rate field in Python version

File:

• issm/trunk-jpl/src/m/classes/slr.py (modified) (3 diffs)

issm/trunk-jpl/src/m/classes/slr.py

@@ -68 +68 @@
         string = "%s\n%s" % (string, fielddisplay(self, 'ocean_area_scaling', 'correction for model representation of ocean area [default: No correction]'))
         string = "%s\n%s" % (string, fielddisplay(self, 'steric_rate', 'rate of steric ocean expansion [mm / yr]'))
+        string = "%s\n%s" % (string, fielddisplay(self, 'hydro_rate', 'rate of hydrological expansion [mm / yr]'))
         string = "%s\n%s" % (string, fielddisplay(self, 'Ngia', 'rate of viscous (GIA) geoid expansion (in mm / yr)'))
         string = "%s\n%s" % (string, fielddisplay(self, 'Ugia', 'rate of viscous (GIA) bedrock uplift (in mm / yr)'))
@@ -110 +111 @@
         #steric:
         self.steric_rate = 0
+        self.hydro_rate = 0
         #how many time steps we skip before we run SLR solver during transient
         self.geodetic_run_frequency = 1
@@ -144 +146 @@
         md = checkfield(md, 'fieldname', 'slr.geodetic_run_frequency', 'size', [1, 1], '>=', 1)
         md = checkfield(md, 'fieldname', 'slr.steric_rate', 'NaN', 1, 'Inf', 1, 'size', [md.mesh.numberofvertices])
+        md = checkfield(md, 'fieldname', 'slr.hydro_rate', 'NaN', 1, 'Inf', 1, 'size', [md.mesh.numberofvertices])
         md = checkfield(md, 'fieldname', 'slr.degacc', 'size', [1, 1], '>=', 1e-10)
         md = checkfield(md, 'fieldname', 'slr.requested_outputs', 'stringrow', 1)