source: issm/oecreview/Archive/21724-22754/ISSM-22237-22238.diff

Index: ../trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- ../trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp    (revision 22237)
+++ ../trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp    (revision 22238)
@@ -1227,6 +1227,8 @@
        IssmDouble re_bot;
        IssmDouble gdn_bot;
        IssmDouble gsp_bot;
+       IssmDouble EI_bot;
+       IssmDouble EW_bot;
        bool        top=false;
     
        IssmDouble* Zcum=NULL;
@@ -1650,6 +1652,8 @@
                re_bot=re[n-1];
                gdn_bot=gdn[n-1];
                gsp_bot=gsp[n-1];
+               EI_bot=EI[n-1];
+               EW_bot=EW[n-1];
         
                dz_add=dz_bot;
         
@@ -1662,6 +1666,8 @@
                newcell(&re,re_bot,top,n);
                newcell(&gdn,gdn_bot,top,n);
                newcell(&gsp,gsp_bot,top,n);
+               newcell(&EI,EI_bot,top,n);
+               newcell(&EW,EW_bot,top,n);
                n=n+1;
        }
        else if (Ztot > zMax){
@@ -1685,6 +1691,8 @@
                celldelete(&re,n,D,D_size);
                celldelete(&gdn,n,D,D_size);
                celldelete(&gsp,n,D,D_size);
+               celldelete(&EI,n,D,D_size);
+               celldelete(&EW,n,D,D_size);
                n=D_size;
                xDelete<int>(D);