#
# Choas.gnu by Mathew Peet
# Based on Plot template of Ananth
#
# Use 'gnuplot choas.gnu > out.ps'
#
#

set term postscript enh "Helvetica" 12

set nokey

#
# This sets the border of the graph itself 
# The lines x and y axis, x2 and y2 axis 
# thickness can be changed like this
#

set border -1 lw 1


#### to set the size of the graph

# set size ratio 1.5

#### set line style 

set ls 1 lt 1 lw 2 pt 1 ps 1 
set ls 2 lt 2 lw 2 pt 1 ps 1 
set ls 3 lt 1 lw 2 pt 1 ps 1 
set ls 4 lt 2 lw 2 pt 1 ps 1

#### increase error bars tics

#set bar 2
#show bar

#### put minor tics

#set mxtics 4    
set nomxtics
set nomytics

# Removing the label beside the tics for each axis
# %g is normal, %e sets exponential format
# %0.1e means to 1 decimal place in exponential format
#

#set format y " "
#set format x " "

#### to put grid lines and also on the minor axis tics

#set grid mxtics 2 lw .1

### rotate the xlabel

#set xtics rotate 4 

### set title atrequired place ###

#set title "g -> d" 9,-3  "Symbol, 30" 




set yrange[0:1]
#set xrange[3:33]  
#set xtics 0.2 
#set logscale x   
#set xlabel "r"  
#set ylabel "x"


#
# To put labels at desired places 
#
# set label 1 "WQ" at graph 0,-0.09 rotate 
# set label 2 "100 ^{o}C\n15 min" at graph 0.14,-0.15 rotate  
# set label 3 "200 ^{o}C\n15 min" at graph 0.32,-0.15 rotate 
# set label 4 "250 ^{o}C\n15 min" at graph 0.47,-0.15 rotate 
# set label 5 "350 ^{o}C\n1 h" at graph 0.65,-0.15 rotate 
# set label 6 "350 ^{o}C\n2 h" at graph 0.83,-0.15 rotate  
# set label 7 "500 ^{o}C\n3 h" at graph 0.98,-0.15 rotate 


#
# Using title here will override the default title by
# gnuplot which is the data file itself
#



set pointsize 0.05

set multiplot
set origin 0,0
set size 1,1
set xlabel "r"
set ylabel "x"
plot 'out' using ($3<300 ? 1/0: $1):2 title "Populations"

set size 0.25,0.3
set origin 0.1,0.66
set xlabel "i"
set ylabel "x"
set title "Steady State"
plot 'out' using ($3>20 ? 1/0: $3):($1>1.4011 ? 1/0:$1<1.4 ? 1/0: $2) title "Converging to Steady State" w l

set size 0.25,0.3
set origin 0.30,0.66
set xlabel "i"
set ylabel ""
set title "Period 2
plot 'out' using ($3>20 ? 1/0: $3):($1>3.30011 ? 1/0:$1<3.30 ? 1/0: $2) title "Converging to Period 2" w l

set size 0.25,0.3
set origin 0.5,0.66
set xlabel "i"
set ylabel ""
set title "Period 4"
plot 'out' using ($3>20 ? 1/0: $3):($1>3.50011 ? 1/0:$1<3.5 ? 1/0: $2) title "Converging to Period 4" w l

set size 0.25,0.3
set origin 0.1,0.34
set xlabel "i"
set ylabel "x"
set title "Choatic"
plot 'out' using ($3>20 ? 1/0: $3):($1>3.80011 ? 1/0:$1<3.8 ? 1/0:  $2) title "Becoming Choatic" w l