...
In the macro first we call flexpart_total_column() to compute the "tcmd" fields with units of "kg m**-2".:
| Code Block | ||
|---|---|---|
| ||
dIn="result_fwd/" inFile=dIn & "conc_s001.grib" #Compute the total column integrated mass g=flexpart_total_column(source: inFile, param: "mdc") |
Next, we define the contouring. The "tcmd" fields are automatically scaled into "ng have the units of "kg m**-2" for contouring (see here for details) but with the current value range it is better to use "g m**-2" units would better fit for contouring. To achieve it we simply multiply the "tcmd" fieldset with 1000:
...
| Code Block | ||||||
|---|---|---|---|---|---|---|
| ||||||
cont_list=[0.00001,0.0001,0.0005,0.001,0.002,0.005,0.01,0.05]
#Define contour shading
conc_shade = mcont(
legend : "on",
contour : "off",
contour_level_selection_type : "level_list",
contour_level_list : cont_list,
contour_label : "off",
contour_shade : "on",
contour_shade_method : "area_fill",
contour_shade_max_level_colour : "red",
contour_shade_min_level_colour : "RGB(0.14,0.37,0.86)",
contour_shade_colour_direction : "clockwise",
contour_method : "linear",
grib_scaling_of_derived_fields : "off"
)
|
...
Next, we build the title with flexpart_build_title(). Please note that we need to explicitly specify the plotting units!
...