Page History

...

dIn="result_fwd_conc/"
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,-1)

Next, we define the contouring definition. The units we used here are is 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.

...

Computing and plotting total column mass

Using the simulation results we will compute and plot the total column integrated mass. The macro to visualise the concentration fields on a given level use is 'plot_total.mv'.

First, we define the level (8000 m) and the parameter ("mdc") we want to plot. Then we call call mvl_flexpart_readtotal_hlcolumn() to filter the data into a fieldsetperform the computations. The result is a fieldset containing "tcmd" fields with units of "kg m**-2".

dIn="result_fwd_conc/"
inFile=dIn  & "conc_s001.grib"
lev=8000
par="mdc"
 
#Read fields on the given height level
#Compute the total column integrated mass
g=mvl_flexpart_readtotal_hlcolumn(inFile,par,lev,-1,-1)

Next, we define the contouring. The "tcmd" fields are autocratically scaled into "ng m**-2" for contouring 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*1000

Then specify the contour levels

cont_list=[0.0001,0.001,0.005,0.01,0.02,0.05,0.1,0.5]

{code}