set term pdf enhanced color font "helvetica,20" linewidth 1 size 10,6
set output "si_coolingrate.pdf"
set key bottom

set xlabel "Cooling rate (K/day) [positive means heating]"
set ylabel "Altitude (km)"
plot \
     "./cooling_rate.txt" u 3:(($1+$2)/2000) title "total" w lines lt -1 lw 1 \
    ,"./gas_ground_coolingrate.txt" u 3:(($1+$2)/2000) title "gas-ground" w lines lt 0 lw 1 \
    ,"./gas_space_coolingrate.txt" u 3:(($1+$2)/2000) title "gas-space" w lines lt 1 lw 1 \
    ,"./gas_gas_coolingrate.txt" u 3:(($1+$2)/2000) title "gas-gas" w lines lt 2 lw 1 \

quit
