Ignore:
Timestamp:
Mar 3, 2015, 2:41:13 PM (10 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes -r2186:2216 into testing branch

Location:
LMDZ5/branches/testing
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/testing

  • LMDZ5/branches/testing/libf/phylmd/cva_driver.F90

    r2160 r2220  
    77                      u1, v1, tra1, &
    88                      p1, ph1, &
    9                       Ale1, Alp1, &
     9                      Ale1, Alp1, omega1, &
    1010                      sig1feed1, sig2feed1, wght1, &
    1111                      iflag1, ft1, fq1, fu1, fv1, ftra1, &
    … …  
    2424                      da1, phi1, mp1, phi21, d1a1, dam1, sigij1, wghti1, & ! RomP, RL
    2525                      clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, &        ! RomP, RL
    26                       wdtrainA1, wdtrainM1)                                ! RomP
     26                      wdtrainA1, wdtrainM1, qtc1, sigt1, tau_cld_cv, &
     27                      coefw_cld_cv)                                        ! RomP, AJ
    2728! **************************************************************
    2829! *
    … …  
    5556! iflag_ice_thermo Integer        Input        accounting for ice thermodynamics (0/1)
    5657! iflag_clos    Integer        Input        version of closure (0/1)
     58! tau_cld_cv    Real           Input        characteristic time of dissipation of mixing fluxes
     59! coefw_cld_cv  Real           Input        coefficient for updraft velocity in convection
    5760! ok_conserv_q  Logical        Input        when true corrections for water conservation are swtiched on
    5861! delt          Real           Input        time step
    … …  
    119122! phi1          Real           Output     used in tracer transport (cvltr)
    120123! mp1           Real           Output     used in tracer transport (cvltr)
    121                                          
     124! qtc1          Real           Output     specific humidity in convection
     125! sigt1         Real           Output     surface fraction in adiabatic updrafts                                         
    122126! phi21         Real           Output     used in tracer transport (cvltr)
    123127                                         
    … …  
    163167  INTEGER iflag_clos
    164168  LOGICAL ok_conserv_q
     169  REAL tau_cld_cv
     170  REAL coefw_cld_cv
    165171  REAL delt
    166172  REAL t1(len, nd)
    … …  
    178184  REAL Ale1(len)
    179185  REAL Alp1(len)
     186  REAL omega1(len,nd)
    180187  REAL sig1feed1 ! pressure at lower bound of feeding layer
    181188  REAL sig2feed1 ! pressure at upper bound of feeding layer
    … …  
    225232  REAL asupmaxmin1(len)
    226233  INTEGER lalim_conv(len)
     234  REAL qtc1(len, nd)         ! cld
     235  REAL sigt1(len, nd)        ! cld
     236
    227237! RomP >>>
    228238  REAL wdtrainA1(len, nd), wdtrainM1(len, nd)
    … …  
    455465  REAL supmax(nloc, klev)
    456466  REAL Ale(nloc), Alp(nloc), coef_clos(nloc)
     467  REAL omega(nloc,klev)
    457468  REAL sigd(nloc)
    458469! real mp(nloc,klev), qp(nloc,klev), up(nloc,klev), vp(nloc,klev)
    … …  
    487498  REAL hnk(nloc), unk(nloc), vnk(nloc)
    488499
     500  REAL qtc(nloc, klev)         ! cld
     501  REAL sigt(nloc, klev)        ! cld
     502 
    489503! RomP >>>
    490504  REAL wdtrainA(nloc, klev), wdtrainM(nloc, klev)
    … …  
    593607
    594608! RomP >>>
     609  sigt1(:, :) = 0.
     610  qtc1(:, :) = 0.
    595611  wdtrainA1(:, :) = 0.
    596612  wdtrainM1(:, :) = 0.
    … …  
    776792                         h1_wake, lv1_wake, lf1_wake, cpn1_wake, tv1_wake, &
    777793                         sig1, w01, ptop21, &
    778                          Ale1, Alp1, &
     794                         Ale1, Alp1, omega1, &
    779795                         iflag, nk, icb, icbs, &
    780796                         plcl, tnk, qnk, gznk, hnk, unk, vnk, &
    … …  
    786802                         h_wake, lv_wake, lf_wake, cpn_wake, tv_wake, &
    787803                         sig, w0, ptop2, &
    788                          Ale, Alp)
     804                         Ale, Alp, omega)
    789805
    790806! print*,'tv ',tv
    … …  
    877893        CALL cv3p1_closure(nloc, ncum, nd, icb, inb, &         ! na->nd
    878894                           pbase, plcl, p, ph, tv, tvp, buoy, &
    879                            supmax, ok_inhib, Ale, Alp, &
     895                           supmax, ok_inhib, Ale, Alp, omega, &
    880896                           sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
    881897                           Plim1, plim2, asupmax, supmax0, &
    … …  
    978994                     cbmf, upwd, dnwd, dnwd0, ma, mip, &
    979995                     tls, tps, qcondc, wd, &
    980                      ftd, fqd)
     996                     ftd, fqd, qnk, qtc, sigt, tau_cld_cv, coefw_cld_cv)
    981997    END IF
    982998
    … …  
    10321048                           clw, elij, evap, ep, epmlmMm, eplaMm, &       ! RomP
    10331049                           wdtrainA, wdtrainM, &                         ! RomP
     1050                           qtc, sigt, &
    10341051                           iflag1, kbas1, ktop1, &
    10351052                           precip1, cbmf1, plcl1, plfc1, wbeff1, sig1, w01, ptop21, &
    … …  
    10431060                           da1, phi1, mp1, phi21, d1a1, dam1, sigij1,  & ! RomP
    10441061                           clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & ! RomP
    1045                            wdtrainA1, wdtrainM1)                         ! RomP
     1062                           wdtrainA1, wdtrainM1,                       & ! RomP
     1063                           qtc1, sigt1)
    10461064    END IF
    10471065
Note: See TracChangeset for help on using the changeset viewer.