Download Polyphemus 1.8 User's Guide - GForge

2.7. CHECKING RESULTS
41
Please note that 4DVar assimilation scheme only works with RACM chemistry module.
As for drivers, BaseDriver is the simplest and the most used of them. The other drivers
available are:
• PlumeDriver: for Gaussian plume model (with or without aerosol species).
• PuffDriver: for Gaussian puff model (with or without aerosol species).
• Data assimilation drivers: OptimalInterpolationDriver (optimal interpolation), EnKFDriver
(ensemble Kalman filter), RRSQRTDriver (reduced rank square root filter), FourDimVarDriver
(4D-Var), to be associated with Polair3DAssimConc model.
• Drivers for the verification of adjoint coding in variational assimilation: AdjointDriver,
GradientDriver, Gradient4DVarDriver.
2.7
Checking Results
It is highly recommended to check the fields generated by Polyphemus programs: meteorological
fields, deposition velocities, output concentrations, . . .
First, you should check the size of the binary files. Be it with preprocessing or processing
programs, results are saved as floating point numbers with single precision. As we will express
file size in bytes, a factor 4 will appear below.
2.7.1
Checking the output file size of preprocessing programs
In what follows:
• take Nz , Ny , Nx and Nt from the section [domain] of the configuration files;
• Ndays is the number of preprocessed days and
• Ntotal is the total number of preprocessing time steps: Ntotal = Nt × Ndays
Here are some indications on what to check depending on the program.
Ground Data
As there are 14 categories in the GLCF description, the size of the binary file obtained from
luc-glcf should be Nx × Ny × 4 × 14. As of luc-usgs, it should be Nx × Ny × 4 × 24.
Meteorological fields
For programs meteo, attenuation, Kz, MM5-meteo, WRF-meteo and Kz TM, related file size :
• for 2D fields is Ntotal × Ny × Nx × 4;
• for 3D fields is Ntotal × Nz × Ny × Nx × 4;
• except for the fields Kz and Kz TM, Ntotal × (Nz + 1) × Ny × Nx × 4;
• for the field MeridionalWind, Ntotal × Nz × (Ny + 1) × Nx × 4;
• and for the field ZonalWind, Ntotal × Nz × Ny × (Nx + 1) × 4.