Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 24735)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 24736)
@@ -937,6 +937,6 @@
 
 		//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
@@ -947,5 +947,5 @@
 
 		// 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]