#!/usr/bin/env gnuplot

if (!exists("ntraj")) load 'corr_fcn_basics.gnuplot'

trajstep=2
lastframe=converged/2.5

if ((!exists("ntraj")) && (purpose == 'mov')) set output 'convergence.gif'
if (!exists("ntraj")) ntraj=1


p '../../frenkel/metrop_rmsd.data' u 2:($3+$4):($3-$4) every ::len*(ntraj-1)::len*ntraj-2 w filledcurves lt 3 noti, \
'mod_frenkel_rmsd.data' u 2:($3+$4):($3-$4) every ::len*(ntraj-1)::len*ntraj-2 w filledcurves lt 4 noti , \
'mod_frenkel_rmsd.data' u 2:3 every ::len*(ntraj-1)::len*ntraj-2 w l lt 2 lw 2 ti 'Frenkel (sorted)', \
'../../frenkel/metrop_rmsd.data' u 2:3 every ::len*(ntraj-1)::len*ntraj-2 w l lt 1 lw 2 ti 'Metropolis', \
'mod_frenkel_rmsd.data' u 2:3 every ::len*(converged-1)::len*converged-2 w l lt -1 lw 2 ti 'Converged (Frenkel, 100k trajs)'

ntraj = ntraj+trajstep
if (ntraj <= lastframe) reread