Index: /issm/trunk-jpl/src/m/classes/SMBsemic.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBsemic.m	(revision 28024)
+++ /issm/trunk-jpl/src/m/classes/SMBsemic.m	(revision 28025)
@@ -12,7 +12,8 @@
 		dailywindspeed		= NaN;
 		dailypressure		= NaN;
-		dailyairdensity		= NaN;
+		dailyairdensity   = NaN;
 		dailyairhumidity	= NaN;
 		dailytemperature	= NaN;
+
 		Tamp              = NaN;
 		mask              = NaN;
@@ -246,4 +247,5 @@
 				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]');
Index: /issm/trunk-jpl/src/m/classes/SMBsemic.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBsemic.py	(revision 28024)
+++ /issm/trunk-jpl/src/m/classes/SMBsemic.py	(revision 28025)
@@ -25,4 +25,5 @@
         self.dailyairhumidity = np.nan
         self.dailytemperature = np.nan
+
         self.Tamp = np.nan
         self.mask = np.nan
@@ -115,5 +116,5 @@
         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,5 +173,5 @@
         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,4 +296,5 @@
             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')