;+
; NAME:
;     starmap_plot_2019_ii
; PURPOSE:
;     plot star map during 07-02-2019 solar eclipse
; CATEGORY:
; CALLING SEQUENCE:
; INPUTS: 
; OPTIONAL (KEYWORD) INPUT PARAMETERS:
; OUTPUTS:  
; COMMON BLOCKS:
; SIDE EFFECTS:
; RESTRICTIONS:   
; MODIFICATION HISTORY:
;     (HH) Jul 07 2011
;-


thick = 1 & xthick = 1 & ythick = 1 & charthick = 1
if !d.name eq 'PS' then !p.font = 0 else !p.font = -1
if !D.name eq 'PS' then begin
  thick=5 & xthick=5 & ythick=5 & charthick = 3
  device,/encapsulate,bits=8,file='starmap_2019.eps'
  device,/color
endif

infil = 'starlist_07_02_2019'
stars = rd_tfile(infil,17,1)
mag = float(stars[7,*])
ra = float(stars[10,*])
dec = float(stars[11,*])
name = string(stars[15,*])+string(' ')+string(stars[16,*])
;
; the sunc variables are SUN's true RA and Dec
;
sunc1 = [15.*6.7674559, 23.014800]
sunc2 = [15.*6.7713669, 23.010490]
sunc3 = [15.*6.7752798, 23.006155]
;
; the srad variables are SUN's semidiameters
;
srad1 = 943.88407/3600.
srad2 = 943.88337/3600.
srad3 = 943.88267/3600.
;
xx = cos(2.*!pi*findgen(1001)/1000)
yy = sin(2.*!pi*findgen(1001)/1000)
plot,ra,dec,psym=3,/iso,xr=[sunc3[0]-2.5,sunc3[0]+2.5],$
  yr=[sunc3[1]-2.5,sunc3[1]+2.5],xst=1,yst=1,$
  xtit='Right Ascension, degrees', ytit='Declination, degrees',$
  thick=thick, xthick=xthick, ythick=ythick, charthick=charthick
nstars = n_elements(stars[0,*])
for ii = 0, nstars-1 do oplot, [ra[ii],ra[ii]], [dec[ii],dec[ii]], $
  psym = sym(1), symsiz= (9-mag[ii])/3, thick=thick
oplot,sunc3[0]+srad3*xx,sunc3[1]+srad3*yy, thick=thick
oplot,sunc2[0]+srad2*xx,sunc2[1]+srad2*yy, thick=thick
oplot,sunc1[0]+srad1*xx,sunc1[1]+srad1*yy, thick=thick
for jj = 100,104 do oplot,[jj,jj],[20,26],lines=1,thick=thick
for jj = 21,25 do oplot,[99,105],[jj,jj],lines=1,thick=thick

xyouts,/data,102.9,21.85,'36 Gem',charsi=charsi
xyouts,/data,103.0,23.7,'87 Gem',charsi=charsi
xyouts,/data,99.8,22.1,'HR 2436',charsi=charsi
xyouts,/data,101.6,23.4,'HD 48805',charsi=charsi
xyouts,/data,100.4,22.9,'HD48548',charsi=charsi
xyouts,/data,102.35,23.1,'HD 49381',charsi=charsi
xyouts,/data,102.0,21.55,'HD 49120',charsi=charsi
xyouts,/data,101.4,24.7,'HD 48640',charsi=charsi
xyouts,/data,101.8,22.05,'HD 48994',charsi=charsi
xyouts,/data,102.35,22.4,'HD 49406',charsi=charsi
xyouts,/data,100.05,23.75,'HD 47525',charsi=charsi

if !D.name eq 'PS' then device,/close 
end