Index: /issm/trunk-jpl/test/NightlyRun/test2001.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2001.py	(revision 14923)
+++ /issm/trunk-jpl/test/NightlyRun/test2001.py	(revision 14924)
@@ -10,13 +10,36 @@
 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,\
+		]
+