Changeset 22979

Timestamp:

07/19/18 15:56:14 (7 years ago)

Author:

Eric.Larour

Message:

CHG: new logic to account for fractional groundedness. Tends to spread the load if you don't take care of
putting the grounding line in the right position. Use expcontourtolevelset to do this.

File:

• issm/trunk-jpl/test/NightlyRun/test2002.m (modified) (3 diffs)

issm/trunk-jpl/test/NightlyRun/test2002.m

@@ -41 +41 @@
 md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1); 
 
+%make sure that the elements that have loads are fully grounded:
+pos=find(md.slr.deltathickness);
+md.mask.groundedice_levelset(md.mesh.elements(pos,:))=1;
+
 %make sure wherever there is an ice load, that the mask is set to ice: 
-pos=find(md.slr.deltathickness); md.mask.ice_levelset(md.mesh.elements(pos,:))=-1;
+pos=find(md.slr.deltathickness); 
+md.mask.ice_levelset(md.mesh.elements(pos,:))=-1;
 % }}}
 
@@ -70 +75 @@
 md.slr.reltol=NaN;
 md.slr.abstol=1e-3;
+md.slr.geodetic=1;
 
 % max number of iteration reverted back to 10 (i.e., the original default value)
@@ -75 +81 @@
 
 %eustatic run: 
-md.slr.rigid=0; md.slr.elastic=0;
+md.slr.rigid=0; md.slr.elastic=0;md.slr.rotation=0;
 md=solve(md,'Sealevelrise');
 Seustatic=md.results.SealevelriseSolution.Sealevel;
 
 %eustatic + rigid run: 
-md.slr.rigid=1; md.slr.elastic=0;
+md.slr.rigid=1; md.slr.elastic=0;md.slr.rotation=0;
 md=solve(md,'Sealevelrise');
 Srigid=md.results.SealevelriseSolution.Sealevel;
 
 %eustatic + rigid + elastic run: 
-md.slr.rigid=1; md.slr.elastic=1;
+md.slr.rigid=1; md.slr.elastic=1;md.slr.rotation=0;
 md=solve(md,'Sealevelrise');
 Selastic=md.results.SealevelriseSolution.Sealevel;