Index: /issm/workshop/2012/Talks/12_SeaRISE/Code/plotparam.m
===================================================================
--- /issm/workshop/2012/Talks/12_SeaRISE/Code/plotparam.m	(revision 14074)
+++ /issm/workshop/2012/Talks/12_SeaRISE/Code/plotparam.m	(revision 14075)
@@ -1,2 +1,2 @@
->> plotmodel(md,'data',md.initialization.thickness);
+>> plotmodel(md,'data',md.geometry.thickness);
 >> plotmodel(md,'data',md.initialization.vel,'caxis',[1e-1 1e4],'log',10);
Index: /issm/workshop/2012/Talks/12_SeaRISE/runme.m
===================================================================
--- /issm/workshop/2012/Talks/12_SeaRISE/runme.m	(revision 14074)
+++ /issm/workshop/2012/Talks/12_SeaRISE/runme.m	(revision 14075)
@@ -22,5 +22,4 @@
 	%Mesh Greenland
 	md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8);
-	plotmodel (md,'data','mesh');
 	
 	%convert x,y coordinates (Polar stereo) to lat/lon
@@ -28,4 +27,6 @@
 
 	save ./Models/Greenland.Mesh_generation md;
+
+	plotmodel (md,'data','mesh');
 end
 
@@ -50,5 +51,4 @@
 	md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
 	md.inversion.maxiter_per_step=5*ones(md.inversion.nsteps,1);
-	md.verbose=verbose('solution',true,'control',true);
 
 	%Cost functions
@@ -66,6 +66,5 @@
 
 	%Additional parameters
-	md.diagnostic.restol=0.01;
-	md.diagnostic.reltol=0.1;
+	md.diagnostic.restol=0.01; md.diagnostic.reltol=0.1; 
 	md.diagnostic.abstol=NaN;
 
@@ -91,5 +90,6 @@
 	y1  = ncread(ncdata,'y1');
 	smb = ncread(ncdata,'smb');
-	smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0)*1000/md.materials.rho_ice;
+	smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0);
+	smb = smb*md.materials.rho_freshwater/md.materials.rho_ice;
 	smb = [smb smb smb-1.0];
 	md.surfaceforcings.mass_balance = [smb;1 10 20];