if __name__ == '__main__':
from argparse import ArgumentParser
p = ArgumentParser(
description='demonstration of reading a RINEX 2 Observation file')
p.add_argument('obsfn',
help='pathname of RINEX Observation file',
type=str)
p.add_argument('--profile',
help='debug of code via profiling',
action='store_true')
p.add_argument(
'-m',
'--maxchunk',
help='how many chunks to read from start of file (default: whole file)',
type=int,
default=None)
p = p.parse_args()
if not p.profile:
data = demorinex(p.obsfn, p.maxchunk)
data = estimateTEC(data)
plotdata(data)
show()
else:
import cProfile
from pstats import Stats
proffn = 'rinprof_old.pstats'
cProfile.run('demorinex(p.obsfn)', proffn)
p = Stats(proffn)
p.sort_stats('time', 'cumulative').print_stats(50)
"""
f1 = 1575.42e6 #[hz]
f2 = 1227.60e6 #[hz]
k = 80.62 #[m^3/s^2]
tic = time()
data.ix[:,:,'TECslp'] = 2*(f1*f2)**2 / (k*(f1**2-f2**2)) * (data.ix[:,:,'P2'] -data.ix[:,:,'P1'])
print('computed TEC for {} time steps in {} milliseconds.'.format(data.shape[1],(time()-tic)*1000))
return data
if __name__ == '__main__':
from argparse import ArgumentParser
p = ArgumentParser(description='demonstration of reading a RINEX 2 Observation file')
p.add_argument('obsfn',help='pathname of RINEX Observation file',type=str)
p.add_argument('--profile',help='debug of code via profiling',action='store_true')
p.add_argument('-m','--maxchunk',help='how many chunks to read from start of file (default: whole file)',type=int,default=None)
p = p.parse_args()
if not p.profile:
data = demorinex(p.obsfn, p.maxchunk)
data = estimateTEC(data)
plotdata(data)
show()
else:
import cProfile
from pstats import Stats
proffn = 'rinprof_old.pstats'
cProfile.run('demorinex(p.obsfn)',proffn)
p = Stats(proffn)
p.sort_stats('time','cumulative').print_stats(50)
