Download source and data
Cross Section Parallel Wind Component with Map Example
#Metview Macro
# **************************** LICENSE START ***********************************
#
# Copyright 2019 ECMWF. This software is distributed under the terms
# of the Apache License version 2.0. In applying this license, ECMWF does not
# waive the privileges and immunities granted to it by virtue of its status as
# an Intergovernmental Organization or submit itself to any jurisdiction.
#
# ***************************** LICENSE END ************************************
#
# read pressure level data - superstorm Sandy
g = read("sandy_pl.grib")
# read wind fields and z500
f_uv = read(data:g , param: ["u","v"])
f_z = read(data:g , param: "z", levelist: 500)
# ---------------------------
# Cross section
# ----------------------------
# define cross section line
line = [30.30,-79.83,36.95,-63.92]
# define cross section view wind projected to the cross section plane
xs = mxsectview(
bottom_level : 1000,
top_level : 100,
line : line,
wind_parallel : "on",
wind_perpendicular : "off",
wind_intensity : "off",
vertical_scaling : "log"
)
# define contour shading for normal wind component
wind_shade = mcont(
legend : "on",
contour : "off",
contour_level_selection_type : "interval",
contour_shade_max_level : 32,
contour_shade_min_level : -32,
contour_interval : 4,
contour_label : "off",
contour_shade : "on",
contour_shade_colour_method : "palette",
contour_shade_method : "area_fill",
contour_shade_palette_name : "m_brown_purple_16"
)
# ---------------------------------------------
# Map view with z500 and cross section line
# ---------------------------------------------
# define coastlines
coast = mcoast(
map_coastline_colour : "RGB(0.4449,0.4414,0.4414)",
map_coastline_resolution : "low",
map_coastline_land_shade : "on",
map_coastline_land_shade_colour : "RGB(0.5333,0.5333,0.5333)",
map_coastline_sea_shade : "on",
map_coastline_sea_shade_colour : "RGB(0.7765,0.8177,0.8941)",
map_boundaries : "on",
map_boundaries_colour : "mustard",
map_boundaries_thickness : 2,
map_grid_colour : "RGB(0.2627,0.2627,0.2627)",
map_grid_latitude_increment : 5,
map_grid_longitude_increment : 5
)
# define map view
map = geoview(
map_area_definition : "corners",
area : [25,-84,40,-60],
coastlines : coast,
subpage_y_lenght : 75
)
# define contoring for z500
cont_z = mcont(
contour_line_thickness : 2,
contour_line_colour : "black",
contour_highlight : "off",
contour_level_selection_type : "interval",
contour_interval : 4,
grib_scaling_of_derived_fields : "on"
)
# plot title
map_title = mtext(text_lines : "Z 500 hpa")
# cross section line visualiser
vis_line = input_visualiser(
input_plot_type : "geo_points",
input_longitude_values : [line[2],line[4]],
input_latitude_values : [line[1],line[3]]
)
graph_line = mgraph(
graph_line_colour : "red",
graph_line_thickness : 4
)
# ---------------------------------------------
# Define layout
# --------------------------------------------
map_page = plot_page(
bottom: 30,
view : map
)
xs_page = plot_page(
top : 30,
view : xs
)
dw = plot_superpage(
pages : [map_page, xs_page]
)
# ---------------------------------------------
# Generate plot
# --------------------------------------------
# define the output plot file
setoutput(pdf_output(output_name : 'cross_section_wind_parallel'))
# generate plot
plot(dw[1], f_z, cont_z, vis_line, graph_line, map_title)
plot(dw[2], f_uv, wind_shade)
Cross Section Parallel Wind Component with Map Example
# **************************** LICENSE START ***********************************
#
# Copyright 2019 ECMWF. This software is distributed under the terms
# of the Apache License version 2.0. In applying this license, ECMWF does not
# waive the privileges and immunities granted to it by virtue of its status as
# an Intergovernmental Organization or submit itself to any jurisdiction.
#
# ***************************** LICENSE END ************************************
#
import metview as mv
# read pressure level data - superstorm Sandy
g = mv.read("sandy_pl.grib")
# read wind fields and z500
f_uv = mv.read(data=g , param=["u","v"])
f_z = mv.read(data=g , param="z", levelist=500)
# ---------------------------
# Cross section
# ----------------------------
# define cross section line
line = [30.30,-79.83,36.95,-63.92]
# define cross section view wind projected to the cross section plane
xs = mv.mxsectview(
bottom_level = 1000,
top_level = 100,
line = line,
wind_parallel = "on",
wind_perpendicular = "off",
wind_intensity = "off",
vertical_scaling = "log"
)
# define contour shading for normal wind component
wind_shade = mv.mcont(
legend = "on",
contour = "off",
contour_level_selection_type = "interval",
contour_shade_max_level = 32,
contour_shade_min_level = -32,
contour_interval = 4,
contour_label = "off",
contour_shade = "on",
contour_shade_colour_method = "palette",
contour_shade_method = "area_fill",
contour_shade_palette_name = "m_brown_purple_16"
)
# ---------------------------------------------
# Map view with z500 and cross section line
# ---------------------------------------------
# define coastlines
coast = mv.mcoast(
map_coastline_colour = "RGB(0.4449,0.4414,0.4414)",
map_coastline_resolution = "low",
map_coastline_land_shade = "on",
map_coastline_land_shade_colour = "RGB(0.5333,0.5333,0.5333)",
map_coastline_sea_shade = "on",
map_coastline_sea_shade_colour = "RGB(0.7765,0.8177,0.8941)",
map_boundaries = "on",
map_boundaries_colour = "mustard",
map_boundaries_thickness = 2,
map_grid_colour = "RGB(0.2627,0.2627,0.2627)",
map_grid_latitude_increment = 5,
map_grid_longitude_increment = 5
)
# define map view
map = mv.geoview(
map_area_definition = "corners",
area = [25,-84,40,-60],
coastlines = coast,
subpage_y_lenght = 75
)
# define contoring for z500
cont_z = mv.mcont(
contour_line_thickness = 2,
contour_line_colour = "black",
contour_highlight = "off",
contour_level_selection_type = "interval",
contour_interval = 4,
grib_scaling_of_derived_fields = "on"
)
# plot title
map_title = mv.mtext(text_lines = "Z 500 hpa")
# cross section line visualiser
vis_line = mv.input_visualiser(
input_plot_type = "geo_points",
input_longitude_values = [line[1],line[3]],
input_latitude_values = [line[0],line[2]]
)
graph_line = mv.mgraph(
graph_line_colour = "red",
graph_line_thickness = 4
)
# ---------------------------------------------
# Define layout
# --------------------------------------------
map_page = mv.plot_page(
bottom= 30,
view = map
)
xs_page = mv.plot_page(
top = 30,
view = xs
)
dw = mv.plot_superpage(
pages = [map_page, xs_page]
)
# ---------------------------------------------
# Generate plot
# --------------------------------------------
# define the output plot file
mv.setoutput(mv.pdf_output(output_name = 'cross_section_wind_parallel'))
# generate plot
mv.plot(dw[0], f_z, cont_z, vis_line, graph_line, map_title,
dw[1], f_uv, wind_shade)