#!/bin/bash
#
# -- Get the inductance from the energies
#
 PROG=`basename $0`
# echo Start $PROG
 typ_d=_A
 typ_s=all


# ==========   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   
help_doc
 }
# ==========   Help routine end

  if [ $# -gt 0 ]; then

    set -- `getopt s:d:  $*`
    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;;
        --)  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/OUTPOI.TXT
      if ! test -f $fnam ; then
        echo Missing file $fnam
      else
        awk -v cur=$cur '
          $1=="ENERGY" {aind=2*$2/cur/cur; printf"%7.1f  %7.2f\n",cur,aind}
         ' $fnam 
      fi
    fi
  done
  exit