Index: /issm/trunk/examples/Jakobshavn/Jks.par
===================================================================
--- /issm/trunk/examples/Jakobshavn/Jks.par	(revision 16453)
+++ /issm/trunk/examples/Jakobshavn/Jks.par	(revision 16454)
@@ -44,6 +44,5 @@
 disp('   Construct basal friction parameters');
 md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
-pos=find(md.mask.vertexonfloatingice);
-md.friction.coefficient(pos)=0; %no friction applied on floating ice
+md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
 md.friction.p=ones(md.mesh.numberofelements,1);
 md.friction.q=ones(md.mesh.numberofelements,1);
@@ -58,4 +57,3 @@
 disp('   Set other boundary conditions');
 md=SetMarineIceSheetBC(md,'./Front.exp');
-md.diagnostic.icefront(:,end)=1; %Apply water pressure at the ice front
 md.basalforcings.melting_rate=zeros(md.mesh.numberofvertices,1);
Index: /issm/trunk/examples/Jakobshavn/runme.m
===================================================================
--- /issm/trunk/examples/Jakobshavn/runme.m	(revision 16453)
+++ /issm/trunk/examples/Jakobshavn/runme.m	(revision 16454)
@@ -42,5 +42,5 @@
 	md=loadmodel('JksPar.mdl');
 
-	md=setflowequation(md,'macayeal','all');
+	md=setflowequation(md,'SSA','all');
 
 	%Control general
@@ -64,13 +64,11 @@
 
 	%Additional parameters
-	md.diagnostic.restol=0.01;
-	md.diagnostic.reltol=0.1;
-	md.diagnostic.abstol=NaN;
+	md.stressbalance.restol=0.01;
+	md.stressbalance.reltol=0.1;
+	md.stressbalance.abstol=NaN;
 
 	%Go solve
 	md.cluster=generic('name',oshostname,'np',4);
-	md.solver=addoptions(md.solver,NoneAnalysisEnum,asmoptions);
-	md.solver=addoptions(md.solver,DiagnosticVertAnalysisEnum,jacobiasmoptions);
-	md=solve(md,DiagnosticSolutionEnum);
+	md=solve(md,StressbalanceSolutionEnum);
 	
 	save JksControl.mdl md
@@ -83,9 +81,9 @@
 		'FontSize#all',12,...
 		'data',md.inversion.vel_obs,'title','Observed velocity',...
-		'data',md.results.DiagnosticSolution.Vel,'title','Modeled Velocity',...
+		'data',md.results.StressbalanceSolution.Vel,'title','Modeled Velocity',...
 		'colorbar#1','off','colorbar#2','on','colorbartitle#2','[m/yr]',...
 		'caxis#1-2',[0,7000],...
 		'data',md.geometry.bed,'title','Bed elevation',...
-		'data',md.results.DiagnosticSolution.FrictionCoefficient,...
+		'data',md.results.StressbalanceSolution.FrictionCoefficient,...
 		'title','Friction Coefficient',...
 		'colorbar#3','on','colorbartitle#3','[m]','colorbar#4','on');
Index: /issm/trunk/examples/SquareIceShelf/Square.par
===================================================================
--- /issm/trunk/examples/SquareIceShelf/Square.par	(revision 16453)
+++ /issm/trunk/examples/SquareIceShelf/Square.par	(revision 16454)
@@ -11,7 +11,6 @@
 
 disp('      creating drag');
-pos=find(md.mask.vertexonfloatingice);
 md.friction.coefficient=200*ones(md.mesh.numberofvertices,1);
-md.friction.coefficient(pos)=0;
+md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
 md.friction.p=ones(md.mesh.numberofelements,1);
 md.friction.q=ones(md.mesh.numberofelements,1);
Index: /issm/trunk/examples/SquareIceShelf/runme.m
===================================================================
--- /issm/trunk/examples/SquareIceShelf/runme.m	(revision 16453)
+++ /issm/trunk/examples/SquareIceShelf/runme.m	(revision 16454)
@@ -3,5 +3,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'Square.par');
-md=setflowequation(md,'macayeal','all');
+md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname,'np',2);
-md=solve(md,DiagnosticSolutionEnum);
+md=solve(md,StressbalanceSolutionEnum);