#bin/bash
ast=`echo $1 | tr A-Z a-z`
JD=$2
n=$3
st=$4
obs=$5
H=$6
pV=$7
eta=$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 ((NF==4) && ($4 == ast)) print $2; if ((NF==3) && ($3 == ast)) print $1}' ~/Documents/thermops_marco1/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

#from Pierre Vernazza SOPHIA Oct 2, 2012 email
#awk -v D=$D -v pV=$pV -v eta=$eta '/^ *[0-9]/ {print D,pV,eta,$1,$2,$3}' .geom | tmflux -l 6.6 20 4.5 | awk '{print $0}' > .flx
#awk -v D=$D -v pV=$pV -v eta=$eta '/^ *[0-9]/ {print D,pV,eta,$1,$2,$3}' .geom | tmflux -l 19.7 40 15.1 | awk '{print $0}' > .flx
awk -v D=$D -v pV=$pV -v eta=$eta '/^ *[0-9]/ {print D,pV,eta,$1,$2,$3}' .geom | tmflux -l 6. 20. 4 | awk '{print $0}' > .flx
#awk -v D=$D -v pV=$pV -v eta=$eta '/^ *[0-9]/ {print D,pV,eta,$1,$2,$3}' .geom | tmflux -l 3 13 2 -V 1.5 | awk '{print $0}' > .flx
#awk -v D=$D -v pV=$pV -v eta=$eta '/^ *[0-9]/ {print D,pV,eta,$1,$2,$3}' .geom | tmflux -l 19.7 35 15.1 | awk '{print $0}' > .flx
#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 -v ast=$ast '{print ast, $1, $2, $3}' .datetime > .geom2
paste .geom2 .flx
rm .eph .flx .datetime .geom .geom2
