common comp,r,x1,y1,t,cx,cy,ax,ay,fo,sa,clx,cly,sax,say,ashx,ashy,tx,ty
common mpsf,nx,ny,np,psfm,wp,xp,yp,gp
common kluge,nxf,nyf
common done,done,crve,iscr
common numb,n
common filen,fname,logfil
common prtt,wave,mfc
common iter,titer
common spher,csph
common camer,camera,lam
common batch, inter
common wfpc, obscsize

logfil=99
reads,'image6',psf,nx,ny
nxo=nx & nyo=ny
fname='image6temp'             
openw,logfil,'image6.log'
psf=psf(25:nx-25,25:ny-25)
nx=n_elements(psf(*,0)) & ny=n_elements(psf(0,*))
np=1 &ip=1
gp=[0]
xp=[400] &yp=[400]
wp=[1.]
psfm=fltarr(nx,ny,1)
psfm(0,0,0)=psf
done=0 & crve=1 & iscr=0
n=5*nx+3
print,'n,nx',n,nx
titer=-1
csph=0.
camera='FCFULL'
lam=.1200
obscsize=0.33
a=[0.,0.,-1.157,0.,0.,0.,0.,-0.5]*0.547/lam
fun,x,a,f
reads,fname+'sim',image6,nx,ny
out=fltarr(nxo,nyo)
out(nxo/2-nx/2,nyo/2-ny/2)=image6
writes,'image6sim',out
close,logfil
stop
end