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

Parameters:

  • 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) –

gnssrefl.daily_avg.fbias_daily_avg(station)

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

Parameters:

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

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

worker code for daily_avg_cl.py

Parameters:

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

  • extension (str) – folder extension - usually ‘’

  • 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 (str) – subdirectory for output files

Returns:

  • 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)

write out the daily average RH values

Parameters:

  • obstimes (datetime object) –

  • 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

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 removed it

Parameters:

  • 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 –

Returns:

tvall

Return type:

??