Index: /issm/trunk-jpl/test/NightlyRun/test2002.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 26067)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 26068)
@@ -7,7 +7,17 @@
 %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,4 +34,12 @@
 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,4 +51,5 @@
 %Materials: 
 md.materials=materials('hydro');
+
 
 %Miscellaneous