1 |
#!/bin/sh |
2 |
case $# in |
3 |
0|1) echo "Usage: `basename $0` input_fil outp_sufix [field_name_list]" |
4 |
echo " from ASCII Diagnostics-Stats file (=input_fil) and for each field" |
5 |
echo " (in field_name_list), extract the associated numerical values and" |
6 |
echo " put them in a separated file: stDiag_{VAR}.{outp_sufix}" |
7 |
exit 9 ;; |
8 |
2) listV='T' ;; |
9 |
*) ;; |
10 |
esac |
11 |
|
12 |
inpFil=$1 |
13 |
#outFil=$2 |
14 |
sufx=$2 |
15 |
tmpFil='TTT.'$$ |
16 |
if test $# = '2' |
17 |
then listV='T' |
18 |
else |
19 |
shift; shift; listV=$* |
20 |
fi |
21 |
echo $inpFil $sufx $listV $tmpFil |
22 |
grep '^# Fields ' $inpFil |
23 |
|
24 |
for vv in $listV |
25 |
do |
26 |
outFil='stDiag_'$vv.$sufx |
27 |
case $vv in |
28 |
'Eta') var='ETAN';; |
29 |
'Et2') var='ETANSQ';; |
30 |
'T') var='THETA';; |
31 |
'S') var='SALT' ;; |
32 |
'U') var='UVEL' ;; |
33 |
'V') var='VVEL' ;; |
34 |
'W') var='WVEL' ;; |
35 |
'U2') var='UVELSQ' ;; |
36 |
'V2') var='VVELSQ' ;; |
37 |
*) var=$vv ;; |
38 |
esac |
39 |
#echo 'grep -m 1' "^ field : ${var} " $inpFil |
40 |
#grep -m 1 "^ field : ${var} " $inpFil |
41 |
nLev=`grep -m 1 "^ field : ${var} " $inpFil | sed 's/.*=//'` |
42 |
if test ${nLev}'xx' = 'xx' |
43 |
then |
44 |
echo 'WARNING variable:' $var 'not found in file' $inpFil |
45 |
else |
46 |
echo 'extract variable:' $var ', outFil=' $outFil |
47 |
#echo 'nLev=' $nLev |
48 |
#- echo set of sed instruction to file sedFil : |
49 |
sedFil=$tmpFil.$vv |
50 |
rm -f $sedFil |
51 |
echo "/^ field : $var /{" > $sedFil |
52 |
if test $nLev = 1 ; then k=0 ; else k=-1 ; fi |
53 |
while [ $k -le $nLev ] |
54 |
do |
55 |
echo "N" >> $sedFil |
56 |
k=`expr $k + 1` |
57 |
done |
58 |
echo "p" >> $sedFil |
59 |
echo "}" >> $sedFil |
60 |
#- extract records of variable "var" using sed command & sedFil |
61 |
sed -n -f $sedFil $inpFil > $tmpFil |
62 |
head -1 $tmpFil |
63 |
sed '/^ k /d' $tmpFil | sed '/^ field /d' > $outFil |
64 |
rm -f $sedFil $tmpFil |
65 |
fi |
66 |
done |
67 |
exit |