...
In the macro first we define the level (8000 m) and the parameter ("feflux" = eastward flux) we want to plot. Then we call the Macro Library Function flexpart_filter() to extract the data.
...
Next, we define the contouring definition. The units we use here are ng m**-2 s**-1 because the native flux units (kg m**-2 s**-1) are automatically scaled by the plotting library (see details about the this scaling for various FLEXPART GRIB fields here.
| Code Block | ||||||
|---|---|---|---|---|---|---|
| ||||||
#The contour levels
cont_list=[1,10,50,100,150,200,250,500,750,1000,2000,4000]
#Define contour shading
flux_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"
) |
...
| Anchor | ||||
|---|---|---|---|---|
|
Finally we define the map view with the map:
| 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 ) |
...