Changeset 23899

Timestamp:

04/26/19 21:44:15 (6 years ago)

Author:

seroussi

Message:

CHG: finished ice only part of the test

File:

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

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

@@ -167 +167 @@
         %Grab lat,long from MITgcm: 
         long=readbin('run/XG.data',[3 200]);
-        long=[long long(:,end)]; long=[long; -105.0*ones(1,size(long,2))];
+        long=[long long(:,end)]; long=[long; -105.1250*ones(1,size(long,2))];
         lat=readbin('run/YG.data',[3 200]);
         lat=[lat -73.8832*ones(size(lat,1),1)]; lat=[lat; lat(end,:)];
@@ -256 +256 @@
         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 = [0*ones(md.mesh.numberofvertices,1); 1];
+        md.smb.mass_balance = 0*ones(md.mesh.numberofvertices,1);
 
         %Flow law 
@@ -307 +307 @@
         %timestepping: 
         md.timestepping.final_time=100;
-        md.timestepping.time_step=0.25;
+        md.timestepping.time_step=0.5;
         md.transient.isgroundingline=0;
         md.transient.isthermal=0;
@@ -313 +313 @@
         md.groundingline.melt_interpolation='NoMeltOnPartiallyFloating';
         md.groundingline.friction_interpolation='SubelementFriction2';
+        md.masstransport.stabilization=0;
 
         md.cluster=generic('name',oshostname(),'np',2);