include file: prognosticate.inc 1,22

fqr=0
fqv=0
frc=0
fccn=0
fth=0
fu=0
fv=0
fug=0
fnc=0
fvg=0
fsmoke=0
fnr=0
ftke=0
fqc=0


kqr=>kh
kqv=>kh
krc=>kh
kccn=>kh
kth=>kh
ku=>kh
kv=>kh
kug=>kh
knc=>kh
kvg=>kh
ksmoke=>kh
knr=>kh
kqc=>kh


dqrdt=0
dqvdt=0
drcdt=0
dccndt=0
dthdt=0
dudt=0
dvdt=0
dugdt=0
dncdt=0
dvgdt=0
dsmokedt=0
dnrdt=0
dtkedt=0
dqcdt=0


call ad_flux(fth,th)
call ad_flux(fu,u)
call ad_flux(fv,v)
call ad_flux(fug,ug)
call ad_flux(fvg,vg)


call ad_flux(fqr,qr,qr_sed)
call ad_flux(fqv,qv,qv_sed)
call ad_flux(frc,rc,rc_sed)
call ad_flux(fccn,ccn,ccn_sed)
call ad_flux(fnc,nc,nc_sed)
call ad_flux(fsmoke,smoke,smoke_sed)
call ad_flux(fnr,nr,nr_sed)
call ad_flux(fqc,qc,qc_sed)


call ad_flux_f(ftke,tke)


if (adcorr.eq.1) call mono_corr(fqr,qr)
if (adcorr.eq.1) call mono_corr(fqv,qv)
if (adcorr.eq.1) call mono_corr(frc,rc)
if (adcorr.eq.1) call mono_corr(fccn,ccn)
if (adcorr.eq.1) call mono_corr(fnc,nc)
if (adcorr.eq.1) call mono_corr(fsmoke,smoke)
if (adcorr.eq.1) call mono_corr(fnr,nr)
if (adcorr.eq.1) call mono_corr(fqc,qc)