Index: /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 25696)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 25697)
@@ -129,4 +129,5 @@
 	iomodel->FetchDataToInput(inputs,elements,"md.geometry.base",BaseEnum);
 	iomodel->FetchDataToInput(inputs,elements,"md.geometry.bed",BedEnum);
+	iomodel->FetchDataToInput(inputs,elements,"md.solidearth.sealevel",SealevelEnum,0); /*Needed for friction*/
 	iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.geothermalflux",BasalforcingsGeothermalfluxEnum);
 	iomodel->FetchDataToInput(inputs,elements,"md.basalforcings.groundedice_melting_rate",BasalforcingsGroundediceMeltingRateEnum);
Index: /issm/trunk-jpl/src/c/classes/Loads/Moulin.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Moulin.cpp	(revision 25696)
+++ /issm/trunk-jpl/src/c/classes/Loads/Moulin.cpp	(revision 25697)
@@ -6,5 +6,5 @@
 /*{{{*/
 #ifdef HAVE_CONFIG_H
-	#include <config.h>
+#include <config.h>
 #else
 #error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
@@ -367,5 +367,7 @@
 
 	/*Load vector*/
-	Ke->values[0] = +Am/(rho_water*g)/dt;
+	if(dt>0){
+		Ke->values[0] = +Am/(rho_water*g)/dt;
+	}
 
 	/*Clean up and return*/
@@ -393,5 +395,8 @@
 	element->GetInputValue(&moulin_load,node,HydrologyMoulinInputEnum);
 
-	pe->values[0] = moulin_load +Am/(rho_water*g) * phi_old/dt;
+	pe->values[0] = moulin_load;
+	if(dt>0.){
+		pe->values[0] += Am/(rho_water*g) * phi_old/dt;
+	}
 
 	/*Clean up and return*/
@@ -447,5 +452,5 @@
 	/*Clean up and return*/
 	return pe;
- }
+}
 /*}}}*/
 ElementVector* Moulin::CreatePVectorHydrologyDCEfficient(void){/*{{{*/