Index: /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 26945)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 26946)
@@ -384,5 +384,6 @@
 
    	meltrate = 1/latentheat*(G+frictionheat+rho_water*g*conductivity*(dh[0]*dh[0]+dh[1]*dh[1])-PMPheat);
-		_assert_(meltrate>0.);
+//		_assert_(meltrate>0.);
+      //  if(meltrate<0.) meltrate = 0.;
 
                   for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*
@@ -632,5 +633,6 @@
 
 		meltrate = 1/latentheat*(G+frictionheat+rho_water*g*conductivity*(dh[0]*dh[0]+dh[1]*dh[1])-PMPheat);
-		_assert_(meltrate>0.);
+//		_assert_(meltrate>0.);
+//		if(meltrate<0.) meltrate = 0.;
 
 		newgap += gauss->weight*Jdet*(gap+dt*(