Changeset 20755

Timestamp:

06/16/16 12:39:16 (9 years ago)

Author:

seroussi

Message:

CHG: finished melting sensitivity run

Location:

issm/trunk-jpl/examples/MeltingSensitivity

Files:

• Help.m (modified) (1 diff)
• runme.m (modified) (4 diffs)

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

@@ -1 @@
-if any(steps==6)
-        % Load Model
-        md = loadmodel('./Models/PIG.Control_drag');
-        md.inversion.iscontrol=0;
+if step==4
+        %Load model
+        md = loadmodel('./Models/PIG_Transient');
 
-        disp('   Extruding mesh')
-        number_of_layers=3;
-        md=extrude(md,number_of_layers,0.9);
+        %Change external forcing basal melting rate and surface mass balance)
+        md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
+        md.basalforcings.floatingice_melting_rate=25*ones(md.mesh.numberofvertices,1);
+        md.smb.mass_balance=2*md.smb.mass_balance;
 
-        disp('   Using HO Ice Flow Model')
-        md=setflowequation(md, 'HO', 'all');
+        %Define time steps and time span of the simulation
+        md.timestepping.time_step=0.1;
+        md.timestepping.final_time=10;
 
-        % Solve
-        md=solve(md,StressbalanceSolutionEnum);
+        %Request additional outputs
+        md.transient.requested_outputs={'default','IceVolume','IceVolumeAboveFloatation'}
 
-        % Save Model
-        save ./Models/PIG.ModelHO md;
+        %Solve
+        md=solve(md,TransientSolutionEnum);
+
+        %Save model
+        save ./Models/PIG_SMB md;
 end

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

@@ -3 +3 @@
 if step==1 %Transient Run #1
 
-        md = loadmodel('../Pig/Models/PIG.Control_drag');       
+        md = loadmodel('../Pig/Models/PIG_Control_drag');       
 
         md.inversion.iscontrol=0;
@@ -23 +23 @@
 
         % Save model
-        save ./Models/PIG.Transient md;
+        save ./Models/PIG_Transient md;
 end
 
 if step==2 %High Melt #2
-        md = loadmodel('./Models/PIG.Transient');       
+        md = loadmodel('./Models/PIG_Transient');       
 
         md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
@@ -41 +41 @@
 end
 
-if step==3 %High surface mass balance #3
-        md = loadmodel('./Models/PIG.Transient');       
-
-        md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
-        md.basalforcings.floatingice_melting_rate=25*ones(md.mesh.numberofvertices,1);
-
-        md.smb.mass_balance=0*md.smb.mass_balance;
-        
-        md.timestepping.time_step=0.1;
-        md.timestepping.final_time=10;
-
-        md.transient.requested_outputs={'default','IceVolume','IceVolumeAboveFloatation'}
-
-        md=solve(md,TransientSolutionEnum);
-
-        save ./Models/PIG.HighMelt md;
-end
-
-if step==4 %Ice Front retreat
-        md = loadmodel('./Models/PIG.Transient');       
+if step==3 %Ice Front retreat
+        md = loadmodel('./Models/PIG_Transient');       
 
         md2=extract(md,'~FrontRetreat.exp');
@@ -75 +57 @@
         md2=solve(md2,TransientSolutionEnum);
 
-        save ./Models/PIG.FrontRetreat md2;
+        save ./Models/PIG_FrontRetreat md2;
 end
+
+if step==4 %High surface mass balance #3
+        %Load model
+
+        %Change external forcing basal melting rate and surface mass balance)
+        
+        %Refine time steps and time span of the simulation
+
+        %Request additional outputs
+
+        %Solve
+
+        %Save model
+
+end