#!/bin/csh
####################################
# theTab = table name
# Argument is P=period
# Vgraph is "gif" "ps" or "txt"
####################################

#################################################################
# Get the Arguments: tab and Period
#################################################################
set Star = "GJ 3473"

if ("$1" =~ P=*) then
    set P = `echo "$1" | cut -c3-` ; shift
    if (("$P" == 0) && ("$1" == "-x" )) then	### REMOVE Limits 0/1
	if (("$2" == "0") && ("$3" == "1" )) then 
	    shift; shift; shift
	endif
    endif
else
    set P = ""
endif
if ("$P" == "") set P = 0 


#################################################################
# Vgraph not defined ==> just the plot introduction
#################################################################
if (! $?Vgraph) then
    if ($P != 0) then
        set plotx = "-x 0 1"
    else
        set plotx = ""
    endif
    set ploty = "" # "-y - - -"       # Reverse Y axis
    # These values are superseded by what's contained in $argv

    echo1 "\section*{RV of \object{$Star} (P=$P)}"
    cgigraph -img -ps:"Postscript Figure" -txt:"Data as a Table" -form -htx \
	P=$P -graph $plotx $ploty $argv:q
    exit 0
endif

################################################################
# Actual gif or ps plot (Vgraph contains gif or plot)
################################################################

set  plotarg = \
  "-f 0.0325 -h 0.75 -w 0.75 -u 0.125 -r 0.15 --title-font-size 0.05 -m -1 -S 3"

if ($P != 0) then
    set X = "Phase"
else
    set X = "BJD-2458000"
endif
set Y = 'V\sbr\eb [km/s]'

#fcat catalog.dat | gawk -v P=$P "/^$Star /"'{\
asu -data <<====EOQ | gawk -v P=$P 'BEGIN{FS="\t"}\
   /^$/{print; next}\
   /---/ {next}\
  {\
    jd = $1; \
    if (P>0) phase = (jd/P) - int(jd/P) ; \
    else phase = jd-2458000; \
    printf("%.7f %s %s\n", phase, $2,$3) \
  }' \
| graph -T $Vgraph -C -I e -L "$Star" --pen-colors 1=red3:2=red \
  $plotarg -X $X:q -Y $Y:q $argv:q
-source=$theCat/$theTab
-out=HJD,HRV,e_HRV
====EOQ
exit 0