Changeset 19006

Timestamp:

01/13/15 14:46:37 (10 years ago)

Author:

bdef

Message:

BUG:fixing epl thickness 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*min(0.0,(eplhead[i]-bed[i]))));
+                                EPL_N=gravity*((rho_ice*ice_thickness[i])-(rho_water*max(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*/
@@ -598 +598 @@
                                 /*And proceed to the real thing*/
                                 thickness[i] = old_thickness[i]*(1+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2-2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n)));
+
+                                printf("old: %e, new: %e, grad: %e \n",old_thickness[i],thickness[i],EPLgrad2);
                                 
                                 /*Take care of otherthikening*/