Index: /issm/trunk-jpl/examples/Greenland/runme.m
===================================================================
--- /issm/trunk-jpl/examples/Greenland/runme.m	(revision 20531)
+++ /issm/trunk-jpl/examples/Greenland/runme.m	(revision 20532)
@@ -1,4 +1,4 @@
 clear all;
-steps=[1 2 3 4 5];
+steps=[1];
 
 %Location of SeaRISE dataset
@@ -146,5 +146,5 @@
 	lat  = ncread(ncbox,'lat');
 	lon  = ncread(ncbox,'lon');
-	smbbox = ncread(ncbox,'SMB');
+	smbbox = ncread(ncbox,'MassFlux');
 	[x1 y1]=ll2xy(lat,lon,+1,45,70);
 
Index: /issm/trunk-jpl/examples/ISMIP/CheatyRunme.m
===================================================================
--- /issm/trunk-jpl/examples/ISMIP/CheatyRunme.m	(revision 20531)
+++ /issm/trunk-jpl/examples/ISMIP/CheatyRunme.m	(revision 20532)
@@ -3,5 +3,5 @@
 %step 8 is specific to ISMIPF
 
-steps=[1:7];
+steps=[1];
 
 % parameter file to be used, choose between CheatyIsmipA.par or CheatyIsmipF.par
Index: /issm/trunk-jpl/examples/IceBridge/Greenland.par
===================================================================
--- /issm/trunk-jpl/examples/IceBridge/Greenland.par	(revision 20531)
+++ /issm/trunk-jpl/examples/IceBridge/Greenland.par	(revision 20532)
@@ -62,6 +62,6 @@
 
 disp('   Interpolating surface mass balance');
-md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
-md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
+md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
+md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
 
 disp('   Construct basal friction parameters');
Index: /issm/trunk-jpl/examples/IceBridge/Greenland_cheatsheet.par
===================================================================
--- /issm/trunk-jpl/examples/IceBridge/Greenland_cheatsheet.par	(revision 20531)
+++ /issm/trunk-jpl/examples/IceBridge/Greenland_cheatsheet.par	(revision 20532)
@@ -66,6 +66,6 @@
 
 disp('   Interpolating surface mass balance');
-md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
-md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
+md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
+md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
 
 disp('   Construct basal friction parameters');
Index: /issm/trunk-jpl/examples/IceBridge/Greenland_noOIB.par
===================================================================
--- /issm/trunk-jpl/examples/IceBridge/Greenland_noOIB.par	(revision 20531)
+++ /issm/trunk-jpl/examples/IceBridge/Greenland_noOIB.par	(revision 20532)
@@ -40,6 +40,6 @@
 
 disp('   Interpolating surface mass balance');
-md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
-md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
+md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
+md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
 
 disp('   Construct basal friction parameters');
Index: /issm/trunk-jpl/examples/IceBridge/runme.m
===================================================================
--- /issm/trunk-jpl/examples/IceBridge/runme.m	(revision 20531)
+++ /issm/trunk-jpl/examples/IceBridge/runme.m	(revision 20532)
@@ -46,5 +46,5 @@
 
 	md = setmask(md,'','');
-	md = parameterize(md,'./Greenland_cheatsheet.par');
+	md = parameterize(md,'./Greenland.par');
 	%md = parameterize(md,'./Greenland_noOIB.par');
 	md = setflowequation(md,'SSA','all');
@@ -108,5 +108,5 @@
 	smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0)*1000/md.materials.rho_ice;
 	smb = [smb smb smb-1.0];
-	md.surfaceforcings.mass_balance = [smb;1 10 20];
+	md.smb.mass_balance = [smb;1 10 20];
 
 	%Set transient options, run for 20 years, saving every year
@@ -117,5 +117,5 @@
 	%Additional options
 	md.inversion.iscontrol=0;
-	md.transient.requested_outputs={'IceVolume','TotalSmb','SurfaceforcingsMassBalance'};
+	md.transient.requested_outputs={'IceVolume','TotalSmb','SmbMassBalance'};
 	md.verbose=verbose('solution',true,'module',true,'convergence',true);
 
@@ -134,5 +134,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(end).SurfaceforcingsMassBalance, ...
+		'data', md.results.TransientSolution(end).SmbMassBalance, ...
 		'title', 'Surface mass balance (m/y)','gridded',1, ...
 		'data',md.results.TransientSolution(end).Thickness,...
@@ -146,5 +146,5 @@
 	%Plot surface mass balance
 	surfmb=[]; for i=1:100; surfmb=[surfmb ...
-		md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end
+		md.results.TransientSolution(i).SmbMassBalance]; end
 	subplot(3,1,1); plot([0.2:0.2:20],mean(surfmb)); title('Mean Surface mass balance');
 
Index: /issm/trunk-jpl/examples/IceflowModels/EISMINT.par
===================================================================
--- /issm/trunk-jpl/examples/IceflowModels/EISMINT.par	(revision 20531)
+++ /issm/trunk-jpl/examples/IceflowModels/EISMINT.par	(revision 20532)
@@ -29,5 +29,5 @@
 sb=10^-2/1000.; %m/yr/m
 rel=450.*1000.; %m
-md.surfaceforcings.mass_balance=min(smb_max,sb*(rel-radius));
+md.smb.mass_balance=min(smb_max,sb*(rel-radius));
 
 disp('      creating velocities');
Index: /issm/trunk-jpl/examples/Jakobshavn/Jks.par
===================================================================
--- /issm/trunk-jpl/examples/Jakobshavn/Jks.par	(revision 20531)
+++ /issm/trunk-jpl/examples/Jakobshavn/Jks.par	(revision 20532)
@@ -38,6 +38,6 @@
 
 disp('   Interpolating surface mass balance');
-md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
-md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
+md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
+md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
 
 disp('   Construct basal friction parameters');
Index: /issm/trunk-jpl/examples/Jakobshavn/runme.m
===================================================================
--- /issm/trunk-jpl/examples/Jakobshavn/runme.m	(revision 20531)
+++ /issm/trunk-jpl/examples/Jakobshavn/runme.m	(revision 20532)
@@ -1,3 +1,3 @@
-steps=[4];
+steps=[1];
 
 if any(steps==1)
Index: /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m
===================================================================
--- /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m	(revision 20531)
+++ /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m	(revision 20532)
@@ -1,4 +1,4 @@
 %PIG Uncertainty Quantification Application
-steps=[1:3]; 
+steps=[1]; 
 
 if any(steps==1) 
@@ -140,8 +140,17 @@
 
 	md.stressbalance.restol=10^-5; %tighten tolerances for UQ analyses
+
+	%Turn off verbosity
+	md.verbose=verbose(0);
 	
 	%solve
 	md.qmu.isdakota=1; md.inversion.iscontrol=0;
 	md.cluster=generic('name',oshostname,'np',4);
+
+	%Dakota runs in parallel with a master/slave configuration.
+	% At least 2 cpu's are needed to run the UQ
+	md.qmu.params.evaluation_scheduling='master';
+	md.qmu.params.processors_per_evaluation=md.cluster.np-1;
+
 	md=solve(md,StressbalanceSolutionEnum,'overwrite','y');
 
@@ -217,4 +226,10 @@
 	md.inversion.iscontrol=0;
 	md.cluster=generic('name',oshostname,'np',2);
+
+	%Dakota runs in parallel with a master/slave configuration.
+	% At least 2 cpu's are needed to run the UQ
+	md.qmu.params.evaluation_scheduling='master';
+	md.qmu.params.processors_per_evaluation=md.cluster.np-1
+
 	md.verbose=verbose('qmu',true);
 	md=solve(md,StressbalanceSolutionEnum,'overwrite','y');