Index: /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 26794)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 26795)
@@ -472,5 +472,5 @@
 
 	/*Skip if water or ice shelf element*/
-	if(element->IsAllFloating()) return;
+	//if(element->IsAllFloating()) return;
 
 	/*Intermediaries */
@@ -543,4 +543,9 @@
 		if(h_new[iv]<AEPS) h_new[iv] = AEPS;
 	}
+	
+	/*Force floating ice to have zero sheet thickness*/
+	if(element->IsAllFloating()){
+				for(int iv=0;iv<numvertices;iv++) h_new[iv] = 0.;
+				}
 
 	element->AddInput(HydrologySheetThicknessEnum,h_new,P1Enum);
@@ -561,5 +566,5 @@
 
 	/*Skip if water or ice shelf element*/
-	if(element->IsAllFloating()) return;
+	//if(element->IsAllFloating()) return;
 
 	/*Intermediary*/
@@ -606,4 +611,9 @@
 	element->AddInput(EffectivePressureEnum,N,element->FiniteElement());
 
+   /*Force floating ice to have zero effective pressure*/
+   if(element->IsAllFloating()){
+            for(int iv=0;iv<numnodes;iv++) N[iv] = 0.;
+            }
+	
 	/*Clean up and return*/
 	delete gauss;