source: issm/oecreview/Archive/20545-21336/ISSM-20751-20752.diff@ 21337

Index: ../trunk-jpl/examples/MeltingSensitivity/runme.m
===================================================================
--- ../trunk-jpl/examples/MeltingSensitivity/runme.m    (revision 20751)
+++ ../trunk-jpl/examples/MeltingSensitivity/runme.m    (revision 20752)
@@ -1,8 +1,8 @@
-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;
        md.transient.ismasstransport=1;
@@ -11,8 +11,9 @@
        md.transient.ismovingfront=0;
        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;
        md.timestepping.final_time=10;
@@ -43,9 +44,9 @@
        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);
+       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;
        md.timestepping.final_time=10;
@@ -60,14 +61,18 @@
 if step==4 %Ice Front retreat
        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