Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18265)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18266)
@@ -72,5 +72,4 @@
 		IssmDouble* b_pos       = xNew<IssmDouble>(numvertices); // Hs-SMB relation parameter
 		IssmDouble* b_neg       = xNew<IssmDouble>(numvertices); // Hs-SMB relation paremeter
-		IssmDouble* h           = xNew<IssmDouble>(numvertices); // ice thickness (m)
 		IssmDouble* s           = xNew<IssmDouble>(numvertices); // surface elevation (m)
 		IssmDouble* smb         = xNew<IssmDouble>(numvertices);
@@ -82,6 +81,5 @@
 		element->GetInputListOnVertices(b_neg,SurfaceforcingsBNegEnum);
 
-		/*Recover surface elevatio at vertices: */
-		element->GetInputListOnVertices(h,ThicknessEnum);
+		/*Recover surface elevation at vertices: */
 		element->GetInputListOnVertices(s,SurfaceEnum);
 
@@ -98,5 +96,5 @@
 				smb[v]=Smbref[v]+b_neg[v]*(s[v]-Href[v]);
 			}
-			smb[v]=smb[v]/rho_ice;      // SMB in m/y ice
+			smb[v]=smb[v]/1000*rho_water/rho_ice;      // SMB in m/y ice
 		}  //end of the loop over the vertices
 
@@ -107,5 +105,4 @@
 		xDelete<IssmDouble>(b_pos);
 		xDelete<IssmDouble>(b_neg);
-		xDelete<IssmDouble>(h);
 		xDelete<IssmDouble>(s);
 		xDelete<IssmDouble>(smb);