Changeset 14826

Timestamp:

05/01/13 09:03:48 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: updated SquareSheetConstrained.py and added GIA stuff

Location:

issm/trunk-jpl/test/Par

Files:

• SquareSheetConstrained.par (modified) (1 diff)
• SquareSheetConstrained.py (modified) (1 diff)

issm/trunk-jpl/test/Par/SquareSheetConstrained.par

@@ -50 +50 @@
 
 %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:

issm/trunk-jpl/test/Par/SquareSheetConstrained.py

@@ -59 +59 @@
 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)