#bin/bash
ast=`echo $1 | tr A-Z a-z`
JD=$2
n=$3
st=$4
obs=$5
pV=$6
eta=$7
base=$8
#echo $ast
#simpsDpv=`awk -v ast=$ast 'tolower($2) ~ tolower(ast) {print $6,$4}' ~/Documents/DATA/iras/SIMPS.txt`
#H=`awk -v ast=$ast '{if ($3 == ast) print $2}' ~/codemdb/eph/hdes`
H=`awk -v ast=$ast '{if ((NF==3) && ($3 == ast)) print $2; if ((NF==2) && ($2 == ast)) print $1}' ~/dev/thermops/eph/hdes`
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
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

awk -v D=$D -v pV=$pV -v eta=$eta '/^ *[0-9]/ {print D,pV,eta,$1,$2,$3}' .geom | tmflux -l 8 13 2 | \
    awk '{printf("%4.1f %5.2f %4.1f %5.2f %4.1f %5.2f\n",$1,$2,$3,$4,$5,$6)}' > .flx
awk -v D=$D -v base=$base '/^ *[0-9]/ {print D/0.72/$2, base,  9.0}' .geom | theta2vis -l 8 13 2 >.visibil

#awk '{if (NF==15) printf("%10s%02d  %02d %02d %05.2f  %02d %02d %04.1f\t%4.1f\t%5.1f\t%5.1f\n",$1,$2,$3,$4,$5,$6,$7,$8,$10,$11,$12)}' .eph
#awk '/^ *[0-9]/ {print $1$2}' .eph > .eph2
echo H=$H D=$D"(km)" pV=$pV eta=$eta "--" r Delta pha sun_elon lambda1 flux1 lambda2 flux2 ...
awk '{print $1, $2, $3}' .datetime > .geom2
paste .geom2 .flx .visibil
rm .eph .flx .datetime .geom .geom2 .visibil
