Index: /issm/trunk-jpl/examples/MeltingSensitivity/runme.m
===================================================================
--- /issm/trunk-jpl/examples/MeltingSensitivity/runme.m	(revision 20751)
+++ /issm/trunk-jpl/examples/MeltingSensitivity/runme.m	(revision 20752)
@@ -1,7 +1,7 @@
-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,6 +12,7 @@
 	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,7 +45,7 @@
 
 	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,13 +62,17 @@
 	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