Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/UpdateElementsHydrologyDCEfficient.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/UpdateElementsHydrologyDCEfficient.cpp	(revision 14934)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/UpdateElementsHydrologyDCEfficient.cpp	(revision 14935)
@@ -50,5 +50,5 @@
 	iomodel->FetchDataToInput(elements,MaskElementonwaterEnum);
 	iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
-	iomodel->FetchDataToInput(elements,SedimentHeadEnum);
+	iomodel->FetchDataToInput(elements,EplHeadEnum);
 
 	/*Free data: */
Index: /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp	(revision 14934)
+++ /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp	(revision 14935)
@@ -21,4 +21,5 @@
 	bool       save_results;
 	bool       modify_loads=true;
+	bool       isefficientlayer;
 	IssmDouble starttime,final_time;
 	IssmDouble time,dt;
@@ -79,4 +80,5 @@
 
 		else if (hydrology_model==HydrologydcEnum){
+			femmodel->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
 			if(VerboseSolution()) _pprintLine_("   computing water head");
 			solver_hydro_nonlinear(femmodel);
@@ -85,5 +87,7 @@
 				InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SedimentHeadEnum);
 				InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SedimentHeadResidualEnum);
-
+				if(isefficientlayer){
+					InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,EplHeadEnum);
+				}
 				/*unload results*/
 				if(VerboseSolution()) _pprintLine_("   saving temporary results");