root / examples / Test / Dep_tscan_m @ 11
Historique | Voir | Annoter | Télécharger (2,84 ko)
| 1 | 1 | pfleura2 | #!/bin/bash |
|---|---|---|---|
| 2 | 1 | pfleura2 | export LANG=C |
| 3 | 11 | pfleura2 | if [ $# -lt 1 ]; then |
| 4 | 1 | pfleura2 | echo "Use: $0 File.out [PES] [T]" |
| 5 | 1 | pfleura2 | exit |
| 6 | 1 | pfleura2 | fi |
| 7 | 1 | pfleura2 | Fout=$1 |
| 8 | 1 | pfleura2 | ItMax=`grep MAXCYC $Fout | tail -1 | awk '{print $3}'`
|
| 9 | 1 | pfleura2 | Nom=`grep PATHNAME $Fout | tail -1 | awk '{print $3}'`
|
| 10 | 1 | pfleura2 | PES=$2 |
| 11 | 1 | pfleura2 | Tgt=$3 |
| 12 | 1 | pfleura2 | |
| 13 | 1 | pfleura2 | WDir=`dirname $0` |
| 14 | 1 | pfleura2 | |
| 15 | 1 | pfleura2 | NomGplot=${Nom}_m.gplot
|
| 16 | 1 | pfleura2 | yr="[-0.50:-0.42]" |
| 17 | 1 | pfleura2 | |
| 18 | 1 | pfleura2 | if [ ! -s list ]; then |
| 19 | 1 | pfleura2 | echo "File list is missing" |
| 20 | 1 | pfleura2 | exit |
| 21 | 1 | pfleura2 | fi |
| 22 | 1 | pfleura2 | awk -f $WDir/Dep_tgt.awk $Fout |
| 23 | 1 | pfleura2 | cat <<EOF > $NomGplot |
| 24 | 1 | pfleura2 | #!/usr/bin/gnuplot -persist |
| 25 | 1 | pfleura2 | set pointsize 2 |
| 26 | 1 | pfleura2 | unset key |
| 27 | 1 | pfleura2 | EOF |
| 28 | 11 | pfleura2 | let ItM=ItMax-2 |
| 29 | 1 | pfleura2 | for i in `seq 0 $ItM` |
| 30 | 1 | pfleura2 | do |
| 31 | 1 | pfleura2 | xyz2scan ${Nom}_cart.${i} > /dev/null
|
| 32 | 1 | pfleura2 | mv Scan.dat ${Nom}_${i}.dat
|
| 33 | 1 | pfleura2 | cat <<EOF >> $NomGplot |
| 34 | 1 | pfleura2 | set multiplot |
| 35 | 1 | pfleura2 | set origin 0.,0. |
| 36 | 11 | pfleura2 | set size 0.5,1. |
| 37 | 1 | pfleura2 | clear |
| 38 | 1 | pfleura2 | set xrange [1.:2.5] |
| 39 | 1 | pfleura2 | set yrange [-10:190] |
| 40 | 11 | pfleura2 | plot "HCN_irc.dat" u 1:3 w l ls 1 |
| 41 | 1 | pfleura2 | EOF |
| 42 | 1 | pfleura2 | if [ "$PES" != "" ]; then |
| 43 | 11 | pfleura2 | echo "replot \"HCN_PES_2D.dat\" u 1:2 t \"\" w l ls 2, \"HCN_PES_2D.dat\" u 1:(-\$2) t \"\" w l ls 2,\"HCN_PES_2D.dat\" u 1:(360-\$2) t \"\" w l ls 2" >> $NomGplot |
| 44 | 1 | pfleura2 | fi |
| 45 | 1 | pfleura2 | if [ -s ${Nom}_spline.${i} ]; then
|
| 46 | 11 | pfleura2 | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p ls 5" >> $NomGplot
|
| 47 | 11 | pfleura2 | echo "replot \"${Nom}_spline.${i}\" u 2:4 w l ls 5" >> $NomGplot
|
| 48 | 1 | pfleura2 | elif [ -s ${Nom}_dbgtgt_${i}.dat ]; then
|
| 49 | 11 | pfleura2 | echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp 5" >> $NomGplot
|
| 50 | 11 | pfleura2 | echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 2:(\$4*180./pi) w l ls 4" >> $NomGplot
|
| 51 | 1 | pfleura2 | else |
| 52 | 11 | pfleura2 | echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp ls 5" >> $NomGplot
|
| 53 | 1 | pfleura2 | fi |
| 54 | 1 | pfleura2 | if [ "$Tgt" != "" ]; then |
| 55 | 11 | pfleura2 | echo "replot \"${Nom}_tgt_${i}.dat\" u 1:3 w l ls 4" >> $NomGplot
|
| 56 | 1 | pfleura2 | fi |
| 57 | 1 | pfleura2 | cat <<EOF >> $NomGplot |
| 58 | 1 | pfleura2 | set origin .505,0. |
| 59 | 11 | pfleura2 | # set size 0.5,1. |
| 60 | 1 | pfleura2 | clear |
| 61 | 1 | pfleura2 | set yrange $yr |
| 62 | 1 | pfleura2 | set xrange [-10:190] |
| 63 | 1 | pfleura2 | EOF |
| 64 | 11 | pfleura2 | echo "plot \"${Nom}_${i}.dat\" u 3:4 t \"E\" w lp ls 1" >> $NomGplot
|
| 65 | 1 | pfleura2 | echo "unset multiplot" >> $NomGplot |
| 66 | 1 | pfleura2 | echo "pause -1" >> $NomGplot |
| 67 | 1 | pfleura2 | done |
| 68 | 11 | pfleura2 | let i=ItMax-1 |
| 69 | 1 | pfleura2 | xyz2scan ${Nom}_cart.${i} > /dev/null
|
| 70 | 1 | pfleura2 | mv Scan.dat ${Nom}_${i}.dat
|
| 71 | 1 | pfleura2 | cat <<EOF >> $NomGplot |
| 72 | 1 | pfleura2 | set multiplot |
| 73 | 1 | pfleura2 | set origin 0.,0. |
| 74 | 11 | pfleura2 | # set size 0.5,1. |
| 75 | 1 | pfleura2 | clear |
| 76 | 1 | pfleura2 | set xrange [1.:2.5] |
| 77 | 1 | pfleura2 | set yrange [-10:190] |
| 78 | 11 | pfleura2 | plot "HCN_irc.dat" u 1:3 t "IRC" w l ls 1 |
| 79 | 1 | pfleura2 | EOF |
| 80 | 1 | pfleura2 | if [ "$PES" != "" ]; then |
| 81 | 11 | pfleura2 | echo "replot \"HCN_PES_2D.dat\" u 1:2 t \"\" w l ls 2, \"HCN_PES_2D.dat\" u 1:(-\$2) t \"\" w l ls 2,\"HCN_PES_2D.dat\" u 1:(360-\$2) t \"\" w l ls 2" >> $NomGplot |
| 82 | 1 | pfleura2 | fi |
| 83 | 1 | pfleura2 | if [ -s ${Nom}_spline.${i} ]; then
|
| 84 | 11 | pfleura2 | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p ls 5" >> $NomGplot
|
| 85 | 11 | pfleura2 | echo "replot \"${Nom}_spline.${i}\" u 2:4 w l ls 5" >> $NomGplot
|
| 86 | 1 | pfleura2 | elif [ -s ${Nom}_dbgtgt_${i}.dat ]; then
|
| 87 | 11 | pfleura2 | echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp ls 5" >> $NomGplot
|
| 88 | 11 | pfleura2 | echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 2:(\$4*180./pi) w l ls 4" >> $NomGplot
|
| 89 | 1 | pfleura2 | else |
| 90 | 11 | pfleura2 | echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp ls 5" >> $NomGplot
|
| 91 | 1 | pfleura2 | fi |
| 92 | 1 | pfleura2 | cat <<EOF >> $NomGplot |
| 93 | 1 | pfleura2 | set origin .505,0. |
| 94 | 11 | pfleura2 | # set size 0.5,1. |
| 95 | 1 | pfleura2 | clear |
| 96 | 1 | pfleura2 | set yrange $yr |
| 97 | 1 | pfleura2 | set xrange [-10:190] |
| 98 | 1 | pfleura2 | EOF |
| 99 | 11 | pfleura2 | echo "plot \"${Nom}_${i}.dat\" u 3:4 t \"E\" w lp ls 1" >> $NomGplot
|
| 100 | 1 | pfleura2 | echo "unset multiplot" >> $NomGplot |
| 101 | 1 | pfleura2 | echo "pause -1" >> $NomGplot |
| 102 | 1 | pfleura2 | |
| 103 | 1 | pfleura2 | chmod u+x ./$NomGplot |
| 104 | 11 | pfleura2 | echo "./$NomGplot to see the path" |