Commit c74aca2c authored by Dany Dumont's avatar Dany Dumont

Modification du critere de nage active

parent 4e3703cc
......@@ -614,7 +614,7 @@
if ((par(ci).lt.parcrit).and.(cc(p,ci) .gt. pmin)) then
ws(z,ci) = 0.0
else
ws(z,ci) = -1.0*w_zmax*tanh(bertha*(par(ci)-parcrit))
ws(z,ci) = -1.0*w_zmax*tanh(log((par(ci)+1e-6)/parcrit))
end if
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment