Index: ../trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- ../trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 22727)
+++ ../trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 22728)
@@ -820,7 +820,7 @@
 		dT_turb = turb  / TCs;
 
 		// upward longwave contribution
-		ulw = - SB * pow(Ts,4.0)* teValue * dt;
+		ulw = - (SB * pow(Ts,4.0)* teValue) * dt ;
 		dT_ulw = ulw / TCs;
 		
 		// new grid point temperature
@@ -1512,7 +1512,7 @@
 		for(int i=0;i<n;i++)W[i] += dW[i];
 
 		//calculate Rsum:
-		Rsum=cellsum(R,n);
+		Rsum=cellsum(R,n) + flxDn[n];
 
 		// delete all cells with zero mass
 		D_size=0; for(int i=0;i<n;i++)if(m[i]>0.0)D_size++; 
@@ -1539,7 +1539,7 @@
 		xDelete<IssmDouble>(F);
 		xDelete<IssmDouble>(R);
 	}
-    
+
 	//Merging of cells as they are burried under snow.
 	Zcum=xNew<IssmDouble>(n);
 	dzMin2=xNew<IssmDouble>(n);
@@ -1668,8 +1668,8 @@
 				cellsplit(&T, n, j,1.0);
 				cellsplit(&d, n, j,1.0);
 				cellsplit(&a, n, j,1.0);
-				cellsplit(&EI, n, j,1.0);
-				cellsplit(&EW, n, j,1.0);
+				cellsplit(&EI, n, j,.5);
+				cellsplit(&EW, n, j,.5);
 				cellsplit(&re, n, j,1.0);
 				cellsplit(&gdn, n, j,1.0);
 				cellsplit(&gsp, n, j,1.0);