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
35 pfile=$pdir/pair_corr.scr
43 set title 'Pair correlation function'
45 set ylabel 'g(r) [a.u.]'
47 set terminal postscript eps enhanced color solid lw 1 'Helvetica' 14
48 set output '$pdir/pair_corr.eps'
51 echo -en "plot [1.0:6.0] " >> $pfile
55 for i in $pdir/pair_corr-*; do
57 time=`basename $i | awk -F- '{print $2}'`
59 if [ ! -z `echo $4 | grep a` ]; then
61 echo -en ", " >> $pfile
62 echo -en "\"$i\" u 1:2 w l t \"ab $time\"" >> $pfile
66 if [ ! -z `echo $4 | grep b` ]; then
68 echo -en ", " >> $pfile
69 echo -en "\"$i\" u 1:3 w l t \"aa $time\"" >> $pfile
73 if [ ! -z `echo $4 | grep c` ]; then
75 echo -en ", " >> $pfile
76 echo -en "\"$i\" u 1:4 w l t \"bb $time\"" >> $pfile
82 echo -en "\n" >> $pfile