Index: ../trunk-jpl/test/NightlyRun/test236.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test236.m	(revision 18723)
+++ ../trunk-jpl/test/NightlyRun/test236.m	(revision 18724)
@@ -43,6 +43,10 @@
 md.timestepping.time_step=20;
 md.timestepping.final_time=60;
 
+md.surfaceforcings.Pfac=[1;1];
+md.surfaceforcings.Tdiff=[1;1];
+md.surfaceforcings.sealev=[1;1];
+
 % 
 md.transient.requested_outputs={'default','SurfaceforcingsMonthlytemperatures'};
 md=setflowequation(md,'SSA','all');
Index: ../trunk-jpl/test/NightlyRun/test237.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test237.m	(revision 18723)
+++ ../trunk-jpl/test/NightlyRun/test237.m	(revision 18724)
@@ -40,6 +40,10 @@
     md.surfaceforcings.precipitations_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
 end
 
+md.surfaceforcings.Pfac=[1;1];
+md.surfaceforcings.Tdiff=[1;1];
+md.surfaceforcings.sealev=[1;1];
+
 % time steps and resolution
 md.timestepping.time_step=20;
 md.settings.output_frequency=1;
Index: ../trunk-jpl/src/m/classes/SMBpdd.m
===================================================================
--- ../trunk-jpl/src/m/classes/SMBpdd.m	(revision 18723)
+++ ../trunk-jpl/src/m/classes/SMBpdd.m	(revision 18724)
@@ -7,14 +7,14 @@
 	properties (SetAccess=public) 
 		precipitation             = NaN;
 		monthlytemperatures       = NaN;
-		desfac                    = 0.5;
+		desfac                    = 0;
 		s0p                       = 0;
 		s0t                       = 0;
-                rlaps                     = 6.5;
-                rlapslgm                  = 6.5;                
-                Pfac                      = NaN;
-                Tdiff                     = NaN;
-                sealev                    = NaN;
+      rlaps                     = 0;
+      rlapslgm                  = 0;                
+      Pfac                      = NaN;
+      Tdiff                     = NaN;
+      sealev                    = NaN;
 		isdelta18o                = 0;
 		delta18o                  = NaN;
 		delta18o_surface          = NaN;
@@ -67,9 +67,9 @@
 				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);
-                                md = checkfield(md,'fieldname','surfaceforcings.Tdiff','NaN',1);
-                                md = checkfield(md,'fieldname','surfaceforcings.sealev','NaN',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(obj.isdelta18o==0)
 					md = checkfield(md,'fieldname','surfaceforcings.monthlytemperatures','forcing',1,'NaN',1);
 					md = checkfield(md,'fieldname','surfaceforcings.precipitation','forcing',1,'NaN',1);
