set term png

#INPUTS
#=====================
#XC carbon content in mole fraction 
XC = 0.0109
#assuming amount in martensite is the same as bulk (calculated by mucg)

set grid
set ylabel "FT0 (J/mol)"
set xlabel "Temperature (degrees C)"
set key left
set output "FT0_GP.png"
plot \
"data1.txt" using 4:9 title "steps of 1 degrees 0 celcius", \
"data2.txt" using 4:9 title "steps of 20 degrees from 200 celcius", \
-400 title "stored energy bainite", \
-1120-10563*XC+94 title "total stored energy martensite", \
-1120 title "-1120 strain energy martensite"

set term postscript enhanced col
set output "FTO_GP.ps"
set ylabel "FT0 / Jmol^{-1}"
set xlabel "Temperature /^oC"
set nogrid
plot \
"data1.txt" using 4:9 title "steps of 1 degrees 0 celcius", \
"data2.txt" using 4:9 title "steps of 20 degrees from 200 celcius", \
-400 title "stored energy bainite", \
-1120-10563*XC+94 title "total stored energy martensite", \
-1120 title "-1120 strain energy martensite"