Index: /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 26631)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 26632)
@@ -357,4 +357,8 @@
 		frictionheat=alpha2*(vx*vx+vy*vy);
 
+		/* *** TEST take out frictional heat for large gap height *** */
+		if(gap>br)
+		 frictionheat=0.;
+
 		/*Get water and ice pressures*/
 		IssmDouble pressure_ice   = rho_ice*g*thickness;    _assert_(pressure_ice>0.); 
@@ -648,4 +652,8 @@
 		frictionheat=alpha2*(vx*vx+vy*vy);
 
+		/* *** TEST take out frictional heat for large gap height *** */
+		if(gap>br)
+		 frictionheat=0.;
+
 		/*Get water and ice pressures*/
 		IssmDouble pressure_ice   = rho_ice*g*thickness;    _assert_(pressure_ice>0.);