Changeset 15767 for issm/trunk-jpl/test/NightlyRun/test1201.py

Timestamp:

08/09/13 09:05:02 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: Prognostic is now Masstransport

File:

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

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

@@ -18 +18 @@
 
 for stabilization in xrange(1,4):
-        #The goal is to test the prognostic model
+        #The goal is to test the masstransport model
         md=bamg(model(),'domain','../Exp/SquareEISMINT.exp','hmax',3000.)
         md=setmask(md,'all','')
@@ -32 +32 @@
         #Stabilization
         if stabilization==2:
-                md.prognostic.stabilization=0
+                md.masstransport.stabilization=0
         else:
-                md.prognostic.stabilization=stabilization
+                md.masstransport.stabilization=stabilization
 
         #spc thickness
         pos=numpy.nonzero(md.mesh.y>199999.9)[0]
         times=numpy.arange(0,501)
-        md.prognostic.spcthickness=float('NaN')*numpy.ones((md.mesh.numberofvertices+1,numpy.size(times)))
-        md.prognostic.spcthickness[-1,:]=times
-        md.prognostic.spcthickness[pos,:]=numpy.tile(500.+100.*numpy.sin(2.*numpy.pi*times/200.),(numpy.size(pos),1))
+        md.masstransport.spcthickness=float('NaN')*numpy.ones((md.mesh.numberofvertices+1,numpy.size(times)))
+        md.masstransport.spcthickness[-1,:]=times
+        md.masstransport.spcthickness[pos,:]=numpy.tile(500.+100.*numpy.sin(2.*numpy.pi*times/200.),(numpy.size(pos),1))
         if stabilization==3:
-                pos=numpy.nonzero(numpy.isnan(md.prognostic.spcthickness))
-                md.prognostic.spcthickness[pos]=500.    #No NaN for DG
+                pos=numpy.nonzero(numpy.isnan(md.masstransport.spcthickness))
+                md.masstransport.spcthickness[pos]=500.    #No NaN for DG
 
         #solve