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

cd $PBS_O_WORKDIR

#title Debye Spectral Density (JCP110 4828 fig2)

#tag Debye
#tag DebyeOld

#blog <p>After previous failed attempts to reproduce the results of the
#blog bottom line of fig 2 from JCP110 4828, I'm trying to reproduce
#blog <emph>all</emph> the lines from that figure. We'll see if I see
#blog behavior that is qualitatively similar as eta increases.</p>

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

aprun -n 16 ${DYNAMIQ} old_fig2a.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} old_fig2b.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} old_fig2c.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} old_fig2d.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} old_fig2e.inp > /dev/null &

aprun -n 16 ${DYNAMIQ} fig2a.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} fig2b.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} fig2c.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} fig2d.inp > /dev/null &
aprun -n 16 ${DYNAMIQ} fig2e.inp > /dev/null &

wait

gnuplot fig2.gnuplot

researchblog.py fig2.job