# Metview Macro

# **************************** LICENSE START ***********************************
#
# Copyright 2012 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 the BUFR file and extract t2m
obs_area_bufr = read("../obs_area.bufr")

t2_geo = obsfilter(
        output  :  "geopoints",
        data    :  obs_area_bufr
        )


# read the GRIB file of 2T analysis
t2_grib = read("../t2_an.grib")


# Compute the difference
diff = t2_geo - t2_grib

# plotting visdef
diff_symb = msymb(
    symbol_type                          : "marker",
    symbol_table_mode                    : "advanced",
    symbol_advanced_table_selection_type : "list",
    symbol_advanced_table_level_list     : [-1000,-1,1,1000],
    symbol_advanced_table_colour_method  : "list",
    symbol_advanced_table_colour_list    : ["blue","grey","red"]
    )

# plot the data
plot(diff, diff_symb)


