#!/usr/bin/env python
"""
reading PFISR data down to IQ samples
See README.rst for the data types this file handles.
"""
from matplotlib.pyplot import show
#
from isrutils.common import boilerplateapi
from isrutils.switchyard import isrselect
from isrutils.plots import plotsnr
#%% read data
p,isrfn,odir,tlim = boilerplateapi()
spec,freq,snrsamp,azel,isrlla,snrint,snr30int = isrselect(
isrfn,p.beamid,tlim,p.zlim,p.t0,p.acf,p.samples)
#%% plot data
plotsnr(snrint,isrfn,tlim)
show()
#!/usr/bin/env python
"""
reading PFISR data down to IQ samples
See Examples/ for more updated specific code
"""
from isrutils import Path
from isrutils.common import boilerplateapi
from isrutils.looper import simpleloop
#%% user parameters
p = boilerplateapi()
P={
'path':Path(p.isrfn).expanduser(),
'beamid': 64157,
'acf': p.acf,
'vlimacf': (18,45),
'zlim_pl': p.zlim,
'vlim_pl': [72,90],
'flim_pl': [3.5,5.5],
'odir': p.odir,
'vlim': [25, 55],
'zlim': p.zlim,
'tlim': p.tlim,
'verbose': True,}
#%% read and plot
#flist = [x for x in Path(P['path']).expanduser().iterdir() if x.suffix=='.h5']
flist = Path(P['path']).expanduser().glob('*.dt{}.h5'.format(p.rtype))
simpleloop(flist,P)
def main():
overlayisrhist(boilerplateapi())
show()
