#! /bin/bash 
# Last edited on 2018-06-17 11:51:28 by jstolfi

tmins=( "-60.00" "-59.30" "-58.50" "-57.75" "-57.00" "-55.50" "-54.00" "-52.00"  )
expos=( "4.93" "5.00" "5.08" "5.11" "5.13" "5.18" "5.25" )

tmins=( "-55.00" )
expos=( "5.00" )

for tmin in "${tmins[@]}" ; do
  for expo in  "${expos[@]}" ; do 
    cat .data.txt \
      | gawk \
          -v tmin="${tmin}" \
          -v expo="${expo}" \
          '//{t=$1;c=$2; tp=exp(expo*log(t-tmin)); printf "%3d %6.2f %10.4f %3.1f\n", n,c,tp,1.0; n=n+1}' \
      > .data_fit.txt
    linear_fit -terms 2 -writeFormula .coefs.txt < .data_fit.txt > .out_fit.txt 2>.err.txt
    printf "=== %8.2f %6.2f : " "${tmin}" "${expo}" 1>&2
    tail -1 .err.txt
  done
done