Changeset 26739

Timestamp:

12/20/21 12:25:02 (3 years ago)

Author:

seroussi

Message:

CHG: working on ice ocean test

File:

• issm/trunk-jpl/test/NightlyRun/test4003.m (modified) (4 diffs)

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

@@ -4 +4 @@
 %Script control parameters
 steps=1:11;
-steps=[1:3 6:11];
+%steps=[1:3 6:11];
 
 %To download and recompile MITgcm from scratch:
-%!rm -rf ${ISSM_DIR}/test/MITgcm/install
-%!rm -rf ${ISSM_DIR}/test/MITgcm/build/*
-%!rm -rf Models
+!rm -rf ${ISSM_DIR}/test/MITgcm/install
+!rm -rf ${ISSM_DIR}/test/MITgcm/build/*
+!rm -rf Models
 
 %Organizer
@@ -273 +273 @@
 
     %level sets:
-    md.mask.ocean_levelset=iceshelf_mask(:);
-    md.mask.ice_levelset=-ones(md.mesh.numberofvertices,1);
-
-    pos=find(~ice_mask(:) & thickness(:)==0); md.mask.ice_levelset(pos)=1;
+    md.mask.ice_levelset=iceshelf_mask(:);
+    pos=find(md.mask.ice_levelset==1); md.mask.ice_levelset(pos)=-1;
+    pos=find(md.mask.ice_levelset==0); md.mask.ice_levelset(pos)=1;
+    md.mask.ocean_levelset=-ones(md.mesh.numberofvertices,1);
 
     savemodel(org,md);
@@ -324 +324 @@
         pos=find(md.mesh.lat==min(md.mesh.lat));
         md.masstransport.spcthickness(pos)=md.geometry.thickness(pos);
-        md.stressbalance.spcvx(pos)=-800;
-        md.stressbalance.spcvy(pos)=0;
+        md.stressbalance.spcvy(pos)=800;
+        md.stressbalance.spcvx(pos)=0;
 
         %deal with boundaries, excluding icefront:
         pos=find(md.mesh.long==min(md.mesh.long) | md.mesh.long==max(md.mesh.long));
-        md.stressbalance.spcvy(pos)=0;
+        md.stressbalance.spcvx(pos)=0;
 
         point1=find(md.mesh.y==min(md.mesh.y)); point2=find(md.mesh.x==max(md.mesh.x));
@@ -362 +362 @@
         md.groundingline.melt_interpolation='NoMeltOnPartiallyFloating';
         md.groundingline.friction_interpolation='SubelementFriction2';
-        md.masstransport.stabilization=0;
+        md.masstransport.stabilization=1;
 
         md.cluster=generic('name',oshostname(),'np',2);