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()
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();
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();
def specr(u): ''' Pulls a spectral-range from the file u. ''' calib.specr(u, S)