Index: sm/trunk-jpl/test/NightlyRun/test2071.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2071.m	(revision 25290)
+++ 	(revision )
@@ -1,31 +1,0 @@
-%Test Name: GiaIvinsBenchmarksAB2dC1
-% Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.) 
-md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
-md=setmask(md,'','');
-md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
-
-%% indicate what you want to compute 
-md.gia.cross_section_shape=1;    % for square-edged x-section 
-
-%% define loading history 
-md.timestepping.start_time=0.3;     % for t \approx 0 kyr : to get eleastic response!   
-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; md.timestepping.start_time],...
-	];
-
-%% solve for GIA deflection 
-md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
-md=solve(md,'Gia');
-
-%Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-	(md.results.GiaSolution.UGia),...
-	(md.results.GiaSolution.UGiaRate),...
-	};
Index: sm/trunk-jpl/test/NightlyRun/test2071.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2071.py	(revision 25290)
+++ 	(revision )
@@ -1,37 +1,0 @@
-#Test Name: GiaIvinsBenchmarksAB2dC1
-from socket import gethostname
-
-import numpy as np
-
-from model import *
-from parameterize import *
-from setmask import *
-from solve import *
-from triangle import *
-
-
-#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
-md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
-md = setmask(md, '', '')
-md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
-
-#indicate what you want to compute
-md.gia.cross_section_shape = 1  # for square-edged x-section
-
-#define loading history
-md.timestepping.start_time = 0.3  # for t \approx 0 kyr : to get eleastic response!
-md.timestepping.final_time = 2500000  # 2,500 kyr
-md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0),
-                                  np.append(md.geometry.thickness / 2.0, 0.1),
-                                  np.append(md.geometry.thickness, 0.2),
-                                  np.append(md.geometry.thickness, md.timestepping.start_time)]).T
-
-#solve for GIA deflection
-md.cluster = generic('name', gethostname(), 'np', 3)
-md.verbose = verbose('1111111')
-md = solve(md, 'Gia')
-
-#Fields and tolerances to track changes
-field_names = ['UGia', 'UGiaRate']
-field_tolerances = [1e-13, 1e-13]
-field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]
Index: sm/trunk-jpl/test/NightlyRun/test2072.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2072.m	(revision 25290)
+++ 	(revision )
@@ -1,31 +1,0 @@
-%Test Name: GiaIvinsBenchmarksAB2dC2
-% Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.) 
-md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
-md=setmask(md,'','');
-md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
-
-%% indicate what you want to compute 
-md.gia.cross_section_shape=1;    % for square-edged x-section 
-
-%% define loading history 
-md.timestepping.start_time=1000.3;  % for t \approx 1 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; md.timestepping.start_time],...
-	];
-
-%% solve for GIA deflection 
-md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
-md=solve(md,'Gia');
-
-%Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-	(md.results.GiaSolution.UGia),...
-	(md.results.GiaSolution.UGiaRate),...
-	};
Index: sm/trunk-jpl/test/NightlyRun/test2072.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2072.py	(revision 25290)
+++ 	(revision )
@@ -1,37 +1,0 @@
-#Test Name: GiaIvinsBenchmarksAB2dC2
-from socket import gethostname
-
-import numpy as np
-
-from model import *
-from parameterize import *
-from setmask import *
-from solve import *
-from triangle import *
-
-
-#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
-md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
-md = setmask(md, '', '')
-md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
-
-#indicate what you want to compute
-md.gia.cross_section_shape = 1  # for square-edged x - section
-
-#define loading history
-md.timestepping.start_time = 1000.3  # for t \approx 1 kyr
-md.timestepping.final_time = 2500000  # 2, 500 kyr
-md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0),
-                                  np.append(md.geometry.thickness / 2.0, 0.1),
-                                  np.append(md.geometry.thickness, 0.2),
-                                  np.append(md.geometry.thickness, md.timestepping.start_time)]).T
-
-#solve for GIA deflection
-md.cluster = generic('name', gethostname(), 'np', 3)
-md.verbose = verbose('1111111')
-md = solve(md, 'Gia')
-
-#Fields and tolerances to track changes
-field_names = ['UGia', 'UGiaRate']
-field_tolerances = [1e-13, 1e-13]
-field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]
Index: sm/trunk-jpl/test/NightlyRun/test2073.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2073.m	(revision 25290)
+++ 	(revision )
@@ -1,31 +1,0 @@
-%Test Name: GiaIvinsBenchmarksAB2dC3
-% Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.) 
-md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
-md=setmask(md,'','');
-md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
-
-%% indicate what you want to compute 
-md.gia.cross_section_shape=1;    % for square-edged x-section 
-
-%% define loading history 
-md.timestepping.start_time=2400000; % for t \approx \infty 
-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; md.timestepping.start_time],...
-	];
-
-%% solve for GIA deflection 
-md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
-md=solve(md,'Gia');
-
-%Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-	(md.results.GiaSolution.UGia),...
-	(md.results.GiaSolution.UGiaRate),...
-	};
Index: sm/trunk-jpl/test/NightlyRun/test2073.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2073.py	(revision 25290)
+++ 	(revision )
@@ -1,37 +1,0 @@
-#Test Name: GiaIvinsBenchmarksAB2dC3
-from socket import gethostname
-
-import numpy as np
-
-from model import *
-from parameterize import *
-from setmask import *
-from solve import *
-from triangle import *
-
-
-#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
-md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
-md = setmask(md, '', '')
-md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
-
-#indicate what you want to compute
-md.gia.cross_section_shape = 1  # for square-edged x - section
-
-#define loading history
-md.timestepping.start_time = 2400000  # for t \approx \infty
-md.timestepping.final_time = 2500000  # 2, 500 kyr
-md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0),
-                                  np.append(md.geometry.thickness / 2.0, 0.1),
-                                  np.append(md.geometry.thickness, 0.2),
-                                  np.append(md.geometry.thickness, md.timestepping.start_time)]).T
-
-#solve for GIA deflection
-md.cluster = generic('name', gethostname(), 'np', 3)
-md.verbose = verbose('1111111')
-md = solve(md, 'Gia')
-
-#Fields and tolerances to track changes
-field_names = ['UGia', 'UGiaRate']
-field_tolerances = [1e-13, 1e-13]
-field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]
Index: sm/trunk-jpl/test/NightlyRun/test2081.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2081.m	(revision 25290)
+++ 	(revision )
@@ -1,31 +1,0 @@
-%Test Name: GiaIvinsBenchmarksAB2dD1
-% Benchmark experiments (Figure A2d Ivins and James, 1999, Geophys. J. Int.) 
-md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
-md=setmask(md,'','');
-md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
-
-%% indicate what you want to compute 
-md.gia.cross_section_shape=2;    % for elliptical edge  
-
-%% define loading history 
-md.timestepping.start_time=0.3;     % for t \approx 0 kyr : to get eleastic response!   
-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; md.timestepping.start_time],...
-	];
-
-%% solve for GIA deflection 
-md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
-md=solve(md,'Gia');
-
-%Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-	(md.results.GiaSolution.UGia),...
-	(md.results.GiaSolution.UGiaRate),...
-	};
Index: sm/trunk-jpl/test/NightlyRun/test2081.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2081.py	(revision 25290)
+++ 	(revision )
@@ -1,37 +1,0 @@
-#Test Name: GiaIvinsBenchmarksAB2dD1
-from socket import gethostname
-
-import numpy as np
-
-from model import *
-from parameterize import *
-from setmask import *
-from solve import *
-from triangle import *
-
-
-#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
-md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
-md = setmask(md, '', '')
-md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
-
-#indicate what you want to compute
-md.gia.cross_section_shape = 2  # for square-edged x - section
-
-#define loading history
-md.timestepping.start_time = 0.3  # for t \approx 0 kyr : to get eleastic response!
-md.timestepping.final_time = 2500000  # 2, 500 kyr
-md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0),
-                                  np.append(md.geometry.thickness / 2.0, 0.1),
-                                  np.append(md.geometry.thickness, 0.2),
-                                  np.append(md.geometry.thickness, md.timestepping.start_time)]).T
-
-#solve for GIA deflection
-md.cluster = generic('name', gethostname(), 'np', 3)
-md.verbose = verbose('1111111')
-md = solve(md, 'Gia')
-
-#Fields and tolerances to track changes
-field_names = ['UGia', 'UGiaRate']
-field_tolerances = [1e-13, 1e-13]
-field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]
Index: sm/trunk-jpl/test/NightlyRun/test2082.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2082.m	(revision 25290)
+++ 	(revision )
@@ -1,32 +1,0 @@
-%Test Name: GiaIvinsBenchmarksAB2dD2
-% Benchmark experiments (Figure A2d Ivins and James, 1999, Geophys. J. Int.) 
-md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
-md=setmask(md,'','');
-md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
-
-%% indicate what you want to compute 
-md.gia.cross_section_shape=2;    % for elliptical edge 
-
-
-%% define loading history 
-md.timestepping.start_time=1000.3;  % for t \approx 1 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; md.timestepping.start_time],...
-	];
-
-%% solve for GIA deflection 
-md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
-md=solve(md,'Gia');
-
-%Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-	(md.results.GiaSolution.UGia),...
-	(md.results.GiaSolution.UGiaRate),...
-	};
Index: sm/trunk-jpl/test/NightlyRun/test2082.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2082.py	(revision 25290)
+++ 	(revision )
@@ -1,37 +1,0 @@
-#Test Name: GiaIvinsBenchmarksAB2dD2
-from socket import gethostname
-
-import numpy as np
-
-from model import *
-from parameterize import *
-from setmask import *
-from solve import *
-from triangle import *
-
-
-#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
-md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
-md = setmask(md, '', '')
-md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
-
-#indicate what you want to compute
-md.gia.cross_section_shape = 2  # for square-edged x - section
-
-#define loading history
-md.timestepping.start_time = 1000.3  # for t \approx 1 kyr
-md.timestepping.final_time = 2500000  # 2, 500 kyr
-md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0),
-                                  np.append(md.geometry.thickness / 2.0, 0.1),
-                                  np.append(md.geometry.thickness, 0.2),
-                                  np.append(md.geometry.thickness, md.timestepping.start_time)]).T
-
-#solve for GIA deflection
-md.cluster = generic('name', gethostname(), 'np', 3)
-md.verbose = verbose('1111111')
-md = solve(md, 'Gia')
-
-#Fields and tolerances to track changes
-field_names = ['UGia', 'UGiaRate']
-field_tolerances = [1e-13, 1e-13]
-field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]
Index: sm/trunk-jpl/test/NightlyRun/test2083.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2083.m	(revision 25290)
+++ 	(revision )
@@ -1,31 +1,0 @@
-%Test Name: GiaIvinsBenchmarksAB2dD3
-% Benchmark experiments (Figure A2d Ivins and James, 1999, Geophys. J. Int.) 
-md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
-md=setmask(md,'','');
-md=parameterize(md,'../Par/GiaIvinsBenchmarksCD.par');
-
-%% indicate what you want to compute 
-md.gia.cross_section_shape=2;    % for elliptical edge 
-
-%% define loading history 
-md.timestepping.start_time=2400000; % for t \approx \infty 
-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; md.timestepping.start_time],...
-	];
-
-%% solve for GIA deflection 
-md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
-md=solve(md,'Gia');
-
-%Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-	(md.results.GiaSolution.UGia),...
-	(md.results.GiaSolution.UGiaRate),...
-	};
Index: sm/trunk-jpl/test/NightlyRun/test2083.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2083.py	(revision 25290)
+++ 	(revision )
@@ -1,37 +1,0 @@
-#Test Name: GiaIvinsBenchmarksAB2dD3
-from socket import gethostname
-
-import numpy as np
-
-from model import *
-from parameterize import *
-from setmask import *
-from solve import *
-from triangle import *
-
-
-#Benchmark experiments (Figure A2c Ivins and James, 1999, Geophys. J. Int.)
-md = triangle(model(), '../Exp/RoundFrontEISMINT.exp', 200000)
-md = setmask(md, '', '')
-md = parameterize(md, '../Par/GiaIvinsBenchmarksCD.py')
-
-#indicate what you want to compute
-md.gia.cross_section_shape = 2  # for square-edged x - section
-
-#define loading history
-md.timestepping.start_time = 2400000  # for t \approx \infty
-md.timestepping.final_time = 2500000  # 2, 500 kyr
-md.geometry.thickness = np.array([np.append(md.geometry.thickness * 0.0, 0.0),
-                                  np.append(md.geometry.thickness / 2.0, 0.1),
-                                  np.append(md.geometry.thickness, 0.2),
-                                  np.append(md.geometry.thickness, md.timestepping.start_time)]).T
-
-#solve for GIA deflection
-md.cluster = generic('name', gethostname(), 'np', 3)
-md.verbose = verbose('1111111')
-md = solve(md, 'Gia')
-
-#Fields and tolerances to track changes
-field_names = ['UGia', 'UGiaRate']
-field_tolerances = [1e-13, 1e-13]
-field_values = [md.results.GiaSolution.UGia, md.results.GiaSolution.UGiaRate]