pro READ_N2_N3E_INT, aaaajjjd,hd,aaaajjjf,hf,dsec,vmin,vmax,snmin, s,v,th,ph,t,f, $
path, help=help, keep_real_freq=keep_real_freq, source=source, antenna=antenna
if keyword_set(help) then begin
print,'READ_N2_N3E_INT, aaaajjjd,hd,aaaajjjf,hf,dsec,vmin,vmax,snmin, s,v,th,ph,t,f'
return
endif
if strmid(path,0,1,/reverse_offset) ne "/" then path=path+"/"
plist = make_file_list(aaaajjjd,hd,aaaajjjf,hf,level='n2')
if ptr_valid(plist) then begin
dataFileListN2 = *plist
nfiles = n_elements(dataFileListN2)
ptr_free,plist
endif else begin
print,'No N2 data... aborting'
stop
endelse
plist = make_file_list(aaaajjjd,hd,aaaajjjf,hf,level='n3e',source=source)
if ptr_valid(plist) then begin
dataFileListN3e = *plist
nfiles = n_elements(dataFileListN3d)
ptr_free,plist
endif else begin
print,'No N3e data... aborting'
stop
endelse
nDataFile = n_elements(dataFileListN2)
iDHd = aaaajjjd*1.d0+hd/24.d0
iDHf = aaaajjjf*1.d0+hf/24.d0
dsec=dsec*1.d0
td=aj_t97(aaaajjjd)+hd/24.d0
tf=aj_t97(aaaajjjf)+hf/24.d0
nt=(tf-td)*86400.d0/dsec
t=dindgen(nt)*dsec/86400.d0 + td
read_data_binary,path+'bg/bg_'+strmid(path,12,12,/reverse_offset),bg,level='bg'
xf = bg.xf
bg=0b
nf=n_elements(xf)
xs=fltarr(nt,nf) & xth=xs & xph=xs & xv=xs & nsvthph=xs
for iDataFile = 0,nDataFile-1 do begin
strPosTmp = strpos(dataFileListN3e(iDataFile),'n3e/N3e_')+11
iDay = long(strmid(dataFileListN3e(iDataFile),strPosTmp,7))
iHour = fix(strmid(dataFileListN3e(iDataFile),strPosTmp+8,2))
iDH1 = iDay*1.d0+iHour/24.d0
iDH2 = idH1+1.d0/24.
if ((iDH1 ge iDHd) and (iDH2 le iDHf)) then begin
print,dataFileListN2(iDataFile),dataFileListN3e(iDataFile)
read_data_binary,dataFileListN2(iDataFile),data_n2,level='n2'
read_data_binary,dataFileListN3e(iDataFile),data_n3,valid,level='n3e'
if valid then w=where((data_n3.v ge vmin) and (data_n3.v le vmax) and $
(data_n3.sn(0) ge snmin) and (data_n3.sn(1) ge snmin)) $
else w=-1
print,n_elements(w)
if w(0) ne -1 then begin
nn=n_elements(w)
data_n2 = data_n2(data_n3.num)
tt = (data_n2(w).t97 - td)*86400.d0/dsec
ff = data_n2(w).f
ss = data_n3(w).s
vv = data_n3(w).v
tth = data_n3(w).th
pph = data_n3(w).ph
for i=0L,nn-1L do begin
iff=where(abs(xf-ff(i)) le 0.01)
itt=tt(i)
xs(itt,iff)=xs(itt,iff)+ss(i)
xv(itt,iff)=xv(itt,iff)+vv(i)
xth(itt,iff)=xth(itt,iff)+tth(i)
xph(itt,iff)=xph(itt,iff)+pph(i)
nsvthph(itt,iff)=nsvthph(itt,iff)+1.
endfor
endif
endif
endfor
nsvthph=nsvthph > 1
xs=xs/nsvthph
xv=xv/nsvthph
xth=xth/nsvthph
xph=xph/nsvthph
if not keyword_set(keep_real_freq) then begin
fmin=3.7 & fmax=16075.
nff=300
f=10.^(findgen(nff)*(alog10(fmax)-alog10(fmin))/(nff-1.) + alog10(fmin))
s=fltarr(nt,nff)
v=s
th=s
ph=s
for i=0l,nt-1 do begin
s(i,*)=interpol(xs(i,*),xf,f)
v(i,*)=interpol(xv(i,*),xf,f)
th(i,*)=interpol(xth(i,*),xf,f)
ph(i,*)=interpol(xph(i,*),xf,f)
endfor
endif else begin
s = xs
v = xv
th = xth
ph = xph
f = xf
endelse
return
end