The next is a mapping facility. To use it, you'll need to install the "maps" package. Installing a package in R is easy - just type (with R running)
install.packages("maps")and answer queries about mirror sites etc. With linux, this is better done as root.
If you don't want to do this, there's a variable at the top of the GlobTempLSv1.4.r called UseMaps - set it to F and everything but maps will be done.
It provides for each selection that you run a world map with the stations that qualified marked as dots. It's a useful check.
V1.4 also provides plot of number of stations in your sample, by year. This helps to keep from overinterpreting when numbers get small.
V1.4 also corrects an error with the ocean modified data, which had the years mislabelled by a year - it said 1881-2009 when it was actually showing 1880-2008. That made a minor difference to results. Plots look not much different now, but slightly better.
The selection function chooser() is simplified using a R switch statement. Now you just add lines to the list like:
"SST" = tv$urban=" ", # the last comma is needed
This works here because for sea "stations" many fields like urban are blank. As usual, the code is on the repository here, and at the document store.
A summary of this series of posts is here.
An SST exampleHere's a map of the artificial SST stations that are used in the global analyses, and a plot of station numbers over the years. Note the WW1 and WW2 dips
"GtLakes" = tv$urban=" " & getdist(42, -83),
Here's what the metadata shows: