Index: ../trunk-jpl/test/NightlyRun/test236.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test236.py	(revision 18970)
+++ ../trunk-jpl/test/NightlyRun/test236.py	(revision 18971)
@@ -47,6 +47,7 @@
 
 # creating precipitation
 md.surfaceforcings.precipitations_presentday=numpy.zeros((md.mesh.numberofvertices+1,12))
+md.surfaceforcings.precipitations_lgm=numpy.zeros((md.mesh.numberofvertices+1,12))
 for imonth in xrange(0,12):
     md.surfaceforcings.precipitations_presentday[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
     md.surfaceforcings.precipitations_presentday[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
@@ -54,11 +55,11 @@
     md.surfaceforcings.precipitations_lgm[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.)
 
 # Interpolation factors
-md.surfaceforcings.Tdiff[1,1:md.timestepping.final_time]=0.5;
-md.surfaceforcings.sealev[1,1:md.timestepping.final_time]=0.5;
+md.surfaceforcings.Tdiff=0.5*numpy.ones((2,md.timestepping.final_time))
+md.surfaceforcings.sealev=0.5*numpy.ones((2,md.timestepping.final_time))
 # Year of each data point
-md.surfaceforcings.Tdiff[2,1:md.timestepping.final_time]=1:1:md.timestepping.final_time;
-md.surfaceforcings.sealev[2,1:md.timestepping.final_time]=1:1:md.timestepping.final_time;
+md.surfaceforcings.Tdiff[1,1:md.timestepping.final_time]=numpy.arange(1.,md.timestepping.final_time)
+md.surfaceforcings.sealev[1,1:md.timestepping.final_time]=numpy.arange(1.,md.timestepping.final_time)
 
 # time steps and resolution
 md.timestepping.time_step=20.
