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.]'
42 set terminal postscript eps enhanced color solid lw 1 'Helvetica' 14
43 set output '$pdir/pair_corr.eps'
46 echo -en "plot [1.5:3.0] " >> $pfile
50 for i in $pdir/pair_corr-*; do
52 time=`basename $i | awk -F- '{print $2}'`
54 if [ ! -z `echo $4 | grep a` ]; then
56 echo -en ", " >> $pfile
57 echo -en "\"$i\" u 1:2 w l t \"ab $time\"" >> $pfile
61 if [ ! -z `echo $4 | grep b` ]; then
63 echo -en ", " >> $pfile
64 echo -en "\"$i\" u 1:3 w l t \"ab $time\"" >> $pfile
68 if [ ! -z `echo $4 | grep c` ]; then
70 echo -en ", " >> $pfile
71 echo -en "\"$i\" u 1:4 w l t \"ab $time\"" >> $pfile
77 echo -en "\n" >> $pfile