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

cd $PBS_O_WORKDIR

#title Debye Spectral Density (JCP110 4828 fig4)

#tag Debye
#tag DebyeOld

#blog <p>An attempt to reproduce figure 4 also. For all all of these, the 
#blog maximum frequency is 150 times the cutoff frequency. (omega_max = 150
#blog omega_c)</p>

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

aprun -n 8 ${DYNAMIQ}/old_fig4a.inp &
aprun -n 8 ${DYNAMIQ}/old_fig4b.inp &
aprun -n 8 ${DYNAMIQ}/old_fig4c.inp &
aprun -n 8 ${DYNAMIQ}/old_fig4d.inp &
aprun -n 8 ${DYNAMIQ}/old_fig4e.inp &

aprun -n 8 ${DYNAMIQ}/fig4a.inp &
aprun -n 8 ${DYNAMIQ}/fig4b.inp &
aprun -n 8 ${DYNAMIQ}/fig4c.inp &
aprun -n 8 ${DYNAMIQ}/fig4d.inp &
aprun -n 8 ${DYNAMIQ}/fig4e.inp 

wait

gnuplot fig4.gnuplot

researchblog.py fig4.job