Changeset 17149

Timestamp:

01/21/14 15:47:01 (11 years ago)

Author:

jbondzio

Message:

CHG: LSM nightlyrun uses SquareSheetConstrained par file now

File:

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

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

@@ -10 +10 @@
 md=triangle(model(),'../Exp/Square.exp',150000.)
 md=setmask(md,'','')
-md=parameterize(md,'../Par/SquareSheetLevelset.py') 
+md=parameterize(md,'../Par/SquareSheetConstrained.py') 
 md=setflowequation(md,'SSA','all')
 md.cluster=generic('name',oshostname(),'np',3)
+
+# parameters
+md.transient.isstressbalance=True
+md.transient.ismasstransport=False
+md.transient.isthermal=False
+md.transient.isgroundingline=False
+md.transient.isgia=False
+md.transient.islevelset=True
+
+# init levelset function
+ymin=numpy.min(md.mesh.y)
+ymax=numpy.max(md.mesh.y)
+xmin=min(md.mesh.x)
+xmax=max(md.mesh.x)
+xmed=(xmax+xmin)/2
+ymed=(ymax+ymin)/2
+distx=numpy.absolute(md.mesh.x.reshape(-1,1)-xmed)
+disty=numpy.absolute(md.mesh.y.reshape(-1,1)-ymed)
+md.mask.ice_levelset=numpy.maximum(distx,disty)-1.e5
+
 md=solve(md,TransientSolutionEnum())