Download source and data
BUFR Temp 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 TEMP bufr file bd = read("temp.bufr") # define observation plotting - selecting level 250 hPa obsp = mobs( obs_distance_apart : 1, obs_level : 250, obs_size : 0.35, obs_ring_size : 0.4, obs_station_ring : "off" ) # define land-sea shading coast = mcoast( map_coastline_land_shade : "on", map_coastline_land_shade_colour : "grey", map_coastline_sea_shade : "on", map_coastline_sea_shade_colour : "RGB(0.8944,0.9086,0.9330)", map_grid_colour : "charcoal", map_grid_longitude_increment : 10 ) # define map projection view = geoview( map_projection : "polar_stereographic", map_area_definition : "corners", area : [17.44,-73.61,43.74,33.84], map_vertical_longitude : -40, coastlines : coast ) # add title title = mtext(text_lines : "TEMP 250 hpa", text_font_size : 0.4) # define the output plot file setoutput(pdf_output(output_name : 'temp_map')) # generate plot plot(view, bd, obsp, title)
BUFR Temp 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 TEMP bufr file bd = mv.read("temp.bufr") # define observation plotting - selecting level 250 hPa obsp = mv.mobs( obs_distance_apart = 1, obs_level = 250, obs_size = 0.35, obs_ring_size = 0.4, obs_station_ring = "off" ) # define land-sea shading coast = mv.mcoast( map_coastline_land_shade = "on", map_coastline_land_shade_colour = "grey", map_coastline_sea_shade = "on", map_coastline_sea_shade_colour = "RGB(0.8944,0.9086,0.9330)", map_grid_colour = "charcoal", map_grid_longitude_increment = 10 ) # define map projection view = mv.geoview( map_projection = "polar_stereographic", map_area_definition = "corners", area = [17.44,-73.61,43.74,33.84], map_vertical_longitude = -40, coastlines = coast ) # add title title = mv.mtext(text_lines = "TEMP 250 hpa", text_font_size = 0.4) # define the output plot file mv.setoutput(mv.pdf_output(output_name = 'temp_map')) # generate plot mv.plot(view, bd, obsp, title)