...
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'. We will see how this macro works.
First, In the macro we first call mvl_flexpart_total_column() to compute the "tcmd" fields with units of "kg m**-2".
...
Code Block | ||
---|---|---|
| ||
title=mvl_flexpart_title(g,0.3,"g m**-2") |
Finally we define the mapview:
Code Block | ||||||
---|---|---|---|---|---|---|
| ||||||
#Define coastlines
coast_grey = mcoast(
map_coastline_thickness : 2,
map_coastline_land_shade : "on",
map_coastline_land_shade_colour : "grey",
map_coastline_sea_shade : "on",
map_coastline_sea_shade_colour : "RGB(0.89,0.89,0.89)",
map_boundaries : "on",
map_boundaries_colour : "black",
map_grid_latitude_increment : 5,
map_grid_longitude_increment : 5
)
#Define geo view
view = geoview(
map_area_definition : "corners",
area : [40,-25,66,9],
coastlines : coast_grey
) |
view with the map in the same way as above and generate the plot:
Code Block | ||
---|---|---|
| ||
plot(view,g,conc_shade,title) |
...