Page History

...

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

Next, we define the contouring. The "feflux" fields are automatically scaled into "ng m**-2" for contouring (see here for details)  but with the current value range it is better to use "g m**-2" units. To achieve it we simply multiply the "tcmd" fieldset with 1000:

g=g*1E12

The contour definition itself goes like this:

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

...