;readdata,dir='150515_10:22:14',swmavg=32
;mir_save,'150515.swmavg32'
mir_restore,'150515.swmavg32'

; 1. data checking - continuum, spectra, tsys etc.

; checking pointing data

result=dat_filter(s_f,'"integ" lt "10"',/reset)
if result gt 0 then flag,/flag

; checking continuum with plot_continuum

select,source='16',/p,/re
plot_continuum,/preavg
; ant3, phase jump between 4th and 5th cycles, i.e. [303,347]
; check log, ant 3 lost lock around 13:07-13:10

; checking spectra with plot_spectra

select,/p,/re,band='s49'
plot_spectra
; channel 61

select,band='s49',/reset
uti_chanfix,channel=61

; check spectra on s49 again.
; channel 94 could be a line from the source, need to check
; the full res data. Here we remove it just considering the continuum.
select,band='s49',/reset
uti_chanfix,channel=94

select,source='3c279',/p,/re,band='s49'
plot_spectra
; channel 61; flat from 40 to 400 about 1.1 GHz

select,/p,/re,band='s50'
plot_spectra

select,source='3c454.3',/p,/re
plot_spectra
; phase offset 

;select,/p,/re
;uti_avgband, swmch1=40,swmch2=400

; checking tsys with plot_var

select,/p,/re,band='c1'
plot_var

select,/p,/re,band=['c1','s25']
plot_var

; last cycle should be flagged for ant 1 due to subaru blocking >533

select,ant=1,int=[533,583],/reset
flag,/flag

select,ant=3,int=[303,347],/reset
flag,/flag

; 2. apply tsys correction
select,/p,/re
apply_tsys

; 3. bandpass calibration.

; to find out the different channel numbers in the track
print, sp[0:100].nch 

result=dat_filter(s_f,'"nch" eq "512" and "wt" gt "0"',/reset)
pass_cal,cal_type='pha',preavg=4, smoothing=8,tel_bsl='telescope',refant=8,ntrim=4
; 3c279 yes
; 3c454.3 yes

result=dat_filter(s_f,'"nch" eq "256" and "wt" gt "0"',/reset)
pass_cal,cal_type='pha',preavg=2, smoothing=8,tel_bsl='telescope',refant=8,ntrim=4
; 3c279 yes
; 3c454.3 yes

result=dat_filter(s_f,'"nch" eq "32" and "wt" gt "0"',/reset)
pass_cal,cal_type='pha',smoothing=2,tel_bsl='telescope',refant=8,ntrim=2,frames_p=16

select,source='3c454.3',/p,/re
plot_spectra

select,/p,/re
uti_avgband, swmch1=40,swmch2=390

result=dat_filter(s_f,'"nch" eq "512" and "wt" gt "0"',/reset)
pass_cal,cal_type='amp',preavg=4, smoothing=8, ntrim=4
pass_cal,cal_type='amp',preavg=4, smoothing=4, ntrim=35

result=dat_filter(s_f,'"nch" eq "256" and "wt" gt "0"',/reset)
pass_cal,cal_type='amp',preavg=2, smoothing=8, ntrim=2

result=dat_filter(s_f,'"nch" eq "32" and "wt" gt "0"',/reset)
pass_cal,cal_type='amp',smoothing=2, ntrim=1, frames_p=16

select,source='3c454.3',/p,/re
plot_spectra

select,/p,/re
uti_avgband, swmch1=40,swmch2=390

select,/p,/re
mir_save,'150515.swmavg32.ft.pcal'


; 4. flux measurement

;mir_restore,'150515.swmavg32.ft.pcal'
select,/p,/re, band='c1',sideband='l'
gain_cal,cal_type='pha',/connect,tel_bsl='telescope',refant=8,/non_point
; all yes

select,/p,/re,band='c1',sideband='l'
gain_cal,cal_type='amp',x='hours',poly=0,tel_bsl='telescope',refant=8,/preavg,/non_point
; all no
; Neptune yes


select,/p,/re,band='c1',sideband='l'
flux_measure
;Vector average: 
;#   Source   Flags   Nscans  Flux(Jy)   SNR    meantime    REAL       IMAG
;     3c279            120    8.1642     583     10.88      8.1642     -0.0023
;  1625-254             90    0.8361     453     13.26      0.8361      0.0003
;  1626-298             30    0.8240     279     13.27      0.8240      0.0012

; use 1625-254 yes 0.84

; 5. gain calibration.

mir_restore,'150515.swmavg32.ft.pcal'

select,/p,/re
gain_cal,cal_type='pha', x='hours',/connect,/preavg,tel_bsl='telescope',refant=8
; 1625-254 yes 0.84

gain_cal,cal_type='amp',x='hours',smoothing=0.4,/preavg,tel_bsl='telescope',refant=8


select,/p,/re
idl2miriad,dir='150515.gy146.swmavg32.lc.vis',source='gy146',sideband='l',/cont

select,/p,/re,band='c1',source='gy146',sideband='l'
uti_rms
; 1.41 mJy/beam

; real image rms is 1.64 mJy/beam; without swarm 1.97 mJy/beam
; noise improve by 1.20
; with swarm 2.2ghz, noise should improve by sqrt(6.1/3.94)=1.25