Esempio n. 1
0
def inical():
    print "-> INICAL"
    #uvfiles = 'c1.UVF,c2.UVF,t1.UVF,t2.UVF'
    uvfiles = S.uvfiles
    #1 Read in the files and do an initial calibration with calib.infits
    if S.i == 0:
        if S.line.upper() != 'NONE':
            for u in uvfiles.split(','):
                print "SPECR"
                calib.specr(u.replace(S.tag, S.retag), S)

    #2 Run calib.calcals on the calibrator files.
        print "CALCALS"
        calib.calcals([S.cal1, S.cal2])

    #3 Run cal.cal2srcs to copy solutions over to the solutions to the files.
    print "CAL2SRCS"
    calib.cal2srcs(cals=[S.cal1, S.cal2], srcs=[S.src1, S.src2, S.name])

    #4 selfcal

    # Initialize stuff

    # First Dirty Map
    print "0th Loop.\n"
    clean(S)
    S.m4s = S.model
    sc(S)
    S.update(selfcal_options='mfs,amplitude')
    S.save()
Esempio n. 2
0
def inical():
	print "-> INICAL"	
	#uvfiles = 'c1.UVF,c2.UVF,t1.UVF,t2.UVF'
	uvfiles = S.uvfiles;
	#1 Read in the files and do an initial calibration with calib.infits
	if S.i==0:
		if S.line.upper()!='NONE': 
			for u in uvfiles.split(','):
				print "SPECR"
				calib.specr(u.replace(S.tag,S.retag), S)
	
	#2 Run calib.calcals on the calibrator files. 
		print "CALCALS"
		calib.calcals([S.cal1,S.cal2]);
	
	#3 Run cal.cal2srcs to copy solutions over to the solutions to the files.
	print "CAL2SRCS"
	calib.cal2srcs(cals=[S.cal1,S.cal2], srcs=[S.src1,S.src2, S.name]);
	
	#4 selfcal
	
	# Initialize stuff
	
	# First Dirty Map
	print "0th Loop.\n"
	clean(S)
	S.m4s = S.model;
	sc(S);
	S.update(selfcal_options='mfs,amplitude')
	S.save();
Esempio n. 3
0
def inical():
	
	#uvfiles = 'c1.UVF,c2.UVF,t1.UVF,t2.UVF'
	uvfiles = S.uvfiles;
	print uvfiles
	sys.exit(0);
	#1 Read in the files and do an initial calibration with calib.infits
	if S.i==0:
		for u in uvfiles.split(','):
			calib.specr(u.replace('.UVfl','.uv'), S)
	
	#2 Run calib.calcals on the calibrator files. 
		calib.calcals([S.cal1,S.cal2]);
	
	#3 Run cal.cal2srcs to copy solutions over to the solutions to the files.
	calib.cal2srcs(cals=[S.cal1,S.cal2], srcs=[S.src1,S.src2, S.name]);
	
	#4 selfcal
	
	# Initialize stuff
	
	# First Dirty Map
	print "0th Loop....... \n -------"
	gocal(S)
	S.m4s = S.model;
	sc(S);
	S.update(selfcal_options='mfs,amplitude')
	S.save();
Esempio n. 4
0
def specr(u):
    '''
	Pulls a spectral-range from the file u.
	'''
    calib.specr(u, S)
Esempio n. 5
0
def specr(u):
	'''
	Pulls a spectral-range from the file u.
	'''
	calib.specr(u, S)