Changeset 26067

Timestamp:

03/10/21 12:51:48 (4 years ago)

Author:

jdquinn

Message:

CHG: Removed references to solidearth.surfaceload; committing so Eric can advise on solution type

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.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;
+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);
 
 %elastic loading from love numbers:
@@ -34 +24 @@
 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;
 % }}}
 
@@ -51 +33 @@
 %Materials: 
 md.materials=materials('hydro');
-
 
 %Miscellaneous