casa6

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

# Regression Script for simdata of  #

#        a protoplanetary disk         #

​

import os, time

​

# configs are in the repository

l=locals() 

if not l.has_key("repodir"): 

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

​

startTime = time.time()

startProc = time.clock()

​

print '--Running simdata of input672GHz_50pc.image--'

​

my_project="psim2"

my_modelimage="diskmodel.im2"

​

# Clear out results from previous runs.

os.system('rm -rf '+my_project+'.* '+my_modelimage)

tb.clearlocks()

​

print casa['build']

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

#importfits(fitsimage=repodir+"/data/alma/simmos/input50pc_672GHz.fits",imagename=my_modelimage)

​

default("simanalyze")

default("simobserve")

project=my_project

skymodel=my_modelimage

skymodel=repodir+"/data/regression/simdata/input50pc_672GHz.fits"

direction="J2000 18h00m00.031s -22d59m59.6s"

​

complist="star672GHz.cl"

if os.path.exists(complist):

    import shutil

    shutil.rmtree(complist)

cl.done()

cl.addcomponent(dir=direction,flux=0.0003,freq="672GHz")

cl.rename("star672GHz.cl")

cl.done()

​

setpointings=True

mapsize="0.76arcsec"

pointingspacing="0.5arcsec"

integration="10s"

​

obsmode="int"

antennalist=repodir+"/data/alma/simmos/alma.out20.cfg"

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

totaltime="1200s"

​

thermalnoise="tsys-atm"

#tau0=0.

#t_ground=0.

#t_sky=0.

user_pwv=0.5

​

maptype="ALMA2012"

​

verbose=True

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