Index: /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par
===================================================================
--- /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par	(revision 14901)
+++ /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par	(revision 14902)
@@ -18,9 +18,9 @@
 
 %GIA parameters specific to Experiments A  and B
+md.gia.mantle_viscosity=10^21*ones(md.mesh.numberofvertices,1);    %in Pa.s
 md.gia.lithosphere_thickness=100*ones(md.mesh.numberofvertices,1); %in km
 md.materials.lithosphere_shear_modulus=6.7*10^10;                  %in Pa
 md.materials.lithosphere_density=3.36;                             %in g/cm^3
 md.materials.mantle_shear_modulus=1.45*10^11;                      %in Pa
-md.materials.mantle_viscosity=10^21;                               %in Pa.s
 md.materials.mantle_density=3.38;                                  %in g/cm^3
 
Index: /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par
===================================================================
--- /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par	(revision 14901)
+++ /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par	(revision 14902)
@@ -18,9 +18,9 @@
 
 %GIA parameters specific to Experiments A  and B
+md.gia.mantle_viscosity=10^21*ones(md.mesh.numberofvertices,1);    %in Pa.s
 md.gia.lithosphere_thickness=100*ones(md.mesh.numberofvertices,1); %in km
 md.materials.lithosphere_shear_modulus=6.7*10^10;                  %in Pa
 md.materials.lithosphere_density=3.32;                             %in g/cm^3
 md.materials.mantle_shear_modulus=1.45*10^11;                      %in Pa
-md.materials.mantle_viscosity=10^21;                               %in Pa.s
 md.materials.mantle_density=3.34;                                  %in g/cm^3
 
Index: /issm/trunk-jpl/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 14901)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 14902)
@@ -51,8 +51,8 @@
 %GIA: 
 md.gia.lithosphere_thickness=100.*ones(md.mesh.numberofvertices,1); % in km
+md.gia.mantle_viscosity=1.0*10^21*ones(md.mesh.numberofvertices,1); % in Pa.s
 md.materials.lithosphere_shear_modulus=6.7*10^10;                   % in Pa
 md.materials.lithosphere_density=3.32;                              % in g/cm^-3
 md.materials.mantle_shear_modulus=1.45*10^11;                       % in Pa
-md.materials.mantle_viscosity=1.*10^21;                             % in Pa.s
 md.materials.mantle_density=3.34;                                   % in g/cm^-3
 
Index: /issm/trunk-jpl/test/Par/SquareSheetConstrained.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.py	(revision 14901)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.py	(revision 14902)
@@ -61,8 +61,8 @@
 #GIA:
 md.gia.lithosphere_thickness=100.*numpy.ones((md.mesh.numberofvertices,1)); # in km
+md.gia.mantle_viscosity=1.*10**21*numpy.ones((md.mesh.numberofvertices,1)); # in Pa.s
 md.materials.lithosphere_shear_modulus=6.7*10**10;                          # in Pa
 md.materials.lithosphere_density=3.32;                                      # in g/cm^-3
 md.materials.mantle_shear_modulus=1.45*10**11;                              # in Pa
-md.materials.mantle_viscosity=1.*10**21;                                    # in Pa.s
 md.materials.mantle_density=3.34;                                           # in g/cm^-3