--- trunk/Sources/phylmd/Interface_surf/calcul_fluxs.f 2017/03/22 13:40:27 214 +++ trunk/phylmd/Interface_surf/calcul_fluxs.f 2018/04/19 17:54:55 266 @@ -14,7 +14,7 @@ ! L. Fairhead, April 2000 - USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats + USE fcttre, ONLY: foede, foeew use nr_util, only: assert_eq USE suphec_m, ONLY: rcpd, rd, retv, rlstt, rlvtt, rtt USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2 @@ -90,7 +90,7 @@ sl = merge(RLSTT, RLVTT, tsurf < RTT) ! Q - oq = 1. - (beta * coef * peqBcoef * dtime) + oq = 1. - beta * coef * peqBcoef * dtime mq = beta * coef * (peqAcoef - qsat + dq_s_dt * tsurf) / oq nq = beta * coef * (- 1. * dq_s_dt) / oq