Page History

...

In the macro first we call mvl_flexpart_total_column() to compute the "tcmd" fields with units of "kg m**-2".

dIn="result_fwd/"
inFile=dIn  & "conc_s001.grib"

#Compute the total column integrated mass
g=mvl_flexpart_total_column(source: inFile,1)

Next, we define the contouring. The "tcmd" 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:

...

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: "g m**-2")   

Finally we define the mapview:

...