Changeset 18183

Timestamp:

06/25/14 10:44:52 (11 years ago)

Author:

bdef

Message:

BUG: correction to epl gradient and modif in transfer scheme

Location:

issm/trunk-jpl/src/c/analyses

Files:

• HydrologyDCEfficientAnalysis.cpp (modified) (1 diff)
• HydrologyDCInefficientAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp ¶

@@ -559 +559 @@
                                 if(EPL_N<0.0)EPL_N=0.0;
                                 /*Get then the square of the gradient of EPL heads*/
-                                EPLgrad2 = (epl_slopeX[i]*epl_slopeX[i]+epl_slopeY[i]*epl_slopeY[i]);
+                                EPLgrad2 = (epl_slopeX[i]+epl_slopeY[i])*(epl_slopeX[i]+epl_slopeY[i]);
                                 
                                 /*And proceed to the real thing*/

issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp ¶

@@ -604 +604 @@
                 sed_trans_input->GetInputValue(&sediment_transmitivity,gauss);
                 element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
+                
                 transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
+                
                 /* if(epl_head>sed_head){ */
                 /*      h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input); */
@@ -619 +621 @@
                 /*      //transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage); */
                 /* } */
+                
                 break;
         default: