--- trunk/phylmd/Interface_surf/climb_hq_down.f 2018/07/26 16:45:51 298 +++ trunk/phylmd/Interface_surf/climb_hq_down.f 2018/08/02 14:27:11 299 @@ -4,9 +4,9 @@ contains - subroutine climb_hq_down(pkf, cq, dq, ch, dh, paprs, pplay, t, coef, dtime, & - delp, q) + subroutine climb_hq_down(pkf, cq, dq, ch, dh, paprs, pplay, t, coef, delp, q) + use comconst, only: dtphys USE conf_phys_m, ONLY: iflag_pbl USE dimphy, ONLY: klev USE suphec_m, ONLY: rcpd, rd, rg, rkappa @@ -26,8 +26,6 @@ ! Le coefficient d'echange (m**2 / s) multiplie par le cisaillement ! du vent (dV / dz) - REAL, intent(in):: dtime ! intervalle du temps (s) - REAL, intent(in):: delp(:, :) ! (knon, klev) ! epaisseur de couche en pression (Pa) @@ -52,7 +50,7 @@ ! Convertir les coefficients en variables convenables au calcul: forall (k = 2:klev) zx_coef(:, k) = coef(:, k) & / (pplay(:, k - 1) - pplay(:, k)) & - * (paprs(:, k) * 2 / (t(:, k) + t(:, k - 1)) / RD)**2 * dtime * RG**2 + * (paprs(:, k) * 2 / (t(:, k) + t(:, k - 1)) / RD)**2 * dtphys * RG**2 ! Preparer les flux lies aux contre-gardients