#Metview Macro

#  **************************** LICENSE START ***********************************
# 
#  Copyright 2019 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 t field (the only field in the file)
t = read(source : "t850.grb", param : "t")

# compute gradient 
grad = gradient(t)

# scale gradient
grad = grad*1E5

# define the vector structure for plotting 
v = grib_vectors(u_component: grad[1],
                 v_component: grad[2])


# define wind plotting for gradient
wp = mwind(
    wind_thinning_factor                  : 1,
    legend                                : "on",
    wind_advanced_method                  : "on",
    wind_advanced_colour_min_value        : 1,
    wind_advanced_colour_max_level_colour : "red",
    wind_advanced_colour_min_level_colour : "violet",
    wind_advanced_colour_direction        : "clockwise",
    wind_arrow_unit_velocity              : 8,
    wind_arrow_thickness                  : 1,
    wind_arrow_legend_text                : "1E-2 K/km"
    )

# define contouring for temeprature
t_cont = mcont(
    contour_line_colour     : "RGB(0.4,0.4,0.4)",
    contour_line_style      : "dash",
    contour_line_thickness  : 2,
    contour_highlight       : "off",
    contour_level_selection_type : "interval",
    contour_interval             : 2,
    contour_label_height         : 0.25,
    grib_scaling_of_derived_fields: "on"
    )
    
# define caastline    
coast = mcoast(
    map_coastline_land_shade        : "on",
    map_coastline_land_shade_colour : "grey",
    map_coastline_sea_shade         : "on",
    map_coastline_sea_shade_colour  : "RGB(0.8944,0.9086,0.933)",
    map_grid_colour                 : "charcoal",
    map_grid_longitude_increment    : 10
    )

# define map view
view = geoview(
    map_projection         : "polar_stereographic",
    map_area_definition    : "corners",
    area                   : [37.71,-53.33,58.46,15.21],
    map_vertical_longitude : -30,
    coastlines             : coast
    )


# define title    
title = mtext(text_lines: ["Temperature 850 hPa [C]","Temperature Gradient 850 hPa [1E-2 K/km]"])    
 
# define the output plot file
setoutput(pdf_output(output_name : 'gradient_vector'))   
    
# generate plot    
plot(view, t, t_cont, v, wp, title)

#  **************************** LICENSE START ***********************************
# 
#  Copyright 2019 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 t field (the only field in the file)
t = mv.read(source = "t850.grb", param = "t")

# compute gradient 
grad = mv.gradient(t)

# scale gradient
grad = grad * 1E5

# define the vector structure for plotting 
v = mv.grib_vectors(u_component = grad[0],
                    v_component = grad[1])


# define wind plotting for gradient
wp = mv.mwind(
    wind_thinning_factor                  = 1,
    legend                                = "on",
    wind_advanced_method                  = "on",
    wind_advanced_colour_min_value        = 1,
    wind_advanced_colour_max_level_colour = "red",
    wind_advanced_colour_min_level_colour = "violet",
    wind_advanced_colour_direction        = "clockwise",
    wind_arrow_unit_velocity              = 8,
    wind_arrow_thickness                  = 1,
    wind_arrow_legend_text                = "1E-2 K/km"
    )

# define contouring for temeprature
t_cont = mv.mcont(
    contour_line_colour     = "RGB(0.4,0.4,0.4)",
    contour_line_style      = "dash",
    contour_line_thickness  = 2,
    contour_highlight       = "off",
    contour_level_selection_type = "interval",
    contour_interval             = 2,
    contour_label_height         = 0.25,
    grib_scaling_of_derived_fields = "on"
    )
    
# define caastline    
coast = mv.mcoast(
    map_coastline_land_shade        = "on",
    map_coastline_land_shade_colour = "grey",
    map_coastline_sea_shade         = "on",
    map_coastline_sea_shade_colour  = "RGB(0.8944,0.9086,0.933)",
    map_grid_colour                 = "charcoal",
    map_grid_longitude_increment    = 10
    )

# define map view
view = mv.geoview(
    map_projection         = "polar_stereographic",
    map_area_definition    = "corners",
    area                   = [37.71,-53.33,58.46,15.21],
    map_vertical_longitude = -30,
    coastlines             = coast
    )


# define title    
title = mv.mtext(text_lines = ["Temperature 850 hPa [C]","Temperature Gradient 850 hPa [1E-2 K/km]"])    
 
# define the output plot file
mv.setoutput(mv.pdf_output(output_name = 'gradient_vector'))   
    
# generate plot    
mv.plot(view, t, t_cont, v, wp, title)