The ubcone routine can be run many different ways and can produce large and complicated output files. To simplify things, I have written a short bash shell scriptthat runs ubcone. A python script (written in my early days of python-hood) takes this file and writes the positions, magnitudes, and star-galaxy (SG) estimates to a simple ASCII file of ploating point data. I can then build a variety of code (fortran, python, shell) that handle these files as input.
% ubc01 11.2337736 +44.3347747 0.2 19.5
% usno_00.py UBC01_usnob.asc
Arguments for ubc01:
arg1 - Ra center (floating point in HOURS)
arg2 - Dec center (floating point in DEGREES)
arg3 - search radius in degrees (0.2)
arg4 - R2 faint magnitude limit.
Output of ubs01 = UBC01_usnob.asc
First 5 lines (2 standard out) for the usno_00.py run:
% usno_00.py UBC01_usnob.asc > a
% head -5 a
Ra_hrs,Dec_deg,magB1,sgB1,magR1,sgR1,magB2,sgB2,magR2,sgR2, I use all 440 characters here, but it is alright!!!!
11.22171093 +44.1951028 20.950 4 19.730 10 19.740 1 19.270 3
11.22198296 +44.1873444 20.790 8 18.820 4 17.770 1 18.070 1
11.22238593 +44.1928944 20.570 9 19.060 1 19.900 2 18.970 9
11.22489741 +44.1764389 20.970 6 18.050 10 20.190 1 18.410 2
I am referring to output of usno_00.py as a "cdfp" (CoorDinate
Floating Point) file. The first line will always be a comma-delineated
strings that deignate the parameter names for the following lines
of floating point data. The first two floating point values in
each line will always be RA in hours units and DEC in degree units.
I have written the OTW code called cdfp00 that can process the above file and create a new file with various coordinate-related columns. I have also modified the XYP plotting code so that it can plot data from such a file if it has a name like: AAA.cdpf. A sample run is shown below:
===============================================================
#!/bin/bash
# ./BRUN 11.2337736 +44.3347747 0.2 19.5
# Check command line arguments
if [ -z "$1" ]
then
printf " RA_hrs DEC_deg r' Rlim\n"
printf "./BRUN 11.2337736 +44.3347747 0.200 19.5 "
exit
fi
ubc01 $1 $2 $3 $4
usno_00.py UBC01_usnob.asc > temp1.cdfp
cdfp00.sh temp1.cdfp $1 $2
===============================================================
% ls
BRUN* S/
% ./BRUN 11.2337736 +44.3347747 0.4 22 > A.cdfp
% ls
A.cdfp BRUN* junc runner* S/ temp1 UBC01_usnob.asc Usno.file1
NOW, I cun run XYP on the files A.cdfp and temp1.cdfp.