source: issm/oecreview/Archive/19101-20495/ISSM-19277-19278.diff@ 20498

../trunk-jpl/test/NightlyRun/test237.py

@@ -53 +53 @@
     md.surfaceforcings.precipitations_lgm[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5
     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
 
 #
 md.transient.requested_outputs=['default','SurfaceforcingsMonthlytemperatures']