I guess I can safely say most of time is dedicated to making our radar website cooler. I’m pretty lucky. So we currently have three things on the website, the radar (obviously), and then two GFS visualizations. The first is a synoptic scale map of the circulation around South-Africa, mainly mimicking the South-African Weather Service (SAWS) Synoptic Map, but the SAWS is not as lazy as me so my maps are drawn by a computer and not by hand. The second product, and I’m fond of this, is the general circulation of the Southern-Hemisphere. These maps were inspired by the work of Harry von Loon and JJ Taljaard.
Now, here comes the actual reason I’m writing this. The third map that’s currently in experimental phase is a provincial scale map. I wanted to achieve two things; (1) allow people to more clearly what’s happing in their little part of the word, this is not really clear on the synoptic map, and (2) give people values they can understand more accurately and not rely on interpreting a scale, by value I mean a number. The map themselves are made with the same script as the other maps, well more or less the same, and you can read on that here. This post details how to get a value of a Lat/Lon point from GFS data, and in this case minimum and maximum temperature.
First we need to grap the GFS data, this pretty simple with wget. I assumed two things here, rightly or wrongly so, (1) that 0600 gfs product will be an accurate reflection of the day’s minimum temperature and (2) that 1200 will be an accurate reflection of the day’s maximum temperature. My notes in the script should make it obvious what’s being done at each stage, but the most important is to convert the data to netcdf and then make the values human friendly (Kelvin -> Celsius). I use the Apparent Temperature at 1200 as the map background. We’re not trying to make an animation here instead just a one day static map that’s updated every morning.
This gives us three values; tmin06.nc tmax12.nc and aptmp12.nc (the aptmp just serves as the background map, so no need to care too much). To get the values is the next part. Here we use Climate Data Operators (CDO) outputtab function to get a value on, or very close the specified Lat/Lon point. The value also had to be easily plotted on a map using GMT’s pstext function, this means I had to hack bash’s echo function. Here is the script.
So GMT reads a point value to plot via pstext as LAT LON FONTSIZE ANGLE FONTTYPE CAPTION so in this case I needed 27.09 -26.71 12 0 1 LM Potchefstroom in the script, as you see before every CDO command I echo this part with the details as needed. The sed command was to remove decimal places, but at the momement it removes the coordinate decimal places to, which I don’t want. The output looks like this (named nwtmp.dat):
Now the actual plot comes in, using GMT I do:
Here we have our map then
The pixels are a little ugly, the way forward here is to play with grdcontour and also the cpt file. The eventual plan is to make map like this for every province, but it takes time to enter the coordinates, as soon as that’s done it’s easy. Another issue is the .dat file, I need to figure out how to remove the temperature decimal places, but not the Lat/Lon. I think this is possible with AWK. I also plan to add MSLP and CAPE values to the output which could be interesting. It should not be to hard to add anything actually, percentage cloud cover could also be cool.
Ps. If you use the script give credit to GMT and CDO.