Chimie Théorique » OpenPath

#!/bin/bash

if [ $# -lt 4 ]; then

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

 exit

fi

Fout=$1

ItMax=$2

Nom=$3

NGeomF=$4

xcol=1

if [ $# -eq 5 ]; then

xcol=$5

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"

for i in `seq 0 $ItMax`

 do

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

 cat Scan.dat >> $Nom.datl

 echo " " >> $Nom.datl

 echo " " >> $Nom.datl

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 $ItMax`

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 $ItMax`

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 $ItMax`

do

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

done

chmod u+x  ${Nom}_l3.gplot