Wasnt that hard. green is land, blue is sea. new data awaits you. old data taken down For Zeke: Green is landcoast cells, blue is seacoast. The cells are AREA weighted per cce request but not Fractiionally weighted. That is if a cell is 90% land and 10% water, this weighting is NOT applied. I 
UPDATE: cce asked for a chart of coastal locations broken out by their SST component and their Land componement. That’s a little tricky, but once you figure out the mask logic its just addition and such.There are TWO masks: One mask contains ALL zeros, except for coast cells which are fraction of LAND.The other mask 
Replaces 4.0, minor changes made here and there. totally replaces 4.0 as the testing move forward. see 4.1 Run: DownloadAll. Run: DataSetUp. SSTTest reruns the SST analsyis and compares at HADSST2 Land1900Present: reruns the old version of program with a bench against CRU land figures. Next version will integrate more of raster to the analysis, 
UPDATE: Ron Broberg has a more definitive explanation of the difference which indicates that 5sig issue is not the main cause of the difference. See his exposition here CRU, it appears, trims out station data when it lies outside 5 sigma. Well, for certain years where there was actually record cold weather that leads to 
UPDATE: Ron Broberg has a more definitive explanation of the difference which indicates that 5sig issue is not the main cause of the difference. See his exposition here A short update. I’m in the process of integration the Land Analysis and the SST analysis into one application. The principle task in front of me is integrating 
Now that SST has been processed with raster we can return to the land processing and put that on the same footing. To get there I’ll rebuild from scratch making some minor adjustments along the way. As we go foward the amount of code we have to write should become less and less because we 
Reader P was kind enough to make us a new color map so I promptly played around with it and other parameters. Need to figure out how to drop the labels and ticks on the “map” map.axes() is no help. In anycase, I had a day long struggle with my R set up, its all 
Update: new zip, correcting bug found by Steve McIntyre: if(!file.exists(HadSST2ncdf)) downloadHADSST2() if(!file.exists(HadSST2ncdf)) downloadHadSST2() issue pending with another line as well. Checking raster versions. I’ve also, added some code into “downloadHadSST2″ that corrects for the “NA” problem with HadSST. (currently commented out). There is an issue with “ncdf” handling CF standards, which has been addressed in 
Zero Inflated Models and Generalized Linear Mixed Models with R.
Zuur, Saveliev, Ieno (2012).
