# Graph of July insolation at 65 N, CO_2 and temperature anomalies # of the next 800 ka; gnuplot script by Jan Hollan # Data sources: # file -800-800ka.txt has been computed online at... This is based on # orbit91 is from: # Berger, A., 1992, Orbital Variations and Insolation Database. # IGBP PAGES/World Data Center-A for Paleoclimatology # Data Contribution Series # 92-007. NOAA/NGDC Paleoclimatology Program, # Boulder CO, USA. See # ftp://ftp.ncdc.noaa.gov/pub/data/paleo/insolation/ # and its file orbit91 (info on it is in readme_insolation.txt) # ORIGINAL REFERENCE: # Berger A. and Loutre M.F., 1991. # Insolation values for the climate of the last 10 million years. # Quaternary Sciences Review, Vol. 10 No. 4, pp. 297-317, 1991. # The very script to plot them, with Czech titles: set encoding iso_8859_2 set term pdf color enhanced font "Arial-Medium" set out "next800ka_cz65.pdf" set grid set title \ "Oslunění atmosféry v půli července na 65° pro příštích 800 tisíc let" set ylabel \ "ozářenost / 2 W/m^2 + 50, koncentrace / ppm, \n \ teplotní odchylka × (8/K) + 225" set xlabel "stáří / ka" plot [800:0] [140:370] \ "-800-800.txt" using (-$1):(50 + $6/2) \ with lines lt 1 lw 3 title 'ozářenost, půlka července, 65° sev. šířky',\ "t.txt" using ($3/1000):(225+$5*8) with lines lt 2 \ title 'teplotní odchylka v Antarktidě' ,\ "comp800.txt" using ($1/1000):2 with lines lt 3 \ title 'CO_2' quit