#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 ************************************
#
# ------------------------------------------------------------------
# Tags: ODB, Map
# Title: ODB satellite radiance map
# Description: Demonstrates how to plot AMSU-A radiances from ODB
# ------------------------------------------------------------------
# define channel
channel=5 #5,6
# read ODB
amsua_odb = read("amsua.odb")
# filter channel and define plot data
vis = odb_visualiser(
odb_where : "vertco_reference_1@body = " & channel,
odb_data : amsua_odb
)
# find out value range
f = odb_filter(odb_query : "select obsvalue@body as val where vertco_reference_1@body = " & channel,
odb_data : amsua_odb)
val = values(f,"val")
v1 = minvalue(val)
v2 = maxvalue(val)
# symbol plotting (adjusted to the actual min and max)
symbol_fixed = msymb(
legend : "on",
symbol_type : "marker",
symbol_table_mode : "advanced",
symbol_advanced_table_selection_type : "interval",
symbol_advanced_table_min_value : v1,
symbol_advanced_table_max_value : v2,
symbol_advanced_table_interval : 5,
symbol_advanced_table_max_level_colour : "red",
symbol_advanced_table_min_level_colour : "blue",
symbol_advanced_table_colour_direction : "clockwise",
symbol_advanced_table_marker_list : 3,
symbol_advanced_table_height_list : 0.15
)
# coastlines
coast = mcoast(
map_coastline_colour : "RGB(0.5,0.5,0.5)",
map_coastline_resolution : "low",
map_coastline_land_shade : "on",
map_coastline_land_shade_colour : "RGB(0.21,0.21,0.21)",
map_coastline_sea_shade : "on",
map_coastline_sea_shade_colour : "RGB(0.67,0.69,0.69)",
map_grid_colour : "RGB(0.38,0.37,0.37)"
)
# legend with histogram
legend = mlegend(
legend_display_type : "histogram"
)
# title
title = mtext(
text_line_count : 2,
text_line_1 : "Date: 2010-12-22 Time: 00",
text_line_2 : "Sat: METOP-A Sensor: AMSU-A Channel: " & channel & " Param: Tb [K]"
)
# define the output plot file
setoutput(pdf_output(output_name : 'odb_radiance'))
# plot
plot(coast, vis, symbol_fixed, 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
# define channel
channel = 5 #5,6
# read ODB
amsua_odb = mv.read("amsua.odb")
# filter channel and define plot data
vis = mv.odb_visualiser(
odb_where = "vertco_reference_1@body = " + str(channel),
odb_data = amsua_odb
)
# find out value range
f = mv.odb_filter(odb_query = "select obsvalue@body as val where vertco_reference_1@body = " + str(channel),
odb_data = amsua_odb)
val = mv.values(f, "val")
v1 = mv.minvalue(val)
v2 = mv.maxvalue(val)
# symbol plotting (adjusted to the actual min and max)
symbol_fixed = mv.msymb(
legend = "on",
symbol_type = "marker",
symbol_table_mode = "advanced",
symbol_advanced_table_selection_type = "interval",
symbol_advanced_table_min_value = v1,
symbol_advanced_table_max_value = v2,
symbol_advanced_table_interval = 5,
symbol_advanced_table_max_level_colour = "red",
symbol_advanced_table_min_level_colour = "blue",
symbol_advanced_table_colour_direction = "clockwise",
symbol_advanced_table_marker_list = 3,
symbol_advanced_table_height_list = 0.15
)
# coastlines
coast = mv.mcoast(
map_coastline_colour = "RGB(0.5,0.5,0.5)",
map_coastline_resolution = "low",
map_coastline_land_shade = "on",
map_coastline_land_shade_colour = "RGB(0.21,0.21,0.21)",
map_coastline_sea_shade = "on",
map_coastline_sea_shade_colour = "RGB(0.67,0.69,0.69)",
map_grid_colour = "RGB(0.38,0.37,0.37)"
)
# legend with histogram
legend = mv.mlegend(
legend_display_type = "histogram"
)
# title
title = mv.mtext(
text_line_count = 2,
text_line_1 = "Date: 2010-12-22 Time: 00",
text_line_2 = "Sat: METOP-A Sensor: AMSU-A Channel: " + str(channel) + " Param: Tb [K]"
)
# define the output plot file
mv.setoutput(mv.pdf_output(output_name = 'odb_radiance'))
# plot
mv.plot(coast, vis, symbol_fixed, legend, title)
|
|
|