Revision 11 examples/Test/Dep_tscan_m

Dep_tscan_m (revision 11)
1 1
#!/bin/bash
2 2
export LANG=C
3
if [ $# -lt 3 ]; then
3
if [ $# -lt 1 ]; then
4 4
 echo "Use: $0 File.out [PES] [T]"
5 5
 exit
6 6
fi
......
25 25
 set pointsize 2
26 26
 unset key
27 27
EOF
28
let ItM=ItMax-1
28
let ItM=ItMax-2
29 29
for i in `seq 0 $ItM`
30 30
  do
31 31
 xyz2scan ${Nom}_cart.${i} > /dev/null
......
33 33
 cat <<EOF >>  $NomGplot
34 34
 set multiplot
35 35
 set origin 0.,0.
36
 set size 0.5
36
 set size 0.5,1.
37 37
 clear
38 38
 set xrange [1.:2.5]
39 39
 set yrange [-10:190]
40
 plot "HCN_irc.dat" u 1:3 w l
40
 plot "HCN_irc.dat" u 1:3 w l ls 1
41 41
EOF
42 42
 if [ "$PES" != "" ]; then
43
  echo "replot \"HCN_PES_2D.dat\" u 1:2 w l 2, \"HCN_PES_2D.dat\" u 1:(-\$2) w l 2,\"HCN_PES_2D.dat\" u 1:(360-\$2) w l 2" >> $NomGplot
43
  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 44
 fi
45 45
 if [ -s ${Nom}_spline.${i} ]; then
46
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
47
 echo "replot \"${Nom}_spline.${i}\" u 2:4 w l" >> $NomGplot
46
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w p ls 5" >> $NomGplot
47
 echo "replot \"${Nom}_spline.${i}\" u 2:4 w l ls 5" >> $NomGplot
48 48
 elif [  -s ${Nom}_dbgtgt_${i}.dat ]; then
49
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
50
 echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 2:(\$4*180./pi) w l" >> $NomGplot
49
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp 5" >> $NomGplot
50
 echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 2:(\$4*180./pi) w l ls 4" >> $NomGplot
51 51
 else
52
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp 5" >> $NomGplot
52
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp ls 5" >> $NomGplot
53 53
fi
54 54
 if [ "$Tgt" != "" ]; then
55
  echo "replot \"${Nom}_tgt_${i}.dat\" u 1:3 w l" >> $NomGplot
55
  echo "replot \"${Nom}_tgt_${i}.dat\" u 1:3 w l ls 4" >> $NomGplot
56 56
 fi
57 57
 cat <<EOF >>  $NomGplot
58 58
 set origin .505,0.
59
# set size 0.5
59
# set size 0.5,1.
60 60
 clear
61 61
 set yrange $yr
62 62
 set xrange [-10:190]
63 63
EOF
64
 echo "plot \"${Nom}_${i}.dat\" u 3:4 t \"E\" w lp" >> $NomGplot
64
 echo "plot \"${Nom}_${i}.dat\" u 3:4 t \"E\" w lp ls 1" >> $NomGplot
65 65
 echo "unset multiplot" >> $NomGplot
66 66
 echo "pause -1" >> $NomGplot
67 67
 done
68
let i=ItMax
68
let i=ItMax-1
69 69
 xyz2scan ${Nom}_cart.${i} > /dev/null
70 70
 mv Scan.dat ${Nom}_${i}.dat
71 71
 cat <<EOF >>  $NomGplot
72 72
 set multiplot
73 73
 set origin 0.,0.
74
# set size 0.5
74
# set size 0.5,1.
75 75
 clear
76 76
 set xrange [1.:2.5]
77 77
 set yrange [-10:190]
78
 plot "HCN_irc.dat" u 1:3 w l
78
 plot "HCN_irc.dat" u 1:3 t "IRC" w l ls 1
79 79
EOF
80 80
 if [ "$PES" != "" ]; then
81
  echo "replot \"HCN_PES_2D.dat\" u 1:2 w l 2, \"HCN_PES_2D.dat\" u 1:(-\$2) w l 2,\"HCN_PES_2D.dat\" u 1:(360-\$2) w l 2" >> $NomGplot
81
  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 82
 fi
83 83
 if [ -s ${Nom}_spline.${i} ]; then
84
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
85
 echo "replot \"${Nom}_spline.${i}\" u 2:4 w l" >> $NomGplot
84
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w p ls 5" >> $NomGplot
85
 echo "replot \"${Nom}_spline.${i}\" u 2:4 w l ls 5" >> $NomGplot
86 86
 elif [  -s ${Nom}_dbgtgt_${i}.dat ]; then
87
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
88
 echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 2:(\$4*180./pi) w l" >> $NomGplot
87
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp ls 5" >> $NomGplot
88
 echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 2:(\$4*180./pi) w l ls 4" >> $NomGplot
89 89
 else
90
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp 5" >> $NomGplot
90
 echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp ls 5" >> $NomGplot
91 91
fi
92 92
 cat <<EOF >>  $NomGplot
93 93
 set origin .505,0.
94
# set size  0.5
94
# set size  0.5,1.
95 95
 clear
96 96
 set yrange $yr
97 97
 set xrange [-10:190]
98 98
EOF
99
 echo "plot \"${Nom}_${i}.dat\" u 3:4 t \"E\" w lp" >> $NomGplot
99
 echo "plot \"${Nom}_${i}.dat\" u 3:4 t \"E\" w lp ls 1" >> $NomGplot
100 100
 echo "unset multiplot" >> $NomGplot
101 101
 echo "pause -1" >> $NomGplot
102 102

  
103 103
chmod u+x ./$NomGplot
104
echo "./$NomGplot      to see the path"
104
echo "./$NomGplot      to see the path"

Also available in: Unified diff