Index: /issm/trunk-jpl/examples/Greenland/Greenland.par
===================================================================
--- /issm/trunk-jpl/examples/Greenland/Greenland.par	(revision 19627)
+++ /issm/trunk-jpl/examples/Greenland/Greenland.par	(revision 19628)
@@ -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/Greenland/runme.m
===================================================================
--- /issm/trunk-jpl/examples/Greenland/runme.m	(revision 19627)
+++ /issm/trunk-jpl/examples/Greenland/runme.m	(revision 19628)
@@ -92,5 +92,5 @@
 	smb = smb*md.materials.rho_freshwater/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 timestep
@@ -102,5 +102,5 @@
 	md.inversion.iscontrol=0;
 	md.transient.requested_outputs={'IceVolume','TotalSmb', ...
-		                             'SurfaceforcingsMassBalance'};
+		                             'SmbMassBalance'};
 	md.verbose=verbose('solution',true,'module',true,'convergence',true);
 
@@ -124,5 +124,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');
 
@@ -181,5 +181,5 @@
 	smb_mo = InterpFromMeshToMesh2d(index,x1(:),y1(:),smbmean(:),xi,yi);
 	smb = smb_mo*12/1000*md.materials.rho_freshwater/md.materials.rho_ice;
-	md.surfaceforcings.mass_balance = [smb;1 ];
+	md.smb.mass_balance = [smb;1 ];
 
 	%Set transient options, run for 20 years, saving every timestep
@@ -191,5 +191,5 @@
 	md.inversion.iscontrol=0;
 	md.transient.requested_outputs={'IceVolume','TotalSmb', ...
-		'SurfaceforcingsMassBalance'};
+		'SmbMassBalance'};
 	md.verbose=verbose('solution',true,'module',true);
 
@@ -272,5 +272,5 @@
 		end
 	end
-	md.surfaceforcings.mass_balance = ...
+	md.smb.mass_balance = ...
 		[smb*12/1000*md.materials.rho_freshwater/md.materials.rho_ice; ...
 		[1/24:1/12:length(years_of_simulation)]];
@@ -284,5 +284,5 @@
 	md.inversion.iscontrol=0;
 	md.transient.requested_outputs={'IceVolume','TotalSmb', ...
-		'SurfaceforcingsMassBalance'};
+		'SmbMassBalance'};
 	md.verbose=verbose('solution',true,'module',true);