GISTEMP: There's a light burning in the fireplace
Herein lie some sample runs with the new GISTEMP code, comparing and contrasting my results with the web published, and comparing and contrasting the public GISS v2.inv file for metadata with my self-generated v2.giss.inv metadata file.
Downloaded the GISTEMP source
Unzipped the file.
Changed directory to GISTEMP_sources/STEP0/input_files
Downloaded GHCN v2.mean
Copy 9641C_201005_F52 to USHCNv2.avg
Downloaded new Antarctic data files
(use this handy-dandy script)
Code and Script
The scripts seem cleaner than the last time I looked at them. However there are two changes that still need to be made for most Linux boxen. First, change the “tail +100” to “tail -100” in the do_comb_step0.sh script. Second, search all the scripts for calls to executables and other scripts with out an explicit directory path. In other words, change calls like “zonav.exe” to “./zonav.exe”
The code is cleaner too. No need to monkey with PApars.f to remove the infinite loop for those of us with lower precision tools. The logic to handle the new ‘night light flags’ or retain the old behavior (US = nightlights, ROW = Rural/Urban) is in the text_to_binary.f code.
Don’t forget to compile the two python modules. I see that the older unused modules were removed. And speaking of compiling, don’t forget to set and export a variable FC to define your fortran compiler for the ‘jit’ compiling that occurs in the scripts. In my case: export FC=f95
You will have to create the “STEPx\temp_files” directories. In STEP 2, you will have to create the “work_files” directory.
GISS WEB -v- GISS LINUX
Comparing the GISS web release of GLB.Ts.txt with the GLB.Ts.GHCN.CL.txt generated on my gfortan Ubuntu 9.10 Sun Vbox using the public source code.
In all of these comparisons, red is the color of the first listed data source. Sorry about the odd scale on the ‘diff’ charts. The perils of sw reuse (these are some older perl scripts I really should rewrite in R)
GISS Old Rural/Urban -v- New Nightlights
Here we compare the linux-compiled GISS code using the traditional setting (US=nightlights, ROW = rural/urban)
GISS Rural -v- Whiteboard Rural
Here we use the older rural/urban method for PApars and compare the GISS v2.inv values as delivered in the source code with the newly generated Whiteboard v2.giss.inv values.
GISS Nightlights -v- Whiteboard Nightlights
Here we use the newer nightlights method for PApars and compare the GISS v2.inv values with the newly generated Whiteboard v2.giss.inv values.
There is very good matching between the ‘web’ version of the GISTEMP land record and the one that I generated – until we get past 1980. I suspect that this difference arises from different versions of the USHCN file being used.
There is relatively strong divergence between the GISS v2.inv files and my self-generated v2.giss.inv files in the early years of the record. I suspect that this is due to sensitivity to the definition of rural/urban or dark/bright when relatively few stations are available.
Overall, I am satisfied that my self-generated rural/urban and dark/bright fields are doing the same job as the GISS/GHCN generated values.