4 # calculate pair correlation functions for a whole simulation output
5 # frank.zirkelbach@physik.uni-augsburg.de
9 echo "processing $1 ..."
10 ./pair_correlation_calc $1 $2
11 trg=`echo $1 | sed 's%s-%pair_corr-%' | sed 's%.save%%'`
12 mv pair_corr_func.txt $trg
17 for file in $1/*.save; do
28 if [ "$3" = "g" ]; then
31 pfile=$pdir/pair_corr.scr
39 set title 'Pair correlation function'
41 set ylabel 'g(r) [a.u.]'
43 set terminal postscript eps enhanced color solid lw 1 'Helvetica' 14
44 set output '$pdir/pair_corr.eps'
47 echo -en "plot [1.5:6.0] " >> $pfile
51 for i in $pdir/pair_corr-*; do
53 time=`basename $i | awk -F- '{print $2}'`
55 if [ ! -z `echo $4 | grep a` ]; then
57 echo -en ", " >> $pfile
58 echo -en "\"$i\" u 1:2 w l t \"ab $time\"" >> $pfile
62 if [ ! -z `echo $4 | grep b` ]; then
64 echo -en ", " >> $pfile
65 echo -en "\"$i\" u 1:3 w l t \"aa $time\"" >> $pfile
69 if [ ! -z `echo $4 | grep c` ]; then
71 echo -en ", " >> $pfile
72 echo -en "\"$i\" u 1:4 w l t \"bb $time\"" >> $pfile
78 echo -en "\n" >> $pfile