Index: /issm/workshop/2012/Talks/13_IceBridge/runme.m
===================================================================
--- /issm/workshop/2012/Talks/13_IceBridge/runme.m	(revision 14122)
+++ /issm/workshop/2012/Talks/13_IceBridge/runme.m	(revision 14123)
@@ -1,4 +1,4 @@
 clear all;
-steps=[1:5];
+steps=[1];
 
 %Location of SeaRISE dataset
@@ -21,7 +21,6 @@
 
 	%Mesh greenland without refinement in Jak basin
-	%md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8);
-	%plotmodel (md,'data','mesh');
-	%return
+	md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8);
+	return
 	
 	%Refine mesh in the region of Jakobshavn (resolution = 1000 m)
@@ -30,9 +29,6 @@
 		'./Exp_Par/Jak_outline.exp','node',1);
 	hmaxVertices(find(in))=1000;
-	
-	%Adapt the mesh to minimize error in velocity interpolation
 	md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,...
 		'err',8,'hmaxVertices',hmaxVertices);
-	plotmodel (md,'data','mesh');
 
 	%convert x,y coordinates (Polar stereo) to lat/lon
@@ -47,8 +43,8 @@
 
 	md = setmask(md,'','');
-	md = parameterize(md,'./Exp_Par/Greenland_solution.par');
+	md = parameterize(md,'./Exp_Par/Greenland.par');
 	md = setflowequation(md,'macayeal','all');
 
-	save ./Models/Greenland.Parameterization md; 
+	save ./Models/Greenland.Parameterization2 md; 
 end
 
@@ -134,5 +130,5 @@
 	plotmodel(md,'data',md.results.TransientSolution(end).Vel,'caxis',[1e-1 6000],...
 		'log', 10, 'title', 'Velocity (m/y)','gridded',1, ...
-		'data', md.results.TransientSolution(1).SurfaceforcingsMassBalance, ...
+		'data', md.results.TransientSolution(end).SurfaceforcingsMassBalance, ...
 		'title', 'Surface mass balance (m/y)','gridded',1, ...
 		'data',md.results.TransientSolution(end).Thickness,...