Index: /issm/trunk-jpl/src/m/classes/SMBd18opdd.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBd18opdd.m	(revision 21709)
+++ /issm/trunk-jpl/src/m/classes/SMBd18opdd.m	(revision 21710)
@@ -12,4 +12,5 @@
 		rlapslgm                  = 0; 
 		dpermil                   = 0; 
+		f                         = 0;
 		Tdiff                     = NaN;
 		sealev                    = NaN;
@@ -63,4 +64,5 @@
 		  self.rlapslgm   = 6.5;
 		  self.dpermil    = 2.4;
+		  self.f          = 0.169;
                   
 		end % }}}
@@ -78,4 +80,5 @@
 					md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
 					md = checkfield(md,'fieldname','smb.dpermil','>=',0,'numel',1);
+				   md = checkfield(md,'fieldname','smb.f','>=',0,'numel',1);
 				end
 			end
@@ -96,4 +99,5 @@
 				fielddisplay(self,'delta18o','delta18o [per mil], required if pdd is activated and d18opd activated');  
 				fielddisplay(self,'dpermil','degree per mil, required if d18opd is activated');                            
+			   fielddisplay(self,'f','precip/temperature scaling factor, required if d18opd is activated');
 			end
 			fielddisplay(self,'requested_outputs','additional outputs requested');
@@ -123,4 +127,5 @@
 				WriteData(fid,prefix,'object',self,'class','smb','fieldname','delta18o','format','DoubleMat','mattype',1,'timeserieslength',2,'yts',md.constants.yts);
 				WriteData(fid,prefix,'object',self,'class','smb','fieldname','dpermil','format','Double');
+			   WriteData(fid,prefix,'object',self,'class','smb','fieldname','f','format','Double');
 			end
 			
Index: /issm/trunk-jpl/src/m/classes/SMBd18opdd.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBd18opdd.py	(revision 21709)
+++ /issm/trunk-jpl/src/m/classes/SMBd18opdd.py	(revision 21710)
@@ -20,4 +20,5 @@
 		self.rlapslgm                  = 0.
 		self.dpermil                   = 0.
+		self.f                         = 0.
 		self.Tdiff                     = float('NaN')
 		self.sealev                    = float('NaN')
@@ -83,5 +84,5 @@
 		self.rlapslgm   = 6.5
 		self.dpermil    = 2.4
-
+      self.f          = 0.169
 		return self
 	#}}}
@@ -100,5 +101,5 @@
 				md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'Inf',1,'size',[2,np.nan],'singletimeseries',1)
 				md = checkfield(md,'fieldname','smb.dpermil','>=',0,'numel',[1])
-		
+		      md = checkfield(md,'fieldname','smb.f','>=',0,'numel',[1])
 		md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)
 
@@ -126,5 +127,5 @@
 			WriteData(fid,prefix,'object',self,'class','smb','fieldname','delta18o','format','DoubleMat','mattype',1,'timeserieslength',2,'yts',md.constants.yts)
 			WriteData(fid,prefix,'object',self,'class','smb','fieldname','dpermil','format','Double')
-			
+			WriteData(fid,prefix,'object',self,'class','smb','fieldname','f','format','Double')
 		#process requested outputs
 		outputs = self.requested_outputs