Objectives

• Set-up a cartesian projection for the Cross section
• Configure the axis
• Load a Netcdf file, and understand the information needed by Magics.
• Apply a shading
• Position the legend box
• Add a text.

You will need to download

• The Xsection Netcdf Data
• and later ... the full solution

from Magics.macro import *
#setting the output 
output = output(    
                output_formats = ['png'],
                output_name = "xsect_step1",
                output_name_first_page_number = "off"
        )
# Setting the cartesian view
projection = mmap(
        subpage_map_projection='cartesian',
        subpage_x_axis_type='geoline',
        subpage_y_axis_type='logarithmic',
        subpage_x_min_latitude=50.,
        subpage_x_max_latitude=30.,
        subpage_x_min_longitude=-90.,
        subpage_x_max_longitude=-60.,
        subpage_y_min=1020.,
        subpage_y_max=200.,
        )
# Vertical axis
vertical = maxis(
        axis_orientation='vertical',
        )
# Horizontal axis
horizontal = maxis(
        axis_orientation='horizontal',
        )
        
plot(output, projection, horizontal, vertical)