Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.stsci.edu/~mperrin/software/sources/allskyplot.pro
Дата изменения: Sat Aug 11 03:00:40 2007
Дата индексирования: Sat Mar 1 15:44:16 2014
Кодировка:

Поисковые слова: http www.badastronomy.com phpbb index.php
;+
; NAME: allskyplot
; PURPOSE:
; plot some points on the whole sky
;
; INPUTS:
; KEYWORDS:
; /usecolor read another column from the database and use this to set
; the color for each object. Use this for plotting survey
; completion status!
; color= use this specific color
; if both /usecolor and color= are set, then /usecolor takes
; priority.
;
; /overplot Just overplot the stars
; /ps plot to a postscript file
; /label Write the star names next to the stars as labels.
; declimit= Declination limit. Default = -30 (Lick)
; OUTPUTS:
;
; HISTORY:
; Began 2004-08-26 17:34:39 by Marshall Perrin
; 2005-07-11 /overplot option added.
; 2005-07-12 color handling improved
;-

PRO allskyplot,filename,_extra=_extra,label=label,usecolors=usecolors,ps=ps,$
overplot=overplot,color=color,declimit=declimit,$
database=database

; NOTE: use /color to plot survey status from haebes.lst

if ~(keyword_set(declimit)) then declimit=-30


if keyword_set(database) then begin

common mysql, SQLhandle
mysqlcheck,SQLhandle

q = "select name,ra2000,dec2000 from sources where ra2000 != '' and dec2000 != '';"
mysqlquery,SQLhandle,q,name,ra2000,dec2000

ra = ten_string(ra2000)
dec = ten_string(dec2000)
colorarr =intarr(n_elements(name))+!p.color
thickarr = intarr(n_elements(name)) + 1

endif else begin
if ~(keyword_set(usecolors)) then begin
readcol,filename,name,rh,rm,rs,dd,dm,ds,format="(A,F,F,F,F,F,F)"
if ~(keyword_set(color)) then $
colorarr =intarr(n_elements(name))+!p.color $
else colorarr = replicate(color,n_elements(name))
endif else begin
readcol,filename,name,rh,rm,rs,dd,dm,ds,color1,color2,format="(A,F,F,F,F,F,F,F,F)"
colorarr = intarr(n_elements(name)) + !p.color
thickarr = intarr(n_elements(name)) + 1
if usecolors eq 1 then color =color1 else color=color2
wg = where(color ge 2)
wr = where(color eq 1)
ww = where(color eq 0)
colorarr[wg] = getcolor('green')
colorarr[wr] = getcolor('red')
colorarr[ww] = getcolor('white')
thickarr[wg] = 2
thickarr[wr] = 2
thickarr[ww] = 1
endelse
ra = tenv(rh,rm,rs)
dec = tenv(dd,dm,ds)

endelse

if ~(keyword_set(overplot)) then $
map_set,/mollweide,/horizon,/noborder

; plot constellations
if ~(keyword_set(overplot)) then $
constellations,color=getcolor('blue')

; and the galactic plane
galeql = findgen(360)
galeqb = fltarr(360)
glactc,galeqra,galeqdec,2000.0,galeql,galeqb,2

if ~(keyword_set(overplot)) then $
oplot,galeqra*15,galeqdec,color=getcolor('red')

; now plot the stars
for i=0L,n_elements(name)-1 do begin
oplot,[ra[i]*15],[dec[i]],psym=1,color=colorarr[i],thick=thickarr[i],symsize=thickarr[i]
print,name[i]
;stop
endfor

; and possibly label them
if (keyword_set(label)) then begin
for i=0L,n_elements(name)-1 do begin
xyouts,ra[i]*15,dec[i],name[i],_extra=_extra,color=colorarr[i]
;print,name[i]

endfor
endif

; plot the negative Dec limit for LICK
horiz_ra = findgen(360)
horiz_dec = fltarr(360)+declimit
if ~(keyword_set(overplot)) then $
oplot,horiz_ra,horiz_dec,color=getcolor('yellow')

stop

end