source: util/src/CalculateDensity.sh@ 3f0c46

#!/bin/sh
#
# calculates the density of a molecule contained in a espack config file

if [ -z $1 ]; then
        echo "Usage: $0 <config file>"
        echo "We assume, that there is also an xyz file to the config file."
        exit 1
else
        config=$1
fi

mass=`/home/heber/workspace/ESPACK/bin/CalculateMass.sh ${config}.xyz | awk -F" " {'print $3'}`
echo "Mass is $mass atomicmassunit"
boxlengths=`grep ^Box --after-context=3 ${config} | tail -n 3`
volume=`echo $boxlengths | awk -F" " {'print $1*$3*$6'}`
if [ `grep IsAngstroem ${config} | awk -F" " {'print $2'}` == "1" ]; then
        unit="angstrom^3"
else
        unit="atomiclength^3"
fi
echo "Volume is $volume $unit"
density=`units "$mass atomicmassunit / $volume $unit" "g/cm^3" | grep \* | awk -F" " {'print $2'}`
echo "Density is $density g/cm^3"

exit 0