Index: ../trunk-jpl/src/m/classes/SMBpdd.py
===================================================================
--- ../trunk-jpl/src/m/classes/SMBpdd.py	(revision 18727)
+++ ../trunk-jpl/src/m/classes/SMBpdd.py	(revision 18728)
@@ -18,6 +18,12 @@
 		self.monthlytemperatures       = float('NaN')
 		self.desfac                    = 0.
 		self.s0p                       = 0.
+		self.s0t                       = 0.
+		self.rlaps                     = 0.
+		self.rlapslgm                  = 0.
+		self.Pfac                      = float('NaN')
+		self.Tdiff                     = float('NaN')
+		self.sealev                    = float('NaN')
 		self.isdelta18o                = 0
 		self.delta18o                  = float('NaN')
 		self.delta18o_surface          = float('NaN')
@@ -35,6 +41,12 @@
 		string="%s\n%s"%(string,fielddisplay(self,'desfac','desertification elevation factor (between 0 and 1, default is 0.5) [m]'))
 		string="%s\n%s"%(string,fielddisplay(self,'isdelta18o','is temperature and precipitation delta18o parametrisation activated (0 or 1, default is 0)'))
 		string="%s\n%s"%(string,fielddisplay(self,'s0p','should be set to elevation from precip source (between 0 and a few 1000s m, default is 0) [m]'))
+		string="%s\n%s"%(string,fielddisplay(self,'s0t','should be set to elevation from temperature source (between 0 and a few 1000s m, default is 0) [m]'))
+		string="%s\n%s"%(string,fielddisplay(self,'rlaps','present day lapse rate [degree/km]'))
+		string="%s\n%s"%(string,fielddisplay(self,'rlapslgm','LGM lapse rate [degree/km]'))
+		string="%s\n%s"%(string,fielddisplay(self,'Pfac','time interpolation parameter for precipitation, 1D (year)'))
+		string="%s\n%s"%(string,fielddisplay(self,'Tdiff','time interpolation parameter for temperature, 1D (year)'))
+		string="%s\n%s"%(string,fielddisplay(self,'sealev','sea level [m], 1D(year)'))
 		string="%s\n%s"%(string,fielddisplay(self,'monthlytemperatures','monthly surface temperatures [K], required if pdd is activated and delta18o not activated'))
 		string="%s\n%s"%(string,fielddisplay(self,'temperatures_presentday','monthly present day surface temperatures [K], required if pdd is activated and delta18o activated'))
 		string="%s\n%s"%(string,fielddisplay(self,'temperatures_lgm','monthly LGM surface temperatures [K], required if pdd is activated and delta18o activated'))
@@ -67,6 +79,9 @@
 		self.isdelta18o = 0
 		self.desfac     = 0.5
 		self.s0p        = 0.
+		self.s0t        = 0.
+		self.rlaps      = 6.5 
+		self.rlapslgm   = 6.5
 
 		return self
 	#}}}
@@ -75,6 +90,12 @@
 		if MasstransportAnalysisEnum() in analyses:
 			md = checkfield(md,'fieldname','surfaceforcings.desfac','<=',1,'numel',[1]);
 			md = checkfield(md,'fieldname','surfaceforcings.s0p','>=',0,'numel',[1]);
+			md = checkfield(md,'fieldname','surfaceforcings.s0t','>=',0,'numel',[1]);
+			md = checkfield(md,'fieldname','surfaceforcings.rlaps','>=',0,'numel',[1]);
+			md = checkfield(md,'fieldname','surfaceforcings.rlapslgm','>=',0,'numel',[1]);
+			md = checkfield(md,'fieldname','surfaceforcings.Pfac','NaN',1,'size',[2,NaN]);
+			md = checkfield(md,'fieldname','surfaceforcings.Tdiff','NaN',1,'size',[2,NaN]);
+			md = checkfield(md,'fieldname','surfaceforcings.sealev','NaN',1,'size',[2,NaN]);
 			if not self.isdelta18o:
 				md = checkfield(md,'fieldname','surfaceforcings.monthlytemperatures','forcing',1,'NaN',1)
 				md = checkfield(md,'fieldname','surfaceforcings.precipitation','forcing',1,'NaN',1)
@@ -97,6 +118,12 @@
 		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','isdelta18o','format','Boolean')
 		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','desfac','format','Double');
 		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','s0p','format','Double');
+		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','s0t','format','Double');
+		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','rlaps','format','Double');
+		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','rlapslgm','format','Double');
+		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','Pfac','format','DoubleMat','mattype',1);
+		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','Tdiff','format','DoubleMat','mattype',1);
+		WriteData(fid,'object',self,'class','surfaceforcings','fieldname','sealev','format','DoubleMat','mattype',1);
 		if self.isdelta18o:
 			WriteData(fid,'object',self,'class','surfaceforcings','fieldname','temperatures_presentday','format','DoubleMat','mattype',1)
 			WriteData(fid,'object',self,'class','surfaceforcings','fieldname','temperatures_lgm','format','DoubleMat','mattype',1)
Index: ../trunk-jpl/src/m/classes/SMBpdd.m
===================================================================
--- ../trunk-jpl/src/m/classes/SMBpdd.m	(revision 18727)
+++ ../trunk-jpl/src/m/classes/SMBpdd.m	(revision 18728)
@@ -95,7 +95,7 @@
 			fielddisplay(obj,'Pfac','time interpolation parameter for precipitation, 1D(year)');
 			fielddisplay(obj,'Tdiff','time interpolation parameter for temperature, 1D(year)');
 			fielddisplay(obj,'sealev','sea level [m], 1D(year)');
-                        fielddisplay(obj,'monthlytemperatures','monthly surface temperatures [K], required if pdd is activated and delta18o not activated');
+			fielddisplay(obj,'monthlytemperatures','monthly surface temperatures [K], required if pdd is activated and delta18o not activated');
 			fielddisplay(obj,'precipitation','surface precipitation [m/yr water eq]');
 			fielddisplay(obj,'temperatures_presentday','monthly present day surface temperatures [K], required if pdd is activated and delta18o activated');
 			fielddisplay(obj,'temperatures_lgm','monthly LGM surface temperatures [K], required if pdd is activated and delta18o activated');
@@ -116,9 +116,9 @@
 			WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','s0t','format','Double');
 			WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','rlaps','format','Double');
 			WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','rlapslgm','format','Double');
-                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','Pfac','format','DoubleMat','mattype',1);
-                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','Tdiff','format','DoubleMat','mattype',1);
-                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','sealev','format','DoubleMat','mattype',1);
+			WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','Pfac','format','DoubleMat','mattype',1);
+			WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','Tdiff','format','DoubleMat','mattype',1);
+			WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','sealev','format','DoubleMat','mattype',1);
 			WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','isdelta18o','format','Boolean');
 			if obj.isdelta18o
 				WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','temperatures_presentday','format','DoubleMat','mattype',1);
