Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Info

The macros described in this tutorial will write PostScript, PDF and PNG output files to the 'Figures' directory in the 'openifs_2018' folder. 

An animated gif is also produced is the same folder.

To save any other images during these exercises for discussion later, you can either use:

"Export" button in Metview's display window under the 'File' menu to save to PNG image format. This will also allow animations to be saved into postscript.

or use the following ksnapshot command to take a 'snapshot' of the screen :

Code Block
languagebash
titleCommand for screen snapshot
ksnapshot

and save it to a file.

If you want to create animations from other images, save the figures as postscript and then use the convert command:

...

Panel
borderColorlightgrey
borderStylesolid
titleMetview icons in Analysis folder

an_1x1.mv : this plots horizontal maps of parameters from the ECMWF analyses overlaid on one plot.

an_2x2.mv : this plots horizontal maps of parameters from the ECMWF analyses four plots to a page (two by two).

TODO: check macros have not changed

Right-click on the 'an_1x1.mv' icon and select the 'Visualise' menu item (see figure)

After a pause, this will generate a map showing mean-sea-level pressure (MSLP).

Use the an_1x1.mv or the an_2x2.mv macros to compare the ECMWF analyses with the satellite images.

Drag and drop the mv_track.mv icon onto the map to add the track of Hurricane Nadine.

...

Panel
titlePlot mslp and z500

With the edit window open, find the line that defines 'plot1':

Code Block
languagebash
#Define plot list (min 1- max 4)
plot1=["mslp"]                      # use square brackets when overlaying multiple fields per plot

Change this line to:

Code Block
plot1=["z500.s","mslp"]

The '.s' means plot the 500hPa geopotential as a shaded plot instead of using contours (this style is not available for all fields). Make sure 'mslp' is second to plot the contours on top of the z500 shaded colour map.

Click the play button and then animate the map that appears

Use the an_1x1.mv (or the an_2x2.mv) to plot fields of your choice and compare with the satellite images.

Panel
titlePlot PV at 320K

Change the value of "plot1" again to animate the PV at 320K (similar to Figure 13 in Pantillon et al).

Code Block
plot1=["pv320K"]

You might add the mslp or z500 fields to this plot e.g.

Code Block
plot1=["z500.s","pv320K","mslp"]

Note that the fields are plotted in the order specified in the list!

...

Exercise 1 looked at the synoptic development up to the 20-Sept-2012. This exercise looks at the ECMWF HRES forecast from this date and how the IFS model developed the interaction between Hurricane Nadine and the cut-off low.

Enter the folder 'HRES_forecast' in the 'openifs_2018' folder to begin.

...

Data is provided for a single 5 day forecast starting from 20th Sept 2012, as used in the paper by Pantillon et al. referred to at the start of this tutorial.

Warning

TODO: this may change if we get the data to day 8 in the forecasts.

Data is Data is provided at the same resolution as the operational model, in order to give the best representation of the Hurricane and cut-off low iterations. This may mean that some plotting will be slow.

...

The parameters available in the analyses are also available in the forecast data.

Available plot types

Warning

TODO: check macros have not changed


Panel

For this exercise, you will use the metview icons in the folder 'HRES_forecast' shown above.

hres_1x1.mv & hres_2x2.mv    : these work in a similar way to the same icons used in the previous task where parameters from a single lead time can be plotted either in a single frame or 4 frame per page.
hres_xs.mv
                                 : this plots a vertical cross section and can be used to compare the vertical structure of Hurricane Nadine and the cut-off low.

Task 1: Synoptic development (day 0-5)

Warning

TODO: include data to day 8 (if possible). Then have task1 to look at day 0-4, and include another task which looks at tracking the systems to day 8 in the forecast & overlay the manual_track.csv file with the mv_track.mv icon.

day 0 - 8 : track & synoptic development

MSLP, Z500, T500. key plots.

x-section

CPS diagram (NCL script -- not on VM)

From 20th to 25th September, Nadine had lost hurricane intensity and was slowly moving over the North Atlantic. A short-wave trough separated into a cutoff and approached from the north. Nadine interacted with a cold front associated with the cutoff on 21st September.Study the interaction of the systems and their subsequent development.

MSLP and geopotential height

First plot the geopotential at 500hPa and MSLP (as in figure 1 in Pantillon et al.) using the hres_1x1.mv macro script to see the large scale circulation and track the development of the systems:

Code Block
plot1=["z500.s","mslp"]

Temperature and mid-troposphere patterns

Next use the Using either the hres_1x1.mv or hres_2x2.mv macro plot some of these other maps to study the synoptic development. Note for example that there should not be a strong horizontal temperature gradient around Nadine (why?) and the winds should be stronger for Nadine than the cutoff.

...

borderColorlightblue
bgColorlightlightblue
titleBGColorlightblue
titleIsobaric maps

to plot additional fields. Right click and select 'Edit' and change the plot list to:

Code Block
plot1=["t2","mslp"]
plot2=["mslp","wind10"]
plot3=["t500","z500"]
plot4=["t850","z850"]
  • "t2", "mslp"            :  2m-temperature and MSLP - low level signature of Nadine and temperature.
  • "mslp", "wind10"

...

  •     : MSLP + 10m winds  -   interesting for Nadine's tracking and primary circulation.
  • "

...

  • t500","

...

  • z500"

...

  •         : Geopotential + temperature at 500hPa  -  large scale patterns, mid-troposphere position of warm Nadine and the cold Atlantic cutoff.
  • "t850", "z850"

...

  •        : Geopotential + temperature at 850hPa  -  lower level conditions, detection of fronts

...


Using either the hres_1x1.mv or hres_2x2.mv macro plot some of these other maps to study the synoptic development. Note for example that there should not be a strong horizontal temperature gradient around Nadine (why?) and the winds should be stronger for Nadine than the cutoff.

Panel
titleOther suggested isobaric maps
  • "vo850", "mslp"           : vorticity at 850hPa and MSLP  -  low level signature of Nadine and disturbance associated with the cutoff low.
  • "r700", "mslp",            : MSLP + relative humidity at 700hPa  -  with mid-level humidity of the systems
  • "pv320K", "mslp"        : 320K potential vorticity (PV) + MSLP  -   upper level conditions, upper level jet and the cutoff signature in PV, interaction between Nadine and the

It would be instructive to view a cross-section through Nadine and the cutoff. Choose from the following (note the cross-section macro hres_xs.mv uses slightly different names for the parameters)

Panel
bgColorlightlightblue
titleBGColorlightblue
titleVertical cross-sections
  • "pt", "pv"              : Potential temperature + potential vorticity to characterize the cold core and warm core structures of Hurricane Nadine and the cut-off low.
  • "rwind850", "ww700"                 : Humidity + vertical motion :   another view of the cold core and warm core structures of Hurricane Nadine and the cut-off low.
  • "pv", "w" ("r")       : Potential vorticity + vertical velocity (+ relative humidity) :  a classical cross-section to see if a PV anomaly is accompanied with vertical motion or not.

You do not need to plot all of these maps. Focus on what interests you and prepare some plots for group discussion.

  •       : Winds at 850hPa + vertical velocity at 700hPa (+MSLP) : focus on moist and warm air in the lower levels and associated vertical motion.


Warning

TODO:

If we get data for day 5-8, include another task here to look at vertical structure and track storm.

Overlay the manual_track.csv file on top of the mslp & z500 map.

This will lead to Fred's discussion of the CPS (cyclone phase space) diagram.

Fred/Etienne - Any other tasks / questions? I'm not sure what questions you want to ask.

Vertical structure

Warning

TODO:this becomes part of new task 2 if we get days 5-8 data.

This part of the task is to look at cross-sections through Hurricane Nadine and the cut-off low.

Right-click on the icon 'hres_xs.mv' icon, select 'Edit' and push the play Image Addedbutton.

The plot shows the potential vorticity (PV), wind vectors projected onto the plane of the cross-section and potential temperature drawn approximately through the centre of the Hurricane and the cut-off low. The red line on the map of MSLP shows the location of the cross-section.

Panel
borderColorred
borderStylesolid

Q. Look at the PV field, how do the vertical structures of Nadine and the cut-off low differ?

Changing forecast time

Cross-section data is only available every 24hrs.

This means the 'steps' value in the macros is only valid for the times:  [2012-09-20 00:00], [2012-09-21 00:00], [2012-09-22 00:00], [2012-09-23 00:00], [2012-09-24 00:00], [2012-09-25 00:00]

To change the date/time of the plot, edit the macro and change the line:

Code Block
steps=[2012-09-22 00:00]
Changing fields

A reduced number of fields is available for cross-sections compared to the isobaric maps: temperature (t), potential temperature (pt), relative humidity (r), potential vorticity (pv), vertical velocity (w), wind-speed (speed; sqrt(u*u+v*v)) and wind vectors (wind3).

Choose from the following (note the cross-section macro hres_xs.mv uses slightly different names for the parameters).

Panel
titleSuggestions for other vertical cross-sections
  • "pt", "pv"              : Potential temperature + potential vorticity to characterize the cold core and warm core structures of Hurricane Nadine and the cut-off low.
  • "r", "w"                 : Humidity + vertical motion :   another view of the cold core and warm core structures of Hurricane Nadine and the cut-off low.
  • "pv", "w" ("r")       : Potential vorticity + vertical velocity (+ relative humidity) :  a classical cross-section to see if a PV anomaly is accompanied with vertical motion or not.


Changing cross-section location
Code Block
#Cross section line [ South, West, North, East ]
line = [30,-29,45,-15]

The cross-section location (red line) can be changed in this macro by defining the end points of the line as shown above.

Remember that if the forecast time is changed, the storm centres will move and the cross-section line will need to be repositioned to follow specific features. This is not computed automatically, but must be changed by altering the coordinates above.

For help on how to save images, see the beginning of this tutorial.

Panel
borderColorred
borderStylesolid
titleQuestions to consider

Q. How strongly does Nadine appear to interact with the cutoff?
Q. What would you forecast to happen to Nadine?
Q. What is the fate of the cutoff low and what synoptic conditions did it create over France?

Task 2: Precipitation over France

Warning

TODO: if we get data to day 8, this becomes task 3.


Choose a hres macro to use, plot the total precipitation (parameter: tp), near surface wind field (parameter: wind10), relative humidity (parameter: r) (and any other parameters of interest).

...