#!/usr/bin/env gnuplot

set term png transparent size 960, 720  \
        xffffff x000000 x404040  \
        xff0000 x0000ff x25bc24 x8e00ba x5f3e1a

set output '/dev/null'
set title '700 Bath Modes, by omega_max'
plot 'fig3_N700_max20wc_corr50000.data' w l ti 'max = 20wc'
rep 'fig3_N700_max60wc_corr50000.data' w l ti 'max = 60wc'
rep 'fig3_N700_max100wc_corr50000.data' w l ti 'max = 100wc'
rep 'fig3_N700_max10wc_corr50000.data' w l ti 'max = 10wc'
rep 'fig3_N700_max200wc_corr50000.data' w l ti 'max = 200wc'
set output 'fig3_N700_50ktraj.png'
replot

set output '/dev/null'
set title '1000 Bath Modes, by omega_max'
plot 'fig3_N1000_max20wc_corr50000.data' w l ti 'max = 20wc'
rep 'fig3_N1000_max60wc_corr50000.data' w l ti 'max = 60wc'
rep 'fig3_N1000_max100wc_corr50000.data' w l ti 'max = 100wc'
rep 'fig3_N1000_max10wc_corr50000.data' w l ti 'max = 10wc'
rep 'fig3_N1000_max200wc_corr50000.data' w l ti 'max = 200wc'
set output 'fig3_N1000_50ktraj.png'
replot

set output '/dev/null'
set title '200 Bath Modes, by omega_max'
plot 'fig3_N200_max20wc_corr50000.data' w l ti 'max = 20wc'
rep 'fig3_N200_max60wc_corr50000.data' w l ti 'max = 60wc'
rep 'fig3_N200_max100wc_corr50000.data' w l ti 'max = 100wc'
rep 'fig3_N200_max10wc_corr50000.data' w l ti 'max = 10wc'
rep 'fig3_N200_max200wc_corr50000.data' w l ti 'max = 200wc'
set output 'fig3_N200_50ktraj.png'
replot

set output '/dev/null'
set title '500 Bath Modes, by omega_max'
plot 'fig3_N500_max20wc_corr50000.data' w l ti 'max = 20wc'
rep 'fig3_N500_max60wc_corr50000.data' w l ti 'max = 60wc'
rep 'fig3_N500_max100wc_corr50000.data' w l ti 'max = 100wc'
rep 'fig3_N500_max10wc_corr50000.data' w l ti 'max = 10wc'
rep 'fig3_N500_max200wc_corr50000.data' w l ti 'max = 200wc'
set output 'fig3_N500_50ktraj.png'
replot