Révision 4 utils/AnaPath

AnaPath (revision 4)
1 1
#!/bin/bash
2
if [ $# -lt 4 ]; then
3
 echo "Use: $0 File.out MaxCyc PathName NGeomF [x column]"
2
if [ $# -lt 1 ]; then
3
 echo "Use: $0 File.out [x column]"
4 4
 exit
5 5
fi
6 6

  
7 7
Fout=$1
8
ItMax=$2
9
Nom=$3
10
NGeomF=$4
11 8

  
9
ItMax=`grep MAXCYC $Fout | tail -1 | awk '{print $3}'`
10
Nom=`grep PATHNAME $Fout | tail -1 | awk '{print $3}'`
11
NGeomF=`grep NGEOMF $Fout | tail -1 | awk '{print $3}'`
12

  
12 13
xcol=1
13
if [ $# -eq 5 ]; then
14
xcol=$5
14
if [ $# -eq 2 ]; then
15
xcol=$2
15 16
fi
16 17

  
17 18
export LANG=C
......
41 42
fi
42 43
echo "Creating $Nom.datl"
43 44

  
45
ItDone=-1
44 46
for i in `seq 0 $ItMax`
45 47
 do
46
  xyz2path  ${Nom}${Ext}.$i 
47
 cat Scan.dat >> $Nom.datl
48
 echo " " >> $Nom.datl
49
 echo " " >> $Nom.datl
48
 if [ -s ${Nom}${Ext}.$i ]; then
49
     xyz2path  ${Nom}${Ext}.$i 
50
     cat Scan.dat >> $Nom.datl
51
     echo " " >> $Nom.datl
52
     echo " " >> $Nom.datl
53
     let ItDone=ItDone+1
54
 fi
50 55
done 
51 56

  
52 57
icol=`tail -1 Scan.dat | wc -w`
......
58 63
EOF
59 64

  
60 65

  
61
for i in `seq 1 $ItMax`
66
for i in `seq 1 $ItDone`
62 67
do
63 68
 echo "plot \"$Nom.datl\" i 0 u xcol:$icol w lp " >> ${Nom}_l.gplot
64 69
 echo "replot \"$Nom.datl\" i $i u xcol:$icol w lp " >> ${Nom}_l.gplot
......
76 81
EOF
77 82

  
78 83
echo "plot \"$Nom.datl\" i 0 u xcol:$icol w lp " >>  ${Nom}_l2.gplot
79
for i in `seq 1 $ItMax`
84
for i in `seq 1 $ItDone`
80 85
do
81 86
 echo "replot \"$Nom.datl\" i $i u xcol:$icol w lp " >>  ${Nom}_l2.gplot
82 87
done
......
89 94
EOF
90 95

  
91 96
 echo "plot \"$Nom.datl\" i 0 u 0:$icol w lp " >>  ${Nom}_l3.gplot
92
for i in `seq 1 $ItMax`
97
for i in `seq 1 $ItDone`
93 98
do
94 99
 echo "replot \"$Nom.datl\" i $i u 0:$icol w lp " >>  ${Nom}_l3.gplot
95 100
done

Formats disponibles : Unified diff