#Metview Macro

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

# read synop bufr
bf = read("synop.bufr")

# extract station ids. The problem with SYNOP station id is that it is the
# combination of two bufr parameters: blockNumber abd stationNumber. 
# We use NCOLS geopoints as an output beacuse the combined station id is automatically
# extracted as a column. So the actual parameter set for the filter is not relevant
# at all (it just has to exist for the messages).
gpt = obsfilter(
    output    : "ncols",
    parameter : "year",
    data      : bf, 
    location_filter: "area",
    area: [45,-8,55,5]
    )

# read station ids into a list
id_lst = stnids(gpt)

# define station id text plotting
sym_txt = msymb(
    legend                : "off",
    symbol_type           : "text",
    symbol_colour         : "black",
    symbol_text_list      : id_lst,
    symbol_text_position  : "top",
    symbol_text_font_size : 0.5,
    symbol_text_blanking: "on"
    )

#  define station location plotting
sym_loc = msymb(
    legend              : "off",
    symbol_type         : "marker",
    symbol_colour       : "coral",
    symbol_height       : 0.4,
    symbol_marker_index : 15
    )

# define coastlines
coast = mcoast(
    map_coastline_colour            : "charcoal",
    map_coastline_thickness         : 2,
    map_coastline_land_shade        : "on",
    map_coastline_land_shade_colour : "RGB(0.4669,0.6782,0.5373)",
    map_coastline_sea_shade         : "on",
    map_coastline_sea_shade_colour  : "RGB(0.7475,0.8504,0.9466)",
    map_boundaries                  : "on",
    map_boundaries_colour           : "rust",
    map_boundaries_thickness        : 2,
    map_grid_colour                 : "charcoal",
    map_grid_latitude_increment     : 2,
    map_grid_longitude_increment    : 2
    )

# define map view
view = geoview(
    map_area_definition : "corners",
    area                : [49,-2.5,54,4],
    coastlines          : coast
    )

# define title
title = mtext(
        text_line_1: "SYNOP Station Identifiers", 
        text_font_size: 0.5
        )

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

# generate plot
plot(view, gpt, sym_txt, sym_loc, title)

"""
Plotting Station Ids from SYNOP BUFR
"""

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

import metview as mv

# read synop bufr
bf = mv.read("synop.bufr")

# extract station ids. The problem with SYNOP station id is that it is the
# combination of two bufr parameters= blockNumber abd stationNumber.
# We use NCOLS geopoints as an output beacuse the combined station id is automatically
# extracted as a column. So the actual parameter set for the filter is not relevant
# at all (it just has to exist for the messages).
gpt = mv.obsfilter(
    output="ncols",
    parameter="year",
    data=bf,
    location_filter="area",
    area=[45, -8, 55, 5],
)

# read station ids into a list
id_lst = mv.stnids(gpt)

# define station id text plotting
sym_txt = mv.msymb(
    legend="off",
    symbol_type="text",
    symbol_colour="black",
    symbol_text_list=id_lst,
    symbol_text_position="top",
    symbol_text_font_size=0.5,
    symbol_text_blanking="on",
)

#  define station location plotting
sym_loc = mv.msymb(
    legend="off",
    symbol_type="marker",
    symbol_colour="coral",
    symbol_height=0.4,
    symbol_marker_index=15,
)

# define coastlines
coast = mv.mcoast(
    map_coastline_colour="charcoal",
    map_coastline_thickness=2,
    map_coastline_land_shade="on",
    map_coastline_land_shade_colour="RGB(0.4669,0.6782,0.5373)",
    map_coastline_sea_shade="on",
    map_coastline_sea_shade_colour="RGB(0.7475,0.8504,0.9466)",
    map_boundaries="on",
    map_boundaries_colour="rust",
    map_boundaries_thickness=2,
    map_grid_colour="charcoal",
    map_grid_latitude_increment=2,
    map_grid_longitude_increment=2,
    map_layer_mode="background",
)

# define map view
view = mv.geoview(
    map_area_definition="corners", area=[49, -2.5, 54, 4], coastlines=coast
)

# define title
title = mv.mtext(text_line_1="SYNOP Station Identifiers", text_font_size=0.5)

# define the output plot file
mv.setoutput(mv.pdf_output(output_name="bufr_synop_station_id"))

# generate plot
mv.plot(view, gpt, sym_txt, sym_loc, title)