Page History

...

dIn="result_fwd/"
inFile=dIn  & "conc_s001.grib"
lev=8000
par="mdc"
 
#Read fields on the given height level
g=flexpart_filter(source: inFile,
                  param: par,
                  levType: "hl", 
                  level: lev)

Next, we define the contouring. The "mdc" fields have the units of "kg m**-3" but with the current value range "ng m**-3" units would better fit for contouring. To achieve it we simply multiply the "mdc" fieldset with 10¹²:

g=g*1E12

The contour definition itself goes like this:

Next, we define the contouring. The units we use here are ng m**-3 because for parameter "mdc" the native units (kg m**-3) are automatically scaled by the plotting library (see details about the this scaling for various FLEXPART GRIB fields here),

...