#  Template v2d file for DiFX correlation of rk18cj

#  Ground + space
#  Both C,K bands in same VEX setup

vex = rk18cj.vex.obs

singleScan = True
singleSetup = True
tweakIntTime = True

# antennas = GT, IR, GB   # stronges fringes GT-GB K-band, GT-IR C-band according to ASC test dataset info
# antennas = GB, GT, YS  # both bands
# antennas = GT,GB,MC,TR # K-band
antennas = GB, GT

########################################################
# For Radioastron and launching calcif3 without
# a background calc server running
DelayPolyInterval = 60
delayServerType = CALC_9_1_RA_Server
delayServerHandlerType = Pipe
########################################################

ANTENNA GB {
  filelist=filelist.gbt
  phaseCalInt = 0
}
ANTENNA YS { filelist=filelist.ys format=VDIF/8032 phaseCalInt=0 }
ANTENNA NT { filelist=filelist.nt format=VDIF/8032 phaseCalInt=0 }
ANTENNA IR { filelist=filelist.ir format=VDIF/8032 phaseCalInt=0 }
ANTENNA MC { filelist=filelist.mc format=VDIF/8032 phaseCalInt=0 }
ANTENNA TR { filelist=filelist.tr phaseCalInt=0 }

ANTENNA GT {
    filelist=filelist.ra
    ephemObject=RADIOASTRON
    ephemType=RUSSCF
    ephemFile=/Exps2/rk18cj/v1/RA171019_0400_v02.scf.2  # copy from M. Lisakov
    ephemDeltaT=1.0

    ## <filelist_RA_v2d.txt> from RDF->Mark5B conversion:
    SC_time_type=GroundReception

    ## Closed-loop
    #--- scan 1, K-band   prior is -3.3 usec
    #SC_GS_clock_break=start@2017y292d09h00m00s/sync@2017y292d09h00m00s/clockfudge@14.92
    SC_GS_clock_break=start@2017y292d09h00m00s/sync@2017y292d09h00m00s/clockfudge@-2.3
    #-- scan 2, K-band
    # SC_GS_clock_break=start@2017y292d09h15m00s/sync@2017y292d09h15m00s/clockfudge@14.9
    SC_GS_clock_break=start@2017y292d09h15m00s/sync@2017y292d09h15m00s/clockfudge@-2.28

    ## Rubidium
    # -- scan 3, K-band   prior is -3.30469 usec
    SC_GS_clock_break=start@2017y292d09h30m00s/sync@2017y292d09h30m00s/clockfudge@-2.340  # -2.50 us attempt, dly pos near zero, resid 0.160
    # -- scan 4, K-band   prior is -3.25000 usec
    SC_GS_clock_break=start@2017y292d09h45m00s/sync@2017y292d09h45m00s/clockfudge@-2.3125

    ## from wiki; https://deki.mpifr-bonn.mpg.de/01VLBI_Group/Projects/Space_VLBI/RadioAstron/DiFX-RA_correlation
    ## GBT tracking station
    SC_rec_delay=0.0
    GS_Name=Gt
    GS_difxName=Gt
    GS_calcName=Gt
    GS_X= 882879.84013
    GS_Y= -4924482.30801
    GS_Z= 3944130.68519
    GS_dX= -0.01514
    GS_dY= -0.00047
    GS_dZ= -0.00159
    GS_pos_epoch= 51544.5
    GS_axisType=AZEL
    GS_axisOffset0=0.0000
    GS_axisOffset1=0
    GS_axisOffset2=0
    GS_clockEpoch=2017y292d07h56m00s
    GS_clock0=0.0
    GS_clock1=0.0
    SC_pos_offsetEpoch=2017y292d07h56m00s

    ## Rubidium mode
    # "For the Rubidium mode correlation, a constant Rate offset should be introduced: d_rate = -3 * 10e-10 sec/sec 
    #  (sign may be oposite, check which works for us. )"
    #
    # deltaClockRate = [<if0 usec/sec>, <if1 usec/sec>, ...]
    # deltaClockRate = 0,0,+3e-10,+3e-10 # job 3 K resid -17.2 us rate -0.000615 Hz
    # deltaClockRate = 0,0,-3e-10,-3e-10 # job 3 K resid -17.2 us rate -0.000629 Hz
    # deltaClockRate = 0,0,-3e-4,-3e-4   # job 3 K no fringe
    # deltaClockRate = 0,0,+3e-4,+3e-4   # job 3 K fringe but very high rate 6.656 Hz (6.656/22236e6 = 3e-10)
}

SETUP defaultSetup
{
  doPolar = True

  # Fringe rotation: 0=post-F, 1=linear (default), 2=quadratic
  #  90sec scan 2  post-F    : SNR 5.85
  #  90sec scan 2  linear    : SNR 9.48
  #  90sec scan 2  quadratic : SNR 9.48 ie no improvement
  fringeRotOrder = 2

  # freqId:  printDiFXInput.py rk18cj_2.input 
  # fq 0,1 = 22236.000000 L,R
  # fq 2,3 = 4836.000000  L,R
  freqId = 0,1

# Fourfit:
#tInt = 0.512
#nChan = 2048

#tInt = 0.064
#tInt = 0.032
#nChan = 1024

# ASC reference
tInt = 0.128  # docs say 0.125 but looks like 0.128 in ASC FITS
#nChan = 64

# nChan = 32768
fftSpecRes = 0.0078125  # corresponding FFT time is 128 usec
outputSpecRes = 0.25      # ASC: 64ch
#outputSpecRes = 0.125     # 128 ch
#outputSpecRes = 0.0625   # 256 ch
#outputSpecRes = 0.03125
#outputSpecRes = 0.015625

  # For radioastron
#  guardNS = 2048
guardNS = 16384
  maxNSBetweenACAvg = 2000000
}

RULE defaultRule
{
  setup = defaultSetup
  # mode = ra1cm2
  # mode = ra6cm2
  mode = ra6cm1cm
}

# Sources (pointing centers) with recorded data but no offset pointing centers:
SOURCE FAKERA { }
SOURCE 0716+714 { }

# EOPs from /cluster/difx/usno_finals.erp last updated at 2019-01-07 10:32:05 (+0000)
# oper@frontend: v1> getEOP.py --local 2017y292d08h52m00s
# Reading EOP data from  /cluster/difx/usno_finals.erp
# Reading leap second data from /cluster/difx/ut1ls.dat
EOP 58043 { xPole=0.203830 yPole=0.276930 tai_utc=37 ut1_utc=0.299930 }
EOP 58044 { xPole=0.202610 yPole=0.275470 tai_utc=37 ut1_utc=0.298567 }
##EOP 58045 { xPole=0.201300 yPole=0.273810 tai_utc=37 ut1_utc=0.297179 } # DiFX
EOP 58045 { xPole=0.20123 yPole=0.27372 tai_utc=37 ut1_utc=0.2971859 } # ASC reference
EOP 58046 { xPole=0.200010 yPole=0.272210 tai_utc=37 ut1_utc=0.295805 }
EOP 58047 { xPole=0.198720 yPole=0.270790 tai_utc=37 ut1_utc=0.294534 }