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

../trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp

@@ -820 +820 @@
                 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 +1512 @@
                 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 +1539 @@
                 xDelete<IssmDouble>(F);
                 xDelete<IssmDouble>(R);
         }
-    
+
         //Merging of cells as they are burried under snow.
         Zcum=xNew<IssmDouble>(n);
         dzMin2=xNew<IssmDouble>(n);
@@ -1668 +1668 @@
                                 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);