--- trunk/Sources/phylmd/Interface_surf/calcul_fluxs.f 2015/09/29 19:48:59 171 +++ trunk/phylmd/Interface_surf/calcul_fluxs.f 2018/04/19 17:54:55 266 @@ -6,47 +6,50 @@ SUBROUTINE calcul_fluxs(dtime, tsurf, p1lay, cal, beta, coef1lay, ps, & qsurf, radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, petAcoef, & - peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, fluxsens, & + peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, flux_t, & dflux_s, dflux_l) - ! Cette routine calcule les fluxs en h et q à l'interface et une + ! Cette routine calcule les flux en h et q à l'interface et une ! température de surface. - ! L. Fairhead April 2000 + ! L. Fairhead, April 2000 - USE abort_gcm_m, ONLY: abort_gcm - USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats, thermcep - USE indicesol, ONLY: is_ter + USE fcttre, ONLY: foede, foeew use nr_util, only: assert_eq - USE suphec_m, ONLY: rcpd, rd, retv, rkappa, rlstt, rlvtt, rtt + USE suphec_m, ONLY: rcpd, rd, retv, rlstt, rlvtt, rtt USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2 real, intent(IN):: dtime - real, intent(IN):: tsurf(:) ! (knon) temperature de surface - real, intent(IN):: p1lay(:) ! (knon) pression 1er niveau (milieu de couche) + real, intent(IN):: tsurf(:) ! (knon) température de surface + + real, intent(IN):: p1lay(:) ! (knon) + ! pression première couche (milieu de couche) + real, intent(IN):: cal(:) ! (knon) capacité calorifique du sol - real, intent(IN):: beta(:) ! (knon) evap reelle + real, intent(IN):: beta(:) ! (knon) évaporation réelle real, intent(IN):: coef1lay(:) ! (knon) coefficient d'échange real, intent(IN):: ps(:) ! (knon) pression au sol - real, intent(OUT):: qsurf(:) ! (knon) humidite de l'air au dessus du sol - real, intent(IN):: radsol(:) ! (knon) rayonnement net au sol (LW + SW) + real, intent(OUT):: qsurf(:) ! (knon) humidité de l'air au-dessus du sol + + real, intent(IN):: radsol(:) ! (knon) + ! rayonnement net au sol (longwave + shortwave) real, intent(IN):: dif_grnd(:) ! (knon) - ! coefficient diffusion vers le sol profond + ! coefficient de diffusion vers le sol profond real, intent(IN):: t1lay(:), q1lay(:), u1lay(:), v1lay(:) ! (knon) real, intent(IN):: petAcoef(:), peqAcoef(:) ! (knon) - ! coefficients A de la résolution de la couche limite pour t et q + ! coefficients A de la résolution de la couche limite pour T et q real, intent(IN):: petBcoef(:), peqBcoef(:) ! (knon) - ! coeff. B de la resolution de la CL pour t et q + ! coefficients B de la résolution de la couche limite pour t et q real, intent(OUT):: tsurf_new(:) ! (knon) température au sol real, intent(OUT):: evap(:) ! (knon) - real, intent(OUT):: fluxlat(:), fluxsens(:) ! (knon) - ! flux de chaleur latente et sensible + real, intent(OUT):: fluxlat(:), flux_t(:) ! (knon) + ! flux de chaleurs latente et sensible real, intent(OUT):: dflux_s(:), dflux_l(:) ! (knon) ! dérivées des flux de chaleurs sensible et latente par rapport à @@ -55,13 +58,12 @@ ! Local: integer i integer knon ! nombre de points a traiter - real, dimension(size(ps)):: mh, oh, mq, nq, oq - real, dimension(size(ps)):: dq_s_dt, coef - real qsat(size(ps)) ! qsat en kg/kg - real sl(size(ps)) ! chaleur latente d'evaporation ou de sublimation + real, dimension(size(ps)):: mh, oh, mq, nq, oq, dq_s_dt, coef ! (knon) + real qsat(size(ps)) ! (knon) mass fraction + real sl(size(ps)) ! (knon) chaleur latente d'évaporation ou de sublimation logical delta real zcor - real, parameter:: t_grnd = 271.35, t_coup = 273.15 + real, parameter:: t_grnd = 271.35 !--------------------------------------------------------------------- @@ -69,38 +71,26 @@ size(coef1lay), size(ps), size(qsurf), size(radsol), size(dif_grnd), & size(t1lay), size(q1lay), size(u1lay), size(v1lay), size(petAcoef), & size(peqAcoef), size(petBcoef), size(peqBcoef), size(tsurf_new), & - size(evap), size(fluxlat), size(fluxsens), size(dflux_s), & + size(evap), size(fluxlat), size(flux_t), size(dflux_s), & size(dflux_l)/), "calcul_fluxs knon") - ! Traitement humidite du sol + ! Traitement de l'humidité du sol - IF (thermcep) THEN - DO i = 1, knon - delta = rtt >= tsurf(i) - qsat(i) = MIN(0.5, r2es * FOEEW(tsurf(i), delta) / ps(i)) - zcor = 1. / (1. - retv * qsat(i)) - qsat(i) = qsat(i) * zcor - dq_s_dt(i) = RCPD * FOEDE(tsurf(i), delta, merge(R5IES * RLSTT, & - R5LES * RLVTT, delta) / RCPD / (1. + RVTMP2 * q1lay(i)), & - qsat(i), zcor) / RLVTT - ENDDO - ELSE - DO i = 1, knon - IF (tsurf(i) < t_coup) THEN - qsat(i) = qsats(tsurf(i)) / ps(i) - dq_s_dt(i) = RCPD * dqsats(tsurf(i), qsat(i)) / RLVTT - ELSE - qsat(i) = qsatl(tsurf(i)) / ps(i) - dq_s_dt(i) = RCPD * dqsatl(tsurf(i), qsat(i)) / RLVTT - ENDIF - ENDDO - ENDIF + DO i = 1, knon + delta = rtt >= tsurf(i) + qsat(i) = MIN(0.5, r2es * FOEEW(tsurf(i), delta) / ps(i)) + zcor = 1. / (1. - retv * qsat(i)) + qsat(i) = qsat(i) * zcor + dq_s_dt(i) = RCPD * FOEDE(tsurf(i), delta, merge(R5IES * RLSTT, & + R5LES * RLVTT, delta) / RCPD / (1. + RVTMP2 * q1lay(i)), & + qsat(i), zcor) / RLVTT + ENDDO coef = coef1lay * (1. + SQRT(u1lay**2 + v1lay**2)) * p1lay / (RD * t1lay) 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 @@ -116,7 +106,7 @@ evap = - mq - nq * tsurf_new fluxlat = - evap * sl - fluxsens = mh + dflux_s * tsurf_new + flux_t = mh + dflux_s * tsurf_new dflux_l = sl * nq ! Nouvelle valeur de l'humidité au dessus du sol :