#!/bin/csh -f
'rm' -f wfn.gnu
set fl=GAwfn$1
echo $fl
set fl='"'$fl'"'
echo $fl
echo ' set xlabel "r (bohr)" ' > wfn.gnu
set g=`pwd`
set g=$g:t
echo $g
set ttt='Radial wavefunctions for '$g' with l='$1
echo $ttt
echo 'set title ' '"' $ttt '"' >>wfn.gnu
echo 'set terminal postscript color' >>wfn.gnu
echo 'set output ' '"''GAwfn'$1'.ps''"' >> wfn.gnu
echo 'plot 0,' $fl 'using 1:2 w l,' $fl 'using 1:3 w l,' $fl 'using 1:4 w l,' $fl 'using 1:5 w l,'$fl 'using 1:6 w l,'$fl 'using 1:7 w l,'$fl 'using 1:8 w l,'$fl 'using 1:9 w l' >>wfn.gnu 
echo 'set terminal x11'  >>wfn.gnu
echo 'plot 0,' $fl 'using 1:2 w l,' $fl 'using 1:3 w l,' $fl 'using 1:4 w l,' $fl 'using 1:5 w l,'$fl 'using 1:6 w l,'$fl 'using 1:7 w l,'$fl 'using 1:8 w l,'$fl 'using 1:9 w l' >>wfn.gnu 
echo ' ' >>wfn.gnu
echo 'pause -1 "******Press return to continue******"' >> wfn.gnu
gnuplot wfn.gnu