#Metview Macro

#Read odb
amsua_2e_odb = read("AMSUA.odb")

#Perform the query - the output is an ODB
filter_dep = odb_filter(
		odb_query	:	"select " &
					"  lat@hdr as lat, " &
					"  lon@hdr as lon, " &
					"  fg_depar@body as fg_dep, " &
					"  an_depar@body as an_dep " &
					"where " &
					"  vertco_reference_1@body = 5 ",
		odb_data	:	amsua_2e_odb
		)


#Acces values from the result	
lat = values(filter_dep,"lat")
lon = values(filter_dep,"lon")
fg_dep = values(filter_dep,"fg_dep")
an_dep = values(filter_dep,"an_dep") 

#Compute the analysis increment as the difference of the analysis departure
#and first guess depature
num=count(lat)
incr = vector(num)
incr = an_dep - fg_dep

#Build a geopoints object for visualisation
geo=create_geo(count(lat),"xyv")
geo=set_latitudes(geo,lat)
geo=set_longitudes(geo,lon)
geo=set_values(geo,incr)

#Define symbol plotting icons for different value ranges
sym_small = msymb(
		symbol_type	:	"marker",
		symbol_table_mode	:	"advanced",
		symbol_advanced_table_selection_type	:	"interval",
		symbol_advanced_table_min_value	:	-0.1,
		symbol_advanced_table_max_value	:	0.1,
		symbol_advanced_table_reference_level	:	-0.1,
		symbol_advanced_table_interval	:	0.2,
		symbol_advanced_table_max_level_colour	:	"green",
		symbol_advanced_table_min_level_colour	:	"green",
		symbol_advanced_table_marker_list	:	3,
		symbol_advanced_table_height_list	:	0.05,
		legend: "on"
		)
		
sym_pos = msymb(
		symbol_type	:	"marker",
		symbol_table_mode	:	"advanced",
		symbol_advanced_table_selection_type	:	"interval",
		symbol_advanced_table_min_value	:	0.1,
		symbol_advanced_table_max_value	:	0.5,
		symbol_advanced_table_interval	:	0.1,
		symbol_advanced_table_max_level_colour	:	"red",
		symbol_advanced_table_min_level_colour	:	"yellow",
		symbol_advanced_table_colour_direction	:	"clockwise",
		legend: "on"
		)

sym_neg = msymb(
		symbol_type	:	"marker",
		symbol_table_mode	:	"advanced",
		symbol_advanced_table_selection_type	:	"interval",
		symbol_advanced_table_min_value	:	-0.5,
		symbol_advanced_table_max_value	:	-0.1,
		symbol_advanced_table_interval	:	0.1,
		symbol_advanced_table_max_level_colour	:	"cyan",
		symbol_advanced_table_min_level_colour	:	"blue",
		symbol_advanced_table_colour_direction	:	"clockwise",
		legend: "on")
		
#Define title
title = mtext(
		text_line_count	:	2,
		text_line_1	:	"Date: 2010-12-22 Time: 00",
		text_line_2	:	"Sat: METOP-A  Sensor: AMSU-A  Channel: 5 Param: an-fg [K]"
		)

#Visualise
plot(geo,sym_neg,sym_small,sym_pos,title)