casa6

import os

import time

import regression_utility as regutl

import shutil

import sys

import pdb

​

THISHOME  = "ic2233_regression_data/";

IMSIZE    = 2048;

#IMSIZE    = 4096;

CELLSIZE  = "4arcsec";

STOKES    = "IV";

CHSTART   = 69;

NCHAN     = 11;

#CHSTART = 64;

#NCHAN =  1;

SPW   ="*:69~80";

FIELD  = "0";

PASTEP    = 360.0;

NITER     = 6000;

#IMAGE     = "imIC2233.pbm0";

DECONVOLVER = "hogbom";

FTMACHINE = "awproject"

PSTERM = False;

IMAGE     = "imIC2233."+FTMACHINE;

REUSEREPOS = False;

INTERACTIVE = False;

TIMERANGE = "";"*+0:30:0";

EPS       = 1e-5;  # Logical "zero"

​

​

#--------------------------------------------------------------

#

def ic2233_reg():

    MSFILE    = THISHOME  + "ic2233.lband.ms";

    MYIMAGE   = THISHOME  + IMAGE;

​

    #os.system("rm -rf "+MYIMAGE+"*");

    if (REUSEREPOS==False):

        os.system("rm -rf "+THISHOME);

        os.mkdir(THISHOME);

        REPOSNAME = os.environ.get('CASAPATH').split()[0]

        PREFIX    = REPOSNAME + "/data/regression/ic2233/";

        FITSFILE  = PREFIX    + "ic2233.lband.fits";

        ms.fromfits(fitsfile=FITSFILE,msfile=MSFILE);

        ms.done();

​

    tclean(vis=MSFILE,selectdata=True,field=FIELD,spw=SPW,timerange=TIMERANGE,uvrange="",antenna="",scan="",observation="",

           intent="",datacolumn="corrected",imagename=MYIMAGE,imsize=IMSIZE,cell=CELLSIZE,phasecenter="",stokes=STOKES,projection="SIN",

           startmodel="",specmode="mfs",reffreq="",nchan=-1,start="",width="",outframe="LSRK",veltype="radio",restfreq=[],interpolation="linear",

           gridder=FTMACHINE,facets=1,chanchunks=1,wprojplanes=1,vptable="",aterm=True,psterm=PSTERM,wbawp=True,conjbeams=False,

           cfcache=MYIMAGE+".cf",computepastep=360.0,rotatepastep=5.0,pblimit=0.05,normtype="flatnoise",deconvolver=DECONVOLVER,scales=[],

           nterms=2,smallscalebias=0.6,restoration=True,restoringbeam=[],pbcor=False,outlierfile="",weighting="natural",robust=0.5,npixels=0,

           uvtaper=[],niter=NITER,gain=0.1,threshold=0.0,cycleniter=-1,cyclefactor=1.0,minpsffraction=0.05,maxpsffraction=0.8,interactive=INTERACTIVE,

           usemask="user",mask="",pbmask=0.0,restart=True,savemodel="none",

           calcres=True,calcpsf=True,parallel=False);

​

#    tclean(vis=MSFILE,selectdata=True,field="0",spw="*:69~80",timerange="",uvrange="",antenna="",scan="",observation="",intent="",datacolumn="corrected",imagename=MYIMAGE,imsize=2048,cell="4arcsec",phasecenter="",stokes="IV",projection="SIN",startmodel="",specmode="mfs",reffreq="",nchan=-1,start="",width="",outframe="LSRK",veltype="radio",restfreq=[],interpolation="linear",gridder="awproject",facets=1,chanchunks=1,wprojplanes=1,vptable="",aterm=True,psterm=False,wbawp=True,conjbeams=True,cfcache="junk_newVR.cf",computepastep=360.0,rotatepastep=5.0,pblimit=0.2,normtype="flatnoise",deconvolver="hogbom",scales=[],nterms=2,smallscalebias=0.6,restoration=True,restoringbeam=[],pbcor=False,outlierfile="",weighting="natural",robust=0.5,npixels=0,uvtaper=[],niter=6000,gain=0.1,threshold=0.0,cycleniter=-1,cyclefactor=1.0,minpsffraction=0.05,maxpsffraction=0.8,interactive=True,usemask="user",mask="",pbmask=0.0,maskthreshold="",maskresolution="",nmask=0,autoadjust=False,restart=True,savemodel="none",calcres=True,calcpsf=True,parallel=False);

​

​

​

    # im.open(MSFILE,usescratch=True);