For preparations and running the simulation needed for this tutorial click here ...
To start these tutorial please enter folder 'forward'.
The total column integrated mas is not an output of FLEXPART so we need to compute it. The macro to do the computation and plot the resulting fields is 'plot_total.mv'.
In the macro first we call 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=flexpart_total_column(source: inFile)
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:
g=g*1000
The contour definition itself goes like this:
Please note that by multiplying our fieldset by 1000 it became a derived fieldset. Since we do not want the automatic contour scaling to happen we need to set the contouring parameter grib_scaling_of_derived_fields
to "off".
Next, we build the title with flexpart_build_title(). Please note that we need to explicitly specify the plotting units!
title=flexpart_build_title(data: g,fontsize: 0.3,units: "g m**-2")
Finally we define the mapview:
and generate the plot:
plot(view,g,conc_shade,title)
Having run the macro we will get a plot like this (after navigating to step 39h):