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)