1 | #!/bin/sh |
---|
2 | # |
---|
3 | #Script sh d'extraction de fichier SCA des sequances Diodes a Bruit (DAB) |
---|
4 | #et Data des cycles ON/OFF |
---|
5 | #Bash script that extract from SCA files the DAB/DATA phases in ON/OFF cycles |
---|
6 | #It is supposed that the sequence is: DAB1 DATA1 DAB2 DATA2 |
---|
7 | #and that the DAB1 and DAB2 have equal number of subsequences |
---|
8 | # Auteur: JE Campagne (LAL Orsay) |
---|
9 | # Creation: 11 Avril 2011 |
---|
10 | # Revision: 12 Avril 2011 |
---|
11 | # o 12/04/11: mv date format to %Y-%m-%d %T for further treatment by 'date' |
---|
12 | # system function |
---|
13 | |
---|
14 | #Basic tuning of system fuction used. |
---|
15 | AWK=/bin/awk |
---|
16 | SED=/bin/sed |
---|
17 | GREP=/bin/grep |
---|
18 | WC=/usr/bin/wc |
---|
19 | RM=/bin/rm |
---|
20 | CAT=/bin/cat |
---|
21 | TOUCH=/bin/touch |
---|
22 | DATE=/bin/date |
---|
23 | ECHO=/bin/echo |
---|
24 | LS=/bin/ls |
---|
25 | MKDIR=/bin/mkdir |
---|
26 | DefaultIFS=$' \t\n' |
---|
27 | IFS=$DefaultIFS |
---|
28 | |
---|
29 | |
---|
30 | scriptName="`basename $0`" |
---|
31 | $ECHO "Processing script ${scriptName} at `date`" |
---|
32 | |
---|
33 | #temporary files to synchronize scripts |
---|
34 | tmppublicpath=${TMPPUBLICPATH} |
---|
35 | $LS -l ${tmppublicpath} > /dev/null |
---|
36 | |
---|
37 | #Path where the job will do temporary IO |
---|
38 | . ${SCRIPTPATH}/set_iojobpath.sh |
---|
39 | iojobpath=$(set_iojobpath) |
---|
40 | cd ${iojobpath} |
---|
41 | |
---|
42 | |
---|
43 | if [ ! $# = 1 ] |
---|
44 | then |
---|
45 | $ECHO "FATAL (${scriptName}): usage: ./sca.sh <file>" |
---|
46 | exit 1 |
---|
47 | fi |
---|
48 | |
---|
49 | |
---|
50 | fullfile=$1 |
---|
51 | if [ ! -e $fullfile ] |
---|
52 | then |
---|
53 | $ECHO "FATAL (${scriptName}): $fullfile not found: FATAL" |
---|
54 | exit 1 |
---|
55 | fi |
---|
56 | |
---|
57 | #The output file would be in "fullprocfile" |
---|
58 | filename=$(basename $fullfile) |
---|
59 | #extension=${filename##*.} |
---|
60 | filename=${filename%.*} |
---|
61 | extension=".sum" |
---|
62 | fullprocfile=${filename}${extension} |
---|
63 | #create empty output file |
---|
64 | $RM -f $fullprocfile |
---|
65 | $TOUCH $fullprocfile |
---|
66 | |
---|
67 | |
---|
68 | #extract the PERI lines valid for DAB and DATA |
---|
69 | OUT1=./tmp1.$$ |
---|
70 | $AWK ' |
---|
71 | $1=="PERI" { |
---|
72 | print "CYCLE",$6,$2,$3,$7,$8 |
---|
73 | } |
---|
74 | ' $fullfile > $OUT1 |
---|
75 | |
---|
76 | |
---|
77 | #Sequence of the cycle numbers (perform a unique filtering) |
---|
78 | tableau=( `cat $OUT1 | $AWK '{print $2}'` ) |
---|
79 | IFS=' |
---|
80 | ' |
---|
81 | tableau=( $( printf "%s\n" "${tableau[@]}" | awk 'x[$0]++ == 0' ) ) |
---|
82 | IFS=$DefaultIFS |
---|
83 | |
---|
84 | #scan each cycle |
---|
85 | OUT2=./tmp2.$$ |
---|
86 | for i in ${tableau[@]} |
---|
87 | do |
---|
88 | #extract the DAB calibration phases. Hypothesis: 2 identical sequences |
---|
89 | CALIB=./tmp_cycle_${i}_calib.$$ |
---|
90 | $AWK -v cycle=$i '$2==cycle && $5=="-1" {print $0} ' $OUT1 > $CALIB |
---|
91 | nlines=`$WC $CALIB | $AWK '{print $1}'` |
---|
92 | start1=1 |
---|
93 | let "end1=nlines/2" |
---|
94 | let "start2=end1+1" |
---|
95 | end2=$nlines |
---|
96 | |
---|
97 | #date Start/End of DAB1 |
---|
98 | $SED -e "${start1}!d" $CALIB | $AWK -v cycle=$i 'BEGIN {ORS=","} { |
---|
99 | split($4,b,"."); |
---|
100 | print cycle "," $3,b[1]}' >> ${OUT2} |
---|
101 | $SED -e "${end1}!d" $CALIB | $AWK 'BEGIN {ORS=","} { |
---|
102 | split($4,b,"."); |
---|
103 | print $3,b[1]}' >> ${OUT2} |
---|
104 | |
---|
105 | |
---|
106 | #date Start/End of DATA1: use date system function to add the duration to |
---|
107 | #compute the End,cf. End = Start + duration |
---|
108 | $AWK -v datesys=$DATE -v cycle=$i '$2==cycle && ( $5=="1" ) { |
---|
109 | split($3,a,"/"); |
---|
110 | split($4,b,"."); |
---|
111 | dateStart=a[2]"/"a[1]"/"a[3] " " b[1]; |
---|
112 | duree=$6 " sec"; |
---|
113 | cmd=sprintf("%s --date \"%s %s\" \"+%%d/%%m/%%y %%T\" ",datesys,dateStart,duree); |
---|
114 | (cmd | getline dateEnd); |
---|
115 | close(cmd); |
---|
116 | ORS=","; print a[1]"/"a[2]"/"a[3] " " b[1]; |
---|
117 | ORS=","; print dateEnd; |
---|
118 | } ' $OUT1 >> ${OUT2} |
---|
119 | |
---|
120 | #same as DAB1 for DAB2 |
---|
121 | $SED -e "${start2}!d" $CALIB | $AWK 'BEGIN {ORS=","} { |
---|
122 | split($4,b,"."); |
---|
123 | print $3,b[1]}' >> ${OUT2} |
---|
124 | $SED -e "${end2}!d" $CALIB | $AWK 'BEGIN {ORS=","} { |
---|
125 | split($4,b,"."); |
---|
126 | print $3,b[1]}' >> ${OUT2} |
---|
127 | #same as DATA1 for DATA2 |
---|
128 | $AWK -v datesys=$DATE -v cycle=$i '$2==cycle && ( $5=="2" ) { |
---|
129 | split($3,a,"/"); |
---|
130 | split($4,b,".") |
---|
131 | dateStart=a[2]"/"a[1]"/"a[3] " " b[1]; |
---|
132 | duree=$6 " sec"; |
---|
133 | cmd=sprintf("%s --date \"%s %s\" \"+%%d/%%m/%%y %%T\" ",datesys,dateStart,duree); |
---|
134 | (cmd | getline dateEnd); |
---|
135 | close(cmd); |
---|
136 | ORS=","; print a[1]"/"a[2]"/"a[3] " " b[1]; |
---|
137 | ORS="\n"; print dateEnd; |
---|
138 | } ' $OUT1 >> ${OUT2} |
---|
139 | #clean |
---|
140 | $RM $CALIB |
---|
141 | done |
---|
142 | |
---|
143 | #transform all the date/time into %Y-%m-%d %T |
---|
144 | $AWK -v datesys=$DATE ' |
---|
145 | BEGIN { FS=","; OFS="," } |
---|
146 | { |
---|
147 | for (i=2;i<=NF;i++) { |
---|
148 | split($i,a,"/"); |
---|
149 | newdate=a[2]"/"a[1]"/"a[3]; |
---|
150 | cmd=sprintf("%s --date \"%s\" \"+%%Y-%%m-%%d %%T\" ",datesys,newdate); |
---|
151 | (cmd | getline newdate); |
---|
152 | close(cmd); |
---|
153 | $i=newdate; |
---|
154 | } |
---|
155 | print $0; |
---|
156 | }' ${OUT2} > ${fullprocfile} |
---|
157 | |
---|
158 | |
---|
159 | #clean |
---|
160 | $RM $OUT1 |
---|
161 | $RM $OUT2 |
---|
162 | |
---|
163 | exit 0 |
---|
164 | |
---|
165 | |
---|
166 | |
---|