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)