Index: ../trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp =================================================================== --- ../trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp (revision 26496) +++ ../trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp (revision 26497) @@ -403,7 +403,7 @@ -beta*sqrt(vx*vx+vy*vy) +ieb +storage*head_old/dt - +meltrate*gap/rho_ice*(dgapxx+dgapyy) + +pow(1,-5)*(dgapxx+dgapyy) )*basis[i]; @@ -659,8 +659,8 @@ meltrate/rho_ice -A*pow(fabs(pressure_ice-pressure_water),n-1)*(pressure_ice-pressure_water)*lc +beta*sqrt(vx*vx+vy*vy) - +(meltrate*gap/rho_ice*(dgapxx+dgapyy)) - )); + +pow(1,-5)*(dgapxx+dgapyy)) + ); /*TEST with linear creep rate*/