Page History

...

In the macro first we define the level (8000 m) and the parameter ("mdc") 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="mdc"
 
#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 used 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.

...

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**-3") 

Finally we define the view with the map:

...