Page History

...

In the macro first we define the level (8000 m) and the parameter ("feflux" = eastward flux) we want to plot. Then we call the Macro Library Function mvl_flexpart_read_hlfilter() to extract the data.

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

Next, we define the contouring definition. 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.

...

Next, we build the title with mvlflexpart_flexpartbuild_title(). Please note that we need to explicitly specify the plotting units!

title=mvlflexpart_flexpartbuild_title(data: g,
    fontsize: 0.3, 
    units: "ng m**-2 s**-1")  

Finally we define the view with the map:

...