# 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)
Metview's documentation is now on readthedocs!
Please go to https://metview.readthedocs.io/en/latest/
Overview
Community Forums
Content Tools