Just like Example1 I use a script called BRUN to gather ubcone data (around NGC 4618). Then I use the python code make_circles.py to prepare regions file (fk5 format) for overplotting in an image (here I would just manually run ds9 on the image, then manually load the regions file). This example worked quite well THE FIRTS TIME!

% cat BRUN
#!/bin/bash

#Name_RC3      RA_2000DEC_2000    Type     Input_name
#NGC  4618    12:41:33.00 +41:09:03.0  .SBT9..   N4618     
#dssget 190.387497  41.150833 10.0 dss poss1_red
# RA_hrs = 12.6924998

# ./BRUN 12.6924998 41.150833 0.2 19.5
# Check command line arguments
if [ -z "$1" ]
then
    printf "         RA_hrs    DEC_deg    r'  Rlim\n"
    printf "./BRUN 12.6924998 41.150833 0.200 19.5 "
    exit
fi

ubc01 $1 $2 $3 $4  
usno_00.py UBC01_usnob.asc > temp1
cdfp00.sh temp1 $1 $2

 

% ./BRUN 12.6924998 41.150833 0.200 19.5
% make_circles.py reg_for_py > test.reg