Index: /issm/trunk-jpl/test/NightlyRun/test2061.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2061.m	(revision 14973)
+++ /issm/trunk-jpl/test/NightlyRun/test2061.m	(revision 14973)
@@ -0,0 +1,38 @@
+% Benchmark experiments (Figure A2b Ivins and James, 1999, Geophys. J. Int.) 
+md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/GiaBenchmarksAB.par');
+
+%% indicate what you want to compute 
+md.gia.output_rates=2;           % just want "w" solution 
+md.gia.cross_section_shape=1;    % for square-edged x-section 
+
+%% define loading history 
+md.timestepping.start_time=2002100; % after 2 kr of deglaciation 
+md.timestepping.final_time=2500000; %2,500 kyr
+% In order to run Ivins1999 benchmarks, for now, we do things manually as follows: 
+% 1. numtimes is fixed to be five 
+% 2. final_time given above is hard-coded in "GiaDeflectionCorex.cpp"
+% 3. evaluation time = start_time 
+% 4. define new load times, according to the chosen evaluation time (care that general load history is not altered). This is really important, because of #1 and also becasue of the fact that futre loading times are not recognized.
+% 5. Ice thickness assoiated with evaluation time is actually the same at final_time
+md.geometry.thickness=[...
+	[md.geometry.thickness*0.0; 0.0],...
+	[md.geometry.thickness; 1000],...
+	[md.geometry.thickness; 2000000],...
+	[md.geometry.thickness*0.0; 2000100],...
+	[md.geometry.thickness*0.0; 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','GiadWdt'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	(md.results.GiaSolution.GiaW),...
+	(md.results.GiaSolution.GiadWdt),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test2062.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2062.m	(revision 14973)
+++ /issm/trunk-jpl/test/NightlyRun/test2062.m	(revision 14973)
@@ -0,0 +1,38 @@
+% Benchmark experiments (Figure A2b Ivins and James, 1999, Geophys. J. Int.) 
+md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/GiaBenchmarksAB.par');
+
+%% indicate what you want to compute 
+md.gia.output_rates=2;           % just want "w" solution 
+md.gia.cross_section_shape=1;    % for square-edged x-section 
+
+%% define loading history 
+md.timestepping.start_time=2005100; % after 5 kr of deglaciation 
+md.timestepping.final_time=2500000; %2,500 kyr
+% In order to run Ivins1999 benchmarks, for now, we do things manually as follows: 
+% 1. numtimes is fixed to be five 
+% 2. final_time given above is hard-coded in "GiaDeflectionCorex.cpp"
+% 3. evaluation time = start_time 
+% 4. define new load times, according to the chosen evaluation time (care that general load history is not altered). This is really important, because of #1 and also becasue of the fact that futre loading times are not recognized.
+% 5. Ice thickness assoiated with evaluation time is actually the same at final_time
+md.geometry.thickness=[...
+	[md.geometry.thickness*0.0; 0.0],...
+	[md.geometry.thickness; 1000],...
+	[md.geometry.thickness; 2000000],...
+	[md.geometry.thickness*0.0; 2000100],...
+	[md.geometry.thickness*0.0; 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','GiadWdt'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	(md.results.GiaSolution.GiaW),...
+	(md.results.GiaSolution.GiadWdt),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test2063.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2063.m	(revision 14973)
+++ /issm/trunk-jpl/test/NightlyRun/test2063.m	(revision 14973)
@@ -0,0 +1,38 @@
+% Benchmark experiments (Figure A2b Ivins and James, 1999, Geophys. J. Int.) 
+md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/GiaBenchmarksAB.par');
+
+%% indicate what you want to compute 
+md.gia.output_rates=2;           % just want "w" solution 
+md.gia.cross_section_shape=1;    % for square-edged x-section 
+
+%% define loading history 
+md.timestepping.start_time=2010100; % after 10 kyr of deglaciation 
+md.timestepping.final_time=2500000; %2,500 kyr
+% In order to run Ivins1999 benchmarks, for now, we do things manually as follows: 
+% 1. numtimes is fixed to be five 
+% 2. final_time given above is hard-coded in "GiaDeflectionCorex.cpp"
+% 3. evaluation time = start_time 
+% 4. define new load times, according to the chosen evaluation time (care that general load history is not altered). This is really important, because of #1 and also becasue of the fact that futre loading times are not recognized.
+% 5. Ice thickness assoiated with evaluation time is actually the same at final_time
+md.geometry.thickness=[...
+	[md.geometry.thickness*0.0; 0.0],...
+	[md.geometry.thickness; 1000],...
+	[md.geometry.thickness; 2000000],...
+	[md.geometry.thickness*0.0; 2000100],...
+	[md.geometry.thickness*0.0; 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','GiadWdt'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	(md.results.GiaSolution.GiaW),...
+	(md.results.GiaSolution.GiadWdt),...
+	};