Changeset 14924

Timestamp:

05/06/13 12:35:10 (12 years ago)

Author:

adhikari

Message:

CHG: updated with GiadWdt field_value

File:

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

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

@@ -10 +10 @@
 from MatlabFuncs import *
 
-md=triangle(model(),'../Exp/Square.exp',50000.)
+#Define a model 
+md=model()
+md=triangle(md,'../Exp/Square.exp',100000.)
 md=setmask(md,'','')
-md=parameterize(md,'../Par/SquareSheetConstrained.py')
-md=setflowequation(md,'macayeal','all')
+md=parameterize(md,'../Par/SquareSheetConstrained.par')
+
+#Indicate what you want to compute 
+md.gia.output_rates=1           # just want "w" solution 
+md.gia.cross_section_shape=1    # for square-edged x-section 
+
+#Define loading history (see test2001.m for the description)
+md.timestepping.start_time=2400000 # 2,400 kyr
+md.timestepping.final_time=2500000 # 2,500 kyr
+md.geometry.thickness=[\
+        [md.geometry.thickness*0.0; 0.0],\
+        [md.geometry.thickness/2.0; 0.1],\
+        [md.geometry.thickness; 0.2],\
+        [md.geometry.thickness; 1.0],\
+        [md.geometry.thickness; md.timestepping.start_time],\
+        ]
+
+#Solve for GIA deflection 
 md.cluster=generic('name',oshostname(),'np',3)
+md.verbose=verbose('1111111')
 md=solve(md,GiaSolutionEnum())
 
 #Fields and tolerances to track changes
-field_names     =['GiaW']
-field_tolerances=[1e-13]
-field_values    =[md.results.GiaSolution.GiaW]
+field_names     =['GiaW','GiadWdt']
+field_tolerances=[1e-13,1e-13]
+field_values    =[\
+                md.results.GiaSolution.GiaW,\
+                md.results.GiaSolution.GiadWdt,\
+                ]
+