Index: /issm/trunk-jpl/test/NightlyRun/test237.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 19277)
+++ /issm/trunk-jpl/test/NightlyRun/test237.py	(revision 19278)
@@ -54,13 +54,17 @@
     md.surfaceforcings.precipitations_lgm[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
 
-for iint in xrange(0,(md.timestepping.final_time/md.timestepping.time_step)+2)
+fsize=int(md.timestepping.final_time/md.timestepping.time_step)+2
+md.surfaceforcings.Pfac=numpy.zeros((2,fsize))
+md.surfaceforcings.Tdiff=numpy.zeros((2,fsize))
+md.surfaceforcings.sealev=numpy.zeros((2,fsize))
+for iint in xrange(0,fsize):
     # Interpolation factors
-	 md.surfaceforcings.Pfac[1,iint]=0.15*(iint+1);
-	 md.surfaceforcings.Tdiff[1,iint]=0.15*(iint+1);
-	 md.surfaceforcings.sealev[1,iint]=0.15*(iint+1);
+	 md.surfaceforcings.Pfac[0,iint]=0.15*(iint+1)
+	 md.surfaceforcings.Tdiff[0,iint]=0.15*(iint+1)
+	 md.surfaceforcings.sealev[0,iint]=0.15*(iint+1)
 	 # Year of each data point
-	 md.surfaceforcings.Pfac[2,iint]=(iint)*20;
-	 md.surfaceforcings.Tdiff[2,iint]=(iint)*20;
-	 md.surfaceforcings.sealev[2,iint]=(iint)*20;
+	 md.surfaceforcings.Pfac[1,iint]=(float(iint))*20
+	 md.surfaceforcings.Tdiff[1,iint]=(float(iint))*20
+	 md.surfaceforcings.sealev[1,iint]=(float(iint))*20
 
 #