#!/bin/sh
######################################################
#  From FITS, create a correct table
######################################################

rm -f sp.tmp ; touch sp.tmp
for k in o p ; do
    f=sp/NH2D$k*
    #echo2 "#---Convert $f"
    fits2a -d, -c SCAN,OBJECT,CRVAL2,CRVAL3,DATE-OBS,UT,IMAGFREQ,VELOCITY $f \
    | awk -F, -v k=$k '{ date=substr($5,1,10); t=int($6+0.49);
       h=int(t/3600); t -= h*3600; m=int(t/60); s=t-m*60;
       printf "%s %3d|%s|%08.4f %+08.4f|%sT%02d:%02d:%02d|%8.4f %5.1f\n",
               k, $1,$2,$3,     $4,   date,  h,  m,   s, $7/1.e9,$8/1000 }' \
    >> sp.tmp
done

echo2 "#===Merge spectra"
cat sp.tmp | awk '
  /^o/{ o++; r=substr($0,2); rec[o]=r; orec[r]=o; next }
  /^p/{ p++; r=substr($0,2);           prec[r]=p; next }
    END{ for(i=1; i<=o; i++) printf "%s o%3d p%3d\n", rec[i], i, prec[rec[i]]}'\
> spectra.dat