Changeset 27773 for issm/trunk-jpl/src/m/classes/SMBsemic.py

Timestamp:

05/21/23 20:03:36 (23 months ago)

Author:

inwoo

Message:

CHG: Update SEMIC to enable the output request of SmbRefreeze value. Also, synchronize the SMBsemic class for MATLAB and PYTHON.

File:

• issm/trunk-jpl/src/m/classes/SMBsemic.py (modified) (4 diffs)

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

@@ -113 +113 @@
             s += '{}\n'.format(fielddisplay(self,'albl','background albedo for bare land (default: 0.07)'))
         # albedo_scheme - 0: none, 1: slater, 2: isba, 3: denby, 4: alex.
-        if self.albedo_scheme == 1:
-            s += '\n\tSEMIC snow albedo parameters for Slater et al, (1998).\n'
+        if self.albedo_scheme == 0:
+            s += '\n\tSEMIC snow albedo parameter of None.\n'
+                                s += '\t   albedo of snow is updated from albedo snow max (alb_smax).\n'
+            s += '\t   alb_snow = abl_smax \n '
+        elif self.albedo_scheme == 1:
+            s += '\n\tSEMIC snow albedo parameters of Slater et al, (1998).\n'
             s += '\t   alb = alb_smax - (alb_smax - alb_smin)*tm^(3.0)\n'
             s += '\t   tm  = 1 (tsurf > 273.15 K)\n'
@@ -123 +127 @@
             s += '{}\n'.format(fielddisplay(self, 'tmax', 'maxmium temperature for which albedo decline become effective. This value should be fixed. (default: 273.15 K)[unit: K])'))
         elif self.albedo_scheme == 2:
-            s += '\n\tSEMIC snow albedo parameters for ISBA.? where is citation?\n'
+            s += '\n\tSEMIC snow albedo parameters of ISBA.? where is citation?\n'
             s += '{}\n'.format(fielddisplay(self, 'mcrit', 'critical melt rate (default: 6e-8) [unit: m/sec]'))
             s += '{}\n'.format(fielddisplay(self, 'wcrit', 'critical liquid water content (default: 15) [unit: kg/m2]'))
@@ -129 +133 @@
             s += '{}\n'.format(fielddisplay(self, 'tau_f', 'wet albedo decline [unit: 1/day]'))
         elif self.albedo_scheme == 3:
-            s += '\n\tSEMIC snow albedo parameters for Denby et al. (2002 Tellus)\n'
+            s += '\n\tSEMIC snow albedo parameters of Denby et al. (2002 Tellus)\n'
             s += '{}\n'.format(fielddisplay(self,'mcrit','critical melt rate (defaut: 6e-8) [unit: m/sec]'))
         elif self.albedo_scheme == 4:
-            s += '\n\tSEMIC snow albedo parameters for Alex.?\n'
+            s += '\n\tSEMIC snow albedo parameters of Alex.?\n'
             s += '{}\n'.format(fielddisplay(self,'afac','[unit: ?]'))
             s += '{}\n'.format(fielddisplay(self,'tmid','[unit: ?]'))
@@ -167 +171 @@
     def outputlists(self, md):  # {{{
         if self.ismethod:
-            list = ['SmbMassBalance', 'SmbMassBalanceSnow', 'SmbMelt', 'SmbAccumulation', 'SmbHIce', 'SmbHSnow', 'SmbAlbedo', 'SmbAlbedoSnow', 'TemperatureSEMIC']
+            list = ['SmbMassBalance', 'SmbMassBalanceSnow', 'SmbMassBalanceIce',
+                  'SmbMassBalanceSecmi', 'SmbMelt', 'SmbRefreeze','SmbAccumulation',
+                  'SmbHIce', 'SmbHSnow', 'SmbAlbedo', 'SmbAlbedoSnow', 'TemperatureSEMIC',
+                  'SmbSemicQmr']
         else:
             list = ['SmbMassBalance']