Contents 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	#!/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