Index: /issm/trunk-jpl/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 14825)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 14826)
@@ -50,10 +50,10 @@
 
 %GIA: 
-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
+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=1.*10^21;                             % in Pa.s
+md.materials.mantle_density=3.34;                                   % in g/cm^-3
 
 %Boundary conditions:
Index: /issm/trunk-jpl/test/Par/SquareSheetConstrained.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.py	(revision 14825)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.py	(revision 14826)
@@ -59,4 +59,12 @@
 md.timestepping.final_time=3.
 
+#GIA:
+md.gia.lithosphere_thickness=100.*numpy.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=1.*10**21;                                    # in Pa.s
+md.materials.mantle_density=3.34;                                           # in g/cm^-3
+
 #Boundary conditions:
 md=SetIceSheetBC(md)