<reponame>b8raoult/magics
# (C) Copyright 1996-2016 ECMWF.
# 
# This software is licensed under the terms of the Apache Licence Version 2.0
# which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 
# In applying this licence, ECMWF does not waive the privileges and immunities 
# granted to it by virtue of its status as an intergovernmental organisation nor
# does it submit to any jurisdiction.


#importing Magics module
from Magics.macro import *

#Example reference
ref = 'wind1'

#Setting of the output file name
output = output(output_formats= ['png'],
                output_name_first_page_number= 'off',
                output_name= ref)

#Setting the coordinates of the geographical area
australia = mmap(subpage_upper_right_longitude= 180.,
                 subpage_upper_right_latitude= -5.,
                 subpage_lower_left_longitude= 105.,
                 subpage_clipping='on',
                 subpage_frame_colour = "evergreen",
                             subpage_frame_line_style = "dot",
                                         subpage_frame_thickness = 5,

                 subpage_map_projection= 'cylindrical',
                 subpage_lower_left_latitude= -55.)


#Background Coastlines 
background = mcoast( map_coastline_sea_shade_colour= 'white',
                     map_coastline_land_shade_colour= 'cream',
                     map_grid= 'off',
                     map_coastline_land_shade= 'on',
                     map_coastline_sea_shade= 'on',
                     map_label= 'off',
                     map_coastline_colour= 'tan')

#Foreground Coastlines
foreground = mcoast( map_grid= 'on',
                 map_grid_colour = 'tan',
                 map_label= 'off',
                 map_coastline_colour= 'tan',
                 map_coastline_land_shade= 'off',
                 map_coastline_sea_shade= 'off')

#Import the z500 data 
msl =  mgrib( grib_input_file_name = './msl.grb',
                grib_id= 'z500')

#Define the simple contouring for msl
msl_contour = mcont(
		legend= 'off', 
		contour_level_selection_type= 'interval',
                contour_interval= 5.,
                contour_line_colour= 'black',
                contour_line_thickness= 1,
                contour_hilo= 'on',
                contour_hilo_quality= 'high',
                contour_hi_colour= 'black',
                contour_lo_colour= 'black',
                contour_hilo_height= 0.25,
                contour_label= 'off',
                contour_highlight_colour= 'black',
                contour_highlight_thickness= 2,
                )

#Import the  wind speed at 200hPa speed200 
speed200 =  mgrib( grib_input_file_name = './speed200.grb',
                grib_id= 't850')

#Define the shading for the wind speed
speed200_contour = mcont(
    legend= 'on',
    contour_level_selection_type= 'level_list', 
    contour_level_list= [30., 40., 50., 60., 70., 80., 90., 100.], 
    contour_shade= 'on', 
    contour_shade_max_level_colour= 'evergreen', 
    contour_shade_min_level_colour= 'yellow',
    contour_shade_method= 'area_fill', 
    contour_reference_level= 0., 
    contour_highlight= 'off', 
    contour_hilo= 'on', 
    contour_hilo_format= '(F3.0)', 
    contour_hilo_height= 0.2, 
    contour_hilo_type= 'number', 
    contour_hilo_suppress_radius= 30., 
    contour_hi_min_value= 15., 
    contour_label= 'off' 
         )


#Import the  wind  at 200hPa uv200 
uv200 =  mgrib( grib_input_file_name = './uv200.grb',
                grib_id= 'uv200')

uv200_wind = mwind(
                legend= 'on',
                wind_field_type= 'arrows',
                wind_arrow_thickness= 1,
                wind_thinning_factor= 3.,
                wind_arrow_colour= 'gold' )

title = mtext(text_lines = ["<font size='1'>Wind arrows with legend</font>",
	                    "",
	                    "<grib_info id='z500' key='name' format='%s'/>"],
	      text_justification = 'left',
	      text_font_size = 0.5,
	      text_colour = 'charcoal')

#add a legend
legend = mlegend(legend= 'on',
           legend_text_colour= 'black',
           legend_box_mode= 'positional',
           legend_box_x_position= 27.5,
           legend_box_y_position= 4.,
           legend_box_x_length= 2.,
           legend_box_y_length= 12.,
           legend_border= 'off',
           legend_border_colour= 'black',
           legend_box_blanking= 'on',
           legend_display_type= 'continuous',
           legend_title = "on",
	   legend_title_text= "Wind speed at 200 hPa",
	   legend_text_font_size = "0.5")

#To the plot
plot(output, australia, background,speed200, speed200_contour,uv200, uv200_wind,msl, msl_contour,foreground, title, legend)

#For documentation only
tofortran(ref, output, australia, background,speed200, speed200_contour,uv200, uv200_wind,msl, msl_contour,foreground, title, legend)
tomv4(ref,uv200_wind)
tohtml(ref,uv200_wind)
