#!/bin/bash
#
# -- Make the grid
#
 PROG=`basename $0`
# echo Start $PROG
 typ_d=_A
 typ_s=all
# sf7_f=sf7_wide_grid_1.txt
 sf7_f=sf7_01.txt

# ==========   Help routine start
 help_dis () {
  cat <<help_doc
 $PROG - a script to produce the field maps on the grid
 Call parameters:
 $PROG [-d directories] [-s type of data]
    where:
           -d appendix for the directory names                 default: $typ_d   
           -s type of sumulation (all -all coils on - normal)  default: $typ_s   
           -f sf7 input file                                   default: $sf7_f   
help_doc
 }
# ==========   Help routine end

  if [ $# -gt 0 ]; then

    set -- `getopt s:d: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;;
        -f)  sf7_f=$2 ; shift 2;;
        --)  shift; file=$*; break;;
       esac
    done
  fi

  if ! test -f $sf7_f ; then
    echo Missing file $sf7_f : try to get it from the upper directory
    if ! test -f ../$sf7_f ; then
       echo Missing file $sf7
       help_dis
       exit 2
    else
      cp -a ../$sf7_f .
      echo File $sf7 is copied from the upper directory
    fi
  fi

  drhome=`pwd`
  sf7_ext=`echo $sf7_f | sed s/sf7//g | sed s/.txt// | sed s/_//`
  
  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
      cp -a $sf7_f $dr1/
      cd $dr1
      echo $dr
      if test -f OUTSF7.TXT; then
        mv OUTSF7.TXT tmpold.txt
      fi
      wine ~/.wine/LANL/SF7.EXE $sf7_f *.T35
      if test -f OUTSF7.TXT; then
        mv OUTSF7.TXT OUTSF7_$sf7_ext.TXT
        if test -f tmpold.txt ; then 
          rm -f tmpold.txt
        fi
      fi
      cd $drhome 
    fi
  done
  exit