gnssrefl.daily_avg module

gnssrefl.daily_avg.daily_avg_stat_plots(obstimes, meanRH, meanAmp, station, txtdir, tv, ngps, nglo, ngal, nbei, test)

plots of results for the daily avg code

  • obstimes (datetime object) –

  • meanRH (numpy array) – daily averaged Reflector Height values in meters

  • meanAmp (numpy array) – daily average RH amplitude

  • station (str) – 4 character station name

  • txtdir (str) – directory for the results

  • tv – is the variable of daily results

  • ngps (numpy array) – number of gps satellites each day

  • nglo (numpy array) – number of glonass satellites each day

  • ngal (numpy array) – number of galileo satellites each day

  • nbei (numpy array) – number of beidou satellites each day

  • test (bool) –


reads QC-RH values and the daily averages computes residuals and estimate the frequency bias for all available frequencies which is printed to the screen


station (str) – station name - 4char - lowercase

gnssrefl.daily_avg.quick_raw(alldatafile2, xdir, station, subdir)

quick plot of the raw RH data. No QC

  • alldatafile2 (str) – name of the raw file to be read

  • xdir (str) – code environ variable (I think)

  • station (str) – 4 ch station name

  • subdir (str) – subdirectory name for results in xdir/Files

gnssrefl.daily_avg.readin_plot_daily(station, extension, year1, year2, fr, alldatafile, csvformat, howBig, ReqTracks, azim1, azim2, test, subdir, plot_limits)

worker code for

It reads in RH files created by gnssir. Applies median filter and saves average results for further analysis

if there is only one RH on a given day - there is no median value and thus nothing will be saved for that day.

  • station (str) – station name, 4 ch, lowercase

  • extension (str) – folder extension - usually empty string

  • year1 (integer) – first year

  • year2 (integer) – last year

  • fr (integer) – 0 for all frequencies. otherwise, it must be a legal frequency (101 for Glonass L1)

  • alldatafile (str) – name of the output filename

  • csvformat (boolean) – whether you want output as csv format

  • howBig (float) – criterion for the median filter, i.e. how far in meters can a RH be from the median for that day, in meters

  • ReqTracks (integer) – is the number of retrievals required per day

  • azim1 (integer) – minimum azimuth, degrees

  • azim2 (integer) – maximum azimuth, degrees

  • test (bool) –

  • subdir (bool) – subdirectory for output files

  • subdir – whether plot limits for the median filter are shown


  • tv (numpy array) – with these values [year, doy, meanRHtoday, len(rh), month, day, stdRH, averageAmplitude] len(rh) is the number of RH on a given day stdRH is the standard deviation of the RH values (meters) averageAmplitude is in volts/volts

  • obstimes (list of datetime objects) – observation times

gnssrefl.daily_avg.write_out_RH_file(obstimes, tv, outfile, csvformat, station, extension)

write out the daily average RH values

  • obstimes (datetime object) – time of observation

  • tv (numpy array) – content of a LSP results file

  • outfile (string) – full name of output file

  • csvformat (boolean) – true if you want csv format output

  • station (str) – 4 ch station name

  • extension (str, optional) – analysis extension name

gnssrefl.daily_avg.write_out_all(allrh, csvformat, NG, yr, doy, d, good, gazim, gfreq, gsat, gamp, gpeak2noise, gutcTime, tvall)

writing out all the RH retrievals to a single file: file ID is allrh) tvall had everything in it, but it was slowing everything down, so i do not do anything with it.

  • allrh (fileID for writing) –

  • csvformat (bool) – whether you are writing to csv file

  • NG (int) – number of lines of results

  • yr (int) – year

  • doy (int) – day of year

  • d (datetime object) –

  • good (float) – reflector height - I think

  • gazim (numpy array of floats) – azimuths

  • gfreq (numpy array of int) – frequencies

  • gsat (numpy array of int) – satellite numbers

  • gamp (numpy array of floats) – amplitudes of periodograms

  • gpeak2noise (numpy array of floats) – peak 2 noise for periodograms

  • gutcTime (numpy array of floats) – time of day in hours

  • tvall



Return type: