Index: /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 27362)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 27363)
@@ -418,4 +418,7 @@
 		b_input->GetInputValue(&b,gauss);
 		H_input->GetInputValue(&H,gauss);
+
+		/*Hard code B*/
+		B = Cuffey(273.15-2);
 
 		/*Get norm of gradient of hydraulic potential and make sure it is >0*/
@@ -523,4 +526,7 @@
 		b_input->GetInputValue(&b,gauss);
 		H_input->GetInputValue(&H,gauss);
+		
+		/*Hard code B*/
+		B = Cuffey(273.15-2);
 
 		/*Get basal velocity*/
@@ -674,5 +680,5 @@
 		element->GetInputListOnNodes(&mask[0],MaskOceanLevelsetEnum);
 		element->GetInputListOnNodes(&bed[0],BaseEnum);
-		element->GetInputListOnNodes(&thickness[0],BaseEnum);
+		element->GetInputListOnNodes(&thickness[0],ThicknessEnum);
 		element->GetInputListOnNodes(&ls_active[0],HydrologyMaskNodeActivationEnum);
 
@@ -779,4 +785,7 @@
 		iceLS_input->GetInputValue(&iceLS,gauss);
 
+		/*Hard code B*/
+		B = Cuffey(273.15-2);
+
 		/*Set sheet thickness to zero if floating or no ice*/
 		if(oceanLS<0. || iceLS>0.){
@@ -812,9 +821,4 @@
 	}
 	}
-
-	/*Force floating ice to have zero sheet thickness*/
-	/*if(!element->IsAllGrounded() || !element->IsIceInElement()){
-		for(int iv=0;iv<numvertices;iv++) h_new[iv] = 0.;
-	}*/
 
 	element->AddInput(HydrologySheetThicknessEnum,h_new,P1Enum);
@@ -878,5 +882,5 @@
 		p_i = rho_ice*g*H;
 
-		/*Copmute overburden potential*/
+		/*Compute overburden potential*/
 		phi_0 = phi_m + p_i;