# Gnuplot script file for plotting Results #{ #cd 'C:\Users\r_col\Dropbox\University_of_Manchester\PhD\Chaotic_Flow\Notes\Chaotic Flow\Intermediate_Da\Figures\Data' #load 'PLOT_ActiveLinks.plt' ######################################################################################## #} ################################## General Instruction ################################# #{ reset set key opaque autotitle columnheader set term qt size 956,640 position -7,0 font ',14' unset grid set palette defined ( 0 "dark-orange", 2.5 "orange", 5 "yellow", 7.5 "green", 15 "blue", 20 "red", 25 "black") #set palette defined ( -10 "magenta", -5 "dark-violet", 5 "blue", 10 "cyan", 15 "green", 20 "yellow", 25 "orange", 30 "red", 40 "black") set cbrange [0:25] unset colorbox d ='.dat' arrows=0; Bd="GlobalBd_"; set style line 1 lw 3 dt 1 lc rgb "#CC0000" dB="LocaldB_"; set style line 2 lw 3 dt 2 lc rgb "#0088FF" bD="LocalbD_"; set style line 3 lw 3 dt 2 lc rgb "#7700EE" Db="GlobalDb_"; set style line 4 lw 3 dt 1 lc rgb "#00AA22" BD="DualBD_"; set style line 5 lw 3 dt 1 lc rgb "#00FF00" DB="DualDB_"; set style line 6 lw 3 dt 1 lc rgb "#FF00FF" poR="Random" poL="Lattice" ######################################################################################## #} disorder = 0 ; actlinks = 1 ; Unconditional = 1 ; Fixed_Wildtype = 0 ; Fixed_Mutant = 0 ; damkohler = 0 ; PT = bD ; processes = 1 ; Da = 5 ; if(disorder){ set ylabel "Disorder parameter (rho=(L_{A}/L{T})/(2i(1-i))"; ymin=0.01; ymax=1.3; set ytics 0.1 set logscale y lx=0.045 set label "Mixed" at screen lx,0.82 font ',11' rotate by 90 set arrow from screen lx,0.89 to screen lx,0.96 filled size screen 0.015,15 set label "Clustered by species" at screen lx,0.2 font ',11' rotate by 90 set arrow from screen lx,0.19 to screen lx,0.11 filled size screen 0.015,15 if(Unconditional){col=10;} if(Fixed_Wildtype){col=17;} if(Fixed_Mutant){col=24;} } if(actlinks){ set ylabel "Fraction of active links (L_{Act} /L_{Tot})" offset 0.5,0 ymin=0; ymax=0.08; set ytics 0.02 unset logscale y unset label unset arrow if(Unconditional){col=5;} if(Fixed_Wildtype){col=12;} if(Fixed_Mutant){col=19;} } set xlabel "Evolutionary Sweeps (tau=t*Da)" set logscale x set format x "10 ^{%.1T}" f1='./ActiveLinks/NetOrder_' f2='pshear_Da=1.00e' f3='_taumax=5.00e-1_Ro=0.110' d='.dat' R='_Random'; C='_Connected'; if (damkohler) { set key bottom left maxrows 4 width 1.5 set label PT at screen 0.5,0.92 font ',18' set cbrange [-2:12] plot [0.1:1e4][ymin:ymax]\ for[i=2:4:1] f1.PT.f2.i.f3.C.d u ($2):(column(col)) title 'Da=1e'.i.'_{C}' w l lw 3 dt 1 lc palette cb 2*i ,\ for[i=2:4:1] f1.PT.f2.i.f3.R.d u ($2):(column(col)) title 'Da=1e'.i.'_{R}' w l lw 2 dt 0 lc palette cb 2*i+1 ,\ } if (processes) { set key bottom left maxrows 4 width 1.5 set label 'Da=10 ^{'.Da.'}' at screen 0.5,0.92 font ',18' set cbrange [-2:12] plot [0.1:1e4][ymin:ymax]\ f1.Bd.f2.Da.f3.C.d every 4 u ($2):(column(col)) title 'Bd_{C}' w l ls 1 ,\ f1.Db.f2.Da.f3.C.d every 4 u ($2):(column(col)) title 'Db_{C}' w l ls 4 ,\ f1.bD.f2.Da.f3.C.d every 4 u ($2):(column(col)) title 'bD_{C}' w l ls 3 ,\ f1.dB.f2.Da.f3.C.d every 4 u ($2):(column(col)) title 'dB_{C}' w l ls 2 ,\ f1.Bd.f2.Da.f3.R.d every 4 u ($2):(column(col)) title 'Bd_{R}' w l ls 1 lw 2 dt 0 ,\ f1.Db.f2.Da.f3.R.d every 4 u ($2):(column(col)) title 'Db_{R}' w l ls 4 lw 2 dt 0 ,\ f1.bD.f2.Da.f3.R.d every 4 u ($2):(column(col)) title 'bD_{R}' w l ls 3 lw 2 dt 0 ,\ f1.dB.f2.Da.f3.R.d every 4 u ($2):(column(col)) title 'dB_{R}' w l ls 2 lw 2 dt 0 } #pause 2 #reread