![](/wiki/download/attachments/29328261/coastlines2-thumb.png)
# Metview Macro # **************************** LICENSE START *********************************** # # Copyright 2014 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: City,River,Coastline,Boundary # Title: Boundaries, Cities and Rivers # Description: Demonstrates how to show main cities, rivers, borders # and admnistrative boundaries. # See also: # https://software.ecmwf.int/wiki/display/MAGP/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)