/[lmdze]/trunk/Sources/phylmd/physiq.f

Diff of /trunk/Sources/phylmd/physiq.f

Parent Directory | Revision Log | View Patch Patch

-revision 205 by guez,
Tue Jun 21 15:16:03 2016 UTC
+revision 224 by guez,
Fri Apr 28 13:40:59 2017 UTC
 Line 20 
 contains
      use calltherm_m, only: calltherm
      USE clesphys, ONLY: cdhmax, cdmmax, ecrit_ins, ksta, ksta_ter, ok_kzmin, &
           ok_instan
-     USE clesphys2, ONLY: cycle_diurne, conv_emanuel, nbapp_rad, new_oliq, &
+     USE clesphys2, ONLY: conv_emanuel, nbapp_rad, new_oliq, ok_orodr, ok_orolf
-          ok_orodr, ok_orolf
      USE clmain_m, ONLY: clmain
      use clouds_gno_m, only: clouds_gno
      use comconst, only: dtphys
      USE comgeomphy, ONLY: airephy
      USE concvl_m, ONLY: concvl
-     USE conf_gcm_m, ONLY: offline, day_step, iphysiq, lmt_pas
+     USE conf_gcm_m, ONLY: lmt_pas
      USE conf_phys_m, ONLY: conf_phys
      use conflx_m, only: conflx
      USE ctherm, ONLY: iflag_thermals, nsplit_thermals
-Line 37 
 contains
+Line 36 
 contains
      USE dimsoil, ONLY: nsoilmx
      use drag_noro_m, only: drag_noro
      use dynetat0_m, only: day_ref, annee_ref
-     USE fcttre, ONLY: foeew, qsatl, qsats, thermcep
+     USE fcttre, ONLY: foeew
      use fisrtilp_m, only: fisrtilp
      USE hgardfou_m, ONLY: hgardfou
      USE histsync_m, ONLY: histsync
-Line 54 
 contains
+Line 53 
 contains
      USE phyetat0_m, ONLY: phyetat0, rlat, rlon
      USE phyredem_m, ONLY: phyredem
      USE phyredem0_m, ONLY: phyredem0
-     USE phystokenc_m, ONLY: phystokenc
      USE phytrac_m, ONLY: phytrac
      use radlwsw_m, only: radlwsw
      use yoegwd, only: sugwd
-     USE suphec_m, ONLY: rcpd, retv, rg, rlvtt, romega, rsigma, rtt
+     USE suphec_m, ONLY: rcpd, retv, rg, rlvtt, romega, rsigma, rtt, rmo3, md
      use time_phylmdz, only: itap, increment_itap
      use transp_m, only: transp
      use transp_lay_m, only: transp_lay
-Line 146 
 contains
+Line 144 
 contains
      ! "physiq".
      REAL, save:: radsol(klon) ! bilan radiatif au sol calcule par code radiatif
      REAL, save:: ftsol(klon, nbsrf) ! skin temperature of surface fraction
      REAL, save:: ftsoil(klon, nsoilmx, nbsrf)
      ! soil temperature of surface fraction
      REAL, save:: fevap(klon, nbsrf) ! evaporation
-     REAL, save:: fluxlat(klon, nbsrf)
+     REAL fluxlat(klon, nbsrf)
      REAL, save:: fqsurf(klon, nbsrf)
      ! humidite de l'air au contact de la surface
-     REAL, save:: qsol(klon)
+     REAL, save:: qsol(klon) ! column-density of water in soil, in kg m-2
-     ! column-density of water in soil, in kg m-2
+     REAL, save:: fsnow(klon, nbsrf) ! \'epaisseur neigeuse
-     REAL, save:: fsnow(klon, nbsrf) ! epaisseur neigeuse
      REAL, save:: falbe(klon, nbsrf) ! albedo visible par type de surface
      ! Param\`etres de l'orographie \`a l'\'echelle sous-maille (OESM) :
-Line 208 
 contains
+Line 203 
 contains
      REAL, save:: pfrac_1nucl(klon, llm)
      ! Produits des coefs lessi nucl (alpha = 1)
-     REAL frac_impa(klon, llm) ! fractions d'aerosols lessivees (impaction)
+     REAL frac_impa(klon, llm) ! fraction d'a\'erosols lessiv\'es (impaction)
      REAL frac_nucl(klon, llm) ! idem (nucleation)
      REAL, save:: rain_fall(klon)
-Line 219 
 contains
+Line 214 
 contains
      REAL rain_tiedtke(klon), snow_tiedtke(klon)
-     REAL evap(klon), devap(klon) ! evaporation and its derivative
+     REAL evap(klon) ! flux d'\'evaporation au sol
-     REAL sens(klon), dsens(klon) ! chaleur sensible et sa derivee
+     real devap(klon) ! derivative of the evaporation flux at the surface
-     REAL, save:: dlw(klon) ! derivee infra rouge
+     REAL sens(klon) ! flux de chaleur sensible au sol
+     real dsens(klon) ! derivee du flux de chaleur sensible au sol
+     REAL, save:: dlw(klon) ! derivative of infra-red flux
      REAL bils(klon) ! bilan de chaleur au sol
-     REAL, save:: fder(klon) ! Derive de flux (sensible et latente)
+     REAL fder(klon) ! Derive de flux (sensible et latente)
      REAL ve(klon) ! integr. verticale du transport meri. de l'energie
      REAL vq(klon) ! integr. verticale du transport meri. de l'eau
      REAL ue(klon) ! integr. verticale du transport zonal de l'energie
-Line 236 
 contains
+Line 233 
 contains
      INTEGER julien
      REAL, save:: pctsrf(klon, nbsrf) ! percentage of surface
-     REAL, save:: albsol(klon) ! albedo du sol total visible
+     REAL, save:: albsol(klon) ! albedo du sol total, visible, moyen par maille
      REAL, SAVE:: wo(klon, llm) ! column density of ozone in a cell, in kDU
+     real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
      real, save:: clwcon(klon, llm), rnebcon(klon, llm)
      real, save:: clwcon0(klon, llm), rnebcon0(klon, llm)
-Line 250 
 contains
+Line 248 
 contains
      REAL cldtau(klon, llm) ! epaisseur optique
      REAL cldemi(klon, llm) ! emissivite infrarouge
-     REAL fluxq(klon, llm, nbsrf) ! flux turbulent d'humidite
+     REAL flux_q(klon, nbsrf) ! flux turbulent d'humidite à la surface
-     REAL fluxt(klon, llm, nbsrf) ! flux turbulent de chaleur
+     REAL flux_t(klon, nbsrf) ! flux turbulent de chaleur à la surface
-     REAL fluxu(klon, llm, nbsrf) ! flux turbulent de vitesse u
+     REAL flux_u(klon, nbsrf) ! flux turbulent de vitesse u à la surface
-     REAL fluxv(klon, llm, nbsrf) ! flux turbulent de vitesse v
+     REAL flux_v(klon, nbsrf) ! flux turbulent de vitesse v à la surface
-     REAL zxfluxt(klon, llm)
-     REAL zxfluxq(klon, llm)
-     REAL zxfluxu(klon, llm)
-     REAL zxfluxv(klon, llm)
      ! Le rayonnement n'est pas calcul\'e tous les pas, il faut donc que
      ! les variables soient r\'emanentes.
-Line 280 
 contains
+Line 273 
 contains
      REAL cldl(klon), cldm(klon), cldh(klon) ! nuages bas, moyen et haut
      REAL cldt(klon), cldq(klon) ! nuage total, eau liquide integree
-     REAL zxqsurf(klon), zxsnow(klon), zxfluxlat(klon)
+     REAL zxfluxlat(klon)
      REAL dist, mu0(klon), fract(klon)
      real longi
      REAL z_avant(klon), z_apres(klon), z_factor(klon)
-Line 289 
 contains
+Line 281 
 contains
      REAL zx_t, zx_qs, zcor
      real zqsat(klon, llm)
      INTEGER i, k, iq, nsrf
-     REAL, PARAMETER:: t_coup = 234.
      REAL zphi(klon, llm)
      ! cf. Anne Mathieu, variables pour la couche limite atmosphérique (hbtm)
-Line 299 
 contains
+Line 290 
 contains
      REAL, SAVE:: capCL(klon, nbsrf) ! CAPE de couche limite
      REAL, SAVE:: oliqCL(klon, nbsrf) ! eau_liqu integree de couche limite
      REAL, SAVE:: cteiCL(klon, nbsrf) ! cloud top instab. crit. couche limite
-     REAL, SAVE:: pblt(klon, nbsrf) ! T a la Hauteur de couche limite
+     REAL, SAVE:: pblt(klon, nbsrf) ! T \`a la hauteur de couche limite
      REAL, SAVE:: therm(klon, nbsrf)
      REAL, SAVE:: trmb1(klon, nbsrf) ! deep_cape
      REAL, SAVE:: trmb2(klon, nbsrf) ! inhibition
-Line 314 
 contains
+Line 305 
 contains
      REAL upwd(klon, llm) ! saturated updraft mass flux
      REAL dnwd(klon, llm) ! saturated downdraft mass flux
-     REAL dnwd0(klon, llm) ! unsaturated downdraft mass flux
      REAL, save:: cape(klon)
      INTEGER iflagctrl(klon) ! flag fonctionnement de convect
-Line 347 
 contains
+Line 337 
 contains
      real rain_lsc(klon)
      REAL, save:: snow_con(klon) ! neige (mm / s)
      real snow_lsc(klon)
-     REAL d_ts(klon, nbsrf)
+     REAL d_ts(klon, nbsrf) ! variation of ftsol
      REAL d_u_vdf(klon, llm), d_v_vdf(klon, llm)
      REAL d_t_vdf(klon, llm), d_q_vdf(klon, llm)
-Line 390 
 contains
+Line 380 
 contains
      REAL uq_lay(klon, llm) ! transport zonal de l'eau a chaque niveau vert.
      real date0
-     REAL ztsol(klon)
+     REAL tsol(klon)
      REAL d_t_ec(klon, llm)
-     ! tendance due \`a la conversion Ec en énergie thermique
+     ! tendance due \`a la conversion d'\'energie cin\'etique en
+     ! énergie thermique
-     REAL ZRCPD
      REAL, save:: t2m(klon, nbsrf), q2m(klon, nbsrf)
      ! temperature and humidity at 2 m
      REAL, save:: u10m(klon, nbsrf), v10m(klon, nbsrf) ! vents a 10 m
-     REAL zt2m(klon), zq2m(klon) ! temp., hum. 2 m moyenne s/ 1 maille
+     REAL zt2m(klon), zq2m(klon) ! température, humidité 2 m moyenne sur 1 maille
-     REAL zu10m(klon), zv10m(klon) ! vents a 10 m moyennes s/1 maille
+     REAL zu10m(klon), zv10m(klon) ! vents a 10 m moyennes sur 1 maille
      ! Aerosol effects:
-     REAL sulfate(klon, llm) ! SO4 aerosol concentration (micro g / m3)
-     REAL, save:: sulfate_pi(klon, llm)
-     ! SO4 aerosol concentration, in \mu g / m3, pre-industrial value
-     REAL cldtaupi(klon, llm)
-     ! cloud optical thickness for pre-industrial aerosols
-     REAL re(klon, llm) ! Cloud droplet effective radius
-     REAL fl(klon, llm) ! denominator of re
-     ! Aerosol optical properties
-     REAL, save:: tau_ae(klon, llm, 2), piz_ae(klon, llm, 2)
-     REAL, save:: cg_ae(klon, llm, 2)
      REAL, save:: topswad(klon), solswad(klon) ! aerosol direct effect
-     REAL, save:: topswai(klon), solswai(klon) ! aerosol indirect effect
      LOGICAL:: ok_ade = .false. ! apply aerosol direct effect
-     LOGICAL:: ok_aie = .false. ! apply aerosol indirect effect
      REAL:: bl95_b0 = 2., bl95_b1 = 0.2
      ! Parameters in equation (D) of Boucher and Lohmann (1995, Tellus
-Line 438 
 contains
+Line 409 
 contains
      integer, save:: ncid_startphy
      namelist /physiq_nml/ fact_cldcon, facttemps, ok_newmicro, iflag_cldcon, &
-          ratqsbas, ratqshaut, ok_ade, ok_aie, bl95_b0, bl95_b1, &
+          ratqsbas, ratqshaut, ok_ade, bl95_b0, bl95_b1, iflag_thermals, &
-          iflag_thermals, nsplit_thermals
+          nsplit_thermals
      !----------------------------------------------------------------
-Line 454 
 contains
+Line 425 
 contains
         q2m = 0.
         ffonte = 0.
         fqcalving = 0.
-        piz_ae = 0.
-        tau_ae = 0.
-        cg_ae = 0.
         rain_con = 0.
         snow_con = 0.
-        topswai = 0.
-        topswad = 0.
-        solswai = 0.
-        solswad = 0.
         d_u_con = 0.
         d_v_con = 0.
         rnebcon0 = 0.
         clwcon0 = 0.
         rnebcon = 0.
         clwcon = 0.
         pblh =0. ! Hauteur de couche limite
         plcl =0. ! Niveau de condensation de la CLA
         capCL =0. ! CAPE de couche limite
         oliqCL =0. ! eau_liqu integree de couche limite
         cteiCL =0. ! cloud top instab. crit. couche limite
-        pblt =0. ! T a la Hauteur de couche limite
+        pblt =0.
         therm =0.
         trmb1 =0. ! deep_cape
         trmb2 =0. ! inhibition
-Line 522 
 contains
+Line 484 
 contains
         ! Initialisation des sorties
-        call ini_histins(dtphys)
+        call ini_histins(dtphys, ok_newmicro)
         CALL ymds2ju(annee_ref, 1, day_ref, 0., date0)
         ! Positionner date0 pour initialisation de ORCHIDEE
         print *, 'physiq date0: ', date0
-Line 538 
 contains
+Line 500 
 contains
      ql_seri = qx(:, :, iliq)
      tr_seri = qx(:, :, 3:nqmx)
-     ztsol = sum(ftsol * pctsrf, dim = 2)
+     tsol = sum(ftsol * pctsrf, dim = 2)
      ! Diagnostic de la tendance dynamique :
      IF (ancien_ok) THEN
-Line 574 
 contains
+Line 536 
 contains
      forall (k = 1: llm) zmasse(:, k) = (paprs(:, k) - paprs(:, k + 1)) / rg
-     ! Prescrire l'ozone :
-     wo = ozonecm(REAL(julien), paprs)
      ! \'Evaporation de l'eau liquide nuageuse :
      DO k = 1, llm
         DO i = 1, klon
-Line 595 
 contains
+Line 554 
 contains
      ! la surface.
      CALL orbite(REAL(julien), longi, dist)
-     IF (cycle_diurne) THEN
+     CALL zenang(longi, time, dtphys * radpas, mu0, fract)
-        CALL zenang(longi, time, dtphys * radpas, mu0, fract)
-     ELSE
-        mu0 = - 999.999
-     ENDIF
-     ! Calcul de l'abedo moyen par maille
      albsol = sum(falbe * pctsrf, dim = 2)
      ! R\'epartition sous maille des flux longwave et shortwave
      ! R\'epartition du longwave par sous-surface lin\'earis\'ee
      forall (nsrf = 1: nbsrf)
-        fsollw(:, nsrf) = sollw + 4. * RSIGMA * ztsol**3 &
+        fsollw(:, nsrf) = sollw + 4. * RSIGMA * tsol**3 &
-             * (ztsol - ftsol(:, nsrf))
+             * (tsol - ftsol(:, nsrf))
         fsolsw(:, nsrf) = solsw * (1. - falbe(:, nsrf)) / (1. - albsol)
      END forall
-     fder = dlw
-     ! Couche limite:
      CALL clmain(dtphys, pctsrf, t_seri, q_seri, u_seri, v_seri, julien, mu0, &
           ftsol, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ftsoil, qsol, &
           paprs, play, fsnow, fqsurf, fevap, falbe, fluxlat, rain_fall, &
-          snow_fall, fsolsw, fsollw, fder, rlat, frugs, agesno, rugoro, &
+          snow_fall, fsolsw, fsollw, frugs, agesno, rugoro, d_t_vdf, d_q_vdf, &
-          d_t_vdf, d_q_vdf, d_u_vdf, d_v_vdf, d_ts, fluxt, fluxq, fluxu, &
+          d_u_vdf, d_v_vdf, d_ts, flux_t, flux_q, flux_u, flux_v, cdragh, &
-          fluxv, cdragh, cdragm, q2, dsens, devap, ycoefh, yu1, yv1, t2m, q2m, &
+          cdragm, q2, dsens, devap, ycoefh, yu1, yv1, t2m, q2m, u10m, v10m, &
-          u10m, v10m, pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, &
+          pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, trmb3, plcl, &
-          trmb3, plcl, fqcalving, ffonte, run_off_lic_0)
+          fqcalving, ffonte, run_off_lic_0)
      ! Incr\'ementation des flux
-     zxfluxt = 0.
+     sens = - sum(flux_t * pctsrf, dim = 2)
-     zxfluxq = 0.
+     evap = - sum(flux_q * pctsrf, dim = 2)
-     zxfluxu = 0.
+     fder = dlw + dsens + devap
-     zxfluxv = 0.
-     DO nsrf = 1, nbsrf
-        DO k = 1, llm
-           DO i = 1, klon
-              zxfluxt(i, k) = zxfluxt(i, k) + fluxt(i, k, nsrf) * pctsrf(i, nsrf)
-              zxfluxq(i, k) = zxfluxq(i, k) + fluxq(i, k, nsrf) * pctsrf(i, nsrf)
-              zxfluxu(i, k) = zxfluxu(i, k) + fluxu(i, k, nsrf) * pctsrf(i, nsrf)
-              zxfluxv(i, k) = zxfluxv(i, k) + fluxv(i, k, nsrf) * pctsrf(i, nsrf)
-           END DO
-        END DO
-     END DO
-     DO i = 1, klon
-        sens(i) = - zxfluxt(i, 1) ! flux de chaleur sensible au sol
-        evap(i) = - zxfluxq(i, 1) ! flux d'\'evaporation au sol
-        fder(i) = dlw(i) + dsens(i) + devap(i)
-     ENDDO
      DO k = 1, llm
         DO i = 1, klon
-Line 659 
 contains
+Line 592 
 contains
      ! Update surface temperature:
-     DO i = 1, klon
-        zxfluxlat(i) = 0.
-        zt2m(i) = 0.
-        zq2m(i) = 0.
-        zu10m(i) = 0.
-        zv10m(i) = 0.
-        zxffonte(i) = 0.
-        zxfqcalving(i) = 0.
-        s_pblh(i) = 0.
-        s_lcl(i) = 0.
-        s_capCL(i) = 0.
-        s_oliqCL(i) = 0.
-        s_cteiCL(i) = 0.
-        s_pblT(i) = 0.
-        s_therm(i) = 0.
-        s_trmb1(i) = 0.
-        s_trmb2(i) = 0.
-        s_trmb3(i) = 0.
-     ENDDO
      call assert(abs(sum(pctsrf, dim = 2) - 1.) <= EPSFRA, 'physiq: pctsrf')
      ftsol = ftsol + d_ts
-     ztsol = sum(ftsol * pctsrf, dim = 2)
+     tsol = sum(ftsol * pctsrf, dim = 2)
-     DO nsrf = 1, nbsrf
+     zxfluxlat = sum(fluxlat * pctsrf, dim = 2)
-        DO i = 1, klon
+     zt2m = sum(t2m * pctsrf, dim = 2)
-           zxfluxlat(i) = zxfluxlat(i) + fluxlat(i, nsrf) * pctsrf(i, nsrf)
+     zq2m = sum(q2m * pctsrf, dim = 2)
+     zu10m = sum(u10m * pctsrf, dim = 2)
-           zt2m(i) = zt2m(i) + t2m(i, nsrf) * pctsrf(i, nsrf)
+     zv10m = sum(v10m * pctsrf, dim = 2)
-           zq2m(i) = zq2m(i) + q2m(i, nsrf) * pctsrf(i, nsrf)
+     zxffonte = sum(ffonte * pctsrf, dim = 2)
-           zu10m(i) = zu10m(i) + u10m(i, nsrf) * pctsrf(i, nsrf)
+     zxfqcalving = sum(fqcalving * pctsrf, dim = 2)
-           zv10m(i) = zv10m(i) + v10m(i, nsrf) * pctsrf(i, nsrf)
+     s_pblh = sum(pblh * pctsrf, dim = 2)
-           zxffonte(i) = zxffonte(i) + ffonte(i, nsrf) * pctsrf(i, nsrf)
+     s_lcl = sum(plcl * pctsrf, dim = 2)
-           zxfqcalving(i) = zxfqcalving(i) + &
+     s_capCL = sum(capCL * pctsrf, dim = 2)
-                fqcalving(i, nsrf) * pctsrf(i, nsrf)
+     s_oliqCL = sum(oliqCL * pctsrf, dim = 2)
-           s_pblh(i) = s_pblh(i) + pblh(i, nsrf) * pctsrf(i, nsrf)
+     s_cteiCL = sum(cteiCL * pctsrf, dim = 2)
-           s_lcl(i) = s_lcl(i) + plcl(i, nsrf) * pctsrf(i, nsrf)
+     s_pblT = sum(pblT * pctsrf, dim = 2)
-           s_capCL(i) = s_capCL(i) + capCL(i, nsrf) * pctsrf(i, nsrf)
+     s_therm = sum(therm * pctsrf, dim = 2)
-           s_oliqCL(i) = s_oliqCL(i) + oliqCL(i, nsrf) * pctsrf(i, nsrf)
+     s_trmb1 = sum(trmb1 * pctsrf, dim = 2)
-           s_cteiCL(i) = s_cteiCL(i) + cteiCL(i, nsrf) * pctsrf(i, nsrf)
+     s_trmb2 = sum(trmb2 * pctsrf, dim = 2)
-           s_pblT(i) = s_pblT(i) + pblT(i, nsrf) * pctsrf(i, nsrf)
+     s_trmb3 = sum(trmb3 * pctsrf, dim = 2)
-           s_therm(i) = s_therm(i) + therm(i, nsrf) * pctsrf(i, nsrf)
-           s_trmb1(i) = s_trmb1(i) + trmb1(i, nsrf) * pctsrf(i, nsrf)
-           s_trmb2(i) = s_trmb2(i) + trmb2(i, nsrf) * pctsrf(i, nsrf)
-           s_trmb3(i) = s_trmb3(i) + trmb3(i, nsrf) * pctsrf(i, nsrf)
-        ENDDO
-     ENDDO
      ! Si une sous-fraction n'existe pas, elle prend la valeur moyenne :
      DO nsrf = 1, nbsrf
         DO i = 1, klon
            IF (pctsrf(i, nsrf) < epsfra) then
-              ftsol(i, nsrf) = ztsol(i)
+              ftsol(i, nsrf) = tsol(i)
               t2m(i, nsrf) = zt2m(i)
               q2m(i, nsrf) = zq2m(i)
               u10m(i, nsrf) = zu10m(i)
-Line 734 
 contains
+Line 638 
 contains
         ENDDO
      ENDDO
-     ! Calculer la dérive du flux infrarouge
+     dlw = - 4. * RSIGMA * tsol**3
-     DO i = 1, klon
-        dlw(i) = - 4. * RSIGMA * ztsol(i)**3
-     ENDDO
      ! Appeler la convection
      if (conv_emanuel) then
         CALL concvl(paprs, play, t_seri, q_seri, u_seri, v_seri, sig1, w01, &
              d_t_con, d_q_con, d_u_con, d_v_con, rain_con, ibas_con, itop_con, &
-             upwd, dnwd, dnwd0, Ma, cape, iflagctrl, qcondc, pmflxr, da, phi, mp)
+             upwd, dnwd, Ma, cape, iflagctrl, qcondc, pmflxr, da, phi, mp)
         snow_con = 0.
         clwcon0 = qcondc
         mfu = upwd + dnwd
-        IF (thermcep) THEN
+        zqsat = MIN(0.5, r2es * FOEEW(t_seri, rtt >= t_seri) / play)
-           zqsat = MIN(0.5, r2es * FOEEW(t_seri, rtt >= t_seri) / play)
+        zqsat = zqsat / (1. - retv * zqsat)
-           zqsat = zqsat / (1. - retv * zqsat)
-        ELSE
-           zqsat = merge(qsats(t_seri), qsatl(t_seri), t_seri < t_coup) / play
-        ENDIF
         ! Properties of convective clouds
         clwcon0 = fact_cldcon * clwcon0
-Line 773 
 contains
+Line 669 
 contains
         conv_t = d_t_dyn + d_t_vdf / dtphys
         z_avant = sum((q_seri + ql_seri) * zmasse, dim=2)
         CALL conflx(dtphys, paprs, play, t_seri(:, llm:1:- 1), &
-             q_seri(:, llm:1:- 1), conv_t, conv_q, zxfluxq(:, 1), omega, &
+             q_seri(:, llm:1:- 1), conv_t, conv_q, - evap, omega, &
              d_t_con, d_q_con, rain_con, snow_con, mfu(:, llm:1:- 1), &
              mfd(:, llm:1:- 1), pen_u, pde_u, pen_d, pde_d, kcbot, kctop, &
              kdtop, pmflxr, pmflxs)
-Line 955 
 contains
+Line 851 
 contains
      DO k = 1, llm
         DO i = 1, klon
            zx_t = t_seri(i, k)
-           IF (thermcep) THEN
+           zx_qs = r2es * FOEEW(zx_t, rtt >= zx_t) / play(i, k)
-              zx_qs = r2es * FOEEW(zx_t, rtt >= zx_t) / play(i, k)
+           zx_qs = MIN(0.5, zx_qs)
-              zx_qs = MIN(0.5, zx_qs)
+           zcor = 1. / (1. - retv * zx_qs)
-              zcor = 1. / (1. - retv * zx_qs)
+           zx_qs = zx_qs * zcor
-              zx_qs = zx_qs * zcor
-           ELSE
-              IF (zx_t < t_coup) THEN
-                 zx_qs = qsats(zx_t) / play(i, k)
-              ELSE
-                 zx_qs = qsatl(zx_t) / play(i, k)
-              ENDIF
-           ENDIF
            zx_rh(i, k) = q_seri(i, k) / zx_qs
            zqsat(i, k) = zx_qs
         ENDDO
      ENDDO
-     ! Introduce the aerosol direct and first indirect radiative forcings:
-     tau_ae = 0.
-     piz_ae = 0.
-     cg_ae = 0.
      ! Param\`etres optiques des nuages et quelques param\`etres pour
      ! diagnostics :
      if (ok_newmicro) then
         CALL newmicro(paprs, play, t_seri, cldliq, cldfra, cldtau, cldemi, &
-             cldh, cldl, cldm, cldt, cldq, flwp, fiwp, flwc, fiwc, ok_aie, &
+             cldh, cldl, cldm, cldt, cldq, flwp, fiwp, flwc, fiwc)
-             sulfate, sulfate_pi, bl95_b0, bl95_b1, cldtaupi, re, fl)
      else
         CALL nuage(paprs, play, t_seri, cldliq, cldfra, cldtau, cldemi, cldh, &
-             cldl, cldm, cldt, cldq, ok_aie, sulfate, sulfate_pi, bl95_b0, &
+             cldl, cldm, cldt, cldq)
-             bl95_b1, cldtaupi, re, fl)
      endif
      IF (MOD(itap - 1, radpas) == 0) THEN
-        ! Appeler le rayonnement mais calculer tout d'abord l'albedo du sol.
+        wo = ozonecm(REAL(julien), paprs)
-        ! Calcul de l'abedo moyen par maille
         albsol = sum(falbe * pctsrf, dim = 2)
+        CALL radlwsw(dist, mu0, fract, paprs, play, tsol, albsol, t_seri, &
-        ! Rayonnement (compatible Arpege-IFS) :
-        CALL radlwsw(dist, mu0, fract, paprs, play, ztsol, albsol, t_seri, &
              q_seri, wo, cldfra, cldemi, cldtau, heat, heat0, cool, cool0, &
              radsol, albpla, topsw, toplw, solsw, sollw, sollwdown, topsw0, &
              toplw0, solsw0, sollw0, lwdn0, lwdn, lwup0, lwup, swdn0, swdn, &
-             swup0, swup, ok_ade, ok_aie, tau_ae, piz_ae, cg_ae, topswad, &
+             swup0, swup, ok_ade, topswad, solswad)
-             solswad, cldtaupi, topswai, solswai)
      ENDIF
      ! Ajouter la tendance des rayonnements (tous les pas)
      DO k = 1, llm
         DO i = 1, klon
            t_seri(i, k) = t_seri(i, k) + (heat(i, k) - cool(i, k)) * dtphys &
-Line 1012 
 contains
+Line 888 
 contains
         ENDDO
      ENDDO
-     ! Calculer l'hydrologie de la surface
-     DO i = 1, klon
-        zxqsurf(i) = 0.
-        zxsnow(i) = 0.
-     ENDDO
-     DO nsrf = 1, nbsrf
-        DO i = 1, klon
-           zxqsurf(i) = zxqsurf(i) + fqsurf(i, nsrf) * pctsrf(i, nsrf)
-           zxsnow(i) = zxsnow(i) + fsnow(i, nsrf) * pctsrf(i, nsrf)
-        ENDDO
-     ENDDO
      ! Calculer le bilan du sol et la d\'erive de temp\'erature (couplage)
      DO i = 1, klon
         bils(i) = radsol(i) - sens(i) + zxfluxlat(i)
      ENDDO
-Line 1106 
 contains
+Line 969 
 contains
           pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, tr_seri, &
           zmasse, ncid_startphy)
-     IF (offline) call phystokenc(dtphys, t, mfu, mfd, pen_u, pde_u, pen_d, &
-          pde_d, fm_therm, entr_therm, ycoefh, yu1, yv1, ftsol, pctsrf, &
-          frac_impa, frac_nucl, pphis, airephy, dtphys)
      ! Calculer le transport de l'eau et de l'energie (diagnostique)
      CALL transp(paprs, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, ue, uq)
-Line 1123 
 contains
+Line 982 
 contains
      ! conversion Ec en énergie thermique
      DO k = 1, llm
         DO i = 1, klon
-           ZRCPD = RCPD * (1. + RVTMP2 * q_seri(i, k))
+           d_t_ec(i, k) = 0.5 / (RCPD * (1. + RVTMP2 * q_seri(i, k))) &
-           d_t_ec(i, k) = 0.5 / ZRCPD &
                 * (u(i, k)**2 + v(i, k)**2 - u_seri(i, k)**2 - v_seri(i, k)**2)
            t_seri(i, k) = t_seri(i, k) + d_t_ec(i, k)
            d_t_ec(i, k) = d_t_ec(i, k) / dtphys
-Line 1175 
 contains
+Line 1033 
 contains
      CALL histwrite_phy("precip", rain_fall + snow_fall)
      CALL histwrite_phy("plul", rain_lsc + snow_lsc)
      CALL histwrite_phy("pluc", rain_con + snow_con)
-     CALL histwrite_phy("tsol", ztsol)
+     CALL histwrite_phy("tsol", tsol)
      CALL histwrite_phy("t2m", zt2m)
      CALL histwrite_phy("q2m", zq2m)
      CALL histwrite_phy("u10m", zu10m)
-Line 1199 
 contains
+Line 1057 
 contains
      DO nsrf = 1, nbsrf
         CALL histwrite_phy("pourc_"//clnsurf(nsrf), pctsrf(:, nsrf) * 100.)
         CALL histwrite_phy("fract_"//clnsurf(nsrf), pctsrf(:, nsrf))
-        CALL histwrite_phy("sens_"//clnsurf(nsrf), fluxt(:, 1, nsrf))
+        CALL histwrite_phy("sens_"//clnsurf(nsrf), flux_t(:, nsrf))
         CALL histwrite_phy("lat_"//clnsurf(nsrf), fluxlat(:, nsrf))
         CALL histwrite_phy("tsol_"//clnsurf(nsrf), ftsol(:, nsrf))
-        CALL histwrite_phy("taux_"//clnsurf(nsrf), fluxu(:, 1, nsrf))
+        CALL histwrite_phy("taux_"//clnsurf(nsrf), flux_u(:, nsrf))
-        CALL histwrite_phy("tauy_"//clnsurf(nsrf), fluxv(:, 1, nsrf))
+        CALL histwrite_phy("tauy_"//clnsurf(nsrf), flux_v(:, nsrf))
         CALL histwrite_phy("rugs_"//clnsurf(nsrf), frugs(:, nsrf))
         CALL histwrite_phy("albe_"//clnsurf(nsrf), falbe(:, nsrf))
      END DO
      CALL histwrite_phy("albs", albsol)
+     CALL histwrite_phy("tro3", wo * dobson_u * 1e3 / zmasse / rmo3 * md)
      CALL histwrite_phy("rugs", zxrugs)
      CALL histwrite_phy("s_pblh", s_pblh)
      CALL histwrite_phy("s_pblt", s_pblt)
-Line 1220 
 contains
+Line 1079 
 contains
      CALL histwrite_phy("s_trmb1", s_trmb1)
      CALL histwrite_phy("s_trmb2", s_trmb2)
      CALL histwrite_phy("s_trmb3", s_trmb3)
-     if (conv_emanuel) CALL histwrite_phy("ptop", ema_pct)
+     if (conv_emanuel) then
+        CALL histwrite_phy("ptop", ema_pct)
+        CALL histwrite_phy("dnwd0", - mp)
+     end if
      CALL histwrite_phy("temp", t_seri)
      CALL histwrite_phy("vitu", u_seri)
      CALL histwrite_phy("vitv", v_seri)
-Line 1229 
 contains
+Line 1093 
 contains
      CALL histwrite_phy("dtvdf", d_t_vdf)
      CALL histwrite_phy("dqvdf", d_q_vdf)
      CALL histwrite_phy("rhum", zx_rh)
+     CALL histwrite_phy("d_t_ec", d_t_ec)
+     CALL histwrite_phy("dtsw0", heat0 / 86400.)
+     CALL histwrite_phy("dtlw0", - cool0 / 86400.)
+     CALL histwrite_phy("msnow", sum(fsnow * pctsrf, dim = 2))
+     call histwrite_phy("qsurf", sum(fqsurf * pctsrf, dim = 2))
      if (ok_instan) call histsync(nid_ins)

 Legend:



Removed from v.205
 


changed lines


 
Added in v.224
 Legend:



Removed from v.205
 


changed lines


 
Added in v.224
-Removed from v.205
+Added in v.224

	ViewVC Help
Powered by ViewVC 1.1.21