Changeset 19003 for issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp

Timestamp:

01/13/15 08:37:22 (10 years ago)

Author:

bdef

Message:

CHG:Corection ine Epl thickness and effective pressure computation

File:

• issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp (modified) (2 diffs)

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

@@ -591 +591 @@
                                 
                                 /*Compute first the effective pressure in the EPL*/
-                                EPL_N=gravity*((rho_ice*ice_thickness[i])-(rho_water*(eplhead[i]-bed[i])));
+                                EPL_N=gravity*((rho_ice*ice_thickness[i])-(rho_water*min(0.0,(eplhead[i]-bed[i]))));
                                 if(EPL_N<0.0)EPL_N=0.0;
                                 /*Get then the square of the gradient of EPL heads*/
-                                EPLgrad2 = (epl_slopeX[i]+epl_slopeY[i])*(epl_slopeX[i]+epl_slopeY[i]);
+                                EPLgrad2 = (epl_slopeX[i]*epl_slopeX[i])+(epl_slopeY[i]*epl_slopeY[i]);
                                 
                                 /*And proceed to the real thing*/
@@ -712 +712 @@
                                         vec_mask->SetValue(basalelement->nodes[i]->Sid(),0.,INS_VAL);
                                         epl_thickness[i]=init_thick;
-                                        //If no more active nodes in the element probably colapse it two
                                 }
                                 else{