Download source and data
Coloured Observation Values 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 BUFR data temp = read('temp.bufr') # filter just the 500 hPa temperature from the obs data (Geopoints) gpt = obsfilter( output : "geopoints", parameter : "airTemperature", level : "single", first_level: 500, data : temp) # convert values to C units gpt = gpt - 273.16 # define coloured value (as number) plotting sym = msymb( legend : "on", symbol_table_mode : "advanced", symbol_format : "(F3.1)", symbol_advanced_table_selection_type : "interval", symbol_advanced_table_interval : 1, symbol_advanced_table_max_level_colour : "red", symbol_advanced_table_min_level_colour : "lavender", symbol_advanced_table_colour_director : "clockwise", symbol_advanced_table_text_font_size : 0.6 ) # shaded land and sea coast = mcoast( map_coastline_land_shade : "on", map_coastline_land_shade_colour : "RGB(0.29,0.29,0.29)", map_coastline_sea_shade : "on", map_coastline_sea_shade_colour : "RGB(0.49,0.49,0.49)", map_coastline_colour : "RGB(0.59,0.59,0.59)", map_grid_latitude_increment : 10, map_grid_longitude_increment : 20, map_grid_colour : "RGB(0.25,0.25,0.25)" ) # adjust the legend legend = mlegend(legend_text_font_size : 0.3) # define title title = mtext(text_lines : "Radiosonde - Temperature at 500 hPa", text_font_size: 0.4) # set the view area view = geoview( map_area_definition : "corners", area : [7,-105,35,-55], coastlines : coast) # define the output plot file setoutput(pdf_output(output_name : 'coloured_obs_values')) # plot the data with the style plot(view, gpt, sym, legend, title)
Coloured Observation Values 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 BUFR data temp = mv.read("temp.bufr") # filter just the 500 hPa temperature from the obs data (Geopoints) gpt = mv.obsfilter( output = "geopoints", parameter = "airTemperature", level = "single", first_level = 500, data = temp) # convert values to C units gpt = gpt - 273.16 # define coloured value (as number) plotting sym = mv.msymb( legend = "on", symbol_table_mode = "advanced", symbol_format = "(F3.1)", symbol_advanced_table_selection_type = "interval", symbol_advanced_table_interval = 1, symbol_advanced_table_max_level_colour = "red", symbol_advanced_table_min_level_colour = "lavender", symbol_advanced_table_colour_director = "clockwise", symbol_advanced_table_text_font_size = 0.6 ) # shaded land and sea coast = mv.mcoast( map_coastline_land_shade = "on", map_coastline_land_shade_colour = "RGB(0.29,0.29,0.29)", map_coastline_sea_shade = "on", map_coastline_sea_shade_colour = "RGB(0.49,0.49,0.49)", map_coastline_colour = "RGB(0.59,0.59,0.59)", map_grid_latitude_increment = 10, map_grid_longitude_increment = 20, map_grid_colour = "RGB(0.25,0.25,0.25)" ) # adjust the legend legend = mv.mlegend(legend_text_font_size = 0.3) # define title title = mv.mtext(text_lines = "Radiosonde - Temperature at 500 hPa", text_font_size = 0.4) # set the view area view = mv.geoview( map_area_definition = "corners", area = [7,-105,35,-55], coastlines = coast) # define the output plot file mv.setoutput(mv.pdf_output(output_name = 'coloured_obs_values')) # plot the data with the style mv.plot(view, gpt, sym, legend, title)