--- trunk/phylmd/cv_driver.f 2014/08/29 13:00:05 103 +++ trunk/Sources/phylmd/cv_driver.f 2016/03/16 11:11:27 182 @@ -14,11 +14,6 @@ ! Several modules corresponding to different physical processes - ! Several versions of convect may be used: - ! - iflag_con = 3: version lmd - ! - iflag_con = 4: version 4.3b - - use clesphys2, only: iflag_con use cv3_compress_m, only: cv3_compress use cv3_feed_m, only: cv3_feed use cv3_mixing_m, only: cv3_mixing @@ -28,8 +23,6 @@ use cv3_uncompress_m, only: cv3_uncompress use cv3_unsat_m, only: cv3_unsat use cv3_yield_m, only: cv3_yield - use cv_feed_m, only: cv_feed - use cv_uncompress_m, only: cv_uncompress USE dimphy, ONLY: klev, klon real, intent(in):: t1(klon, klev) ! temperature @@ -46,7 +39,7 @@ real, intent(out):: fv1(klon, klev) ! v-wind tend real, intent(out):: precip1(klon) ! precipitation - real, intent(out):: VPrecip1(klon, klev+1) + real, intent(out):: VPrecip1(klon, klev + 1) ! vertical profile of precipitation real, intent(inout):: cbmf1(klon) ! cloud base mass flux @@ -60,7 +53,10 @@ real, intent(in):: delt ! time step real Ma1(klon, klev) ! Ma1 Real Output mass flux adiabatic updraft - real, intent(out):: upwd1(klon, klev) ! total upward mass flux (adiab+mixed) + + real, intent(out):: upwd1(klon, klev) + ! total upward mass flux (adiab + mixed) + real, intent(out):: dnwd1(klon, klev) ! saturated downward mass flux (mixed) real, intent(out):: dnwd01(klon, klev) ! unsaturated downward mass flux @@ -74,9 +70,9 @@ real, intent(inout):: da1(klon, klev), phi1(klon, klev, klev) real, intent(inout):: mp1(klon, klev) - ! --- ARGUMENTS + ! ARGUMENTS - ! --- On input: + ! On input: ! t: Array of absolute temperature (K) of dimension KLEV, with first ! index corresponding to lowest model level. Note that this array @@ -104,7 +100,7 @@ ! index corresponding to lowest model level. Must be defined ! at same grid levels as T. - ! ph: Array of pressure (mb) of dimension KLEV+1, with first index + ! ph: Array of pressure (mb) of dimension KLEV + 1, with first index ! corresponding to lowest level. These pressures are defined at ! levels intermediate between those of P, T, Q and QS. The first ! value of PH should be greater than (i.e. at a lower level than) @@ -115,7 +111,7 @@ ! delt: The model time step (sec) between calls to CONVECT - ! --- On Output: + ! On Output: ! iflag: An output integer whose value denotes the following: ! VALUE INTERPRETATION @@ -207,11 +203,11 @@ integer icbs(klon) integer inb(klon), inbis(klon) - real cbmf(klon), plcl(klon), tnk(klon), qnk(klon), gznk(klon) + real plcl(klon), tnk(klon), qnk(klon), gznk(klon) real t(klon, klev), q(klon, klev), qs(klon, klev) real u(klon, klev), v(klon, klev) real gz(klon, klev), h(klon, klev), lv(klon, klev), cpn(klon, klev) - real p(klon, klev), ph(klon, klev+1), tv(klon, klev), tp(klon, klev) + real p(klon, klev), ph(klon, klev + 1), tv(klon, klev), tp(klon, klev) real clw(klon, klev) real dph(klon, klev) real pbase(klon), buoybase(klon), th(klon, klev) @@ -232,32 +228,33 @@ real Ma(klon, klev), mike(klon, klev), tls(klon, klev) real tps(klon, klev), qprime(klon), tprime(klon) real precip(klon) - real VPrecip(klon, klev+1) + real VPrecip(klon, klev + 1) real qcondc(klon, klev) ! cld real wd(klon) ! gust !------------------------------------------------------------------- - ! --- SET CONSTANTS AND PARAMETERS - ! -- set simulation flags: + ! SET CONSTANTS AND PARAMETERS + + ! set simulation flags: ! (common cvflag) CALL cv_flag - ! -- set thermodynamical constants: + ! set thermodynamical constants: ! (common cvthermo) CALL cv_thermo - ! -- set convect parameters + ! set convect parameters ! includes microphysical parameters and parameters that ! control the rate of approach to quasi-equilibrium) ! (common cvparam) - if (iflag_con == 3) CALL cv3_param(klev, delt) + CALL cv3_param(klev, delt) - ! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS + ! INITIALIZE OUTPUT ARRAYS AND PARAMETERS do k = 1, klev do i = 1, klon @@ -285,191 +282,113 @@ iflag1(i) = 0 wd1(i) = 0.0 cape1(i) = 0.0 - VPrecip1(i, klev+1) = 0.0 + VPrecip1(i, klev + 1) = 0.0 end do - if (iflag_con == 3) then - do il = 1, klon - sig1(il, klev) = sig1(il, klev) + 1. - sig1(il, klev) = min(sig1(il, klev), 12.1) - enddo - endif - - ! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY - - if (iflag_con == 3) then - CALL cv3_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, & - gz1, h1, hm1, th1) - else - ! iflag_con == 4 - CALL cv_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, & - gz1, h1, hm1) - endif - - ! --- CONVECTIVE FEED - - if (iflag_con == 3) then - CALL cv3_feed(klon, klev, t1, q1, qs1, p1, ph1, hm1, gz1, nk1, icb1, & - icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na - else - ! iflag_con == 4 - CALL cv_feed(klon, klev, t1, q1, qs1, p1, hm1, gz1, nk1, icb1, icbmax, & - iflag1, tnk1, qnk1, gznk1, plcl1) - endif + do il = 1, klon + sig1(il, klev) = sig1(il, klev) + 1. + sig1(il, klev) = min(sig1(il, klev), 12.1) + enddo + + ! CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY + + CALL cv3_prelim(klon, klev, klev + 1, t1, q1, p1, ph1, lv1, cpn1, tv1, & + gz1, h1, hm1, th1) - ! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part - ! (up through ICB for convect4, up through ICB+1 for convect3) + ! CONVECTIVE FEED + + CALL cv3_feed(klon, klev, t1, q1, qs1, p1, ph1, gz1, nk1, icb1, & + icbmax, iflag1, tnk1, qnk1, gznk1, plcl1) ! klev->na + + ! UNDILUTE (ADIABATIC) UPDRAFT / 1st part + ! (up through ICB for convect4, up through ICB + 1 for convect3) ! Calculates the lifted parcel virtual temperature at nk, the ! actual temperature, and the adiabatic liquid water content. - if (iflag_con == 3) then - CALL cv3_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, & - tp1, tvp1, clw1, icbs1) ! klev->na - else - ! iflag_con == 4 - CALL cv_undilute1(klon, klev, t1, q1, qs1, gz1, p1, nk1, icb1, icbmax, & - tp1, tvp1, clw1) - endif - - ! --- TRIGGERING - - if (iflag_con == 3) then - CALL cv3_trigger(klon, klev, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, & - buoybase1, iflag1, sig1, w01) ! klev->na - else - ! iflag_con == 4 - CALL cv_trigger(klon, klev, icb1, cbmf1, tv1, tvp1, iflag1) - end if + CALL cv3_undilute1(klon, klev, t1, q1, qs1, gz1, plcl1, p1, nk1, icb1, & + tp1, tvp1, clw1, icbs1) ! klev->na + + ! TRIGGERING - ! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY + CALL cv3_trigger(klon, klev, icb1, plcl1, p1, th1, tv1, tvp1, pbase1, & + buoybase1, iflag1, sig1, w01) ! klev->na + + ! Moist convective adjustment is necessary ncum = 0 do i = 1, klon - if(iflag1(i) == 0)then - ncum = ncum+1 + if (iflag1(i) == 0) then + ncum = ncum + 1 idcum(ncum) = i endif end do IF (ncum > 0) THEN - ! --- COMPRESS THE FIELDS + ! COMPRESS THE FIELDS ! (-> vectorization over convective gridpoints) - if (iflag_con == 3) then - CALL cv3_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, & - plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, & - v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, & - sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, & - buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, & - tvp, clw, sig, w0) - else - ! iflag_con == 4 - CALL cv_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, cbmf1, & - plcl1, tnk1, qnk1, gznk1, t1, q1, qs1, u1, v1, gz1, h1, lv1, & - cpn1, p1, ph1, tv1, tp1, tvp1, clw1, iflag, nk, icb, cbmf, & - plcl, tnk, qnk, gznk, t, q, qs, u, v, gz, h, lv, cpn, p, ph, & - tv, tp, tvp, clw, dph) - endif - - ! --- UNDILUTE (ADIABATIC) UPDRAFT / second part : - ! --- FIND THE REST OF THE LIFTED PARCEL TEMPERATURES - ! --- & - ! --- COMPUTE THE PRECIPITATION EFFICIENCIES AND THE - ! --- FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD - ! --- & - ! --- FIND THE LEVEL OF NEUTRAL BUOYANCY - - if (iflag_con == 3) then - CALL cv3_undilute2(klon, ncum, klev, icb, icbs, nk, tnk, qnk, gznk, & - t, q, qs, gz, p, h, tv, lv, pbase, buoybase, plcl, inb, tp, & - tvp, clw, hp, ep, sigp, buoy) !na->klev - else - ! iflag_con == 4 - CALL cv_undilute2(klon, ncum, klev, icb, nk, tnk, qnk, gznk, t, q, & - qs, gz, p, dph, h, tv, lv, inb, inbis, tp, tvp, clw, hp, ep, & - sigp, frac) - endif - - ! --- CLOSURE - - if (iflag_con == 3) then - CALL cv3_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, & - buoy, sig, w0, cape, m) ! na->klev - else - ! iflag_con == 4 - CALL cv_closure(klon, ncum, klev, nk, icb, tv, tvp, p, ph, dph, & - plcl, cpn, iflag, cbmf) - endif - - ! --- MIXING - - if (iflag_con == 3) then - CALL cv3_mixing(klon, ncum, klev, klev, icb, nk, inb, ph, t, q, & - qs, u, v, h, lv, qnk, hp, tv, tvp, ep, clw, m, sig, ment, & - qent, uent, vent, nent, sij, elij, ments, qents) - else - ! iflag_con == 4 - CALL cv_mixing(klon, ncum, klev, icb, nk, inb, inbis, ph, t, q, qs, & - u, v, h, lv, qnk, hp, tv, tvp, ep, clw, cbmf, m, ment, qent, & - uent, vent, nent, sij, elij) - endif - - ! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS - - if (iflag_con == 3) then - CALL cv3_unsat(klon, ncum, klev, klev, icb, inb, t, q, qs, gz, u, & - v, p, ph, th, tv, lv, cpn, ep, sigp, clw, m, ment, elij, delt, & - plcl, mp, qp, up, vp, wt, water, evap, b)! na->klev - else - ! iflag_con == 4 - CALL cv_unsat(klon, ncum, klev, inb, t, q, qs, gz, u, v, p, ph, h, & - lv, ep, sigp, clw, m, ment, elij, iflag, mp, qp, up, vp, wt, & - water, evap) - endif + CALL cv3_compress(klon, klon, ncum, klev, iflag1, nk1, icb1, icbs1, & + plcl1, tnk1, qnk1, gznk1, pbase1, buoybase1, t1, q1, qs1, u1, & + v1, gz1, th1, h1, lv1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, & + sig1, w01, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, pbase, & + buoybase, t, q, qs, u, v, gz, th, h, lv, cpn, p, ph, tv, tp, & + tvp, clw, sig, w0) + + ! UNDILUTE (ADIABATIC) UPDRAFT / second part : + ! FIND THE REST OF THE LIFTED PARCEL TEMPERATURES + ! & + ! COMPUTE THE PRECIPITATION EFFICIENCIES AND THE + ! FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD + ! & + ! FIND THE LEVEL OF NEUTRAL BUOYANCY + + CALL cv3_undilute2(klon, ncum, klev, icb, icbs, nk, tnk, qnk, gznk, & + t, qs, gz, p, h, tv, lv, pbase, buoybase, plcl, inb, tp, & + tvp, clw, hp, ep, sigp, buoy) !na->klev + + ! CLOSURE + + CALL cv3_closure(klon, ncum, klev, icb, inb, pbase, p, ph, tv, & + buoy, sig, w0, cape, m) ! na->klev + + ! MIXING + + CALL cv3_mixing(klon, ncum, klev, klev, icb, nk, inb, t, q, qs, u, & + v, h, lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, & + sij, elij, ments, qents) + + ! UNSATURATED (PRECIPITATING) DOWNDRAFTS + + CALL cv3_unsat(klon, ncum, klev, klev, icb, inb, t, q, qs, gz, u, & + v, p, ph, th, tv, lv, cpn, ep, sigp, clw, m, ment, elij, delt, & + plcl, mp, qp, up, vp, wt, water, evap, b)! na->klev - ! --- YIELD + ! YIELD ! (tendencies, precipitation, variables of interface with other ! processes, etc) - if (iflag_con == 3) then - CALL cv3_yield(klon, ncum, klev, klev, icb, inb, delt, t, q, u, v, & - gz, p, ph, h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, & - wt, water, evap, b, ment, qent, uent, vent, nent, elij, sig, & - tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, upwd, dnwd, & - dnwd0, ma, mike, tls, tps, qcondc, wd)! na->klev - else - ! iflag_con == 4 - CALL cv_yield(klon, ncum, klev, nk, icb, inb, delt, t, q, u, v, gz, & - p, ph, h, hp, lv, cpn, ep, clw, frac, m, mp, qp, up, vp, wt, & - water, evap, ment, qent, uent, vent, nent, elij, tv, tvp, & - iflag, wd, qprime, tprime, precip, cbmf, ft, fq, fu, fv, Ma, & - qcondc) - endif + CALL cv3_yield(klon, ncum, klev, klev, icb, inb, delt, t, q, u, v, & + gz, p, ph, h, hp, lv, cpn, th, ep, clw, m, tp, mp, qp, up, vp, & + wt, water, evap, b, ment, qent, uent, vent, nent, elij, sig, & + tv, tvp, iflag, precip, VPrecip, ft, fq, fu, fv, upwd, dnwd, & + dnwd0, ma, mike, tls, tps, qcondc, wd)! na->klev - ! --- passive tracers + ! passive tracers - if (iflag_con == 3) CALL cv3_tracer(klon, klon, ncum, klev, klev, & - ment, sij, da, phi) + CALL cv3_tracer(klon, ncum, klev, ment, sij, da, phi) - ! --- UNCOMPRESS THE FIELDS + ! UNCOMPRESS THE FIELDS ! set iflag1 = 42 for non convective points do i = 1, klon iflag1(i) = 42 end do - if (iflag_con == 3) then - CALL cv3_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, & - ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, wd, cape, & - da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, & - fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, & - cape1, da1, phi1, mp1) - else - ! iflag_con == 4 - CALL cv_uncompress(idcum(:ncum), iflag, precip, cbmf, ft, fq, fu, & - fv, Ma, qcondc, iflag1, precip1, cbmf1, ft1, fq1, fu1, fv1, & - Ma1, qcondc1) - endif + CALL cv3_uncompress(idcum(:ncum), iflag, precip, VPrecip, sig, w0, & + ft, fq, fu, fv, inb, Ma, upwd, dnwd, dnwd0, qcondc, wd, cape, & + da, phi, mp, iflag1, precip1, VPrecip1, sig1, w01, ft1, fq1, & + fu1, fv1, inb1, Ma1, upwd1, dnwd1, dnwd01, qcondc1, wd1, & + cape1, da1, phi1, mp1) ENDIF ! ncum>0 end SUBROUTINE cv_driver