Publications by Steven Mosher
SST with Raster. Complete
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 comment...
4081 sym 50 img
Map colors
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 fixed. old packages stuck in a user library. So, Someone else asked for a 1 fi...
1109 sym 18 img
Moshtemp4.0
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 are letting raster do things for us. As much as I hate t...
3909 sym 16 img
5 Sigma in CRU
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 some new capabil...
2662 sym 22 img
1936.
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 discrepancies between CRU ...
2417 sym 24 img
MoshTemp 4.1
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, t...
1185 sym 16 img
Coastal
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 contains all zeros except th...
2593 sym R (304 sym/1 pcs) 18 img
CCE
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 think thats right for what he is after Related To ...
746 sym 20 img
Raster “Cover” function
The last bits of raster are coming together and everything is working so we are getting really close here. Today i will cover the raster function ‘cover’ which robert has just updated to handle a brick function. Lets review the situation with a toy example: First, we read in the Land data into a Brick. A brick is object with layers of rasters...
2482 sym 18 img
Handling .Z files
A while back Steve Mcintyre was looking for a way to handle .Z files in R Ron Broberg over at the whiteboard had an approach that steve adopted both for untar and for uncompressing .Z files. While the approach is slick, its somewhat of a hack. Nothing wrong with that, but I wanted something a bit more elegant. Long ago a reader Nicholas created...
3496 sym 16 img