#!/bin/bash
#
# -- Run the POISSON for several different currents
#
 PROG=`basename $0`
# echo Start $PROG
 cur1=50
 cur2=1600
 dcur=50
 drel=0
 d_exe=LANL

#echo GGG=$*
# ==========   Help routine start
 help_dis () {
  cat <<help_doc
 $PROG - a script to run POISSON for various currrents
 Call parameters:
 $PROG [-s current1] [-e current2] [-d delta] [-r relative_increment]
    where:
           -s start current (integer)                                   default: $cur1    
           -e end   current (integer)                                   default: $cur2    
           -d delta - current increment (integer)                       default: $dcur
           -r relative small current increment to calculate derivatives default: $drel
              =0 - skip the small increment calculations
           -x directory LANL for the executables                        default: $d_exe    
help_doc
 }
# ==========   Help routine end

  if [ $# -gt 0 ]; then

    set -- `getopt s:e:d:r:  $*`
    if [ $? != 0 ] ; then
      echo *** Error in getopt . Usage:
      help_dis
      exit 2
    fi
    for i in $* ; do
       case $i in
        -s)  cur1=$2 ; shift 2;;
        -e)  cur2=$2 ; shift 2;;
        -d)  dcur=$2 ; shift 2;;
        -r)  drel=$2 ; shift 2;;
        --)  shift; file=$*; break;;
       esac
    done
  fi
  
  nf=`ls -1 *.am | wc -l`
  if [ $nf -ne 1 ]; then
    echo " *** Error: expect 1 file *.am"
  fi

  fnam1=`ls -1 *.am`
  nam1=`echo $fnam1 | cut -f1 -d\.`
  echo $fnam1 " " $nam1
  cur0=`cat $fnam1 | awk 'FNR==1{cur=$NF;gsub("A","",cur);printf"%s",cur}'`
  echo cur=$cur0 
  nam2=MAG
  fnam2=$nam2".am"

  cur=`expr $cur1 - $dcur`
  while let $((cur+=$dcur));[ $cur -le $cur2 ]; do
    echo $cur
    curn=`echo $cur | awk '{printf"%4.4d",$1}'`
    cura=`echo $cur $drel | awk '{a=int($1*$2*10); printf"%d",a}'`
    vcalc='A'
    [ $cura -gt 0 ] && vcalc='A AD'
#    [ $cura -gt 0 ] && vcalc='AD'
    for vc in $vcalc ; do
      curr=$cur
      if [ $vc == "AD" ]; then
        curr=`echo $cur $cura | awk '{a=$1+$2/10.; printf"%5.1f",a}'`
      fi  
      dr=$curn"_"$vc
      if ! test -d $dr ; then
        mkdir $dr
      fi
      fnam2f=$dr/$fnam2
      pwd
      echo "Start --- cur=" $curr $fnam1 $fnam2f
      ls -alF $fnam1
      echo Current=$curr
      awk -v cur="$curr" \
       ' FNR==1 {cur0=$NF;gsub("A","",cur0);
                 fact=cur/cur0;
                 lin=$0; gsub($NF,cur"A",lin);
                 print lin;
                 print ";  Current used ",cur;
               }
         FNR>1  {
                 if($1~/xjfact/ && cur0 != cur ){
                   print "xjfact="fact"         ; Current scaling factor for "cur"A";
                 }else{
                   print;
                 }
               }' $fnam1 > $fnam2f
      cd $dr
      dr1=all
      if ! test -d $dr1 ; then
        mkdir $dr1
      fi
      cd $dr1
      cp -a ../$fnam2 .
      wine ~/.wine/$d_exe/AUTOMESH.EXE $nam2
      wine ~/.wine/$d_exe/POISSON.EXE $nam2
      cd ../../
    done
  done
  exit