Changeset 19278

Timestamp:

04/09/15 14:26:04 (10 years ago)

Author:

schlegel

Message:

BUG: need to use python indexing

File:

• issm/trunk-jpl/test/NightlyRun/test237.py (modified) (1 diff)

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

@@ -54 +54 @@
     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
 
 #