source: issm/oecreview/Archive/24684-25833/ISSM-24735-24736.diff@ 25834

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

@@ -936 +936 @@
                 // new grid point temperature
 
                 //SW penetrates surface
-                for(int j=0;j<m;j++) T[j] = T[j] + dT_sw[j];
-                T[0] = T[0] + dT_dlw + dT_ulw + dT_turb;
+                if(!isconstrainsurfaceT) for(int j=0;j<m;j++) T[j] = T[j] + dT_sw[j];
+                if(!isconstrainsurfaceT) T[0] = T[0] + dT_dlw + dT_ulw + dT_turb;
 
                 // temperature diffusion
                 for(int j=0;j<m;j++) T0[1+j]=T[j];
@@ -946 +946 @@
                 for(int j=0;j<m;j++) T[j] = (Np[j] * T[j]) + (Nu[j] * Tu[j]) + (Nd[j] * Td[j]);
 
                 // calculate cumulative evaporation (+)/condensation(-)
-                EC = EC + (EC_day/dts)*dt;
+                if(!isconstrainsurfaceT) EC = EC + (EC_day/dts)*dt;
 
                 /* CHECK FOR ENERGY (E) CONSERVATION [UNITS: J]
                 //energy flux across lower boundary (energy supplied by underling ice)