Scroll pdf ignore | ||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
...
Case description
In this exercise we will use Metview to produce the plots shown above:
- the plot at the top shows a map with mean sea level pressure forecast fields overlayed with the track of Hurricane Sandy.
- the plot at the bottom contains a graph chart showing the evolution of the minimum of the mean sea level pressure forecast along the storm track.
You We will need to work both interactively with icons interactively and write Macro code, as well. First, you we will create the two plots independently then align them together in the same page.
Info |
---|
This exercise involves:
|
Preparations
...
Verify that the data are as expected.
We will prepare the plot interactively using icons. Then, at the end, we will put it all together into a macro. Remember to give your icons useful names!
Creating the map plot
Setting the
...
View
With a new Geographical View icon, set up a cylindrical projection with its area defined as
...
- the land coloured in grey
- the sea coloured as #dcf0ff (or your favourite shade of blue)
Plotting the Mean Sea Level Pressure field
Plot the GRIB file sandy_msl.grib into this your view using a new Contouring icon. Plot black isolines with an interval of 5 hPa between them. Animate through the fields to see how the forecast is evolving.
Info |
---|
The fields you visualised were taken from the model run at 2012-10-27 0UTC and containing 12 hourly forecast steps from 0 to 120 hours. |
Plotting the storm track
The storm track data is stored in the CSV file called sandy_track.txttxt. If you open this file you will see that it contains the date, time and geographical coordinates of the track points.
...
- set the plot type to geo points
- select the columns holding storing the latitude and longitude by their index (column indexing start at 1)
- carefully specify the table delimiter and header information by setting
- table_delimiter to whitespace (enter a space character here)
- table_combine_delimiters to on (it means we allow multiple whitespaces as separators)
- table header_row to 0 (it means do not have a header in the file)
Now drag your your Table Visualiser icon the plot to overlay the mean sea level pressure forecast with the track.
Customising the storm track
The storm track in its current form does not look great so you need to customise it with needs some customisation. Use a Graph Plotting icon for it by setting the
- the track line to black and thick
- the track points to be white filled circles (their marker index is 15) with red outline.
Plotting date/time labels onto the track
To finalise the track plot you need to we will add the date/time labels to the track points. This can could be done with a Symbol Plotting icon by specifying the list of labels you we want to plot into the map. Since it would require too much editing you will learn a better (programmatic) way of doing it However, to make it reusable for other datasets as well, we will do it programmatically by using Metview Macro.
Create a new Macro and edit it. First, read the CSV file in with the Table Reader:
Code Block |
---|
tbl = read_table( table_delimiter : " ", table_combine_delimiters : "on", table_header_row : 0, table_filename : "sandy_track.txt" ) |
Info |
---|
As you can see we specified the table delimiters exactly in the same way as we did for the Table Visualiser icon. |
In the code above, Here the object referenced by variable tbl
contains all the columns from the CSV file. Now read the date and time (from the first two columns) into separate vectors:
Code Block |
---|
val_date=values(tbl,1) val_time=values(tbl,2) |
Next, you need to we build the list of labels. Each label is made up from a day and an hour part separated by a slash. We convert the date into a string and then take the last two characters to get the day. Use this loop to construct the list of labels:
Code Block |
---|
labels=nil for i=1 to count(val_date) do labels = labels & [" " &dPart = substring(string(val_date[i]),7,8) & "/"tPart &= val_time[i] label = " " & dPart & "/" & tPart labels = labels & [label] end for |
FinallyNext, define a Symbol Plotting visual definition and return itusing the text mode.
Info |
---|
Symbol Plotting in text mode is used to plot string values to the positions of the dataset it is applied to. The rule is that the first string in the list defined by |
The code you we need to add is like this:
Code Block |
---|
sym = msymb(
symbol_type : "text",
symbol_text_font_colour : "black",
symbol_text_font_size: "0.3",
symbol_text_font_style: "bold",
symbol_text_list : labels
)
|
We finish the macro by returning our Visual Definition.
Code Block |
---|
return [sym] |
Info |
---|
By returning the visual definition |
...
our Macro behaves as if it were a real Symbol Plotting icon. |
Save So save the Macro and drag it into the plot to see the labels appearing along the track.
Creating the graph plot
Setting the View
With a new Cartesian View icon, set up a view to cater for the graph showing the mean sea level pressure values in hPa units by setting
- the x-axis type to date
- the x axis minimum to 2012-10-27 and its maximum to 2012-10-31
- the y-axis label to hPa
- the y-axis minimum value to 940 and its maximum is to 1000
Computing the minimum pressure along the track
Since this task is fairly complex you we will use a Macro for it. The idea goes like this:
- you we read the track points from the CSV file
- we define a lat-lon box around each point
- we read the forecast mean sea level data for the box for the corresponding time
- we compute the minimum of the pressure in the box
- from these minimum values you we can build the curve data to plot.
Create new Macro and edit it. First, read the CSV file in the very same way as before but . However, this time, on top of date and time, you we also need to read latitude and longitude into vectors:
Code Block |
---|
val_latlon=values(tbl,3) val_lonlat=values(tbl,4) |
Next, read in the GRIB file containing the mean sea level forecast:
...
The curve data requires two lists: one for the dates and one for the values. First you need to we initialise these lists:
Code Block |
---|
trVal = nil trDate = nil |
Now the main part of the macro follows: you need to we will loop through the track points and build the curve dataset. You can We will use a loop like this:
Code Block |
---|
for i=1 to count(val_date) do ... your code will go here ... end for |
Within the loop first construct an area of e.g. 10 degrees wide centred on the current track point.
...
Code Block |
---|
p=read( data: g, step: (i-1)*2412, area : wbox ) |
Here we used the fact the forecasts steps are stored in hours units in the GRIB file.
Next, compute the minimum of the field you read in the subarea using the minval()
macro function:
Code Block |
---|
pmin=minvalminvalue(p) |
Finally, build the lists list for the values (scaling Pa units stored in the GRIB to hPa units):
Code Block |
---|
trVal= trVal & [vpmin/100] |
and for the And also build the list of dates:
Code Block |
---|
trDatedt = trDate & [date(val_date[i]) + hour(val_time[i]) trDate = trDate & [dt] |
Having finished the body of the loop the last step in our Macro is to define an Input Visualiser and return it. The code you we need to add is like this:
Code Block |
---|
vis = input_visualiser ( input_x_type : "date", input_date_x_values : trDate, input_y_values : trVal ) return [vis] |
Info |
---|
By returning the |
...
visualiser our Macro behaves as if it were an Input Visualiser icon. |
Now visualise your Cartesian View icon and drag your Macro into it.
Customising the graph
Customise the graph with a Graph Plotting icon by setting the
- the line thicker
- the points to be blue filled circles (their marker index is 15) with a reasonable size.
Creating a title
Define a custom title to as shown in the example plot with a new Text Plotting icon.
Putting it all together
With a new Display Window icon design an A4 portrait layout with two views: your Geographical View icon should go top and your Cartesian View icon into the bottom. Now visualise your this icon and populate the view views with the data.
Extra Work
Try the following if you have time.
Ensuring the title has the correct date and time
There are various ways we can ensure that the title has the date and time according to the actual data. The default title in fact contains the date and time, but in this exercise we want more control over it.
Construct the second line of the title by extracting the date and time from the MSLP field and converting into an appropriate string - do this in the Macro code.
Hints:
- if you have a fieldset variable called
msl_grib
, the following line will extract the date at which the field is valid:msl_date = grib_get_long(msl_grib, 'validityDate')
- we can do something similar for the validity time
- these are extracted as integer numbers, but can be combined into a proper date variable in Macro:
full_date = date(msl_date) + hour(msl_time)
- use the
string()
function to construct a date string similar to the one used in the current title- see String Functions for details of how to use it
- insert this into the
mtext()
function instead of the current title - it is now more robust - if you use data from a different date or time, the title will still be correct
- note that this method will not work directly if you want to generate an animation from different time steps of data
Experiment with different backgrounds
Modify the Coastlines icon, for example:
Adding new curves to the x-y plot
On top of the minimum pressure try to add the maximum and average pressure to the graph plot. Use a different colour to each curve and add a custom legend as well.
Hints:
- first, just try to add your Graph Plotting definition to the Macro. In the end return both the Input Visualiser and the Graph Plotting as a list like this
Code Block |
---|
return [vis,graph] |
If you visualise the Macro your Graph Plotting settings will be directly applied to the resulting curve.
- next, compute the maximum of the pressure (with the
maxvalue()
function) in the loop and store its values in another list. Build an input visualiser out of it (e.g. call itvis_max
). Add a Graph Plotting for it (e.g. call it graph_max) using a different colour. In the end you need to return a longer list like this:
Code Block |
---|
return [vis,graph,vis_max,graph_max] |
- the average pressure curve (with the
average()
function) can be derived in a very similar manner
add a Legend with disjoint mode. Set legend_text_composition to user_text_only and carefully set the legend_user_lines to provide a textual description to each curve in the legend. Add your legend to the back of the list you return from the Macro.
Doing the whole task in Macro
Try to write a single Macro that is doing all the tasks in one go and directly produces the composite plot with the map and graph in the end
...