reset
erase
lweight 1.5
data source_hot.out
scolumn 2
data heat_hot.out
xcolumn 2
sevaluate 3.46e4*(x/s)
expand 1.25
angle 0
axis -5.0 5.0 0.5 2.5 0.15 0.15 0.7 1 0
axis -5.0 5.0 0.5 2.5 0.15 0.85 0.7 0 1
angle 90
axis 300. 700. 100. 100. 0.15 0.15 0.7 2 1
axis 300. 700. 100. 100. 0.85 0.15 0.7 0 0
angle 0
ltype 0
location 0.15 0.85 0.15 0.85
limits -5.0 5.0 300. 700. 
relocate 0. 300. 
draw 0. 700. 
EXPAND 1.5
RELOCATE -0.00 265.
PUTLABEL 2 "HEATING RATE (DEGREES/DAY)"
ANGLE 90
RELOCATE -6. 500.
PUTLABEL 8 "ALTITUDE (km)"
EXPAND 1.0
ANGLE 0

# conversion to degress per day 

data source_hot.out
scolumn 2
data heat_hot.out
xcolumn 2
sevaluate 3.46e4*(x/s)

# hydrocarbon cooling rates

ycolumn 1
xcolumn 7
xevaluate s*x
ltype 1
connect

xcolumn 8
xevaluate s*x
ltype 2
connect

xcolumn 9
xevaluate s*x
ltype 4
connect

# solar gas heating rate

xcolumn 6
ycolumn 4
xevaluate x+y
xevaluate s*x
ycolumn 1
ltype 0
connect


# solar dust heating rate

xcolumn 12 
xevaluate s*x
ltype 0
connect

# RELOCATE -0.099 432.
# PUTLABEL 1 "C\d2H\d2"
# RELOCATE 0.089 424. 
# PUTLABEL 3 "C\d2H\d6"
# RELOCATE 0.1366 633.
# PUTLABEL 3 "SOLAR UV"
# RELOCATE 0.047 654.4149
# PUTLABEL 9 "SOLAR IR"
# RELOCATE -0.076 509.
# PUTLABEL 9 "CH\d4"