# Metview Macro
# **************************** LICENSE START ***********************************
#
# Copyright 2015 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: contour
# Title: Interval Shading Contour
# Description: Demonstrates how to define a shading contour using an interval range.
# See also:
# https://software.ecmwf.int/magics/Contour+examples+contour9
# ------------------------------------------------------------------
# read the input grib file
my_data = read("2m_temperature.grib")
# set up the geographical area
my_view = geoview(map_area_definition : "CORNERS",
map_projection : "CYLINDRICAL",
area : [35,-15,70,35])
# set up the coastlines
my_coast = mcoast(map_coastline_land_shade : "ON",
map_coastline_land_shade_colour : "CREAM")
# set up the contour using an interval range
my_contour = mcont(contour_level_selection_type : "INTERVAL",
contour_interval : 2.00,
contour_shade_method : "AREA_FILL",
contour_shade : "ON",
contour_shade_max_level : 42.00,
contour_shade_min_level : -32.00,
contour_highlight_thickness : 5,
contour_line_thickness : 3,
contour_label : "OFF",
contour_highlight_frequency : 100,
legend : "ON",
contour_hilo : "OFF",
contour : "OFF",
contour_shade_colour_method : "LIST",
contour_shade_colour_list : ["rgb(0.3,0.3,0.3)","rgb(0.4,0.4,0.4)","rgb(0.5,0.5,0.5)", "rgb(0.6,0.6,0.6)","rgb(0.7,0.7,0.7)",
"rgb(0.8,0.8,0.8)","rgb(0.35,0,0.6)", "rgb(0.5,0,0.9)", "rgb(0.6,0.2,1)", "rgb(0.75,0.4,1)",
"rgb(0.85,0.6,1)", "rgb(0,0,0.75)", "rgb(0,0,1)", "rgb(0.2,0.4,1)", "rgb(0.4,0.7,1)",
"rgb(0.6,0.9,1)", "rgb(0,0.55,0.19)","rgb(0.15,0.75,0.1)","rgb(0.5,0.85,0)", "rgb(0.65,0.95,0)",
"rgb(0.8,1,0.2)", "rgb(0.65,0.65,0)","rgb(0.8,0.8,0)", "rgb(0.92,0.92,0)","rgb(1,1,0)",
"rgb(1,1,0.6)", "rgb(0.85,0.45,0)","rgb(1,0.5,0)", "rgb(1,0.62,0)", "rgb(1,0.74,0)",
"rgb(1,0.85,0)", "rgb(0.6,0,0)", "rgb(0.8,0,0)", "rgb(1,0,0)", "rgb(1,0.4,0.4)",
"rgb(1,0.6,0.6)", "rgb(1,0.75,0.75)"])
# set-up the title
my_title = mtext(text_font_size : 0.50,
text_lines : ["<font size='0.7'>2m temperature (Range: -32 .. 42)</font>",
"<magics_title/>",
""],
text_colour : "CHARCOAL")
# define the output media
to_psfile = ps_output
(
output_name : "plot" # extension is added automatically
)
# check the runmode and decide which media to putput the plot to
mode = runmode()
if (mode = "execute") then setoutput(to_psfile)
else if (mode = "batch") then setoutput(to_psfile)
else if (mode = "visualise") then print('screen')
else fail("Only execute, batch and visualise allowed")
end if
# plot the data onto the map
plot(my_view,my_coast,my_data,my_contour,my_title)
