|
|
# Metview Macro
# **************************** LICENSE START ***********************************
#
# Copyright 2020 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 ************************************
# -------------------------------------------------------------------------
# Description: Demonstrates how to make a graph plot with two vertical axes
# -------------------------------------------------------------------------
# We will compute a vertical profile from a point in a data set
# that is on model levels. This will be plotted with model levels
# as the left-hand vertical axis, and computed pressure levels on the
# right-hand vertical axis.
lat_lon_point = [56.9, -2.6]
# read the input model level grib file and compute a profile curve from it
t_data = read(source: "fc_ml.grib", parameter:'t')
model_levels = grib_get_long(t_data, 'level')
iv1 = input_visualiser(
input_x_values : nearest_gridpoint(t_data, lat_lon_point),
input_y_values : model_levels
)
dots = msymb(
symbol_type : "marker",
symbol_outline : "on",
symbol_marker_index : 15
)
# read the LNSP file and compute pressure levels for the input model levels
lnsp_data = read("lnsp.grib")
pressure_fields = pressure(lnsp_data, t_data)
pressure_levels = nearest_gridpoint(pressure_fields, lat_lon_point) / 100
# Construct two Cartesian views; the first will be plotted with data, the
# second will be plotted without data in order to plot just its vertical axis
# set up the main Cartesian view - it has a left-hand vertical axis
primary_vertical_axis = maxis(
axis_orientation : "vertical",
axis_title_text : "model levels",
axis_title_height : 0.5,
axis_tick_label_hight : 0.5
)
c1 = cartesianview(
y_min : maxvalue(vector(model_levels)), # larger level usually at bottom
y_max : minvalue(vector(model_levels)),
x_min : 180, # Kelvin
x_max : 320, # Kelvin
vertical_axis : primary_vertical_axis
)
# set up the secondary Cartesian view with a right-hand vertical axis
secondary_vertical_axis = maxis(
axis_orientation : "vertical",
axis_position : "right",
axis_title_text : "pres levels",
axis_title_height : 0.5,
axis_tick_label_hight : 0.5
)
# switch off the plotting of its horizontal axis
secondary_horizontal_axis = maxis(
axis_line : 'off',
axis_tick : 'off',
axis_tick_label : 'off'
)
c2 = cartesianview(
y_min : maxvalue(vector(pressure_levels)),
y_max : minvalue(vector(pressure_levels)),
vertical_axis : secondary_vertical_axis,
horizontal_axis : secondary_horizontal_axis)
# set up the pages that will use these views
p1 = plot_page(view : c1)
p2 = plot_page(view : c2)
dw = plot_superpage(pages : [p1,p2])
# define the output plot file
setoutput(pdf_output(output_name : 'double_axis_1'))
# plot the data into the first view, then plot the second view empty
plot(dw[1],iv1,dots,dw[2])
|
|
|
# Metview Example
# **************************** LICENSE START ***********************************
#
# Copyright 2020 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
# -------------------------------------------------------------------------
# Description: Demonstrates how to make a graph plot with two vertical axes
# -------------------------------------------------------------------------
# We will compute a vertical profile from a point in a data set
# that is on model levels. This will be plotted with model levels
# as the left-hand vertical axis, and computed pressure levels on the
# right-hand vertical axis.
lat_lon_point = [56.9, -2.6]
# read the input model level grib file and compute a profile curve from it
t_data = mv.read(source="fc_ml.grib", parameter='t')
model_levels = mv.grib_get_long(t_data, 'level')
iv1 = mv.input_visualiser(
input_x_values = mv.nearest_gridpoint(t_data, lat_lon_point),
input_y_values = model_levels
)
dots = mv.msymb(
symbol_type = "marker",
symbol_outline = "on",
symbol_marker_index = 15
)
# read the LNSP file and compute pressure levels for the input model levels
lnsp_data = mv.read("lnsp.grib")
pressure_fields = mv.pressure(lnsp_data, t_data)
pressure_levels = mv.nearest_gridpoint(pressure_fields, lat_lon_point)
pressure_levels = [x / 100 for x in pressure_levels]
# Construct two Cartesian views; the first will be plotted with data, the
# second will be plotted without data in order to plot just its vertical axis
# set up the main Cartesian view - it has a left-hand vertical axis
primary_vertical_axis = mv.maxis(
axis_orientation = "vertical",
axis_title_text = "model levels",
axis_title_height = 0.5,
axis_tick_label_hight = 0.5
)
c1 = mv.cartesianview(
y_min = max(model_levels), # larger level usually at bottom
y_max = min(model_levels),
x_min = 180, # Kelvin
x_max = 320, # Kelvin
vertical_axis = primary_vertical_axis
)
# set up the secondary Cartesian view with a right-hand vertical axis
secondary_vertical_axis = mv.maxis(
axis_orientation = "vertical",
axis_position = "right",
axis_title_text = "pres levels",
axis_title_height = 0.5,
axis_tick_label_hight = 0.5
)
# switch off the plotting of its horizontal axis
secondary_horizontal_axis = mv.maxis(
axis_line = 'off',
axis_tick = 'off',
axis_tick_label = 'off'
)
c2 = mv.cartesianview(
y_min = max(pressure_levels),
y_max = min(pressure_levels),
vertical_axis = secondary_vertical_axis,
horizontal_axis = secondary_horizontal_axis)
# set up the pages that will use these views
p1 = mv.plot_page(view = c1)
p2 = mv.plot_page(view = c2)
dw = mv.plot_superpage(pages = [p1,p2])
# define the output plot file
mv.setoutput(mv.pdf_output(output_name = 'double_axis_1'))
# plot the data into the first view, then plot the second view empty
mv.plot(dw[0],iv1,dots,dw[1])
|
|
|
