...
We will simulate a volcano eruption by releasing of some SO2 from the Icelandic volcano Eyjafjallajökull.
The simulation itself is defined by the 'fwd_conc' FLEXPART Run icon and the 'rel_volcano' FLEXPART Release icon, respectively. Both these are encompassed in a single macro called 'fw_cond.mv'. For simplicity will use this macro to examine the settings in detail.
...
- conc_s001.grib is GRIB file containing the gridded concentration fields
- tr_r1.csv is CSV file containing the plume trajectories
| Info |
|---|
Please note that these are not the original outputs form FLEXTRA but were converted to formats more suitable for use in Metview. For details about the FLEXPART outputs please click here. |
Visualising gridded fields
...
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_hl() to filter extract the data into a fieldsetto plot.
| Code Block |
|---|
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) |
...
| Code Block | ||
|---|---|---|
| ||
#The contour levels
cont_list=[1,10,50,100,150,200,250,500,750,1000,2000,6000]
#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"
) |
Next, we build the title with with mvl_flexpart_title(). Please note that we need to explicitly specify the plotting units!
| Code Block | ||
|---|---|---|
| ||
title=mvl_flexpart_title(g,0.3,"ng m**-3") |
...