casa6

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

#                                                                    #

# Use Case Demo Script for Jupiter 6cm VLA                           #

# Trimmed down from Use Case jupiter6cm_usecase.py                   #

#                                                                    #

# Assumes you have already flagged using jupiter6cm_flagdemo.py      #

# and calibratied using jupiter6cm_caldemo.py                        #

# Will do imaging                                                    #

#                                                                    #

# Updated STM 2008-05-26 (Beta Patch 2.0) use new clean task         #

# Updated STM 2008-06-12 (Beta Patch 2.0) for summer school demo     #

#                                                                    #

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

​

import time

import os

​

# 

#=====================================================================

#

# This script has some interactive commands: scriptmode = True

# if you are running it and want it to stop during interactive parts.

​

scriptmode = True

​

#=====================================================================

#

# Set up some useful variables - these will be set during the script

# also, but if you want to restart the script in the middle here

# they are in one place:

​

prefix='jupiter6cm.demo'

​

msfile = prefix + '.ms'

​

#

# Special prefix for calibration demo output

#

calprefix = prefix + '.cal'

​

srcname = 'JUPITER'

srcsplitms = calprefix + '.' + srcname + '.split.ms'

​

#

#=====================================================================

#

# Special prefix for this imaging demo output

#

imprefix = prefix + '.img'

​

# Clean up old files

os.system('rm -rf '+imprefix+'*')

​

#

#=====================================================================

#

# Imaging parameters

#

# This is D-config VLA 6cm (4.85GHz) obs

# Check the observational status summary

# Primary beam FWHM = 45'/f_GHz = 557"

# Synthesized beam FWHM = 14"

# RMS in 10min (600s) = 0.06 mJy (thats now, but close enough)