# Metview Macro

#  **************************** LICENSE START ***********************************
# 
#  Copyright 2016 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 ************************************

# ------------------------------------------------------------------
# Description: Demonstrates how to show observations from a geopoints
#              file using symbol objects.
# ------------------------------------------------------------------

# read the input geopoints file
my_data = read("airep.geo")

# set up the data symbol and colour
my_symbol = msymb(symbol_type                            : "marker",
                  symbol_marker_index                    : 15,
                  symbol_table_mode                      : "advanced",
                  symbol_advanced_table_selection_type   : "interval",
                  symbol_advanced_table_colour_direction : "clockwise",
                  symbol_advanced_table_max_level_colour : "red",
                  symbol_advanced_table_min_level_colour : "blue",
                  symbol_advanced_table_interval         : 5.00,
                  legend                                 : "on")

# set up the contour so that observed values are mapped into a colour list
my_contour = mcont(contour_grid_value_plot          : "on",
                   contour_grid_value_plot_type     : "both",
                   contour_grid_value_height        : 0.30,
                   contour_grid_value_marker_height : 0.20,
                   contour_level_selection_type     : "level_list",
                   contour_level_list               : [-20,-10,-5,-2.5,-1,-0.5,0,0.5,1,2.5,5,10,20],
                   contour_line_colour              : "grey",
                   contour_line_thickness           : 2,
                   contour_highlight                : "off",
                   contour_label                    : "off",
                   contour_grid_value_lon_frequency : 10,
                   contour_grid_value_lat_frequency : 10,
                   legend                           : "off")

# legend type is histogram
my_legend = mlegend(legend_text_colour  : "navy",
                    legend_display_type : "histogram")

# set up the coastlines
my_coast = mcoast(map_coastline_land_shade        : "on",
                  map_coastline_land_shade_colour : "cream",
                  map_grid                        : "on",
                  map_grid_colour                 : "grey")

# set up the title
my_title = mtext(text_mode           : "positional",
                 text_font_size      : 0.60,
                 text_justification  : "left",
                 text_border         : "off",
                 text_lines          : ["Monitoring of AIREP data"],
                 text_box_x_length   : 20.00,
                 text_box_y_length   : 2.50,
                 text_box_x_position : 1.50,
                 text_box_y_position : 16.50,
                 text_colour         : "navy")


# define the output plot file
setoutput(pdf_output(output_name : 'air_report'))

# plot the data onto the map
plot(my_coast, my_data,my_symbol,my_title,my_legend)