casa6

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

# test_task_simanalyze.py

#

# Copyright (C) 2020

# Associated Universities, Inc. Washington DC, USA.

#

# This script is free software; you can redistribute it and/or modify it

# under the terms of the GNU Library General Public License as published by

# the Free Software Foundation; either version 2 of the License, or (at your

# option) any later version.

#

# This library is distributed in the hope that it will be useful, but WITHOUT

# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public

# License for more details.

#

# https://open-jira.nrao.edu/browse/CAS-3669

#

# Based on the requirements listed in plone found here:

# https://casadocs.readthedocs.io/en/stable/api/tt/casatasks.simulation.simanalyze.html

#

# Test case: requirement

#

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

import os

import unittest

import shutil

​

from casatools import ctsys, image, coordsys, componentlist, quanta

ia = image()

cs = coordsys()

cl = componentlist()

qa = quanta()

from casatasks import tclean, simobserve, simanalyze, casalog

​

# DATA #

configpath_int      = ctsys.resolve('alma/simmos/vla.a.cfg')

imagepath_int       = ctsys.resolve('nrao/VLA/CalModels/3C286_Q.im/')

configpath_both_int = ctsys.resolve('alma/simmos/aca.cycle7.cfg')

configpath_sd       = ctsys.resolve('alma/simmos/aca.tp.cfg')

# Path to MS data

mspath_sd = ctsys.resolve('unittest/simanalyze/refim_twopoints_twochan.ms')

​

logpath = casalog.logfile()

​

int_project = 'sim_interferometric'

sd_project = 'sim_single_dish'

imagepath_sd = 'sd_model.image'

both_project = 'sim_both'

both_component_list = 'J0319+4130.cl'

both_model_image_int = both_component_list[:-3]+'.im'

both_model_image_sd = both_component_list[:-3]+'_rebin.im'

​

####    interferometric    ####

class simanalyze_main_usage_modes_test_int(unittest.TestCase):

​

    def setUp(self):

        simobserve(project=int_project, skymodel=imagepath_int, complist='', setpointings=True, direction=[], 

                   mapsize='', maptype='square', pointingspacing='',caldirection='',calflux='1Jy',obsmode='int', 

                   refdate='2020/02/13', hourangle='transit', totaltime='100s', antennalist=configpath_int,

                   outframe='LSRK', thermalnoise='', leakage=0.0, graphics='none',verbose=False, overwrite=False)

​