casa6

############################################

# Regression Script for simdata2 of a 3d cube #

​

import os, time

import pylab as pl

​

# Clear out results from previous runs.

#os.system('rm -rf testcube2 tc2*')

​

startTime = time.time()

startProc = time.clock()

​

print '--Running simdata of test cube--'

# configs are in the repository

l=locals()

if not l.has_key("repodir"):

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

​

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

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

cfgdir=repodir+"/data/alma/simmos/"

rmtables("testcube2")

importfits(fitsimage=datadir+"testcube.fits",imagename="testcube2")

​

# test

#default("simdata")

#image=False

#project="tc2_simdata"

​

default("simobserve")

project="tc2"

​

skymodel="testcube2"

inbright=".1"

indirection="J2000 19h00m00s -40d00m00s"

incell="0.2arcsec"

incenter="350GHz"

inwidth="0.5MHz"

​

setpointings=False

ptgfile=datadir+"testcube.ptg.txt"

​

obsmode="int"

antennalist=cfgdir+"alma.out01.cfg"

refdate="2012/06/21/03:25:00"

totaltime="7200s"

​

thermalnoise=""

verbose=True

overwrite=True

​

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

if interactive:

    graphics="both"

else:

    graphics="file"

​

inp()

go()

​

endTime = time.time()

endProc = time.clock()

​