Changeset 3761 for LMDZ6/branches


Ignore:
Timestamp:
Jul 15, 2020, 10:14:35 PM (4 years ago)
Author:
adurocher
Message:

Moved compressed part of cv3a_driver to separate subroutine

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3a_driver.f90

    r3760 r3761  
    11module cv3a_driver_mod
     2  implicit none
    23contains
    34
    … …  
    164165    USE print_control_mod, ONLY: prt_level, lunout
    165166    USE cv3a_compress_mod
    166     USE cv3p_mixing_mod, ONLY : cv3p_mixing
    167167    USE profiling_physic_mod
    168168    IMPLICIT NONE
    … …  
    262262! Local (non compressed) arrays
    263263    INTEGER nk1(len)
    264     INTEGER icb1(len)
    265264    INTEGER icbs1(len)
    266265
    … …  
    374373                    p1feed1, p2feed1, wght1, &
    375374                    wghti1, tnk1, thnk1, qnk1, qsnk1, unk1, vnk1, &
    376                     cpnk1, hnk1, nk1, icb1, icbmax, iflag1, gznk1, plcl1)
     375                    cpnk1, hnk1, nk1, kbas1, icbmax, iflag1, gznk1, plcl1)
    377376
    378377! --------------------------------------------------------------------
    … …  
    386385      tvp1(:, :) = 0.
    387386      clw1(:, :) = 0.
    388       CALL cv3_undilute1(len, nd, t1, qs1, gz1, plcl1, p1, icb1, tnk1, qnk1, & ! nd->na
     387      CALL cv3_undilute1(len, nd, t1, qs1, gz1, plcl1, p1, kbas1, tnk1, qnk1, & ! nd->na
    389388                         gznk1, tp1, tvp1, clw1, icbs1)
    390389
    … …  
    394393
    395394      call driver_log('cv3_trigger')
    396       CALL cv3_trigger(len, nd, icb1, plcl1, p1, th1, tv1, tvp1, thnk1, & ! nd->na
     395      CALL cv3_trigger(len, nd, kbas1, plcl1, p1, th1, tv1, tvp1, thnk1, & ! nd->na
    397396                       pbase1, buoybase1, iflag1, sig1, w01)
    398397      call exit_profile("cv3a_uncompressed")
    … …  
    505504      call add_array_i1(cv3a_compress_list, iflag1, iflag)
    506505      call add_array_i1(cv3a_compress_list, nk1, nk)
    507       call add_array_i1(cv3a_compress_list, icb1, icb)
     506      call add_array_i1(cv3a_compress_list, kbas1, icb)
    508507      call add_array_i1(cv3a_compress_list, icbs1, icbs)
    509508      call add_array_r1(cv3a_compress_list, plcl1, plcl)
    … …  
    555554      call exit_profile("cv3a_compress")
    556555
    557       if( timers_first ) then
    558         call enter_profile("cv3p_mixing")
    559         call exit_profile("cv3p_mixing")
    560         call enter_profile("cv3_yield")
    561         call exit_profile("cv3_yield")
    562         call enter_profile("cv3_tracer")       
    563         call exit_profile("cv3_tracer")
    564 
    565         timers_first = .false.
    566       endif
    567 
    568       call enter_profile("cv3a_compressed")
    569       IF (ncum > 0) THEN
    570 
    571 ! -------------------------------------------------------------------
    572 ! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
    573 ! ---   FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
    574 ! ---   &
    575 ! ---   COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
    576 ! ---   FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
    577 ! ---   &
    578 ! ---   FIND THE LEVEL OF NEUTRAL BUOYANCY
    579 ! -------------------------------------------------------------------
    580 
    581         call driver_log('cv3_undilute2')
    582         CALL cv3_undilute2(nloc, ncum, nd, iflag, icb, icbs, nk, &
    583                            tnk, qnk, gznk, hnk, t, q, qs, gz, &
    584                            p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &
    585                            inb, tp, tvp, clw, hp, ep, sigp, buoy, &
    586                            frac_a, frac_s, qpreca, qta)
    587 
    588 ! epmax_cape
    589 ! on recalcule ep et hp
    590         call driver_log('cv3_epmax_cape')
    591         call cv3_epmax_fn_cape(nloc, ncum, nd &
    592                                , ep, hp, icb, inb, clw, nk, t, h, hnk, lv, lf, frac_s &
    593                                , pbase, p, ph, tv, buoy, sig, w0, iflag &
    594                                , epmax_diag)
    595 
    596 ! -------------------------------------------------------------------
    597 ! --- MIXING(1)   (if iflag_mix .ge. 1)
    598 ! -------------------------------------------------------------------
    599         call enter_profile("cv3p_mixing")
    600         IF (iflag_mix >= 1) THEN
    601           CALL zilch(supmax, nloc*nd)
    602           call driver_log('cv3p_mixing')
    603           CALL cv3p_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb, &
    604                            ph, t, q, qs, u, v, tra, h, lv, lf, frac_s, qta, &
    605                            unk, vnk, hp, tv, tvp, ep, clw, sig, &
    606                            ment, qent, hent, uent, vent, nent, &
    607                            sigij, elij, supmax, ments, qents, traent)
    608         ELSE
    609           CALL zilch(supmax, nloc*nd)
    610         END IF
    611         call exit_profile("cv3p_mixing")
    612 
    613 ! -------------------------------------------------------------------
    614 ! --- CLOSURE
    615 ! -------------------------------------------------------------------
    616 
    617         ptop2(:) = 0
    618         coef_clos(:) = 1.
    619         IF (iflag_clos == 0) THEN
    620           call driver_log('cv3_closure')
    621           CALL cv3_closure(nloc, ncum, nd, icb, inb, &
    622                            pbase, p, ph, tv, buoy, &
    623                            sig, w0, cape, m, iflag)
    624         END IF   ! iflag_clos==0
    625 
    626         ok_inhib = iflag_mix == 2
    627 
    628         IF (iflag_clos == 1) PRINT *, ' pas d appel cv3p_closure'
    629 
    630         IF (iflag_clos == 2) THEN
    631           call driver_log('cv3p1_closure')
    632           CALL cv3p1_closure(nloc, ncum, nd, icb, inb, &
    633                              pbase, plcl, p, ph, tv, tvp, buoy, &
    634                              supmax, ok_inhib, Ale, Alp, omega, &
    635                              sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
    636                              Plim1, plim2, asupmax, supmax0, &
    637                              asupmaxmin, cbmf, plfc, wbeff)
    638           if (prt_level >= 10) PRINT *, 'cv3p1_closure-> plfc,wbeff ', plfc(1), wbeff(1)
    639         END IF ! iflag_clos==2
    640 
    641         IF (iflag_clos == 3) THEN
    642           call driver_log('cv3p2_closure')
    643           CALL cv3p2_closure(nloc, ncum, nd, icb, inb, &
    644                              pbase, plcl, p, ph, tv, tvp, buoy, &
    645                              supmax, ok_inhib, Ale, Alp, omega, &
    646                              sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
    647                              Plim1, plim2, asupmax, supmax0, &
    648                              asupmaxmin, cbmf, plfc, wbeff)
    649           if (prt_level >= 10) PRINT *, 'cv3p2_closure-> plfc,wbeff ', plfc(1), wbeff(1)
    650         END IF ! iflag_clos==3
    651 
    652 ! -------------------------------------------------------------------
    653 ! --- MIXING(2)
    654 ! -------------------------------------------------------------------
    655 
    656         IF (iflag_mix == 0) THEN
    657           call driver_log('cv3_mixing')
    658           CALL cv3_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb, &             ! na->nd
    659                           ph, t, q, qs, u, v, tra, h, lv, lf, frac_s, qnk, &
    660                           unk, vnk, hp, tv, tvp, ep, clw, m, sig, &
    661                           ment, qent, uent, vent, nent, sigij, elij, ments, qents, traent)
    662           CALL zilch(hent, nloc*nd*nd)
    663         ELSE
    664           mm(:, :) = m(:, :)
    665           CALL cv3_mixscale(nloc, ncum, nd, ment, mm)
    666           IF (debut) PRINT *, ' cv3_mixscale-> '
    667         END IF
    668 
    669         IF (debut) PRINT *, ' cv_mixing ->'
    670 
    671 ! -------------------------------------------------------------------
    672 ! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
    673 ! -------------------------------------------------------------------
    674         IF (debut) PRINT *, ' cva_driver -> cv3_unsat '
    675 
    676         call driver_log('cv3_unsat')
    677         CALL cv3_unsat(nloc, ncum, nd, nd, ntra, icb, inb, iflag, &
    678                        t_wake, q_wake, qs_wake, gz, u, v, tra, p, ph, &
    679                        th_wake, tv_wake, lv_wake, lf_wake, cpn_wake, &
    680                        ep, sigp, clw, frac_s, qpreca, frac_a, qta, &
    681                        m, ment, elij, delt, plcl, coef_clos, &
    682                        mp, qp, up, vp, trap, wt, water, evap, fondue, ice, &
    683                        faci, b, sigd, &
    684                        wdtrainA, wdtrainS, wdtrainM)
    685         IF (prt_level >= 10) THEN
    686           Print *, 'cva_driver after cv3_unsat:mp , water, ice, evap, fondue '
    687           DO k = 1, nd
    688             write (6, '(i4,5(1x,e13.6))'), &
    689               k, mp(igout, k), water(igout, k), ice(igout, k), &
    690               evap(igout, k), fondue(igout, k)
    691           ENDDO
    692           Print *, 'cva_driver after cv3_unsat: wdtrainA, wdtrainS, wdtrainM '
    693           DO k = 1, nd
    694             write (6, '(i4,3(1x,e13.6))'), &
    695               k, wdtrainA(igout, k), wdtrainS(igout, k), wdtrainM(igout, k)
    696           ENDDO
    697         ENDIF
    698 
    699         IF (debut) PRINT *, 'cv_unsat-> '
    700 ! -------------------------------------------------------------------
    701 ! YIELD
    702 ! (tendencies, precipitation, variables of interface with other processes, etc)
    703 ! -------------------------------------------------------------------
    704 
    705         call driver_log('cv3_yield')
    706         call enter_profile("cv3_yield")
    707         CALL cv3_yield(nloc, ncum, nd, nd, ntra, ok_conserv_q, &
    708                        icb, inb, delt, &
    709                        t, q, t_wake, q_wake, s_wake, u, v, tra, &
    710                        gz, p, ph, h, hp, lv, lf, cpn, th, th_wake, &
    711                        ep, clw, qpreca, m, tp, mp, qp, up, vp, trap, &
    712                        wt, water, ice, evap, fondue, faci, b, sigd, &
    713                        ment, qent, hent, iflag_mix, uent, vent, &
    714                        nent, elij, traent, sig, &
    715                        tv, tvp, wghti, &
    716                        iflag, precip, Vprecip, Vprecipi, ft, fq, fu, fv, ftra, &
    717                        cbmf, upwd, dnwd, dnwd0, ma, mip, &
    718                        qcondc, wd, &
    719                        ftd, fqd, qta, qtc, sigt, tau_cld_cv, coefw_cld_cv)
    720         call exit_profile("cv3_yield")
    721         ! Test conseravtion de l'eau
    722         IF (debut) PRINT *, ' cv3_yield -> fqd(1) = ', fqd(igout, 1)
    723         IF (prt_level >= 10) THEN
    724           Print *, 'cva_driver after cv3_yield:ft(1) , ftd(1) ', &
    725             ft(igout, 1), ftd(igout, 1)
    726           Print *, 'cva_driver after cv3_yield:fq(1) , fqd(1) ', &
    727             fq(igout, 1), fqd(igout, 1)
    728         ENDIF
    729 
    730 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    731 !--- passive tracers
    732 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    733 
    734         call driver_log('cv3_tracer')
    735         call enter_profile("cv3_tracer")
    736         CALL cv3_tracer(nloc, len, ncum, nd, nd, &
    737                         ment, sigij, da, phi, phi2, d1a, dam, &
    738                         ep, vprecip, elij, clw, epmlmMm, eplaMm, &
    739                         icb, inb)
    740         call exit_profile("cv3_tracer")
    741       END IF ! ncum>0
    742       call exit_profile("cv3a_compressed")
     556      call cv3a_driver_compressed(nloc, nd, ntra, &
     557                                  debut, ok_conserv_q, iflag_mix, iflag_clos, prt_level, &
     558                                  ! Input fields
     559                                  delt, tau_cld_cv, coefw_cld_cv, nk, icbs, tnk, qnk, gznk, hnk, unk, vnk, pbase, buoybase, s_wake, Ale, Alp, wghti, th, t, q, qs, t_wake, q_wake, qs_wake, u, v, gz, h, th_wake, lv, lf, cpn, p, tv, tp, h_wake, lv_wake, lf_wake, cpn_wake, tv_wake, omega, ph, &
     560                                  ! Input/output
     561                                  iflag, icb, plcl, sig, w0, tvp, clw, &
     562                                  ! Output fields
     563                                  inb, precip, cbmf, plfc, wbeff, ptop2, sigd, wd, cape, cin, Plim1, plim2, supmax0, asupmaxmin, epmax_diag, ft, fq, fu, fv, ma, mip, upwd, dnwd, dnwd0, qcondc, ftd, fqd, asupmax, da, mp, d1a, dam, qta,  evap, ep, eplaMm, wdtrainA, wdtrainS, wdtrainM, qtc, sigt, vprecip, vprecipi, phi, phi2, sigij, elij, epmlmMm)
    743564
    744565! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    … …  
    826647  END SUBROUTINE
    827648
     649  subroutine cv3a_driver_compressed(nloc, nd, ntra, &
     650                                    debut, ok_conserv_q, iflag_mix, iflag_clos, prt_level, &
     651                                    ! Input fields
     652                                    delt, tau_cld_cv, coefw_cld_cv, nk, icbs, tnk, qnk, gznk, hnk, unk, vnk, pbase, buoybase, s_wake, Ale, Alp, wghti, th, t, q, qs, t_wake, q_wake, qs_wake, u, v, gz, h, th_wake, lv, lf, cpn, p, tv, tp, h_wake, lv_wake, lf_wake, cpn_wake, tv_wake, omega, ph, &
     653                                    ! Input/output
     654                                    iflag, icb, plcl, sig, w0, tvp, clw, &
     655                                    ! Output fields
     656                                    inb, precip, cbmf, plfc, wbeff, ptop2, sigd, wd, cape, cin, Plim1, plim2, supmax0, asupmaxmin, epmax_diag, ft, fq, fu, fv, ma, mip, upwd, dnwd, dnwd0, qcondc, ftd, fqd, asupmax, da, mp, d1a, dam, qta,  evap, ep, eplaMm, wdtrainA, wdtrainS, wdtrainM, qtc, sigt, vprecip, vprecipi, phi, phi2, sigij, elij, epmlmMm)
     657    use profiling_physic_mod
     658    USE cv3p_mixing_mod, ONLY : cv3p_mixing
     659
     660    integer, intent(in) :: nloc, nd, ntra ! Arrays sizes
     661    logical, intent(in) :: debut, ok_conserv_q
     662    integer, intent(in) :: iflag_mix, iflag_clos, prt_level
     663    real, intent(in)    :: delt, tau_cld_cv, coefw_cld_cv
     664
     665    ! Input arrays
     666    integer, dimension(nloc), intent(in) :: nk, icbs
     667    real, dimension(nloc), intent(in) :: tnk, qnk, gznk, hnk, unk, vnk, pbase, buoybase, s_wake, Ale, Alp
     668    real, dimension(nloc, nd), intent(in) :: wghti, th, t, q, qs, t_wake, q_wake, qs_wake, u, v, gz, h, th_wake, lv, lf, cpn, p, tv, tp, h_wake, lv_wake, lf_wake, cpn_wake, tv_wake, omega
     669    real, dimension(nloc, nd + 1), intent(in) :: ph
     670
     671    ! Input/output arrays
     672    integer, dimension(nloc), intent(inout) :: iflag, icb
     673    real, dimension(nloc), intent(inout) :: plcl
     674    real, dimension(nloc, nd), intent(inout) :: sig, w0, tvp, clw
     675
     676    ! Output arrays
     677    integer, dimension(nloc), intent(out) :: inb
     678    real, dimension(nloc), intent(out) :: precip, cbmf, plfc, wbeff, ptop2, sigd, wd, cape, cin, Plim1, plim2, supmax0, asupmaxmin, epmax_diag
     679    real, dimension(nloc, nd), intent(out) :: ft, fq, fu, fv, ma, mip, upwd, dnwd, dnwd0, qcondc, ftd, fqd, asupmax, da, mp, d1a, dam, qta,  evap, ep, eplaMm, wdtrainA, wdtrainS, wdtrainM, qtc, sigt
     680    real, dimension(nloc, nd + 1), intent(out) :: vprecip, vprecipi
     681    real, dimension(nloc, nd, nd), intent(out) :: phi, phi2, sigij, elij, epmlmMm
     682
     683    ! Local fields
     684    INTEGER nent(nloc, nd)
     685    REAL hp(nloc, nd), sigp(nloc, nd), qpreca(nloc, nd)
     686    REAL buoy(nloc, nd)
     687    REAL m(nloc, nd)
     688    REAL mm(nloc, nd)
     689    REAL ment(nloc, nd, nd)
     690    REAL qent(nloc, nd, nd)
     691    REAL hent(nloc, nd, nd)
     692    REAL uent(nloc, nd, nd), vent(nloc, nd, nd)
     693    REAL ments(nloc, nd, nd), qents(nloc, nd, nd)
     694    REAL supmax(nloc, nd)
     695    REAL coef_clos(nloc)
     696    REAL, DIMENSION(nloc, nd)     :: qp, up, vp
     697    REAL, DIMENSION(nloc, nd)     :: wt, water
     698    REAL, DIMENSION(nloc, nd)     :: ice, fondue, b
     699    REAL, DIMENSION(nloc, nd)     :: frac_a, frac_s, faci
     700    REAL tra(nloc, nd, ntra), trap(nloc, nd, ntra)
     701    REAL ftra(nloc, nd, ntra), traent(nloc, nd, nd, ntra)
     702
     703    LOGICAL ok_inhib ! True => possible inhibition of convection by dryness
     704    integer :: k
     705
     706    integer :: ncum ! Number of convective cells
     707    logical, save :: timers_first = .true.
     708    !$omp threadprivate(timers_first)
     709    INTEGER, PARAMETER :: igout = 1
     710
     711    if (timers_first) then
     712      call enter_profile("cv3p_mixing")
     713      call exit_profile("cv3p_mixing")
     714      call enter_profile("cv3_yield")
     715      call exit_profile("cv3_yield")
     716      call enter_profile("cv3_tracer")
     717      call exit_profile("cv3_tracer")
     718
     719      timers_first = .false.
     720    endif
     721
     722    ncum = nloc
     723    call enter_profile("cv3a_compressed")
     724    IF (ncum > 0) THEN
     725
     726! -------------------------------------------------------------------
     727! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
     728! ---   FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
     729! ---   &
     730! ---   COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
     731! ---   FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
     732! ---   &
     733! ---   FIND THE LEVEL OF NEUTRAL BUOYANCY
     734! -------------------------------------------------------------------
     735
     736      call driver_log('cv3_undilute2')
     737      CALL cv3_undilute2(nloc, ncum, nd, iflag, icb, icbs, nk, &
     738                         tnk, qnk, gznk, hnk, t, q, qs, gz, &
     739                         p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &
     740                         inb, tp, tvp, clw, hp, ep, sigp, buoy, &
     741                         frac_a, frac_s, qpreca, qta)
     742
     743! epmax_cape
     744! on recalcule ep et hp
     745      call driver_log('cv3_epmax_cape')
     746      call cv3_epmax_fn_cape(nloc, ncum, nd &
     747                             , ep, hp, icb, inb, clw, nk, t, h, hnk, lv, lf, frac_s &
     748                             , pbase, p, ph, tv, buoy, sig, w0, iflag &
     749                             , epmax_diag)
     750
     751! -------------------------------------------------------------------
     752! --- MIXING(1)   (if iflag_mix .ge. 1)
     753! -------------------------------------------------------------------
     754      call enter_profile("cv3p_mixing")
     755      IF (iflag_mix >= 1) THEN
     756        CALL zilch(supmax, nloc*nd)
     757        call driver_log('cv3p_mixing')
     758        CALL cv3p_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb, &
     759                         ph, t, q, qs, u, v, tra, h, lv, lf, frac_s, qta, &
     760                         unk, vnk, hp, tv, tvp, ep, clw, sig, &
     761                         ment, qent, hent, uent, vent, nent, &
     762                         sigij, elij, supmax, ments, qents, traent)
     763      ELSE
     764        CALL zilch(supmax, nloc*nd)
     765      END IF
     766      call exit_profile("cv3p_mixing")
     767
     768! -------------------------------------------------------------------
     769! --- CLOSURE
     770! -------------------------------------------------------------------
     771
     772      ptop2(:) = 0
     773      coef_clos(:) = 1.
     774      IF (iflag_clos == 0) THEN
     775        call driver_log('cv3_closure')
     776        CALL cv3_closure(nloc, ncum, nd, icb, inb, &
     777                         pbase, p, ph, tv, buoy, &
     778                         sig, w0, cape, m, iflag)
     779      END IF   ! iflag_clos==0
     780
     781      ok_inhib = iflag_mix == 2
     782
     783      IF (iflag_clos == 1) PRINT *, ' pas d appel cv3p_closure'
     784
     785      IF (iflag_clos == 2) THEN
     786        call driver_log('cv3p1_closure')
     787        CALL cv3p1_closure(nloc, ncum, nd, icb, inb, &
     788                           pbase, plcl, p, ph, tv, tvp, buoy, &
     789                           supmax, ok_inhib, Ale, Alp, omega, &
     790                           sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
     791                           Plim1, plim2, asupmax, supmax0, &
     792                           asupmaxmin, cbmf, plfc, wbeff)
     793        if (prt_level >= 10) PRINT *, 'cv3p1_closure-> plfc,wbeff ', plfc(1), wbeff(1)
     794      END IF ! iflag_clos==2
     795
     796      IF (iflag_clos == 3) THEN
     797        call driver_log('cv3p2_closure')
     798        CALL cv3p2_closure(nloc, ncum, nd, icb, inb, &
     799                           pbase, plcl, p, ph, tv, tvp, buoy, &
     800                           supmax, ok_inhib, Ale, Alp, omega, &
     801                           sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
     802                           Plim1, plim2, asupmax, supmax0, &
     803                           asupmaxmin, cbmf, plfc, wbeff)
     804        if (prt_level >= 10) PRINT *, 'cv3p2_closure-> plfc,wbeff ', plfc(1), wbeff(1)
     805      END IF ! iflag_clos==3
     806
     807! -------------------------------------------------------------------
     808! --- MIXING(2)
     809! -------------------------------------------------------------------
     810
     811      IF (iflag_mix == 0) THEN
     812        call driver_log('cv3_mixing')
     813        CALL cv3_mixing(nloc, ncum, nd, nd, ntra, icb, nk, inb, &             ! na->nd
     814                        ph, t, q, qs, u, v, tra, h, lv, lf, frac_s, qnk, &
     815                        unk, vnk, hp, tv, tvp, ep, clw, m, sig, &
     816                        ment, qent, uent, vent, nent, sigij, elij, ments, qents, traent)
     817        CALL zilch(hent, nloc*nd*nd)
     818      ELSE
     819        mm(:, :) = m(:, :)
     820        CALL cv3_mixscale(nloc, ncum, nd, ment, mm)
     821        IF (debut) PRINT *, ' cv3_mixscale-> '
     822      END IF
     823
     824      IF (debut) PRINT *, ' cv_mixing ->'
     825
     826! -------------------------------------------------------------------
     827! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
     828! -------------------------------------------------------------------
     829      IF (debut) PRINT *, ' cva_driver -> cv3_unsat '
     830
     831      call driver_log('cv3_unsat')
     832      CALL cv3_unsat(nloc, ncum, nd, nd, ntra, icb, inb, iflag, &
     833                     t_wake, q_wake, qs_wake, gz, u, v, tra, p, ph, &
     834                     th_wake, tv_wake, lv_wake, lf_wake, cpn_wake, &
     835                     ep, sigp, clw, frac_s, qpreca, frac_a, qta, &
     836                     m, ment, elij, delt, plcl, coef_clos, &
     837                     mp, qp, up, vp, trap, wt, water, evap, fondue, ice, &
     838                     faci, b, sigd, &
     839                     wdtrainA, wdtrainS, wdtrainM)
     840      IF (prt_level >= 10) THEN
     841        Print *, 'cva_driver after cv3_unsat:mp , water, ice, evap, fondue '
     842        DO k = 1, nd
     843          write (6, '(i4,5(1x,e13.6))'), &
     844            k, mp(igout, k), water(igout, k), ice(igout, k), &
     845            evap(igout, k), fondue(igout, k)
     846        ENDDO
     847        Print *, 'cva_driver after cv3_unsat: wdtrainA, wdtrainS, wdtrainM '
     848        DO k = 1, nd
     849          write (6, '(i4,3(1x,e13.6))'), &
     850            k, wdtrainA(igout, k), wdtrainS(igout, k), wdtrainM(igout, k)
     851        ENDDO
     852      ENDIF
     853
     854      IF (debut) PRINT *, 'cv_unsat-> '
     855! -------------------------------------------------------------------
     856! YIELD
     857! (tendencies, precipitation, variables of interface with other processes, etc)
     858! -------------------------------------------------------------------
     859
     860      call driver_log('cv3_yield')
     861      call enter_profile("cv3_yield")
     862      CALL cv3_yield(nloc, ncum, nd, nd, ntra, ok_conserv_q, &
     863                     icb, inb, delt, &
     864                     t, q, t_wake, q_wake, s_wake, u, v, tra, &
     865                     gz, p, ph, h, hp, lv, lf, cpn, th, th_wake, &
     866                     ep, clw, qpreca, m, tp, mp, qp, up, vp, trap, &
     867                     wt, water, ice, evap, fondue, faci, b, sigd, &
     868                     ment, qent, hent, iflag_mix, uent, vent, &
     869                     nent, elij, traent, sig, &
     870                     tv, tvp, wghti, &
     871                     iflag, precip, Vprecip, Vprecipi, ft, fq, fu, fv, ftra, &
     872                     cbmf, upwd, dnwd, dnwd0, ma, mip, &
     873                     qcondc, wd, &
     874                     ftd, fqd, qta, qtc, sigt, tau_cld_cv, coefw_cld_cv)
     875      call exit_profile("cv3_yield")
     876      ! Test conseravtion de l'eau
     877      IF (debut) PRINT *, ' cv3_yield -> fqd(1) = ', fqd(igout, 1)
     878      IF (prt_level >= 10) THEN
     879        Print *, 'cva_driver after cv3_yield:ft(1) , ftd(1) ', &
     880          ft(igout, 1), ftd(igout, 1)
     881        Print *, 'cva_driver after cv3_yield:fq(1) , fqd(1) ', &
     882          fq(igout, 1), fqd(igout, 1)
     883      ENDIF
     884
     885!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
     886!--- passive tracers
     887!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
     888
     889      call driver_log('cv3_tracer')
     890      call enter_profile("cv3_tracer")
     891      CALL cv3_tracer(nloc, -1, ncum, nd, nd, &
     892                      ment, sigij, da, phi, phi2, d1a, dam, &
     893                      ep, vprecip, elij, clw, epmlmMm, eplaMm, &
     894                      icb, inb)
     895      call exit_profile("cv3_tracer")
     896    END IF ! ncum>0
     897    call exit_profile("cv3a_compressed")
     898  end subroutine
     899
    828900  subroutine driver_log(message)
    829901    use print_control_mod, only: prt_level
Note: See TracChangeset for help on using the changeset viewer.