casa6

startTime = time.time()

startProc = time.clock()

​

projname="cltest"

import shutil 

​

if os.path.exists(projname): shutil.rmtree(projname)

​

l=locals() 

if not l.has_key("repodir"): 

    repodir=os.getenv("CASAPATH").split(' ')[0]

​

print casa['build']

print 'I think the data repository is at '+repodir

datadir=repodir+"/data/regression/simdata/"

if os.path.exists("6334.cl"): shutil.rmtree("6334.cl")

shutil.copytree(datadir+"6334.cl","6334.cl")

​

default("simobserve") 

project            =  projname

complist           =  "6334.cl"

compwidth          =  "1.875GHz"

setpointings       =  True

integration        =  "10s"

direction          =  "J2000 17h20m53.2s -35d47m00s"

mapsize            =  "13arcsec"

maptype            =  "ALMA"

pointingspacing    =  ""

#observe            =  True

obsmode            =  "int"

refdate            =  "2014/01/01"

hourangle          =  "transit"

totaltime          =  "7200s"

antennalist        =  "alma.cycle0.extended.cfg"

thermalnoise       =  ""

​

if not l.has_key('interactive'): interactive=False

if interactive:

    graphics="both"

else:

    graphics="file"

​

verbose            =  True

overwrite          =  True

​

inp()

simobserve()

​

​

default("simanalyze") 

project            =  projname

image              =  True

vis                =  "default"

imsize             =  300

imdirection        =  direction

cell               = "0.043294054arcsec"

niter              =  500

threshold          =  "0.1mJy"

analyze            =  True

​

if not l.has_key('interactive'): interactive=False

if interactive:

    graphics="both"