Changeset 20752

Timestamp:

06/16/16 11:55:29 (9 years ago)

Author:

seroussi

Message:

CHG: finishing Pig sensitivity tutorial

File:

• issm/trunk-jpl/examples/MeltingSensitivity/runme.m (modified) (4 diffs)

issm/trunk-jpl/examples/MeltingSensitivity/runme.m

@@ -1 @@
-step=[1];
+step=[4];
 
 if step==1 %Transient Run #1
 
-        md = loadmodel('./Models/PIG.Control_drag');    
+        md = loadmodel('../Pig/Models/PIG.Control_drag');       
 
         md.inversion.iscontrol=0;
@@ -12 +12 @@
         md.transient.isthermal=0;
         
+        pos=find(md.mask.groundedice_levelset<0);
         md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
-        md.basalforcings.floatingice_melting_rate=30*ones(md.mesh.numberofvertices,1);
+        md.basalforcings.floatingice_melting_rate=25*ones(md.mesh.numberofvertices,1);
         
         md.timestepping.time_step=0.1;
@@ -44 +45 @@
 
         md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
-        md.basalforcings.floatingice_melting_rate=30*ones(md.mesh.numberofvertices,1);
+        md.basalforcings.floatingice_melting_rate=25*ones(md.mesh.numberofvertices,1);
 
-        md.smb.mass_balance=2*md.smb.mass_balance;
+        md.smb.mass_balance=0*md.smb.mass_balance;
         
         md.timestepping.time_step=0.1;
@@ -61 +62 @@
         md = loadmodel('./Models/PIG.Transient');       
 
-        md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
-        md.basalforcings.floatingice_melting_rate=30*ones(md.mesh.numberofvertices,1);
+        md2=extract(md,'FrontRetreat.exp');
 
-        md.timestepping.time_step=0.1;
-        md.timestepping.final_time=10;
-        md.transient.requested_outputs={'default','IceVolume','IceVolumeAboveFloatation'}
+        md2=SetMarineIceSheetBC(md2);
 
-        md=solve(md,TransientSolutionEnum);
+        md2.basalforcings.groundedice_melting_rate=zeros(md2.mesh.numberofvertices,1);
+        md2.basalforcings.floatingice_melting_rate=25*ones(md2.mesh.numberofvertices,1);
 
-        save ./Models/PIG.HighMelt md;
+        md2.timestepping.time_step=0.1;
+        md2.timestepping.final_time=10;
+        md2.transient.requested_outputs={'default','IceVolume','IceVolumeAboveFloatation'}
+
+        md2=solve(md2,TransientSolutionEnum);
+
+        save ./Models/PIG.FrontRetreat md2;
 end