Changeset 26800 for issm/trunk-jpl/test/Par/GiaIvinsBenchmarksAB.par

Timestamp:

01/20/22 17:02:07 (3 years ago)

Author:

caronlam

Message:

BUG: Sealevelchange core: load weights did not add up correctly; BUG: Sealevel change core viscous stacks were offset in time; partial load area phi was double counted, appearing both in Green functions and loads; NEW: love core optimization, paralelization supported, support for viscous rotational feedback; NEW:Sea level change core: support for viscous rotational feedback; CHG: grd loads now contain the average kg.m^{-2 over the (sub-)elemental area they apply in, not the average*phi, phi is simply accounted for in the loadarea or Green's function; NEW: nightly runs 2070, 2071, 2072, 2091, 2092 test Spada et al (2011) benchmark cases}

File:

• issm/trunk-jpl/test/Par/GiaIvinsBenchmarksAB.par (modified) (1 diff)

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

@@ -14 +14 @@
 md.geometry.surface=md.geometry.thickness+md.geometry.base; 
 
+
+%GIA parameters specific to Experiments A  and B
+
+md.materials=materials('litho','ice');
+md.materials.numlayers=2;
+md.materials.radius =  [10 6271e3 6371e3]; %100km litosphere, the rest is mantle material
+md.materials.density=  [3.38e3 3.36e3];
+md.materials.lame_mu=  [1.45e11         6.7e10];
+md.materials.viscosity=[1e21            1e40];
+
 %Ice density used for benchmarking, not 917 kg/m^3
 md.materials.rho_ice=1000; %kg m^3
-
-%GIA parameters specific to Experiments A  and B
-md.gia=giaivins();
-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_density=3.38;                                  %in g/cm^3