-- NOTE: NONE OF THESE CHANGES ARE NEEDED ANY LONGER

Until the script generator catches up to the CASA 4.2 release, here are changes you will have to make to your reduction script if you choose to run CASA 4.2

  • Comment out the for loops that use setjy to set to fluxes of your quasars. This will be in the gain calibration step (probably step 15)
#  for phaseCalName in ['J0238+1636']:
#    for i in range(len(fc)):
#      if fc[i].find('Flux density for '+phaseCalName) != -1 and re.search('in SpW=[0-9]+(?: \(ref SpW=[0-9]+\))? is: [0-9]+\.[0-9]+', fc[i]) != None:
#        line = (re.search('in SpW=[0-9]+(?: \(ref SpW=[0-9]+\))? is: [0-9]+\.[0-9]+', fc[i])).group(0)
#        spwId = (line.split('='))[1].split()[0]
#        flux = float((line.split(':'))[1].split()[0])
#        setjy(vis = 'uid___A002_X7310ce_X8cf.ms.split',
#          field = phaseCalName.replace(';','*;').split(';')[0],
#          spw = spwId,
#          fluxdensity = [flux,0,0,0])

  • Add in your own setjy commands to set the fluxes of your quasars (the values should come from the .fluxscale file)
  setjy(vis = 'uid___A002_X74a173_X1fdb.ms.split',
        field = '1',
        spw = '0,1,2,3',
        standard = 'manual',
        fluxdensity = [0.690735,0,0,0],
        spix = -0.550191,
        reffreq = '108.405GHz') 

  • Comment out the call to es.fluxscale2
  #if applyonly != True: es.fluxscale2(caltable = 'uid___A002_X74a173_X1fdb.ms.split.ampli_short_inf', removeOutliers=True, msName='uid___A002_X74a173_X1fdb.ms', writeToFile=True, preavg=10000)

-- ScottSchnee - 2013-12-11
Topic revision: r3 - 2014-05-14, ScottSchnee
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding NRAO Public Wiki? Send feedback