50 real :: rlevmass,summe
52 integer :: iflag, k, kk, kuvz
56 real :: cbmfold, precip, qprime
79 phconv(kuvz) = 0.5*(pconv(kuvz)+pconv(k))
80 dpr(k) = phconv(k) - phconv(kuvz)
81 qsconv(k) =
f_qvsat( pconv(k), tconv(k) )
100 pconv_hpa(k)=pconv(k)/100.
101 phconv_hpa(k)=phconv(k)/100.
103 phconv_hpa(nconvlev+1)=phconv(nconvlev+1)/100.
104 call
convect(nconvlevmax, nconvlev, delt, iflag, &
105 precip, wd, tprime, qprime, cbmf)
110 if (iflag .ne. 1 .and. iflag .ne. 4)
then
118 if (cbmf.le.0..and.cbmfold.le.0.)
then
131 fmassfrac(k,kk) = delt*fmass(k,kk)
132 summe = summe + fmassfrac(k,kk)
134 fmassfrac(k,k)=fmassfrac(k,k) + rlevmass - summe
subroutine calcmatrix_gfs(lconv, delt, cbmf)
subroutine convect(ND, NL, DELT, IFLAG, PRECIP, WD, TPRIME, QPRIME, CBMF)
real function f_qvsat(p, t)