Coefficients of Table 1 in Zeng et al. 
are calculated by the Chou and Neelin weakly nonlinear
scheme  with a reference profile, Tqo3.ref.
One can use a command
run coef.sh, and the coefficients
will be in lw-coef.dat.
The reference values of radiative fluxes are also
calculated. The Harshvardhan longwave radiation (lwdrive.f) is also included
for calculating the reference radiative fluxes and the Green's functions.