#!/bin/bash
#PBS -q regular
#PBS -l walltime=22:00:00
#PBS -l mppwidth=8
#PBS -V

cd $PBS_O_WORKDIR

#title Debye Spin-Boson (Old Discretization)
#tag debye

#blog <p>This is a quick run of the old discretization of the Debye
#blog spin-boson, attempting to reproduce the results of JCP 110 4828, fig
#blog 3, bottom line.</p>

#blog <p>Running with revision 316.</p>

DYNAMIQ=${DYNAMIQ_BASEDIR}/branches/dwhs/spinboson/dynamiq


aprun -n 1 ${DYNAMIQ} old_fig3_N700_max100wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N700_max140wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N700_max200wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N700_max20wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N700_max60wc.inp &

aprun -n 1 ${DYNAMIQ} old_fig3_N100_max100wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N100_max140wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N100_max200wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N100_max20wc.inp 
aprun -n 1 ${DYNAMIQ} old_fig3_N100_max60wc.inp  &

#aprun -n 1 ${DYNAMIQ} old_fig3_N200_max100wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N200_max140wc.inp 
#aprun -n 1 ${DYNAMIQ} old_fig3_N200_max200wc.inp &
#aprun -n 1 ${DYNAMIQ} old_fig3_N200_max20wc.inp &
#aprun -n 1 ${DYNAMIQ} old_fig3_N200_max60wc.inp 

aprun -n 1 ${DYNAMIQ} old_fig3_N500_max100wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N500_max140wc.inp
aprun -n 1 ${DYNAMIQ} old_fig3_N500_max200wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N500_max20wc.inp &
aprun -n 1 ${DYNAMIQ} old_fig3_N500_max60wc.inp

wait

gnuplot old_debye_bywmax.gnuplot
gnuplot old_debye_byNmax.gnuplot
gnuplot old_debye_byTraj.gnuplot

# for researchblog.py : old_debye_test_1000.job

researchblog.py old_debye_test.job