Page History

...

title=mvl_flexpart_title(g,0.3,"ng m**-3")

Finally we define the view with the map:

#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
	)

...

Computing and plotting total column mass

Using the simulation results same FLEXPART output as above we will compute and plot the total column integrated mass. The macro to use is 'plot_total.mv'.

...

Next, we define the contouring. The "tcmd" fields are autocratically 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

Then specify the contour levelsThe contour definition itself goes like this:

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

#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"
    )

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 mvl_flexpart_title(). Please note that we need to explicitly specify the plotting units!

title=mvl_flexpart_title(g,0.3,"g m**-2")

Finally we define the view with the map in the same way as above FLEXPART forward 1

and generate the plot:

plot(view,g,conc_shade,title)

Having run the macro we get a plot like this (after navigating to step 39h):