Changeset 18971

Timestamp:

12/30/14 10:45:42 (10 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed some matlab to python conversion

File:

• issm/trunk-jpl/test/NightlyRun/test236.py (modified) (2 diffs)

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

@@ -48 +48 @@
 # 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
@@ -55 +56 @@
 
 # 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