48 real :: rlevmass,summe
50 integer :: iflag, k, kk, kuvz
54 real :: cbmfold, precip, qprime
79 pconv(k) = (akz(kuvz) + bkz(kuvz)*psconv)
80 phconv(kuvz) = (akm(kuvz) + bkm(kuvz)*psconv)
81 dpr(k) = phconv(k) - phconv(kuvz)
82 qsconv(k) =
f_qvsat( pconv(k), tconv(k) )
103 pconv_hpa(k)=pconv(k)/100.
104 phconv_hpa(k)=phconv(k)/100.
106 phconv_hpa(nconvlev+1)=phconv(nconvlev+1)/100.
107 call
convect(nconvlevmax, nconvlev, delt, iflag, &
108 precip, wd, tprime, qprime, cbmf)
113 if (iflag .ne. 1 .and. iflag .ne. 4)
then
121 if (cbmf.le.0..and.cbmfold.le.0.)
then
134 fmassfrac(k,kk) = delt*fmass(k,kk)
135 summe = summe + fmassfrac(k,kk)
137 fmassfrac(k,k)=fmassfrac(k,k) + rlevmass - summe
subroutine convect(ND, NL, DELT, IFLAG, PRECIP, WD, TPRIME, QPRIME, CBMF)
subroutine calcmatrix_ecmwf(lconv, delt, cbmf)
real function f_qvsat(p, t)