Changeset 26068

Timestamp:

03/10/21 15:12:42 (4 years ago)

Author:

jdquinn

Message:

CHG: Reverting changes as Eric and Lambert are still working on SE tests

File:

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

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

@@ -7 +7 @@
 %parameterize solidearth solution:
 %solidearth loading:  {{{
-md.dsl.global_average_thermosteric_sea_level=[0;0];
-md.dsl.sea_surface_height_above_geoid=zeros(md.mesh.numberofvertices+1,1);
-md.dsl.sea_water_pressure_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
+md.solidearth.surfaceload.icethicknesschange=zeros(md.mesh.numberofelements,1);
+md.solidearth.initialsealevel=zeros(md.mesh.numberofvertices,1);
+md.dsl.global_average_thermosteric_sea_level_change=[0;0];
+md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices+1,1);
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
+%antarctica
+late=sum(md.mesh.lat(md.mesh.elements),2)/3;
+longe=sum(md.mesh.long(md.mesh.elements),2)/3;
+pos=find(late < -80);
+md.solidearth.surfaceload.icethicknesschange(pos)=-100;
+%greenland
+pos=find(late>70 & late<80 & longe>-60 & longe<-30);
+md.solidearth.surfaceload.icethicknesschange(pos)=-100;
 
 %elastic loading from love numbers:
@@ -24 +34 @@
 md.mask.ice_levelset=icemask;
 md.mask.ocean_levelset=-icemask;
+
+%make sure that the elements that have loads are fully grounded:
+pos=find(md.solidearth.surfaceload.icethicknesschange);
+md.mask.ocean_levelset(md.mesh.elements(pos,:))=1;
+
+%make sure wherever there is an ice load, that the mask is set to ice:
+pos=find(md.solidearth.surfaceload.icethicknesschange);
+md.mask.ice_levelset(md.mesh.elements(pos,:))=-1;
 % }}}
 
@@ -33 +51 @@
 %Materials: 
 md.materials=materials('hydro');
+
 
 %Miscellaneous