Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 23376)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 23377)
@@ -349,5 +349,5 @@
 	}
 	if(smb_model==SMBgradientscomponentsEnum){
-		surface_runoff_input = basalelement->GetInput(SmbRunoffEnum);
+		surface_runoff_input = basalelement->GetInput(SmbRunoffEnum); _assert_(surface_runoff_input);
 	}
 
@@ -363,5 +363,6 @@
 		/*Loading term*/
 		basal_melt_input->GetInputValue(&water_load,gauss);
-		surface_runoff_input->GetInputValue(&runoff_value,gauss);
+		if(surface_runoff_input) surface_runoff_input->GetInputValue(&runoff_value,gauss);
+		else                     runoff_value = 0.;
 		scalar = Jdet*gauss->weight*(water_load+runoff_value);
 		//scalar = Jdet*gauss->weight*(water_load)/epl_transmitivity;
Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 23376)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 23377)
@@ -407,5 +407,6 @@
 			if(!active_element){
 				basal_melt_input->GetInputValue(&water_load,gauss);
-				surface_runoff_input->GetInputValue(&runoff_value,gauss);
+				if(surface_runoff_input) surface_runoff_input->GetInputValue(&runoff_value,gauss);
+				else                     runoff_value = 0.;
 				scalar = Jdet*gauss->weight*(water_load+runoff_value);
 				//scalar = Jdet*gauss->weight*(water_load)/sediment_transmitivity;