def cleansim(project, nfld=1, alg='clark',imsize=[512,512],
             cell='0.0015arcsec',
             direction="J2000 18h00m00.0 -23d00m0.00"):
    msfile = project + '.ms'
    im.open(msfile)
    im.selectvis(field=range(nfld), spw=0)
    im.defineimage(nx=imsize[0], ny=imsize[1], cellx=cell, phasecenter=direction)
    im.setvp(dovp=True)
    im.setoptions(ftmachine='mosaic')
    im.make(project+'.model')
    im.clean(algorithm='mf'+alg, model=[project+'.model'], image=[project+'.restored'], residual=[project+'.residual'])
#    im.makeimage(type='psf',image=project+'.psf')
#    beam=im.fitpsf(psf=project+'.psf')
    im.done()