Index: /issm/trunk-jpl/test/NightlyRun/test236.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test236.m	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test236.m	(revision 19308)
@@ -41,4 +41,5 @@
 
 md.initialization.temperature=md.surfaceforcings.temperatures_lgm(1:md.mesh.numberofvertices,1); %*ones(md.mesh.numberofvertices,1);
+md.surfaceforcings = initialize(md.surfaceforcings,md);
 
 % creating precipitation
Index: /issm/trunk-jpl/test/NightlyRun/test236.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test236.py	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test236.py	(revision 19308)
@@ -46,4 +46,5 @@
 
 md.initialization.temperature=md.surfaceforcings.temperatures_lgm[0:md.mesh.numberofvertices,0]    #*ones(md.mesh.numberofvertices,1)
+md.surfaceforcings.initialize(md)
 
 # creating precipitation
Index: /issm/trunk-jpl/test/NightlyRun/test237.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test237.m	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test237.m	(revision 19308)
@@ -34,4 +34,5 @@
 
 md.initialization.temperature=md.surfaceforcings.temperatures_lgm(1:md.mesh.numberofvertices,1); %*ones(md.mesh.numberofvertices,1);
+md.surfaceforcings = initialize(md.surfaceforcings,md);
 
 % creating precipitation
Index: /issm/trunk-jpl/test/NightlyRun/test237.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 19308)
@@ -44,4 +44,5 @@
 
 md.initialization.temperature=md.surfaceforcings.temperatures_lgm[0:md.mesh.numberofvertices,0].reshape(-1,1)   #*ones(md.mesh.numberofvertices,1)
+md.surfaceforcings.initialize(md)
 
 # creating precipitation
Index: /issm/trunk-jpl/test/NightlyRun/test238.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test238.m	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test238.m	(revision 19308)
@@ -37,4 +37,5 @@
     md.surfaceforcings.precipitations_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
 end
+md.surfaceforcings = initialize(md.surfaceforcings,md);
 
 % time steps and resolution
Index: /issm/trunk-jpl/test/NightlyRun/test238.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test238.py	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test238.py	(revision 19308)
@@ -36,4 +36,5 @@
 md.thermal.spctemperature=md.thermal.spctemperature-10
 md.initialization.temperature=md.thermal.spctemperature
+md.surfaceforcings.initialize(md)
 
 # creating precipitation
Index: /issm/trunk-jpl/test/NightlyRun/test239.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test239.m	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test239.m	(revision 19308)
@@ -38,4 +38,5 @@
     md.surfaceforcings.precipitations_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
 end
+md.surfaceforcings = initialize(md.surfaceforcings,md);
 
 % time steps and resolution
Index: /issm/trunk-jpl/test/NightlyRun/test239.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test239.py	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test239.py	(revision 19308)
@@ -36,4 +36,5 @@
 md.thermal.spctemperature=md.thermal.spctemperature-10
 md.initialization.temperature=md.thermal.spctemperature
+md.surfaceforcings.initialize(md)
 
 # creating precipitation
Index: /issm/trunk-jpl/test/NightlyRun/test240.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test240.m	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test240.m	(revision 19308)
@@ -25,9 +25,9 @@
     md.surfaceforcings.temperatures_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
 end
+md.surfaceforcings = initialize(md.surfaceforcings,md);
 
 % creating initialization and spc temperatures initialization and
 % spc
 md.thermal.spctemperature=mean(md.surfaceforcings.temperatures_presentday(1:md.mesh.numberofvertices,1:12),2)-10; %-10*ones(md.mesh.numberofvertices,1);
-
 md.initialization.temperature=md.thermal.spctemperature; %md.surfaceforcings.temperatures_presentday(1:md.mesh.numberofvertices,1);
 
Index: /issm/trunk-jpl/test/NightlyRun/test240.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test240.py	(revision 19307)
+++ /issm/trunk-jpl/test/NightlyRun/test240.py	(revision 19308)
@@ -36,4 +36,5 @@
 md.thermal.spctemperature=md.thermal.spctemperature-10
 md.initialization.temperature=md.thermal.spctemperature
+md.surfaceforcings.initialize(md)
 
 # creating precipitation