Tabs Page |
---|
|
Code Block |
---|
language | py |
---|
title | Wind Coloured By Temperature Example |
---|
|
#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 grib file - contains model level data, one timestep
fs = read(source : "fc_ml.grib")
# define model level to read
level = 96 # around 500 hPa
# read temperature and scale it to C
t = read(data : fs, param : "t", level: level)
t = t - 273.16
# read wind components
u = read(data : fs, param : "u", level: level)
v = read(data : fs, param : "v", level: level)
# define wind plotting - will be coloured by temperature
wp = mwind(
wind_thinning_factor : 1,
legend : "on",
wind_advanced_method : "on",
wind_advanced_colour_parameter : "parameter",
wind_advanced_colour_max_level_colour : "red",
wind_advanced_colour_min_level_colour : "violet",
wind_advanced_colour_direction : "clockwise",
wind_arrow_unit_velocity : 35,
wind_arrow_thickness: 2
)
# define coastlines
coast = mcoast(
map_coastline_land_shade : "on",
map_coastline_land_shade_colour : "RGB(0.8,0.8,0.8)",
map_coastline_sea_shade : "on",
map_coastline_sea_shade_colour : "RGB(0.9,0.9,0.9)",
map_coastline_colour : "RGB(0.2,0.2,0.2)",
map_coastline_resolution : "medium"
)
# define geo view
view = geoview(
coastlines: coast,
map_area_definition : "corners",
area : [30,-30,75,40]
)
# define the vector structure for plotting - wind will be coloured by t
v = grib_vectors(u_component: u,
v_component: v,
colouring_field: t)
# define the output plot file
setoutput(pdf_output(output_name : 'wind_coloured_by_t'))
# generate plot
plot(view,v,wp)
|
|
Tabs Page |
---|
|
Code Block |
---|
language | py |
---|
title | Wind Coloured By Temperature Example |
---|
|
# **************************** 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 grib file - contains model level data, one timestep
fs = mv.read(source = "fc_ml.grib")
# define model level to read
level = 96 # around 500 hPa
# read temperature and scale it to C
t = mv.read(data = fs, param = "t", level = level)
t = t - 273.16
# read wind components
u = mv.read(data = fs, param = "u", level = level)
v = mv.read(data = fs, param = "v", level = level)
# define wind plotting - will be coloured by temperature
wp = mv.mwind(
wind_thinning_factor = 1,
legend = "on",
wind_advanced_method = "on",
wind_advanced_colour_parameter = "parameter",
wind_advanced_colour_max_level_colour = "red",
wind_advanced_colour_min_level_colour = "violet",
wind_advanced_colour_direction = "clockwise",
wind_arrow_unit_velocity = 35,
wind_arrow_thickness = 2
)
# define coastlines
coast = mv.mcoast(
map_coastline_land_shade = "on",
map_coastline_land_shade_colour = "RGB(0.8,0.8,0.8)",
map_coastline_sea_shade = "on",
map_coastline_sea_shade_colour = "RGB(0.9,0.9,0.9)",
map_coastline_colour = "RGB(0.2,0.2,0.2)",
map_coastline_resolution = "medium"
)
# define geo view
view = mv.geoview(
coastlines = coast,
map_area_definition = "corners",
area = [30,-30,75,40]
)
# define the vector structure for plotting - wind will be coloured by t
v = mv.grib_vectors(u_component = u,
v_component = v,
colouring_field = t)
# define the output plot file
mv.setoutput(mv.pdf_output(output_name = 'wind_coloured_by_t'))
# generate plot
mv.plot(view, v, wp)
|
|
|