set term postscript eps enhanced colour # define convert functions RA2POS( hour, min, sec) = hour + min/60.0 + sec/60.0/60.0 DEC2POS( deg, arcmin, arcsec) = deg + (deg < 0 ? -arcmin : arcmin)/60.0 + (deg < 0 ? -arcsec : arcsec)/60.0/60.0 set size square set xlabel "RA" set ylabel "DEC" RA = RA2POS(5,20,0) DEC = DEC2POS(3,0,0) delta = 27 dRA = 24 * delta / 360.0 dDEC = delta set xrange [RA+dRA:RA-dRA] set yrange [DEC-dDEC:DEC+dDEC] set xtics 1 set ytics 3 set output "plot03.eps" plot "crop.dat"\ using (($7 < 1.0) ? RA2POS($1,$2,$3) : 1/0):(DEC2POS($4,$5,$6))\ notitle\ with points lt 1 pt 13 ps 1.2,\ ""\ using ((1.0 <= $7 && $7 < 2.0) ? RA2POS($1,$2,$3) : 1/0):(DEC2POS($4,$5,$6))\ notitle\ with points lt 1 pt 13 ps 1.0,\ ""\ using ((2.0 <= $7 && $7 < 3.0) ? RA2POS($1,$2,$3) : 1/0):(DEC2POS($4,$5,$6))\ notitle\ with points lt 1 pt 13 ps 0.8,\ ""\ using ((3.0 <= $7 && $7 < 4.0) ? RA2POS($1,$2,$3) : 1/0):(DEC2POS($4,$5,$6))\ notitle\ with points lt 1 pt 13 ps 0.6,\ ""\ using ((5.0 <= $7) ? RA2POS($1,$2,$3) : 1/0):(DEC2POS($4,$5,$6))\ notitle\ with points lt 1 pt 13 ps 0.2