povconv

   1 #!/bin/sh
   2
   3 #
   4 # convert moldyn output to parcas/rasmol output
   5 # frank.zirkelbach@physik.uni-augsburg.de
   6 #
   7
   8 #if [ ! -d $1 ]; then
   9 #        echo "no valid directory"
  10 #        exit 1
  11 #fi
  12
  13 TRG="mdlat.in.xyz"
  14
  15 rm -f $TRG
  16
  17 ((frame=1))
  18
  19 #for file in $1/atomic_conf_*.xyz; do
  20
  21 file=$1
  22
  23         count=`grep '# \[P\]' $file | awk '{ print $3 }'`
  24         time=`grep '# \[P\]' $file | awk '{ print $4 }'`
  25         x0=`grep '# \[D\]' $file | head -n 1 | awk '{ print $3 '}`
  26         y0=`grep '# \[D\]' $file | head -n 1 | awk '{ print $4 '}`
  27         z0=`grep '# \[D\]' $file | head -n 1 | awk '{ print $5 '}`
  28         x1=`grep '# \[D\]' $file | tail -n 1 | awk '{ print $3 '}`
  29         y1=`grep '# \[D\]' $file | tail -n 1 | awk '{ print $4 '}`
  30         z1=`grep '# \[D\]' $file | tail -n 1 | awk '{ print $5 '}`
  31         dimx=`echo $x0 $x1 | awk '{ print $2-$1 }'`
  32         dimy=`echo $y0 $y1 | awk '{ print $2-$1 }'`
  33         dimz=`echo $z0 $z1 | awk '{ print $2-$1 }'`
  34
  35         echo "   $count" >> $TRG
  36         echo " Frame number $frame $time fs boxsize $dimx $dimy $dimz" >> $TRG
  37
  38         ((ac=1))
  39
  40         cat $file | grep -v '^#' | while read name x y z color temp; do
  41                 # well, i only use Si and C
  42                 [ "$name" = "Si" ] && type=0
  43                 [ "$name" = "C" ] && type=1
  44                 #echo "$name $x $y $z $type ${ac}" >> $TRG
  45                 echo "$name $x $y $z $type" >> $TRG
  46                 ((ac+=1))
  47         done
  48
  49         ((frame+=1))
  50 #done
  51
  52