• The 4 Columns of the Data plotted separately:
• The Script used to plot:

Crowley 2000 Total Forcing updated April 2001

---------------------------------------------------------------------

World Data Center for Paleoclimatology, Boulder

                                    and

                NOAA Paleoclimatology Program

---------------------------------------------------------------------

NOTE: PLEASE CITE ORIGINAL REFERENCE WHEN USING THIS DATA!!!!!

NAME OF DATA SET: Causes of Climate Change Over the Past 1000 Years

LAST UPDATE: 4/2001 (Addition of file forc-total-4.12.01.txt)

CONTRIBUTORS: Thomas J. Crowley

IGBP PAGES/WDCA CONTRIBUTION SERIES NUMBER: 2000-045

SUGGESTED DATA CITATION: Crowley, T.J., 2000, Causes of Climate Change Over the Past 1000 Years, IGBP PAGES/World Data Center for Paleoclimatology Data Contribution Series #2000-045.   NOAA/NGDC Paleoclimatology Program, Boulder CO, USA.

ORIGINAL REFERENCE: Crowley, T.J., 2000, Causes of Climate Change Over the Past 1000 Years, Science, v.289, pp 270-277, July 14, 2000.

GEOGRAPHIC REGION: Global

PERIOD OF RECORD: 1000-1998 AD

FUNDING SOURCES: National Oceanic and Atmospheric Administration NA96GPO415 (USA),

Department of Energy (USA) through Battelle/Pacific Northwest Laboratory.

DESCRIPTION:

(April 2001 versions of estimates of volcano, solar, greenhouse gas, and tropospheric (1000-1998) total forcing prior to accounting for the planetary albedo (Tab-delimited ASCII format).

April 2001 file forc-total-4.12.01.txt contains estimates of volcano, solar, greenhouse gas, and tropospheric (1000-1998) total forcing, prior to accounting for the planetary albedo.   All units are in W/m**2. hl in volcanic time series refers to the fact that eruptions of unknown origin have been assigned a high latitude (hl) origin. There are "tails" to most of the large eruptions that were determined based on the estimated e-folding time of the aerosols as being about 1 year.

Sol.Be10 refers to the Beryllium 10 measurements of Bard et al. scaled to the Lean et al. changes over the last 400 years.

GHG refers to greenhouse gases.

Aer refers to tropospheric aerosols.

Year Vol.hl.cct Sol.Be10/Lean.splice GHG Trop.aer

1000 0.00 0.06 -0.05 0.00

1001 0.00 0.06 -0.05 0.00

1002 0.00 0.06 -0.05 0.00

1003 0.00 0.06 -0.05 0.00

[...]

1998 -0.07 0.46 2.43 -1.00

Data: ftp://ftp.ncdc.noaa.gov/pub/data/paleo/gcmoutput/crowley2000/forc-total-4_12_01.txt

Paper: http://www.ncdc.noaa.gov/paleo/pubs/crowley.html

The 4 Columns of the Data plotted separately:

The Script used to plot:

#!/bin/sh
IN=crowley2000-forc-total-4_12_01.txt
OUT1=${IN%.txt}-1
OUT2=${IN%.txt}-2
OUT3=${IN%.txt}-3
OUT4=${IN%.txt}-4

cat "$IN" | sed -n '53,$p' | awk '{print $1" "$2}' >"$OUT1"
cat "$IN" | sed -n '53,$p' | awk '{print $1" "$3}' >"$OUT2"
cat "$IN" | sed -n '53,$p' | awk '{print $1" "$4}' >"$OUT3"
cat "$IN" | sed -n '53,$p' | awk '{print $1" "$5}' >"$OUT4"

table2png "$OUT1" 1
table2png "$OUT2" 1
table2png "$OUT3" 1
table2png "$OUT4" 1

table2png:

#!/bin/sh
test -z "$1" && echo "Usage: table2png <datafile> [<lines>]" && echo 0
LINES=0
test -n "$2" && LINES=1
SYMBOL="1 .01"
test $LINES = 0 && SYMBOL="2 0.03"
X=X
Y=Y
DEL=tee
test -n "`head -n1 $1 | tr -d '+-0123456789.eE' | awk '{print $1}'`" && \
      X=`head -n1 $1 | awk '{print $1}'` \
      Y=`head -n1 $1 | awk '{print $2}'` \
      DEL='sed -e 1d'
REG=`regression $1 | awk '{print $11}' | cut -c1-5`

cat "$1" | $DEL | graph \
--display-type png \
--input-format a \
--bitmap-size "640x480" \
--width-of-plot .8 \
--height-of-plot .8 \
--right-shift .1 \
--upward-shift .1 \
--line-mode $LINES \
--font-name "HersheySans" \
--font-size .03 \
--title-font-size .03 \
--top-label "$1" \
--x-label "$X (r²=$REG)" \
--y-label "$Y" \
--symbol $SYMBOL \
--grid-style 2 >"$1.png"

regression.pas compiled with Free-Pascal: regression Signature: regression.sig