Changeset 23886

Timestamp:

04/22/19 20:53:33 (6 years ago)

Author:

seroussi

Message:

CHG: starting to put ice part in test4003

File:

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

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

@@ -269 +269 @@
 
         %friction: 
-        md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
-        pos=find(md.mask.groundedice_levelset<=0);
-        md.friction.coefficient(pos)=0;
+        md.friction.coefficient=0*ones(md.mesh.numberofvertices,1);
+        pos=find(md.mask.groundedice_levelset>0);
+        md.friction.coefficient(pos)=5;
         md.friction.p=ones(md.mesh.numberofelements,1);
         md.friction.q=ones(md.mesh.numberofelements,1);
 
         %temperatures and surface mass balance:
-        md.initialization.temperature=(273.15-20)*ones(md.mesh.numberofvertices,1);
+        md.initialization.temperature=(273.15-22)*ones(md.mesh.numberofvertices,1);
         md.initialization.pressure=md.materials.rho_ice*md.constants.g*(md.geometry.surface-md.geometry.base);
         md.smb.mass_balance = [1*ones(md.mesh.numberofvertices,1); 1];
@@ -302 +302 @@
 
         %get some flux at the ice divide: 
-        pos=find(md.mesh.lat==min(md.mesh.lat));
-        md.stressbalance.spcvy(pos)=200;
+        pos=find(md.mesh.y==min(md.mesh.y));
+        md.masstransport.spcthickness(pos)=md.geometry.thickness(pos);
+        md.stressbalance.vx(pos)=0;
+        md.stressbalance.vy(pos)=1500;
 
         %deal with boundaries, excluding icefront: 
@@ -309 +311 @@
         vertex_on_boundary(md.mesh.segments(:,1:2))=1;
         pos=find(vertex_on_boundary & md.mask.groundedice_levelset<=0);
-        md.stressbalance.spcvx(pos)=md.initialization.vx(pos);
-        md.stressbalance.spcvy(pos)=md.initialization.vy(pos);
-        md.stressbalance.spcvz(pos)=md.initialization.vz(pos);
-        md.masstransport.spcthickness(pos)=md.geometry.thickness(pos);
+        md.stressbalance.spcvx(pos)=0;
 
         md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
@@ -334 +333 @@
         md.timestepping.final_time=final_time;
         md.timestepping.time_step=time_step;
-        md.transient.isgroundingline=1;
+        md.transient.isgroundingline=0;
         md.transient.isthermal=0;
         md.groundingline.migration='SubelementMigration';
-        md.groundingline.melt_interpolation='SubelementMelt2';
+        md.groundingline.melt_interpolation='NoMeltOnPartiallyFloating';
         md.groundingline.friction_interpolation='SubelementFriction2';
 
@@ -449 +448 @@
         md.transient.isthermal=0;
         md.groundingline.migration='SubelementMigration';
-        md.groundingline.melt_interpolation='SubelementMelt2';
+        md.groundingline.melt_interpolation='NoMeltOnPartiallyFloating';
         md.groundingline.friction_interpolation='SubelementFriction2';