#
# n21dk01gKQ.v2d K-band plus Q-band
#
# PI: D.J. Kim (MPIfr)
# K-band 256 ch per 64 MHz IF, 0.25 MHz/ch, 1.0 sec AP
# Q-band 128 ch per 64 MHz IF, 0.50 MHz/ch, 0.5 sec AP
#
# VEX has scan "mode=kband" and "mode=qband" according to the active Effelsberg receiver in each scan,
# while Yebes and KVN always has simultaneous K&Q in both of these "mode=..." modes.
#
# The first correlation attempts used v2d 'freqId=<list>', with one v2d for K band, another v2d for Q band.
#
# For some scans this worked fine. However, as soon as array composition is different in some scan, the order of
# frequency entries in of the .input FREQ table changes wrt other scans, and the whole-track 'freqId=<list>' indices
# turn incorrect.
#
# In the end, must correlate K and Q together with the following config:
# ==> 256 ch per 64 MHz IF, 0.25 MHz/ch, 0.5 sec AP
#
# PI must later in CASA/AIPS split out K and Q bands, and
# ==> apply time averaging to K data from 0.5->1.0s, and
# ==> similarly apply spectral averaging to Q data from 0.25->0.50 MHz/channel.
#

vex = n21dk01g.KQ.vex.obs

dataBufferFactor = 24
visBufferLength = 40
startSeries = 1000
exhaustiveAutocorrs = False
minSubarray = 2
singleScan = True
singleSetup = True

# mjdStart = 2021y124d18h00m00s # for Ef: due to weather was stowed and not on source until ~18 UT
antennas = Ky, Ku, Kt, Ef, Ys

ZOOM zoomsKQ
{
    # K USB
    addZoomFreq = freq@21872.0/bw@64.0/noparent@true
    addZoomFreq = freq@21992.0/bw@64.0/noparent@true
    addZoomFreq = freq@22112.0/bw@64.0/noparent@true
    addZoomFreq = freq@22232.0/bw@64.0/noparent@true
    # K LSB
    addZoomFreq = freq@21808.0/bw@64.0/noparent@true
    addZoomFreq = freq@21928.0/bw@64.0/noparent@true
    addZoomFreq = freq@22048.0/bw@64.0/noparent@true
    addZoomFreq = freq@22168.0/bw@64.0/noparent@true

    # Q USB
    addZoomFreq = freq@42848.0/bw@64.0/noparent@true
    addZoomFreq = freq@42968.0/bw@64.0/noparent@true
    addZoomFreq = freq@43088.0/bw@64.0/noparent@true
    addZoomFreq = freq@43208.0/bw@64.0/noparent@true
    # Q LSB
    addZoomFreq = freq@42784.0/bw@64.0/noparent@true
    addZoomFreq = freq@42904.0/bw@64.0/noparent@true
    addZoomFreq = freq@43024.0/bw@64.0/noparent@true
    addZoomFreq = freq@43144.0/bw@64.0/noparent@true
}

ZOOM zoomsKQEfYs
{
    # K LSB
    addZoomFreq = freq@21808.0/bw@64.0/noparent@true
    addZoomFreq = freq@21928.0/bw@64.0/noparent@true
    addZoomFreq = freq@22048.0/bw@64.0/noparent@true
    addZoomFreq = freq@22168.0/bw@64.0/noparent@true

    # Q LSB
    addZoomFreq = freq@42784.0/bw@64.0/noparent@true
    addZoomFreq = freq@42904.0/bw@64.0/noparent@true
    addZoomFreq = freq@43024.0/bw@64.0/noparent@true
    addZoomFreq = freq@43144.0/bw@64.0/noparent@true
}

DATASTREAM EfK { filelist = filelist.ef.k }
DATASTREAM EfQ { filelist = filelist.ef.q }
ANTENNA Ef {
   zoom = zoomsKQEfYs
   datastreams = EfK, EfQ
   # Q-band has +0.12 usec of extra offset vs K-band:
   freqClockOffs = 0,0,0,0, 0,0,0,0, 0.12,0.12,0.12,0.12, 0.12,0.12,0.12,0.12
}

ANTENNA Ys {
   zoom = zoomsKQEfYs
   filelist = filelist.ys
}

ANTENNA Kt {
   zoom = zoomsKQ
   filelist = filelist.kt
   # Q-band has +0.032 usec of extra offset vs K-band:
   freqClockOffs = 0, 0.032
}
ANTENNA Ku {
   zoom = zoomsKQ
   filelist = filelist.ku
}
ANTENNA Ky {
   zoom = zoomsKQ
   filelist = filelist.ky
}

SETUP setupKQ {
  tInt = 0.5
  outputSpecRes = 0.25
  xmacLength = 0
  strideLength = 0
  guardNS = 2000
  doPolar = True
}

RULE kqband {
  setup = setupKQ
  # final clock adj on:
  #   scan No0043 (difx job 1042) 124-172620 3C371 fringe KVN K+Q to Yebes K+Q
  # scan = No0043
  # scan = No0053  # EfYs K 124-1834
  # scan = No0062  # EfYs Q 124-1930
}

# EOPS from /cluster/difx/usno_finals.erp last updated at 2021-06-07 05:42:10 (+0000)
EOP 59337 { xPole=0.109670 yPole=0.437450 tai_utc=37 ut1_utc=-0.182791 }
EOP 59338 { xPole=0.111440 yPole=0.438470 tai_utc=37 ut1_utc=-0.182530 }
EOP 59339 { xPole=0.113220 yPole=0.439390 tai_utc=37 ut1_utc=-0.182397 }
EOP 59340 { xPole=0.115140 yPole=0.440350 tai_utc=37 ut1_utc=-0.182408 }
EOP 59341 { xPole=0.116970 yPole=0.441510 tai_utc=37 ut1_utc=-0.182571 }