# Script to observe two sources for baseline determination
# Set both zpectsource1 and zpectsource2
# SKYXXX sources are 5' from corresponding source for blank sky
Catalog("""
coordmode=J2000
DR21 20:39:01.2 42:19:45
3C48 01:37:41.3 33:09:35
SKY3C48 01:37:41.3 33:04:35
3C138 05:21:09.9 16:38:22
3C147 05:42:36.1 49:51:07
3C286 13:31:08.3 30:30:33

SMMJ22174+0013 22:17:24.7 +00:12:42
SMMJ22175+0009 22:17:33.02 +00:09:06.0

SMMJ03025+0006 03:02:26.2 +00:06:24

SMMJ00266+1708 00:26:34.1 +17:08:34
SMMJ00265+1710 00:26:31.28 +17:10:09.6

SMMJ02399-0136 02:39:51.9 -01:36:00
SMMJ02400-0134 02:39:57.9 -01:34:45

SMMJ04431+0210 04:43:07.2 +02:10:24
SMMJ04433+0210 04:43:15.04 +02:10:01.3

SMMJ09429+4658 09:42:54.6 +46:58:44
SMMJ09429+4659 09:42:53.4 +46:59:54

SMMJ14011+0252 14:01:04.9 +02:52:24
MMJ154127+6616 15:41:26.9 +66:14:37
SMMJ16368+4057 16:36:50.4 +40:57:34

coordmode=B1950
M17S_peak 18:17:32 -16:13:00
W51 19:21:24.4 14:24:40
OriA 5:32:49.0 -5:25:16
UpThere 22:00:00 +40:00:00
UpThere2 22:30:00 +35:00:00
UpThere3 22:30:00 +5:00:00
UpThere4 22:30:00 -5:00:00
UpThere5 22:30:00 -15:00:00
UpThere6 22:30:00 -30:00:00
""")

# *** Define target sources (always in pairs):

# Bonnie & Clyde:
#zpectsource1="SMMJ00266+1708"
#zpectsource2="SMMJ00265+1710"

# Donny & Marie:
#zpectsource1="SMMJ02399-0136"
#zpectsource2="SMMJ02400-0134"

# Fred & Ginger:
#zpectsource1="SMMJ04431+0210"
#zpectsource2="SMMJ04433+0210"

# Lucy & Ricky:
zpectsource1="SMMJ09429+4658"
zpectsource2="SMMJ09429+4659"

# pair #5
#zpectsource1="SMMJ22174+0013"
#zpectsource2="SMMJ22175+0009"

# *** Set up parameters of subreflector nodding...
# ... duration of a subreflector nod cycle, in seconds (can be modified):
zpNodCycSR=10.0
# ... number of subreflector nod cycles (can be modified):
# ... integration time in units of zpNodCycSR seconds: 
#     6 -> 1 minute with zpNodCycSR = 10
#     24 -> 4 minutes with zpNodCycSR = 10
zpNumNodSR=18
# ... number of Zpectrometer dumps per nod cycle (please don't modify):
zpDumpPerNod=14
# ... duration of a subreflector move, in seconds (do NOT modify):
MoveTimeSR=0.5
# ... duration of a Zpectrometer readout, in seconds (do NOT modify):
zpReadTime=0.01152
# ... number of Zpectrometer readouts per dump (do NOT modify):
zpReadPerDump=round((zpNodCycSR/zpReadTime)/zpDumpPerNod,0)
# ... total number of Zpectrometer dumps (do NOT modify):
zpRepeat=zpNumNodSR*zpDumpPerNod
# ... dump time, in seconds (do NOT modify):
zpExposure=zpReadPerDump*zpReadTime
# ... skew of Zpectrometer dump, in seconds (do NOT modify):
zpWait=round(1000.0*(zpExposure - (zpExposure-MoveTimeSR)/2.0),0)

SetValues("Zpectrometer",{"exposure":zpExposure})
SetValues("Zpectrometer",{"repeat":zpRepeat})
SetValues("Zpectrometer",{"wait":zpWait})

# *** Slew to source:
#Slew(zpectsource)

# *** Execute the fast subreflector nodding command.
# *** MR12 means refers beam offsets to midpoint between beams 1 & 2.
# *** SubNod argument should be half of full subreflector nod cycle.
#Slew(zpectsource1)
Track(zpectsource1,None,zpNodCycSR*zpNumNodSR,"MR12",submotion=SubNod(zpNodCycSR/2.0))
#Slew(zpectsource2)
Track(zpectsource2,None,zpNodCycSR*zpNumNodSR,"MR12",submotion=SubNod(zpNodCycSR/2.0))