PRO plotmod,c2,r2,head,cam,filt,exp
!ORDER = 1
file3=''
gettime,head,date,time
;read,yr,mo,da,hr,min,sec,lat,lon,cam,filt,exp
yr = date(2) & mo=date(1) & da=date(0) & hr=time(0) & min=time(1) & sec=time(2) 
lat=58.233 & lon=256.330 

ut= (hr*60. + min)*60. + sec

stephem,yr,mo,da,hr,min,sec,lat,lon,n,desig,az,za,mag,vis
file2 = STRING(cam,filt,da,cam,filt,exp,26,$
	FORMAT='("~steele/data/aur/c",I1,"/f",I1,"/a",I2.2,2I1,I3.3,".",I3.3)')
;OPENR,unt,file2,/GET_LUN
;REPEAT BEGIN
;  READF,unt,utime,file3
;  file3 = STRTRIM(file3,2)
 ; IF (EOF) THEN BEGIN
 ;    print,'enter new utime' & read,hr,min,sec
 ; ENDIF
;ENDREP UNTIL FIX(ut) EQ FIX(utime)
;  CLOSE,unt
;  FREE_LUN,unt

rdkimg,file2,hb,img1
stop
;rdmeandk,0,60,dk60
img=img1
bimg=img
bright=where(mag le 4.00)
;bright=[bright,where(STRMID(desig,11,3) eq 'UMA')]
;bright=where(STRMID(desig,11,3) eq 'UMA')
baz = az(bright)
bza = za(bright)

;pocarot,baz,bza,x,y,z
camaz = -0.40 * !RADEG
warp =0.645403
stretch = 199.870
ccol = 127.964
crow=128.739
;ccol = 130
;crow=128
numst = N_ELEMENTS(bright)
mag = mag(bright)
desig=desig(bright)
az = az(bright)
za =za(bright)
polrec3d,REPLICATE(1.,numst),za,-az,x,y,z,/degrees
rot_3d,3,x,y,z,camaz,xp,yp,zp,/degrees
rho=stretch*SIN(warp*za*!DTOR)
phi=+ATAN(yp,xp)
cl=2.*(ccol+rho*sin(phi))
rw=2.*(crow-rho*COS(phi))
rw = REPLICATE(512.,numst)-rw
;window,0
;tvscl, bimg<300,0,0

notedge = where(SQRT((c2-2.*ccol)^2. + (r2-2.*crow)^2.) LT 2.*135)
stop
ct=c2 & rt=r2
c2 = c2(notedge)
r2 = r2(notedge)
dd=shift(dist(512,512),256,256)
;read,k
x=-x
openw,unit,'/wesson/user/steele/schee/idl/rowcol.dat',/GET_LUN
printf,unit,2,N_ELEMENTS(c2)
FOR i=0,N_ELEMENTS(c2)-1 DO printf,unit,c2(i),r2(i)
close,unit & FREE_LUN,unit
ps=7
rdcols,'/wesson/user/steele/schee/idl/rowcol.dat',n,col,row
dfile = STRING(cam,filt,da,mo,exp,hr,min,sec,$
	FORMAT='("/wesson/user/steele/schee/idl/stpos.c",I1,".f",I1,".g",I2.2,I2.2,I3.3,".",3I2.2)')
OPENW,unit,dfile,/GET_LUN

FOR i=0,N_ELEMENTS(c2)-1 DO BEGIN

  tvscl,bimg < 6000
  plot,c2(i:i),(-r2(i:i)+512),clip=[0,0,511,511],/device,xmargin=[0,0],ymargin=[0,0],/noerase,$
  symsize=2,position=[0,0,511,511],psym=4,xrange=[0,511],yrange=[0,511],xstyle=5, ystyle=5

  j=-1
  REPEAT BEGIN
  j=j+1
  IF j GT (N_ELEMENTS(cl)-1) THEN j=(N_ELEMENTS(cl)-1)
  print,desig(j),mag(j)
  plot,cl(j:j),rw(j:j),clip=[0,0,511,511],/device,xmargin=[0,0],ymargin=[0,0],/noerase,symsize=2,$
   position=[0,0,511,511],psym=ps,xrange=[0,511],yrange=[0,511],xstyle=5, ystyle=5
  ans = '' & read,ans
  if ans eq 'b' then BEGIN
    j=j-5 
    ps=ps+1
    IF ps EQ 8 THEN ps=1  
  ENDIF
  IF ans EQ 'f' THEN j=j+10
  IF ans EQ 'l' THEN j=j+50
  IF ans EQ 's' THEN j=-1

  ;print,desig(j),mag(j)
  ENDREP UNTIL(ans eq 'y' OR ans EQ 'plannet')
  
  IF ans EQ 'y' THEN BEGIN
     print,desig(j),mag(j)
     printf,unit,strmid(desig(j),0,4),az(j)*!DTOR,za(j)*!DTOR,col(i),row(i),$
     format='(I5,2f12.6,2I4)'

     print,desig(j:j),az(j),za(j)
  ENDIF
ENDFOR
close,unit & FREE_LUN,unit
;contour,dd,/follow,levels=[240,260,280,300],c_labels=[intarr(6)+1],/overplot,$
;clip=[0,0,511,511],/device,position=[0,0,511,511]
end