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)

#  **************************** 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)