root / examples / Test / Zmat / Dep_tscan_p @ 1
Historique | Voir | Annoter | Télécharger (3,13 ko)
| 1 | 1 | equemene | #!/bin/bash |
|---|---|---|---|
| 2 | 1 | equemene | export LANG=C |
| 3 | 1 | equemene | if [ $# -lt 3 ]; then |
| 4 | 1 | equemene | echo "Use: $0 File.out MaxCyc PathName [PES] [T]" |
| 5 | 1 | equemene | exit |
| 6 | 1 | equemene | fi |
| 7 | 1 | equemene | Fout=$1 |
| 8 | 1 | equemene | ItMax=$2 |
| 9 | 1 | equemene | Nom=$3 |
| 10 | 1 | equemene | PES=$4 |
| 11 | 1 | equemene | Tgt=$5 |
| 12 | 1 | equemene | |
| 13 | 1 | equemene | NomGplot=${Nom}_p.gplot
|
| 14 | 1 | equemene | if [ ! -s list ]; then |
| 15 | 1 | equemene | echo "File list is missing" |
| 16 | 1 | equemene | exit |
| 17 | 1 | equemene | fi |
| 18 | 1 | equemene | |
| 19 | 1 | equemene | if [ -s ${Nom}_cart.0 ]; then
|
| 20 | 1 | equemene | Ext=_cart |
| 21 | 1 | equemene | elif [ -s ${Nom}.0 ]; then
|
| 22 | 1 | equemene | Ext="" |
| 23 | 1 | equemene | else |
| 24 | 1 | equemene | echo "Cannot find ${Nom}.0 nor ${Nom}_cart.0 -- ERROR"
|
| 25 | 1 | equemene | exit |
| 26 | 1 | equemene | fi |
| 27 | 1 | equemene | |
| 28 | 1 | equemene | awk -f Dep_tgt.awk $Fout |
| 29 | 1 | equemene | cat <<EOF > $NomGplot |
| 30 | 1 | equemene | #!/usr/bin/gnuplot -persist |
| 31 | 1 | equemene | set term x11 |
| 32 | 1 | equemene | set pointsize 2 |
| 33 | 1 | equemene | set origin 0.,0. |
| 34 | 1 | equemene | set size square 1 |
| 35 | 1 | equemene | EOF |
| 36 | 1 | equemene | let ItM=ItMax-1 |
| 37 | 1 | equemene | for i in `seq 0 $ItM` |
| 38 | 1 | equemene | do |
| 39 | 1 | equemene | xyz2scan ${Nom}${Ext}.${i} > /dev/null
|
| 40 | 1 | equemene | mv Scan.dat ${Nom}_${i}.dat
|
| 41 | 1 | equemene | if [ "$Ext" = "" ]; then |
| 42 | 1 | equemene | # We are analysing a calculation done in cartesian. |
| 43 | 1 | equemene | # We have to analyse the spline file and maybe also the tangent file. |
| 44 | 1 | equemene | if [ -s ${Nom}_spline.${i} ]; then
|
| 45 | 1 | equemene | xyz2scan ${Nom}_spline.${i} > /dev/null
|
| 46 | 1 | equemene | mv Scan.dat ${Nom}_spline_${i}.dat
|
| 47 | 1 | equemene | fi |
| 48 | 1 | equemene | if [ -s ${Nom}_dbgtgt.${i} ]; then
|
| 49 | 1 | equemene | xyz2scan ${Nom}_dbgtgt.${i} > /dev/null
|
| 50 | 1 | equemene | mv Scan.dat ${Nom}_dbgtgt_${i}.dat
|
| 51 | 1 | equemene | fi |
| 52 | 1 | equemene | fi |
| 53 | 1 | equemene | |
| 54 | 1 | equemene | cat <<EOF >> $NomGplot |
| 55 | 1 | equemene | clear |
| 56 | 1 | equemene | set xrange [1.:2.5] |
| 57 | 1 | equemene | set yrange [-10:190] |
| 58 | 1 | equemene | plot "HCN_irc.dat" u 1:3 w l |
| 59 | 1 | equemene | EOF |
| 60 | 1 | equemene | if [ "$PES" != "" ]; then |
| 61 | 1 | equemene | 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 |
| 62 | 1 | equemene | fi |
| 63 | 1 | equemene | if [ "$Ext" = "_cart" ]; then |
| 64 | 1 | equemene | if [ -s ${Nom}_spline.${i} ]; then
|
| 65 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
|
| 66 | 1 | equemene | echo "replot \"${Nom}_spline.${i}\" u 2:4 w l" >> $NomGplot
|
| 67 | 1 | equemene | elif [ -s ${Nom}_dbgtgt.${i} ]; then
|
| 68 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
|
| 69 | 1 | equemene | echo "replot \"${Nom}_dbgtgt.${i}\" u 2:(\$4*180./pi) w l" >> $NomGplot
|
| 70 | 1 | equemene | else |
| 71 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp 5" >> $NomGplot
|
| 72 | 1 | equemene | fi |
| 73 | 1 | equemene | else |
| 74 | 1 | equemene | if [ -s ${Nom}_spline_${i}.dat ]; then
|
| 75 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
|
| 76 | 1 | equemene | echo "replot \"${Nom}_spline_${i}.dat\" u 1:3 w lp" >> $NomGplot
|
| 77 | 1 | equemene | elif [ -s ${Nom}_dbgtgt_${i}.dat ]; then
|
| 78 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
|
| 79 | 1 | equemene | echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 1:3 w l" >> $NomGplot
|
| 80 | 1 | equemene | else |
| 81 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp 5" >> $NomGplot
|
| 82 | 1 | equemene | fi |
| 83 | 1 | equemene | fi |
| 84 | 1 | equemene | if [ "$Tgt" != "" ]; then |
| 85 | 1 | equemene | echo "replot \"${Nom}_tgt_${i}.dat\" u 1:3 w l" >> $NomGplot
|
| 86 | 1 | equemene | fi |
| 87 | 1 | equemene | echo "pause -1" >> $NomGplot |
| 88 | 1 | equemene | done |
| 89 | 1 | equemene | let i=ItMax |
| 90 | 1 | equemene | xyz2scan ${Nom}${Ext}.${i} > /dev/null
|
| 91 | 1 | equemene | mv Scan.dat ${Nom}_${i}.dat
|
| 92 | 1 | equemene | cat <<EOF >> $NomGplot |
| 93 | 1 | equemene | clear |
| 94 | 1 | equemene | set xrange [1.:2.5] |
| 95 | 1 | equemene | set yrange [-10:190] |
| 96 | 1 | equemene | plot "HCN_irc.dat" u 1:3 w l |
| 97 | 1 | equemene | EOF |
| 98 | 1 | equemene | if [ "$PES" != "" ]; then |
| 99 | 1 | equemene | 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 |
| 100 | 1 | equemene | fi |
| 101 | 1 | equemene | if [ -s ${Nom}_spline.${i} ]; then
|
| 102 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
|
| 103 | 1 | equemene | echo "replot \"${Nom}_spline.${i}\" u 2:4 w l" >> $NomGplot
|
| 104 | 1 | equemene | elif [ -s ${Nom}_dbgtgt_${i}.dat ]; then
|
| 105 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w p 5" >> $NomGplot
|
| 106 | 1 | equemene | echo "replot \"${Nom}_dbgtgt_${i}.dat\" u 2:(\$4*180./pi) w l" >> $NomGplot
|
| 107 | 1 | equemene | else |
| 108 | 1 | equemene | echo "replot \"${Nom}_${i}.dat\" u 1:3 w lp 5" >> $NomGplot
|
| 109 | 1 | equemene | fi |
| 110 | 1 | equemene | echo "pause -1" >> $NomGplot |
| 111 | 1 | equemene | chmod u+x $NomGplot |
| 112 | 1 | equemene | |
| 113 | 1 | equemene | echo "./$NomGplot to see the path" |