The output (reference results) consists of the following:
Tips and instructions on reading and using the output files
The output for each CIRC case consists of 3 files:
Contains fluxes at the boundaries of the atmospheric column (downwelling surface and upwelling TOA) at 1 cm-1 resolution from LBLRTM.
Contains broadband flux and heating rate profiles from LBLRTM
Contains fluxes at the boundaries of the atmospheric column (downwelling surface and upwelling TOA) at 1 cm-1 resolution from CHARTS
The SW broadband fluxes are given at the beginning of the CHARTS 1 cm-1 output files. Note that the LW flux and heating rate profiles are inverted (values start at TOA and proceed to surface) in contrast to the input profiles where surface values are placed first.
Here is a sample FORTRAN code to read the contents of these files. The output files themselves can be downloaded from this page.
Description of reference observations and calculations
The longwave and shortwave surface irradiance measurements utilized for the CIRC SGP and NSA cases are 5-minute averages of the irradiance values provided by ARM Best Estimate Radiative Flux VAP. The 5-minute averaging window is centered at the time provided for each case. The observed TOA irradiance values are from the temporally closest CERES measurement for Case 4 (NSA) and the temporally closest GOES broadband fluxes inferred from a radiance to flux and narrowband to broadband conversion algorithm, for cases 3 (SGP) and 7 (PYE). For the remaining SGP cases the weighted average of the two GOES-inferred fluxes bracketing the nominal BBHRP calculation time were used. The observed broadband irradiances can be found in the SW and LW summary tables.
The reference radiative transfer calculations are performed with CHARTS v. 4.04 and LBLRTM v. 11.1 for the SW and LW part of the spectrum, respectively. The continuum model used was MT_CKD_2.0 and the line parameter file was AER_V_2.0. LBLRTM is run using 3 streams (i.e, 3 angles) while CHARTS is run using 16 streams with delta-M scaling and in our configuration provides fluxes only at the top and bottom boundary of the atmospheric column. The surface albedo is resolved at 1 cm-1 (the same as the one provided to the participants) and is linearly interpolated to the wavenumber of the calculation. The cloud optical properties (extinction coefficient, single-scattering albedo and asymmetry parameter) for the two cloudy cases are calculated from the retrieved layer water contents and effective radii using Mie theory and an assumed gamma droplet size distribution. Aerosols are treated similarly to BBHRP, i.e., the spectral aerosol optical thickness is obtained from the Angstrom relationship while g and SSA are spectrally invariant; a Henyey-Greenstein phase function is assumed. Output is averaged at 1 cm-1 intervals. In LBLRTM the cloud absorption optical thickness is also obtained from Mie calculations and cloud scattering is not accounted for (if your model accounts for LW scattering you may want to turn it off, if possible, for a more consistent comparison), as is any aerosol effect. The surface emissivity is set to unity across the LW spectrum .
For questions or problems reading/interpreting the output dataset please contact Lazaros Oreopoulos.