Annotation of /MITgcm_contrib/jmc_script/extract_StD

#!/bin/sh
case $# in
  0|1) echo "Usage: `basename $0` input_fil outp_sufix [field_name_list]"
       echo "   from ASCII Diagnostics-Stats file (=input_fil) and for each field"
       echo "   (in field_name_list), extract the associated numerical values and"
       echo "   put them in a separated file: stDiag_{VAR}.{outp_sufix}"
       exit 9 ;;
    2) listV='T' ;;
    *) ;;
esac

inpFil=$1
#outFil=$2
sufx=$2
tmpFil='TTT.'$$
if test $# = '2'
then listV='T'
else
 shift; shift; listV=$*
fi
echo $inpFil $sufx $listV $tmpFil
grep '^# Fields ' $inpFil

for vv in $listV
do
 outFil='stDiag_'$vv.$sufx
 case $vv in
   'Eta') var='ETAN';;
   'Et2') var='ETANSQ';;
     'T') var='THETA';;
     'S') var='SALT' ;;
     'U') var='UVEL' ;;
     'V') var='VVEL' ;;
     'W') var='WVEL' ;;
    'U2') var='UVELSQ' ;;
    'V2') var='VVELSQ' ;;
       *)  var=$vv ;;
 esac
#echo 'grep -m 1' "^ field : ${var} " $inpFil
#grep -m 1 "^ field : ${var} " $inpFil
 nLev=`grep -m 1 "^ field : ${var} " $inpFil | sed 's/.*=//'`
 if test ${nLev}'xx' = 'xx'
 then
  echo 'WARNING variable:' $var 'not found in file' $inpFil
 else
  echo 'extract variable:' $var ', outFil=' $outFil
#echo 'nLev=' $nLev
#- echo set of sed instruction to file sedFil :
  sedFil=$tmpFil.$vv
  rm -f $sedFil
  echo "/^ field : $var /{" > $sedFil
  if test $nLev = 1 ; then k=0 ; else k=-1 ; fi
  while [ $k -le $nLev ]
  do 
   echo "N" >> $sedFil
   k=`expr $k + 1`
  done
  echo "p" >> $sedFil
  echo "}" >> $sedFil
#- extract records of variable "var" using sed command & sedFil
  sed -n -f $sedFil $inpFil > $tmpFil
  head -1 $tmpFil
  sed '/^ k /d' $tmpFil | sed '/^ field /d' > $outFil 
  rm -f $sedFil $tmpFil
 fi
done
exit
1	jmc	1.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