fout=[]
fwhm=[]
bmout=[]
fws=['1','2','3','4','5','6','7','8','9','10']

for fw in fws:
    
    if not os.path.exists("1mm.d-20."+fw+"asec.cl"):
        cl.done()
        cl.addcomponent(dir="J2000 0h00m00s -20d00m00s", flux=1,freq='300GHz',shape="Gaussian",majoraxis=fw+"arcsec",minoraxis=fw+"arcsec",positionangle="0deg")
        cl.rename("1mm.d-20."+fw+"asec.cl")
    if not os.path.exists("1mm.d-20."+fw+"asec.image"):
        simdata(project="1mm.d-20."+fw+"asec",complist="1mm.d-20."+fw+"asec.cl",antennalist="compact.v2.cfg",image=T,cell="0.2arcsec",imsize=512,niter=10000,threshold="1mJy",graphics="both",overwrite=T,totaltime="6h",mask=[206,206,306,306])
    ia.open("1mm.d-20."+fw+"asec.image")
    fout.append(ia.statistics()['flux'][0])
    fwhm.append(fw)
    bmout.append([ia.restoringbeam()['major']['value'],ia.restoringbeam()['minor']['value']])
    ia.done()


pl.clf()
pl.plot(fwhm,fout)
pl.plot(fwhm,fout,'o')

import pickle
out=open("1mm.d-20.v2.pickle","wb")
data={'fout':fout,'fwhm':fwhm,'bmout':bmout}
pickle.dump(data,out)
out.close()

pl.savefig("1mm.d-20.v2.png")