Against my better judgement I'll continue ...
Perhaps you should read the release notes and operating instructions. I expect if you're capable of reading source code in Fortran you're capable of finding them, yes? I certainly had no trouble, and it seems we disagree about what's in that source code. Do you really want me to post annotated code, and the source files and line numbers it came from? You do realize this will make you look like an idiot, correct?
Yes, please do post some *specific* code details where you believe I am in error.
Please note that *THIS*
HTML:
http://data.giss.nasa.gov/gistemp/sources/
HTML:
http://data.giss.nasa.gov/gistemp/sources/GISTEMP_sources.tar.gz
Is the code I specifically referenced.
This is, as far as I can tell, the earliest version of Hansen's code released. Hansen, according to reports, objected to the release as it needed additional work. From comments within the text files this code and related data sets were updated in 2003 if not later. Funny that he would draw economically devastating conclusions from code that was not, 13 years later, ready for review. I am not interested in the later re-writes that have appeared for the purpose of this discussion.
Here is a little help for you. This is a list of every file included:
./GISTEMP_sources
./GISTEMP_sources/STEP0
./GISTEMP_sources/STEP0/USHCN2v2.f
./GISTEMP_sources/STEP0/antarc_comb.f
./GISTEMP_sources/STEP0/antarc_comb.sh
./GISTEMP_sources/STEP0/antarc_to_v2.sh
./GISTEMP_sources/STEP0/cmb.hohenp.v2.f
./GISTEMP_sources/STEP0/cmb2.ushcn.v2.f
./GISTEMP_sources/STEP0/dif.ushcn.ghcn.2005.f
./GISTEMP_sources/STEP0/do_comb_step0.sh
./GISTEMP_sources/STEP0/dump_old.f
./GISTEMP_sources/STEP0/get_USHCN
./GISTEMP_sources/STEP0/get_offset_noFIL
./GISTEMP_sources/STEP0/hohp_to_v2.f
./GISTEMP_sources/STEP0/input_files
./GISTEMP_sources/STEP0/input_files/Ts.discont.RS.alter.IN
./GISTEMP_sources/STEP0/input_files/Ts.strange.RSU.list.IN
./GISTEMP_sources/STEP0/input_files/antarc1.list
./GISTEMP_sources/STEP0/input_files/antarc1.txt
./GISTEMP_sources/STEP0/input_files/antarc2.list
./GISTEMP_sources/STEP0/input_files/antarc2.txt
./GISTEMP_sources/STEP0/input_files/antarc3.list
./GISTEMP_sources/STEP0/input_files/antarc3.txt
./GISTEMP_sources/STEP0/input_files/combine_pieces_helena.in
./GISTEMP_sources/STEP0/input_files/mcdw.tbl
./GISTEMP_sources/STEP0/input_files/preliminary_manual_steps.txt
./GISTEMP_sources/STEP0/input_files/sumofday.tbl
./GISTEMP_sources/STEP0/input_files/t_hohenpeissenberg_200306.txt_as_received_July17_2003
./GISTEMP_sources/STEP0/input_files/ushcn.tbl
./GISTEMP_sources/STEP0/input_files/v2.inv
./GISTEMP_sources/STEP0/step0_README.txt
./GISTEMP_sources/STEP0/to_next_step
./GISTEMP_sources/STEP0/work_files
./GISTEMP_sources/STEP1
./GISTEMP_sources/STEP1/EXTENSIONS.tar.gz
./GISTEMP_sources/STEP1/PYTHON_README.txt
./GISTEMP_sources/STEP1/alter_discont.py
./GISTEMP_sources/STEP1/bdb_to_text.py
./GISTEMP_sources/STEP1/comb_pieces.py
./GISTEMP_sources/STEP1/comb_records.py
./GISTEMP_sources/STEP1/do_comb_step1.sh
./GISTEMP_sources/STEP1/drop_strange.py
./GISTEMP_sources/STEP1/input_files
./GISTEMP_sources/STEP1/input_files/Ts.discont.RS.alter.IN
./GISTEMP_sources/STEP1/input_files/Ts.strange.RSU.list.IN
./GISTEMP_sources/STEP1/input_files/combine_pieces_helena.in
./GISTEMP_sources/STEP1/input_files/mcdw.tbl
./GISTEMP_sources/STEP1/input_files/sumofday.tbl
./GISTEMP_sources/STEP1/input_files/ushcn.tbl
./GISTEMP_sources/STEP1/input_files/v2.inv
./GISTEMP_sources/STEP1/listStats.py
./GISTEMP_sources/STEP1/to_next_step
./GISTEMP_sources/STEP1/v2_to_bdb.py
./GISTEMP_sources/STEP1/work_files
./GISTEMP_sources/STEP2
./GISTEMP_sources/STEP2/PApars
./GISTEMP_sources/STEP2/PApars.f
./GISTEMP_sources/STEP2/do_comb_step2.sh
./GISTEMP_sources/STEP2/flags.f
./GISTEMP_sources/STEP2/input_files
./GISTEMP_sources/STEP2/input_files/v2.inv
./GISTEMP_sources/STEP2/invnt.f
./GISTEMP_sources/STEP2/padjust
./GISTEMP_sources/STEP2/padjust.f
./GISTEMP_sources/STEP2/split_binary.f
./GISTEMP_sources/STEP2/t2fit.f
./GISTEMP_sources/STEP2/text_to_binary.f
./GISTEMP_sources/STEP2/toANNanom
./GISTEMP_sources/STEP2/toANNanom.f
./GISTEMP_sources/STEP2/to_next_step
./GISTEMP_sources/STEP2/tr2.f
./GISTEMP_sources/STEP2/trim_binary.f
./GISTEMP_sources/STEP2/work_files
./GISTEMP_sources/STEP3
./GISTEMP_sources/STEP3/annzon.f
./GISTEMP_sources/STEP3/do_comb_step3.sh
./GISTEMP_sources/STEP3/input_files
./GISTEMP_sources/STEP3/results
./GISTEMP_sources/STEP3/to.SBBXgrid.f
./GISTEMP_sources/STEP3/to_next_step
./GISTEMP_sources/STEP3/trimSBBX
./GISTEMP_sources/STEP3/trimSBBX.f
./GISTEMP_sources/STEP3/work_files
./GISTEMP_sources/STEP3/zonav
./GISTEMP_sources/STEP3/zonav.f
./GISTEMP_sources/STEP4_5
./GISTEMP_sources/STEP4_5/SBBXotoBX.f
./GISTEMP_sources/STEP4_5/annzon.f
./GISTEMP_sources/STEP4_5/convert.HadR2_mod4.upto15full_yrs.f
./GISTEMP_sources/STEP4_5/convert1.HadR2_mod4.f
./GISTEMP_sources/STEP4_5/do.mult_year.TocnHR2.upd
./GISTEMP_sources/STEP4_5/do_comb_step4.sh
./GISTEMP_sources/STEP4_5/do_comb_step5.sh
./GISTEMP_sources/STEP4_5/input_files
./GISTEMP_sources/STEP4_5/input_files/SBBX_LtSN.LnWE.dat.gz
./GISTEMP_sources/STEP4_5/input_files/oisstv2_mod4.clim.gz
./GISTEMP_sources/STEP4_5/trimSBBX
./GISTEMP_sources/STEP4_5/trimSBBX.f
./GISTEMP_sources/STEP4_5/zonav
./GISTEMP_sources/STEP4_5/zonav.f
./GISTEMP_sources/gistemp.txt
Here is a list of every occurence of "trace" in the tarball:
# grep -n -r -i trace ./GISTEMP_sources/*
./GISTEMP_sources/STEP2/PApars.f:34:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP2/PApars.f:127: TRACE=INFO(8)
./GISTEMP_sources/STEP2/PApars.f:193:C?*** Change data if necessary (e.g. trace flag for precip)
./GISTEMP_sources/STEP2/PApars.f:194:C?PRC IF(RDATA(M).EQ.TRACE) RDATA(M)=0.
./GISTEMP_sources/STEP3/annzon.f:34:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP3/to.SBBXgrid.f:56:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP3/to.SBBXgrid.f:162: TRACE=INFOI(8)
./GISTEMP_sources/STEP3/to.SBBXgrid.f:256:C?*** Change data if necessary (e.g. trace flag for precip)
./GISTEMP_sources/STEP3/to.SBBXgrid.f:257:C?PRC IF(RDATA(M).EQ.TRACE) RDATA(M)=0.
./GISTEMP_sources/STEP3/trimSBBX.f:37:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP3/zonav.f:45:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP3/zonav.f:112: TRACE=INFOI(8)
./GISTEMP_sources/STEP4_5/SBBXotoBX.f:55:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP4_5/SBBXotoBX.f:140: TRACE=INFO(8)
./GISTEMP_sources/STEP4_5/annzon.f:34:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP4_5/trimSBBX.f:37:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP4_5/zonav.f:45:C**** 8 = flag for precipitation trace
./GISTEMP_sources/STEP4_5/zonav.f:112: TRACE=INFOI(8)
So where are "the tracers" you claim appear here ? Please indicate the file and line number specifics. I insist.
Sonny, I was playing with Fortran when you weren't a gleam in your daddy's eye. I suggest you not try to tell your grandma how to suck eggs.
You apparently believe your stale old programming skills a special asset ? It's an idiotic claim IMO. FWIW Fortran wasn't invented until after my birth and I wrote Fortran2 on an IBM1620 in 1965/66. The difference is that I've learned better methods since and you apparently have stuck with 1960 methodology.
Welcome to scientific programming. They don't care if it's pretty. What they care about is the algorithms, and that those algorithms be transparent enough that if they need to update them for a change in the understanding of the physics, they can do so quickly. And that it run efficiently- and it will most likely do that, no matter how nasty the commenting is.
They aren't interested in the basic structure. They're interested in whether they can find what they need to modify if the physics changes, and if it runs fast. Everything else is candy.
What an idiotic claim. That isn't scientific programming, that's just old fashioned, quick & dirty, bad programming. Apparently you haven't learned much about software design since your ancient Fortran days. If anyone tried to pass code like this into a MRI image reconstruction system, or even a graduate thesis they'd be canned or caned. If you want to compare notes I'm quite certain you are the one lacking any recent experience with scientific software development (It's been part of my education and career for 37 years).
Not only do you ignore the release notes, the FAQ, and the manual for the software, you have the temerity to lie about it. Seriously, I hope you're not planning on sticking around- because I'm going to make an idiot of you if you do. It's not very pleasant to watch. I'll need the heavy rubber gloves; I think we can forgo the vaseline, all things considered.
Bring it on, you foulmouthed anal-oriented buffoon. *IF* you to stick to specifics instead of your wandering allusions to these secondary documents, then everyone will see that you are wrong. Show me the code !
Just point to the source code lines in the specified tarball that indicate "tracers" to you. That is all I ask. A file name and a line number will do.
You are cute, I'll give you that- but you've missed something important, which is the integrity to admit when I'm wrong and move on.
I've just given you the opportunity to prove your integrity again. You are clearly wrong - the code indicated has no tracers and microscopic level of diagnostics, and is just poorly structured.
But let's not cloak ignorance in nobility. You went to all the trouble to type in a response, including reasonably correct eq'n to a very simple physics (asteroid) question and got it blatantly wrong. No typo, it's clear that you don't know the material well enough to reason about it, and didn't even read/understand the eq'n you typed. Just like a parrot, as I said ! With a parrot at least we can put a cloth over the cage to have it shut up.
If your reply doesn't contain the source file name from above and line number of your mythical "tracers" then we can all see that you are wrong (not to mention remarkably rude and loud).
-S