Index: ../trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp
===================================================================
--- ../trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18583)
+++ ../trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18584)
@@ -30,19 +30,19 @@
 			PositiveDegreeDayx(femmodel);
 			break;
 		case SMBgradientsEnum:
-			if(VerboseSolution())_printf_("	call smb gradients module\n");
+			if(VerboseSolution())_printf0_("	call smb gradients module\n");
 			SmbGradientsx(femmodel);
 			break;
 		case SMBhenningEnum:
-			if(VerboseSolution())_printf_("  call smb Henning module\n");
+			if(VerboseSolution())_printf0_("  call smb Henning module\n");
 			SmbHenningx(femmodel);
 			break;
 		case SMBcomponentsEnum:
-			if(VerboseSolution())_printf_("  call smb Components module\n");
+			if(VerboseSolution())_printf0_("  call smb Components module\n");
 			SmbComponentsx(femmodel);
 			break;
 		case SMBmeltcomponentsEnum:
-			if(VerboseSolution())_printf_("  call smb Melt Components module\n");
+			if(VerboseSolution())_printf0_("  call smb Melt Components module\n");
 			SmbMeltComponentsx(femmodel);
 			break;
 		default:
@@ -262,10 +262,11 @@
 			  smb = (a + b*z)*(f + g*(z-z_critical) + h*(z-z_critical)*(z-z_critical)) + c;
 			}
             
-            /* Compute smb including anomaly,
-             correct for number of seconds in a year [s/yr]*/
-            smb = smb + anomaly*yts;
-            
+			/* Compute smb including anomaly,
+				correct for number of seconds in a year [s/yr]*/
+			smb = smb/yts + anomaly;
+
+
 			/*Update array accordingly*/
 			smblist[v] = smb;
 
