Index: /issm/trunk-jpl/examples/PigSensitivity/FrontRetreat.exp
===================================================================
--- /issm/trunk-jpl/examples/PigSensitivity/FrontRetreat.exp	(revision 21017)
+++ /issm/trunk-jpl/examples/PigSensitivity/FrontRetreat.exp	(revision 21017)
@@ -0,0 +1,17 @@
+## Name:FrontRetreat
+## Icon:0
+# Points Count Value
+11 1.000000
+# X pos Y pos
+-1626858.6376043619 -327176.4356578548
+-1605190.4936934854 -325911.1425827671
+-1604716.0087903275 -316263.2828852235
+-1607088.4333061168 -312625.5652943464
+-1612465.9288752396 -311676.5954880306
+-1618476.0709819060 -312467.4036599604
+-1624011.7281854146 -313890.8583694341
+-1626858.6376043619 -315788.7979820656
+-1627807.6074106777 -319268.3539385567
+-1628756.5772169936 -322589.7482606619
+-1626858.6376043619 -327176.4356578548
+
Index: /issm/trunk-jpl/examples/PigSensitivity/Help.m
===================================================================
--- /issm/trunk-jpl/examples/PigSensitivity/Help.m	(revision 21017)
+++ /issm/trunk-jpl/examples/PigSensitivity/Help.m	(revision 21017)
@@ -0,0 +1,22 @@
+if step==4
+	%Load model
+	md = loadmodel('./Models/PIG_Transient');
+
+	%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;
+
+	%Define time steps and time span of the simulation
+	md.timestepping.time_step=0.1;
+	md.timestepping.final_time=10;
+
+	%Request additional outputs
+	md.transient.requested_outputs={'default','IceVolume','IceVolumeAboveFloatation'};
+
+	%Solve
+	md=solve(md,TransientSolutionEnum);
+
+	%Save model
+	save ./Models/PIG_SMB md;
+end
Index: /issm/trunk-jpl/examples/PigSensitivity/runme.m
===================================================================
--- /issm/trunk-jpl/examples/PigSensitivity/runme.m	(revision 21017)
+++ /issm/trunk-jpl/examples/PigSensitivity/runme.m	(revision 21017)
@@ -0,0 +1,75 @@
+step=[1];
+
+if step==1 %Transient Run #1
+
+	md = loadmodel('../Pig/Models/PIG_Control_drag');	
+
+	md.inversion.iscontrol=0;
+	md.transient.ismasstransport=1;
+	md.transient.isstressbalance=1;
+	md.transient.isgroundingline=1;
+	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=25*ones(md.mesh.numberofvertices,1);
+	
+	md.timestepping.time_step=0.1;
+	md.timestepping.final_time=10;
+	md.transient.requested_outputs={'default','IceVolume','IceVolumeAboveFloatation'};
+
+	md=solve(md,TransientSolutionEnum);
+
+	% Save model
+	save ./Models/PIG_Transient md;
+end
+
+if step==2 %High Melt #2
+	md = loadmodel('./Models/PIG_Transient');	
+
+	md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
+	md.basalforcings.floatingice_melting_rate=60*ones(md.mesh.numberofvertices,1);
+	
+	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==3 %Ice Front retreat
+	md = loadmodel('./Models/PIG_Transient');	
+
+	md2=extract(md,'~FrontRetreat.exp');
+
+	md2=SetMarineIceSheetBC(md2);
+
+	md2.basalforcings.groundedice_melting_rate=zeros(md2.mesh.numberofvertices,1);
+	md2.basalforcings.floatingice_melting_rate=25*ones(md2.mesh.numberofvertices,1);
+
+	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
+
+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