source: issm/oecreview/Archive/21724-22754/ISSM-22727-22728.diff

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);