Objectives


  • Set-Up a cartesian projection
  • Learn how to pass  arrays as input of curves
  • Learn how to use CSV files as input of curves
  • Add a text.
  • Create a layout

You will need to download


 

 

Profile: create the Cartesian Projection

For our vertical profile, we will need :

  • a horizontal regular axis for a temperature range going from -60oC to 20oC.
  • a vertical logarithmic axis going from 1020 hPa to 10 hPa

Have a look at the Subpage Documentation and at the Axis Plotting Documentation, and try to get the 0oC grid line highlighted.

Parameters to check

 

 

Useful projection parameters

subpage_map_projection

subpage_x_axis_type
subpage_y_axis_type
subpage_x_min
subpage_y_min
subpage_x_max
subpage_y_max
Useful axis parameters
axis_orientation
axis_type
axis_grid

axis_grid_reference

axis_title
Python - set-up of cartesian projection
# importing Magics module
from Magics.macro import *
# define the output
output = output(output_formats=['png'],
        output_name_first_page_number='off',
        output_name="profile_step1")
projection = mmap(
    subpage_map_projection='cartesian',
    subpage_x_axis_type='regular',
    subpage_y_axis_type='logarithmic',   
    subpage_x_min=-70.,
    subpage_x_max=20.,
    subpage_y_min=1020.,
    subpage_y_max=10.,)
#define the vertical axis
vertical = maxis(
    axis_orientation='vertical',
    axis_grid='on',
    axis_type='logarithmic',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid_colour='charcoal',
    axis_grid_line_style='dash',
    axis_title='on',
    axis_title_text='Pressure',
    axis_title_font='arial',
    axis_title_font_style='bold',
    axis_title_height=1.,
    )
# define the horizontal axis
horizontal = maxis(
    axis_orientation='horizontal',
    axis_type='regular',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid='on',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_reference_level=0.,
    axis_grid_reference_line_style='solid',
    axis_grid_reference_thickness=1,
    axis_grid_line_style='dash',)
    
plot(output, projection, vertical, horizontal)

 

 

Profile: Using Arrays as Input of a Graph

Here is the list of levels, and the list of Temperature in Kelvin at each levels.

levels = [1.,2.,3,5,7,10,20,30,50,70,100,150,200,250,300,400,500,600,700,800,850,900,925,950,1000]
kelvin = ([263.118652344,254.822738647,242.517868042,223.301025391,219.254943848,216.710174561,216.507095337, 215.398986816,211.643295288,207.18812561,207.172134399,217.097396851,223.809326172,235.13168335,243.377059937,260.635147095,272.02935791,272.145584106,273.448501587,279.562927246,281.745040894,285.503082275,287.543685913,289.284072876,292.170974731]

We can pass these arrays to Magics using the minput object . Check the Input Data Documentation for more information about the parameters. Note that we want to use Celsius in our graph...

The mgraph object can then be used to  preform the visualisation. All the parameters can be found in the Graph Plotting Page.

Parameters to check

 


 

Useful input parameters

input_x_values

input_y_values
Useful graph parameters
graph_line_colour
graph_line_thickness
Python - Input data for graph
# importing Magics module
from Magics.macro import *
# define the output
output = output(output_formats=['png'],
        output_name_first_page_number='off',
        output_name="profile_step2")
projection = mmap(
    subpage_map_projection='cartesian',
    subpage_x_axis_type='regular',
    subpage_y_axis_type='logarithmic',   
    subpage_x_min=-70.,
    subpage_x_max=20.,
    subpage_y_min=1020.,
    subpage_y_max=10.,)
#define the vertical axis
vertical = maxis(
    axis_orientation='vertical',
    axis_grid='on',
    axis_type='logarithmic',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid_colour='charcoal',
    axis_grid_line_style='dash',
    axis_title='on',
    axis_title_text='Pressure',
    axis_title_font='arial',
    axis_title_font_style='bold',
    axis_title_height=1.,
    )
# define the horizontal axis
horizontal = maxis(
    axis_orientation='horizontal',
    axis_type='regular',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid='on',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_reference_level=0.,
    axis_grid_reference_line_style='solid',
    axis_grid_reference_thickness=1,
    axis_grid_line_style='dash',)
levels = [1.,2.,3,5,7,10,20,30,50,70,100,150,200,250,300,400,500,600,700,800,850,900,925,950,1000]
kelvin = numpy.array([263.118652344,254.822738647,242.517868042,223.301025391,219.254943848,216.710174561,216.507095337,
    215.398986816,211.643295288,207.18812561,207.172134399,217.097396851,
    223.809326172,235.13168335,243.377059937,260.635147095,272.02935791,
    272.145584106,273.448501587,279.562927246,281.745040894,285.503082275,287.543685913,289.284072876,292.170974731]) 
celsius = kelvin -273.16
# Define the input of the graph
data = minput(input_x_values = celsius,
        input_y_values = levels, )
#define the graph action.
graph = mgraph( legend='on' ,
            legend_user_text= 'Tempe', 
            graph_line_colour="navy", 
            graph_line_thickness= 3, )
    
plot(output, projection, vertical, horizontal, data, graph)

 

 

Profile: Using Arrays as Input for wind Plotting

Here is the list of pressure levels, and the list of u and v components for each pressure level.

levels = [1.,2.,3,5,7,10,20,30,50,70,100,150,200,250,300,400,500,600,700,800,850,900,925,950,1000]
u_component = [77.815612793,63.232711792,44.9090881348,28.2789916992,17.401763916,
    11.8307037354,-0.725036621094,-4.75889587402,-7.83796691895,-3.76489257812,
    -6.18388366699,-13.1139526367,-13.05027771,-15.1851959229,-23.5844726562,
    -27.866394043,-33.0458984375,-45.8946380615,-41.540222168,-39.6210327148,
    -37.6231689453,-32.947265625,-29.4743499756,-22.0828552246,-10.0008392334]
v_component = [7.38911437988,5.23355102539,0.636322021484,-3.75567626953,
    0.167388916016,-3.12893676758,-3.60003662109,0.0688171386719,
    0.305236816406,7.06889343262,8.10516357422,11.1141815186,
    -2.29931640625,3.08065795898,7.5965423584,11.0112915039,
    13.6437530518,11.2753143311,11.1552124023,13.1494903564,
    6.01538085938,1.37588500977,-0.0724334716797,-0.953598022461,-0.508377075195]

We can pass these arrays to Magics using the minput object . Check the Input Data Documentation for more information about the parameters.

We want to plot these winds as red flags on the vertical of -20.

The mgraph object can then be used to preform the visualisation. All the parameters can be found in the Graph Plotting Page.

Do not forget to add a legend with some meaningful text!

Parameters to check

 


 

Useful input parameters

input_x_values

input_y_values
input_x_component_values
input_y_component_values
Useful graph parameters
graph_type
graph_flag_colour
graph_flag_length
legend
legend_user_text
Python - Input data and Flags plotting
# importing Magics module
from Magics.macro import *
# define the output
output = output(output_formats=['png'],
        output_name_first_page_number='off',
        output_name="profile_step2")
projection = mmap(
    subpage_map_projection='cartesian',
    subpage_x_axis_type='regular',
    subpage_y_axis_type='logarithmic',   
    subpage_x_min=-70.,
    subpage_x_max=20.,
    subpage_y_min=1020.,
    subpage_y_max=10.,)
#define the vertical axis
vertical = maxis(
    axis_orientation='vertical',
    axis_grid='on',
    axis_type='logarithmic',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid_colour='charcoal',
    axis_grid_line_style='dash',
    axis_title='on',
    axis_title_text='Pressure',
    axis_title_font='arial',
    axis_title_font_style='bold',
    axis_title_height=1.,
    )
# define the horizontal axis
horizontal = maxis(
    axis_orientation='horizontal',
    axis_type='regular',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid='on',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_reference_level=0.,
    axis_grid_reference_line_style='solid',
    axis_grid_reference_thickness=1,
    axis_grid_line_style='dash',)
levels = [1.,2.,3,5,7,10,20,30,50,70,100,150,200,250,300,400,500,600,700,800,850,900,925,950,1000]
kelvin = numpy.array([263.118652344,254.822738647,242.517868042,223.301025391,219.254943848,216.710174561,216.507095337,
    215.398986816,211.643295288,207.18812561,207.172134399,217.097396851,
    223.809326172,235.13168335,243.377059937,260.635147095,272.02935791,
    272.145584106,273.448501587,279.562927246,281.745040894,285.503082275,287.543685913,289.284072876,292.170974731]) 
celsius = kelvin -273.16
# Define the input of the graph
data = minput(input_x_values = celsius,
        input_y_values = levels, )
#define the graph action.
graph = mgraph( legend='on' ,
            legend_user_text= 'Tempe', 
            graph_line_colour="navy", 
            graph_line_thickness= 3, )
#define the input for the wind plotting
wind_position = [-20.] * len(levels)
u_component = [77.815612793,63.232711792,44.9090881348,28.2789916992,17.401763916,
    11.8307037354,-0.725036621094,-4.75889587402,-7.83796691895,-3.76489257812,
    -6.18388366699,-13.1139526367,-13.05027771,-15.1851959229,-23.5844726562,
    -27.866394043,-33.0458984375,-45.8946380615,-41.540222168,-39.6210327148,
    -37.6231689453,-32.947265625,-29.4743499756,-22.0828552246,-10.0008392334]
v_component = [7.38911437988,5.23355102539,0.636322021484,-3.75567626953,
    0.167388916016,-3.12893676758,-3.60003662109,0.0688171386719,
    0.305236816406,7.06889343262,8.10516357422,11.1141815186,
    -2.29931640625,3.08065795898,7.5965423584,11.0112915039,
    13.6437530518,11.2753143311,11.1552124023,13.1494903564,
    6.01538085938,1.37588500977,-0.0724334716797,-0.953598022461,-0.508377075195]
wind = minput(
        input_y_component_values = v_component,
        input_x_values = wind_position,
        input_x_component_values = u_component,
        input_y_values = levels)
#define wind plotting
flags = mgraph(graph_type = "flag",
        graph_flag_colour = "red",
        legend = "on",
        graph_flag_length = 1.00,
        legend_user_text = " Wind ")

 plot(output, projection, vertical, horizontal, data, graph, wind, flags)

 

 

Time Series: Setting the Cartesian Projection for the Mean Sea Level Pressure time series

In this exercise, we will learn how to setup a date Coordinate, and browse the parameters available to style the style a date axis.

We want

  • a horizontal date axis going from the 2012-10-27at  00:00:00 to 2012-11-02 12:00:00
  • a vertical regular axis  going from 950 hPa to 1050 hPa

Have a look at the Subpage Documentation and at the Axis Plotting Documentation,

Parameters to check

 


 

Useful projection parameters

subpage_map_projection

subpage_x_axis_type
subpage_y_axis_type
subpage_x_date_min
subpage_y_min
subpage_x_date_max
subpage_y_max
Useful axis parameters
axis_orientation
axis_type
Python - Date Coordinates system
 # importing Magics module
from Magics.macro import *
# define the output
output = output(output_formats=['png'],
        output_name_first_page_number='off',
        output_name="time_serie_step1")
projection = mmap(
    subpage_map_projection='cartesian',
    subpage_x_axis_type='date',
    subpage_y_axis_type='regular',   
    subpage_x_date_min='2012-10-27 00:00:00',
    subpage_x_date_max='2012-11-02 12:00:00',
    subpage_y_max=1050.,
    subpage_y_min=950.,
    )
# Vertical axis
vertical = maxis(
    axis_orientation='vertical',
    axis_grid='on',
    axis_type='regular',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_reference_level=0.,
    axis_grid_reference_line_style='solid',
    axis_grid_reference_thickness=1,
    axis_grid_line_style='dash',
    axis_title='on',
    axis_title_text='Pressure',
    axis_title_font='arial',
    axis_title_height=0.5,
    )
# Horizontal axis
horizontal = maxis(
    axis_orientation='horizontal',
    axis_type='date',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid='on',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_line_style='dash',
    )
    
plot(output, projection, horizontal, vertical)

 

 

 

Time Series: Using a CSV file as input of a curve

The values for our time series are in CSV (Comma-Separeted-Value) ascii file called msl.csv.

first few lines of msl.csv
1,2012-10-27 00:00:00,1020.64125,0
2,2012-10-27 12:00:00,1017.261875,0.1220703125
3,2012-10-28 00:00:00,1014.27125,-9.58071727163e-21
4,2012-10-28 12:00:00,1011.97375,-2.80979885101e-20

The valid date is defined in the second column, the value of msl in the third.

This information need to be given to Magics . This can be done with the mtable data action,. See the full documentation.

The mgraph action will then be used to define the curve attributes. All the parameters can be found in the Graph Plotting Page.

You can always add a title using the mtext object.

Parameters to check

 


 

Useful table parameters

table_filename

table_variable_identifier_type
table_x_type
table_x_variable
table_y_variable
table_header_row
Useful graph parameters
graph_line_colour
legend
legend_user_text
Python - CSV file
# importing Magics module
from Magics.macro import *
# define the output
output = output(output_formats=['png'],
        output_name_first_page_number='off',
        output_name="time_serie_step2")
projection = mmap(
    subpage_map_projection='cartesian',
    subpage_x_axis_type='date',
    subpage_y_axis_type='regular',   
    subpage_x_date_min='2012-10-27 00:00:00',
    subpage_x_date_max='2012-11-02 12:00:00',
    subpage_y_max=1050.,
    subpage_y_min=950.,
    )
# Vertical axis
vertical = maxis(
    axis_orientation='vertical',
    axis_grid='on',
    axis_type='regular',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_reference_level=0.,
    axis_grid_reference_line_style='solid',
    axis_grid_reference_thickness=1,
    axis_grid_line_style='dash',
    axis_title='on',
    axis_title_text='Pressure',
    axis_title_font='arial',
    axis_title_height=0.5,
    )
# Horizontal axis
horizontal = maxis(
    axis_orientation='horizontal',
    axis_type='date',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid='on',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_line_style='dash',)
#load the csv file, and define how to interpret it
data = mtable(table_filename = "msl.csv",
    table_x_type='date',
    table_variable_identifier_type='index',
    table_x_variable = "2",
    table_y_variable = "3",
    table_header_row = 0
 )
#define the visualisation
graph = mgraph( legend='on' ,
    legend_user_text= '<font colour="navy" size="0.3"> Mean Sea Level</font>', 
    graph_line_colour="navy", 
    graph_line_thickness= 3, )
#add a title
title = mtext(
            text_lines= ['Mean Sea Level plotted as curve...', ],
            text_html= 'true',
            text_justification= 'left',
            text_font_size= 0.8,
            text_colour= 'charcoal',
        )
plot(output, projection, horizontal, vertical, data, graph, title)

 

 

Time Series : Using Automatic Axis and bar plotting

Magics can automatically set-up the coordinates system according to the plotted data.

In this exercise, we will use this functionality to draw the time series of the precipitation.

  • The horizontal date coordinate system is set automatically.
  • The vertical coordinate system is regular going from 0 to 5

The data will be pass to Magics using arrays:

dates = ["2012-10-27 00:00:00","2012-10-27 12:00:00","2012-10-28 00:00:00",
        "2012-10-28 12:00:00","2012-10-29 00:00:00","2012-10-29 12:00:00",
        "2012-10-30 00:00:00","2012-10-30 12:00:00",
        "2012-10-31 00:00:00","2012-10-31 12:00:00",
        "2012-11-01 00:00:00","2012-11-01 12:00:00",
        "2012-11-02 00:00:00","2012-11-02 12:00:00"]
precip = [0.,0.1220703125,0.,0., 0., 3.16429138184,
            2.18200683594,1.75476074219,0,0,0,0,0,0.]

We want to display the information using bar plotting. 

We will perhaps need to check quickly

Parameters to check

 


 

Useful projection parameters

subpage_x_automatic

Useful input parameters

input_x_type

input_date_x_values
input_y_values
input_y2_values
Useful graph parameters
graph_type
graph_bar_line_colour
graph_shade_colour
legend
legend_user_text
Python - Automatic projection and bar plotting
# importing Magics module
from Magics.macro import *
# define the output
output = output(output_formats=['png'],
        output_name_first_page_number='off',
        output_name="time_serie_step3")
# Setting the cartesian view
projection = mmap(
    subpage_map_projection='cartesian',
    subpage_x_axis_type='date',
    subpage_y_axis_type='regular',   
    subpage_x_automatic='on',
    subpage_y_max=5.,
    subpage_y_min=0.,
    )
    
# Vertical axis
vertical = maxis(
    axis_orientation='vertical',
    axis_grid='on',
    axis_type='regular',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_reference_level=0.,
    axis_grid_reference_line_style='solid',
    axis_grid_reference_thickness=1,
    axis_grid_line_style='dash',
    axis_title='on',
    axis_title_text='Amount of rain',
    axis_title_font='arial',
    axis_title_height=0.5,
    )
# Horizontal axis
horizontal = maxis(
    axis_orientation='horizontal',
    axis_type='date',
    axis_tick_label_height=0.4,
    axis_tick_label_colour='charcoal',
    axis_grid='on',
    axis_grid_colour='charcoal',
    axis_grid_thickness=1,
    axis_grid_line_style='dash',
    )
dates = ["2012-10-27 00:00:00","2012-10-27 12:00:00","2012-10-28 00:00:00",
    "2012-10-28 12:00:00","2012-10-29 00:00:00","2012-10-29 12:00:00",
    "2012-10-30 00:00:00","2012-10-30 12:00:00",
    "2012-10-31 00:00:00","2012-10-31 12:00:00",
    "2012-11-01 00:00:00","2012-11-01 12:00:00",
    "2012-11-02 00:00:00","2012-11-02 12:00:00"]
    
precip = [0.,0.1220703125,0.,0., 0., 3.16429138184,
        2.18200683594,1.75476074219,0,0,0,0,0,0.]
#defining the data 
data = minput(
    input_x_type='date',
    input_date_x_values= dates,
    input_y2_values= precip,
    input_y_values= [0.] * len(precip),
     )
#defining the visualisation
bar = mgraph( legend='on' ,
        graph_type ='bar',
        graph_bar_line_colour = "navy",
        graph_shade_colour = "evergreen",
        legend_user_text= '<font colour ="evergreen" size="0.4"> Precipitation </font>', 
        graph_line_colour="navy", 
        graph_bar_width = 3600.00 *6,
        graph_line_thickness= 3, )
#adding a title  
title = mtext(
        text_lines=  ['Precipitation as Bar Plotting..', ],
        text_html= 'true',
        text_justification= 'left',
        text_font_size= 0.8,
        text_colour= 'charcoal',
    )
plot(output, projection, vertical, horizontal, data, bar, title)

 

 

Finally, A bit of layout ...

On our image, we want to create 3 pages.. We will use the pseudo action page to do that..

layout

The position of the page is set with the 4 following parameters

  • page_x_position / page_y_position : to position in cm the bottom left corner of the page in its parent
  • page_x_length / page_y_length : the dimension in cm.

The drawing area ( where the plotting is rendered) is called subpage and can position into the page using the 4 following parameters

  • subpage_x_position / subpage_y_position : to position in cm the bottom left corner of the drawing area (subpage)  in its parent page.
  • subpage_x_length / subpage_y_length : the dimension in cm.

 

Parameters to check

 


 

Useful page parameters

layout

page_x_position
page_y_position
page_x_length
page_y_length
subpage_x_position
subpage_y_position
subpage_x_length
subpage_y_length
page_id_line
page_frame
Python - complex layout
# importing Magics module
from Magics.macro import *
from profile import profile
from timeserie import msl_timeserie
from timeserie import precip_timeserie
#Setting the output
output = output(output_name = 'layout', 
                output_formats = ["png"],
                output_name_first_page_number = "off")
profile_left = page(
        layout='positional',
        page_frame='off', 
        page_x_length=15., 
        page_id_line='off',
        subpage_y_length=17.,
        )
        
graph_bottom = page(
        layout='positional',
        page_frame='off', 
        page_x_length=15., 
        page_y_length= 10., 
        page_id_line='off',
        page_x_position=15., 
        page_y_position=0.5,
        subpage_x_position=1.5,
        subpage_y_position=0.5,
        subpage_y_length=7.,
        )
        
graph_top = page(
        layout='positional',
        page_frame='off', 
        page_x_length=15., 
        page_y_length=10., 
        page_id_line='off',
        page_x_position=15.,
        page_y_position=10.5, 
        subpage_x_position=1.5,
        subpage_y_position=0.5,    
        subpage_y_length=7.,
        )
plot(output, 
    profile_left, profile(title_size=0.6),
    graph_bottom, msl_timeserie(title_size=0.6),
    graph_top, precip_timeserie(title_size=0.6),
    )

Go to the next step...

Go to to the Main Magics Tutorial page...

1 Comment

  1. Delete or uncomment line 9 in curves.py. trajectories.py appears in the next tutorial.

    from trajectories import trajectories