Index: /issm/trunk-jpl/examples/MeltingSensitivity/Help.m
===================================================================
--- /issm/trunk-jpl/examples/MeltingSensitivity/Help.m	(revision 20754)
+++ /issm/trunk-jpl/examples/MeltingSensitivity/Help.m	(revision 20755)
@@ -1,18 +1,22 @@
-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
Index: /issm/trunk-jpl/examples/MeltingSensitivity/runme.m
===================================================================
--- /issm/trunk-jpl/examples/MeltingSensitivity/runme.m	(revision 20754)
+++ /issm/trunk-jpl/examples/MeltingSensitivity/runme.m	(revision 20755)
@@ -3,5 +3,5 @@
 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,9 +23,9 @@
 
 	% 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,24 +41,6 @@
 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,4 +57,19 @@
 	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