Chimie Théorique » OpenPath

#!/bin/bash

if [ $# -lt 1 ]; then

 echo "Use: $0 File.out [x column]"

 exit

fi

Fout=$1

ItMax=`grep MAXCYC $Fout | tail -1 | awk '{print $3}'`

Nom=`grep PATHNAME $Fout | tail -1 | awk '{print $3}'`

NGeomF=`grep NGEOMF $Fout | tail -1 | awk '{print $3}'`

xcol=1

if [ $# -eq 2 ]; then

xcol=$2

fi

export LANG=C

echo "#ItMax=$ItMax"

echo "#Nom=$Nom"

echo "#NGeomF"=$NGeomF

if [ -s ${Nom}_cart.0 ]; then

 Ext=_cart

elif [ -s ${Nom}.0 ]; then

 Ext=""

else

 echo "Cannot find ${Nom}.0 nor ${Nom}_cart.0 -- ERROR"

 exit

fi

 echo "Using files ${Nom}${Ext}.XX"

if [ -s $Nom.datl ]; then

 it=1

 while [ -s $Nom.datl_${it} ]

  do

 let it=it+1

 done

 echo "Moving $Nom.datl into $Nom.datl_${it}"

 mv $Nom.datl $Nom.datl_${it}

fi

echo "Creating $Nom.datl"

ItDone=-1

for i in `seq 0 $ItMax`

 do

 if [ -s ${Nom}${Ext}.$i ]; then

     xyz2path  ${Nom}${Ext}.$i 

     cat Scan.dat >> $Nom.datl

     echo " " >> $Nom.datl

     echo " " >> $Nom.datl

     let ItDone=ItDone+1

 fi

done 

icol=`tail -1 Scan.dat | wc -w`

cat << EOF > ${Nom}_l.gplot

#!/usr/bin/gnuplot -persist

 set pointsize 2

 xcol=$xcol

EOF

for i in `seq 1 $ItDone`

do

 echo "plot \"$Nom.datl\" i 0 u xcol:$icol w lp " >> ${Nom}_l.gplot

 echo "replot \"$Nom.datl\" i $i u xcol:$icol w lp " >> ${Nom}_l.gplot

 echo "pause -1" >> ${Nom}_l.gplot

done

 echo "pause -1" >> ${Nom}_l.gplot

chmod u+x  ${Nom}_l.gplot

cat << EOF > ${Nom}_l2.gplot

#!/usr/bin/gnuplot -persist

 set pointsize 2

 xcol=$xcol

EOF

echo "plot \"$Nom.datl\" i 0 u xcol:$icol w lp " >>  ${Nom}_l2.gplot

for i in `seq 1 $ItDone`

do

 echo "replot \"$Nom.datl\" i $i u xcol:$icol w lp " >>  ${Nom}_l2.gplot

done

chmod u+x  ${Nom}_l2.gplot

cat << EOF > ${Nom}_l3.gplot

#!/usr/bin/gnuplot -persist

 set pointsize 2

EOF

 echo "plot \"$Nom.datl\" i 0 u 0:$icol w lp " >>  ${Nom}_l3.gplot

for i in `seq 1 $ItDone`

do

 echo "replot \"$Nom.datl\" i $i u 0:$icol w lp " >>  ${Nom}_l3.gplot

done

chmod u+x  ${Nom}_l3.gplot