Index: /issm/trunk-jpl/src/m/classes/surfaceforcings.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/surfaceforcings.py	(revision 13685)
+++ /issm/trunk-jpl/src/m/classes/surfaceforcings.py	(revision 13686)
@@ -22,4 +22,6 @@
 		self.isdelta18o = 0
 		self.hc = float('NaN')
+		self.href = float('NaN')
+		self.smbref = float('NaN')
 		self.smb_pos_max = float('NaN')
 		self.smb_pos_min = float('NaN')
@@ -56,4 +58,6 @@
 		string="%s\n%s"%(string,fielddisplay(self,'issmbgradients','is smb gradients method activated (0 or 1, default is 0)'))
 		string="%s\n%s"%(string,fielddisplay(self,'hc',' elevation of intersection between accumulation and ablation regime required if smb gradients is activated'))
+		string="%s\n%s"%(string,fielddisplay(self,'href',' reference elevation from which deviation is used to calculate SMB adjustment in smb gradients method'))
+		string="%s\n%s"%(string,fielddisplay(self,'smbref',' reference smb from which deviation is calculated in smb gradients method'))
 		string="%s\n%s"%(string,fielddisplay(self,'smb_pos_max',' maximum value of positive smb required if smb gradients is activated'))
 		string="%s\n%s"%(string,fielddisplay(self,'smb_pos_min',' minimum value of positive smb required if smb gradients is activated'))
@@ -94,4 +98,6 @@
 			elif self.issmbgradients:
 				md = checkfield(md,'surfaceforcings.hc','forcing',1,'NaN',1)
+				md = checkfield(md,'surfaceforcings.href','forcing',1,'NaN',1)
+				md = checkfield(md,'surfaceforcings.smbref','forcing',1,'NaN',1)
 				md = checkfield(md,'surfaceforcings.smb_pos_max','forcing',1,'NaN',1)
 				md = checkfield(md,'surfaceforcings.smb_pos_min','forcing',1,'NaN',1)
@@ -130,4 +136,6 @@
 		if self.issmbgradients:
 			WriteData(fid,'object',self,'fieldname','hc','format','DoubleMat','mattype',1)
+			WriteData(fid,'object',self,'fieldname','href','format','DoubleMat','mattype',1)
+			WriteData(fid,'object',self,'fieldname','smbref','format','DoubleMat','mattype',1)
 			WriteData(fid,'object',self,'fieldname','smb_pos_max','format','DoubleMat','mattype',1)
 			WriteData(fid,'object',self,'fieldname','smb_pos_min','format','DoubleMat','mattype',1)