;W0B16105B.D0D   21-JUN-1990 20:02:05.21  log  -3.97590 ?? -312H
;W05K3E02B.D0D   22-JUN-1990 11:35:52.01  log  -0.8968  0H
;W0B17105B.D0D   23-JUN-1990 15:50:11.90  log  2.72538  342H
reads,'wfc1psf1',a1,nx1,ny1
reads,'wfc1psf2',b1,nx1,ny1
reads,'wfc1psf3',c1,nx1,ny1
tvwindow,1,xsize=nx1*2,ysize=ny1*2
;
m=(30.*0.547/2.)/15
reads,'W0B16105B.D0h',a,nx,ny,0
tvscl,alog((a>385)-384)
print,'input star center coordinates'
tvrdc,x,y,/dev
print,'x,y=',x,y
a=(a(x-m*nx1/2:x+m*nx1/2,y-m*ny1/2:y+m*ny1/2)-384)>1
a2=congrid(a,nx1,ny1)
;
reads,'W05K3E02B.D0h',a,nx,ny,0
tvscl,alog((a>385)-384)
print,'input star center coordinates'
tvrdc,x,y,/dev
print,'x,y=',x,y
a=(a(x-m*nx1/2:x+m*nx1/2,y-m*ny1/2:y+m*ny1/2)-384)>1
b2=congrid(a,nx1,ny1)
;
reads,'W0B17105B.D0h',a,nx,ny,0   
tvscl,alog((a>385)-384)
print,'input star center coordinates'
tvrdc,x,y,/dev
print,'x,y=',x,y
a=(a(x-m*nx1/2:x+m*nx1/2,y-m*ny1/2:y+m*ny1/2)-384)>1
c2=congrid(a,nx1,ny1)
;
tvscl,a1,0,0
tvscl,a2,0,ny1
;tvscl,b1,nx1,0
;tvscl,b2,nx1,ny1
tvscl,c1,nx1,0
tvscl,c2,nx1,ny1
stop
end