# Metview Macro

rh = read("relative_humidity.grib")
z  = read("geopotential.grib")


field_names      = grib_get_string(rh[1], 'name') & ' / ' & grib_get_string(z[1], 'name')
field_valid_date = grib_get_string(rh[1], 'validityDate')
field_valid_time = grib_get_string(rh[1], 'validityTime')
field_level      = "<grib_info key='level' where='shortName=r'/>"

title_text = field_names & ' ' & field_valid_date & ' ' & field_valid_time & ' at level ' & field_level
print(title_text)

title = mtext(text_line_1: title_text)

rh_shade = mcont(
    legend                         : "on",
    contour                        : "off",
    contour_shade                  : "on",
    contour_shade_method           : "area_fill",
    contour_shade_max_level_colour : "red",
    contour_shade_min_level_colour : "white",
    contour_shade_colour_direction : "clockwise"
    )

plot(rh, rh_shade, z, title)
