#!/bin/bash

pack=$1
ast=`echo $pack | mpc2asteroid | awk '{if ($1>0) print $1; else print $2}' | tr A-Z a-z`
JD=$2
n=$3
st=$4
obs=$5
pV=$6
eta=$7


H=`awk -v ast=$pack '{if ((NF==4)&&($1==ast)&&($2>0)){print $2;found=1};}END{if (found==0)print 100}' ~/Documents/Dropbox/thermops_marco1/eph/hdes`

if [ $H == "100" ] 
then
  H=`awk -v ast=$ast '{if ((NF==4)&&($4==ast)&&($2>0)){print $2;found=1};}END{if (found==0)print 100}' ~/Documents/Dropbox/thermops_marco1/eph/hdes`
fi


D=`echo $H $pV | h2d | awk '{print $3}'`
curl -s "http://vo.imcce.fr/webservices/miriade/ephemcc_query.php?-name=a:$ast&-type=&-ep=$JD&-nbd=$n&-step=$st&-observer=$obs&-mime=text&-tscale=UTC&-theory=INPOP&-teph=1&-tcoor=2" > .eph

#cat .eph
gawk 'BEGIN {FIELDWIDTHS="12 1 12 18 18 18 18 18 8 8 50"} 
     {split($1, date, "-")
      split($3, time, ":")
      if ((date[1]>0) && ($7>0)) printf("%04d-%02d-%02dT%02d:%02d:%02d %8.3f %8.3f %6.1f %6.1f\n",
      date[1],date[2],date[3],time[1],time[2],time[3],$8,$7,$9,$10)}' .eph >.datetime
gawk 'BEGIN {FIELDWIDTHS="12 1 12 18 18 18 18 18 8 8 50"} 
     {split($1, date, "-")
      split($3, time, ":")
      if ((date[1]>0) && ($7>0)) print $8, $7, $9}' .eph >.geom

gawk -v D=$D -v pV=$pV -v eta=$eta '/^ *[0-9]/ {print D,pV,eta,$1,$2,$3}' .geom | tmflux -l 11.6 23 11 | gawk '{print $0}' > .flx
gawk '{printf("W3 %6.3f W4 %6.3f\n", -2.5*log($2/29.0448)/log(10), -2.5*log($4/8.2839)/log(10)) }' .flx  >.mag
#echo H=$H D=$D"(km)" pV=$pV eta=$eta "--" r Delta pha sun_elon lambda1 flux1 lambda2 flux2 ...
echo $ast $H $D $pV $eta > .prop
gawk '{print $1, $2, $3}' .datetime > .geom2
paste .prop .geom2 .mag
#rm .eph .flx .datetime .geom .geom2 .mag .prop
