#Metview Macro

#MARS retrieval
mars_hirs = retrieve(
		type	:	"mfb",
		repres	:	"bu",
		obsgroup	: "hirs",
		time	:	00,
		resol	:	"",
		filter	:	"select lat, lon, obsvalue, vertco_reference_1"
		)
	
#Get observation values
val=values(mars_hirs,"obsvalue@body")

#Compute stats
min_v=minvalue(val)
max_v=maxvalue(val)
mean_v=mean(val)

#Print results
print("min: ",min_v," max: ",max_v," mean: ",mean_v)
