Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 22727)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 22728)
@@ -821,5 +821,5 @@
 
 		// upward longwave contribution
-		ulw = - SB * pow(Ts,4.0)* teValue * dt;
+		ulw = - (SB * pow(Ts,4.0)* teValue) * dt ;
 		dT_ulw = ulw / TCs;
 		
@@ -1513,5 +1513,5 @@
 
 		//calculate Rsum:
-		Rsum=cellsum(R,n);
+		Rsum=cellsum(R,n) + flxDn[n];
 
 		// delete all cells with zero mass
@@ -1540,5 +1540,5 @@
 		xDelete<IssmDouble>(R);
 	}
-    
+
 	//Merging of cells as they are burried under snow.
 	Zcum=xNew<IssmDouble>(n);
@@ -1669,6 +1669,6 @@
 				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);