...
We finish the case study by looking into the predictability of the large scale flow pattern by generating spaghetti plots for 500 hPa geopotential from the same ENS run as we investigated before. In a spaghetti plot each ENS member is rendered into the same map using a single isoline value. The plot we want to generate is shown below (it contains the spaghetti plot for 500 hPa geopotential using the 560 gpm isoline value):
...
Create a new Macro and edit it. Drop your Geographical View and the Coastlines icons into the Macro editor . Once the code is generated and tidied up change the map area to
| Code Block |
|---|
[40,-40,70,20] |
...
In this mcont() we turned contour labels off so that to keep the plot uncluttered and defined only a single contour value.
...
| Code Block |
|---|
plot(your_view,f,cont) |
By default if no title definition is specified Metview adds a title line for each filed in the plot. Since we are about to plot 50 fields this would result in 50 titles in the plot. To avoid having too many titles we use a custom Text Plotting icon:
p, li { white-space: pre-wrap; } title=mtext(text_line_1: "Value: 560 gpm T+<grib_info key='step' where='number=50' /> h" )
If we ran the Macro we would see far too many titles in the plot If you run this macro you will see that the "spaghetti" was properly generated but you have too many titles (one for each layer i.e. 50!). To get rid of prevent the extra titles you from popping up we need to define a Text Plotting inside the loop and add it to your plot() command:
...