# 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: Coastline
# Title: Boundaries, Cities and Rivers
# Description: Demonstrates how to show main cities, rivers, borders
# and admnistrative boundaries.
# See also:
# https://software.ecmwf.int/magics/Coastlines+examples+coastlines2
# ------------------------------------------------------------------
# set up the geographical view
my_view = geoview(map_area_definition : "CORNERS",
map_projection : "CYLINDRICAL",
area : [48.00,-3.00,55.00,10.00])
# set up the coastlines so that cities, rivers, borders and admnistrative boundaries are shown
my_coast = mcoast(map_administrative_boundaries_colour : "ORANGE",
map_boundaries : "ON",
map_coastline_resolution : "HIGH",
map_coastline_land_shade_colour : "CREAM",
map_cities : "ON",
map_grid : "OFF",
map_boundaries_colour : "RED",
map_grid_colour : "TAN",
map_rivers : "ON",
map_administrative_boundaries_countries_list : ["FRA","DEU","GBR"],
map_coastline_land_shade : "ON",
map_coastline_colour : "TAN",
map_administrative_boundaries : "ON")
my_text = mtext(text_font_size : 0.80,
text_lines : ["Administrative boundaries, cities and rivers"],
text_justification : "LEFT",
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 coastlines data onto the map
plot(my_view,my_coast,my_text)