Changeset 28025

Timestamp:

12/01/23 23:36:05 (16 months ago)

Author:

inwoo

Message:

CHG: fix bug in classes/SMBsemic.py and synchronize classes/SMBsemic for python and matlab.

Location:

issm/trunk-jpl/src/m/classes

Files:

• SMBsemic.m (modified) (2 diffs)
• SMBsemic.py (modified) (4 diffs)

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

@@ -12 +12 @@
                 dailywindspeed          = NaN;
                 dailypressure           = NaN;
-                dailyairdensity         = NaN;
+                dailyairdensity   = NaN;
                 dailyairhumidity        = NaN;
                 dailytemperature        = NaN;
+
                 Tamp              = NaN;
                 mask              = NaN;
@@ -246 +247 @@
                                 fielddisplay(self,'hsnow','initial thickness of snow [unit: m]');
                                 fielddisplay(self,'mask','masking for albedo. 0: ocean, 1: land, 2: ice (default: 2)');
+                                fielddisplay(self,'qmr','initial net energy difference between melt and refreeze in SEMIC [unit: W m^{-2}]. This variable can be set with zeros because net energy difference between melt and refreeze is dissipated fast.');
                                 fielddisplay(self,'hcrit','critical snow height for albedo [unit: m]');
                                 fielddisplay(self,'rcrit','critical refreezing height for albedo [no unit]');

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

@@ -25 +25 @@
         self.dailyairhumidity = np.nan
         self.dailytemperature = np.nan
+
         self.Tamp = np.nan
         self.mask = np.nan
@@ -115 +116 @@
         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   albedo of snow is updated from albedo snow max (alb_smax).\n'
             s += '\t   alb_snow = abl_smax \n '
         elif self.albedo_scheme == 1:
@@ -172 +173 @@
         if self.ismethod:
             list = ['default','SmbMassBalance', 'SmbMassBalanceSnow', 'SmbMassBalanceIce',
-                  'SmbMassBalanceSecmi', 'SmbMelt', 'SmbRefreeze','SmbAccumulation',
+                  'SmbMassBalanceSemic', 'SmbMelt', 'SmbRefreeze','SmbAccumulation',
                   'SmbHIce', 'SmbHSnow', 'SmbAlbedo', 'SmbAlbedoSnow', 'TemperatureSEMIC',
                   'SmbSemicQmr', 'TotalSmb', 'TotalSmbMelt', 'TotalSmbRefreeze']
@@ -295 +296 @@
             WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'hice', 'format', 'DoubleMat', 'mattype', 1)
             WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'hsnow', 'format', 'DoubleMat', 'mattype', 1)
+            WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'qmr', 'format', 'DoubleMat', 'mattype', 1)
 
             WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'hcrit', 'format', 'Double')