#!/bin/bash

# chartGaia [-c ra,dec] [-d DFS] [-p PIXEL] file.tsv ... ( or pipe)
#-----------------------------------------------------------------------
CENTER=; DFS="\t"; PIX=1; HELP=;
nfRA=3; nfDec=4;
#-----------------------------------------------------------------------



IF=$1
awk -F"$DFS" -v nfRA=$nfRA -v nfDec=$nfDec -v PIX=$PIX\
	'BEGIN{pi=3.141592653589793}
    	 NR>1{
             ++i; RA[i]=$nfRA; Dec[i]=$nfDec
             cRA+=RA[i];cDec+=Dec[i]
             }
 	END {
 	    n=i; cRA=cRA/n; cDec=cDec/n; 
            for (i=1;i<=n;i++) {

            dRA[i]=(RA[i]-cRA)*3600*cos(Dec[i]*pi/180);
            dDec[i]=(Dec[i]-cDec)*3600;
            theta[i]=sqrt(dRA[i]^2+dDec[i]^2);

            pixdRA[i]=int(dRA[i]/PIX); 
            pixDec[i]=int(dDec[i]/PIX);
        
            (pixX<pixdRA[i])? pixX=pixX:pixX=pixdRA[i];
            (pixY<pixdDec[i])?pixY=pixY:pixY=pixDec[i]
                                 }
            for (i=1;i<=n;i++) {
	    ix[i]=-pixX+pixdRA[i]+1;
            iy[i]=-pixY+pixDec[i]+1;
	    printf ("%0.2f %0.2f %0.1f %d %d %d %d\n",
            dRA[i],dDec[i],theta[i],pixdRA[i],pixDec[i],ix[i],iy[i])
            #print ix[i],iy[i]
            }
            }' $IF
