#!/bin/bash
#
# -- Get the field integrals
#
 PROG=`basename $0`
# echo Start $PROG
 typ_d=_A
 typ_s=all
 rr=0
 z1=0
 z2=350
 fld=4


# ==========   Help routine start
 help_dis () {
  cat <<help_doc
 $PROG - a script to analyse payroll tables
 Call parameters:
 $PROG [-d directories] [-s type of data]
    where:
           -d appendix for the directoy names                  default: $typ_d   
           -s type of sumulation (all -all coils on - normal)  default: $typ_s   
           -r radius                                           default: $rr      
           -z z1                                               default: $z1      
           -Z z2                                               default: $z2      
           -F field (column in the sf7 output) =3 Br, 4 - Bz   default: $fld
              5 -full field
help_doc
 }
# ==========   Help routine end

  if [ $# -gt 0 ]; then

    set -- `getopt s:d:r:z:Z:F:  $*`
    if [ $? != 0 ] ; then
      echo *** Error in getopt . Usage:
      help_dis
      exit 2
    fi
    for i in $* ; do
       case $i in
        -s)  typ_s=$2 ; shift 2;;
        -d)  typ_d=$2 ; shift 2;;
        -r)  rr=$2 ; shift 2;;
        -z)  z1=$2 ; shift 2;;
        -Z)  z2=$2 ; shift 2;;
        -F)  fld=$2 ; shift 2;;
        --)  shift; file=$*; break;;
       esac
    done
  fi
  
  for dr in `ls -al -d *$typ_d | awk 'index($1,"d")==1{print $NF}'`; do
    cur=`echo $dr | cut -f1 -d_` 
    dr1=$dr/$typ_s
    if ! test -d $dr1 ; then
      echo Missing directory $dr1
    else
      fnam=$dr1/OUTSF7.TXT
      if ! test -f $fnam ; then
        echo Missing file $fnam
      else
#        echo $fnam
        awk -v cur=$cur -v rr=$rr -v z1=$z1 -v z2=$z2 -v fld=$fld '
          BEGIN {np=0; ok=0}
          $1=="(cm)" {ok=1} 
          ok>0 {if($1==rr && $2>z1 && $2<z2){np++;s+=$fld}}
          END{printf"%7.1f  %7.1f %7.1f\n",cur,s/np,z2-z1}
         ' $fnam 
      fi
    fi
  done
  exit