def testconfigs(avgsepsf="configs.avgsep"): """ Test ALMA configurations. """ # read avgseps. try: avsfile = open(avgsepsf, 'r') avgseps = map(float, avsfile.readlines()) avsfile.close() except: print "Could not read", avgsepsf return dec = -23 direcshun = "J2000 18h00m00.03s %sd0m0.0s" % dec # cl.addcomponent(dir=direcshun,flux=1.0, freq='672.0GHz') compname = "star672GHz_dec%s.cl" % dec # cl.rename(compname) # cl.done() # for (ttime, integrity) in zip(('30s', '14400s'), ('10s', '300s')): ttime = '14400s' integrity = '300s' for confnum in range(1, 29): confname = "out%02d" % confnum print "confname:", confname pixelsize = 0.1 * avgseps[0] / avgseps[confnum - 1] cell = "%farcsec" % pixelsize antlist = "almaconfig.%s" % confname project = "testalmaconfig_%s_%s_%s_briggs_nn" % (confname, ttime, dec) print "Starting", project almasimmos(project=project, modelimage="", complist=compname, antennalist=antlist, direction=direcshun, nmosx=1, nmosy=1, pointingspacing="0.01arcsec", refdate="2008/06/21/03:25:00", totaltime=ttime, integration=integrity, mode="none", alg="clark", niter=50, nchan=-1, startfreq="668.0GHz", chanwidth="8.0GHz", cell=cell, imsize=[256, 256], stokes="I", weighting="briggs", rmode="norm", robust=0.0, display=False) print "Finished", project