source: issm/oecreview/Archive/14312-15392/ISSM-14934-14935.diff

../trunk-jpl/src/c/solutions/hydrology_core.cpp

@@ -20 +20 @@
         int        output_frequency,hydrology_model;
         bool       save_results;
         bool       modify_loads=true;
+        bool       isefficientlayer;
         IssmDouble starttime,final_time;
         IssmDouble time,dt;
 
@@ -78 +79 @@
                 }
 
                 else if (hydrology_model==HydrologydcEnum){
+                        femmodel->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
                         if(VerboseSolution()) _pprintLine_("   computing water head");
                         solver_hydro_nonlinear(femmodel);
                         if(save_results && ((i+1)%output_frequency==0 || (i+1)==nsteps)){
                                 if(VerboseSolution()) _pprintLine_("   saving results ");
                                 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");
                                 OutputResultsx(femmodel->elements, femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,femmodel->results);

../trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/UpdateElementsHydrologyDCEfficient.cpp

@@ -49 +49 @@
         iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum);
         iomodel->FetchDataToInput(elements,MaskElementonwaterEnum);
         iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
-        iomodel->FetchDataToInput(elements,SedimentHeadEnum);
+        iomodel->FetchDataToInput(elements,EplHeadEnum);
 
         /*Free data: */
         iomodel->DeleteData(1,MeshElementsEnum);