gnssrefl.quickLook_function module

gnssrefl.quickLook_function.goodbad(fname, station, year, doy, h1, h2, PkNoise, reqAmp, freq, e1, e2, hires_figs)

makes a plot that shows “good” and “bad” refletor height retrievals as a function of azimuth

Parameters:

  • fname (str) – filename

  • station (str) – 4 char station name

  • year (int) – full year

  • doy (int) – day of year

  • h1 (float) – minimum reflector height (m)

  • h2 (float) – max reflector height (m)

  • PkNoise (float) – peak 2 noise QC

  • reqAmp (float) – required LSP amplitude

  • freq (int) – frequency

  • e1 (float) – minimum elevation angle (deg)

  • e2 (float) – maximum elevation angle (deg)

  • hires_figs (bool) – whether you want eps instead of png

Return type:

plot is written $REFL_CODE/Files/station/quickLook_summary.png

gnssrefl.quickLook_function.quickLook_function(station, year, doy, snr_type, f, e1, e2, minH, maxH, reqAmp, pele, satsel, PkNoise, fortran, pltscreen, azim1, azim2, ediff, delTmax, hires_figs, **kwargs)

This is the main function to compute spectral characteristics of a SNR file. It takes in all user inputs and calculates reflector heights. It makes two png files to summarize the data.

Parameters:

  • station (str) – station name (4 char)

  • year (int) – full year

  • doy (int) – day of year

  • snr_type (int) – snr file extension (i.e. 99, 66 etc)

  • f (int) – frequency (1, 2, 5), etc

  • e1 (float) – minimum elevation angle in degrees

  • e2 (float) – maximum elevation angle in degrees

  • minH (float) – minimum allowed reflector height in meters

  • maxH (float) – maximum allowed reflector height in meters

  • reqAmp (float) – is LSP amplitude significance criterion

  • pele (list of floats) – is the elevation angle limits for the polynomial removal. units: degrees

  • satsel (int) – satellite number

  • PkNoise (float) – peak to noise ratio for QC

  • fortran (bool) – whether external fortran translator is being explicitly called.

  • pltscreen (bool) – whether you want plots sent to the terminal

  • azim1 (float) – minimum azimuth in degrees

  • azim2 (float) – maximum azimuth in degrees

  • ediff (float) – QC parameter - restricts length of arcs (degrees)

  • delTmax (float) – maximum arc length in minutes

  • hires_figs (bool) – whether to use eps instead of png

gnssrefl.quickLook_function.read_snr_simple(obsfile)

Loads the data from a SNR file into memory

Parameters:

obsfile (string) – name of the SNR file

Returns:

  • allGood (int) – one for all is well, zero is for all is bad

  • sat (numpy array of floats) – satellite numbers

  • ele (numpy array of floats) – elevation angle (deg)

  • azi (numpy array of float) – azimuth (deg)

  • t (numpy array of floats) – seconds of the day (no leap seconds)

  • edot (numpy array of floats) – derivative of elevation angle wrt time deg/sec

  • s1 (numpy array of floats) – L1 SNR

  • s2 (numpy array of floats) – L2 SNR

  • s5 (numpy array of floats) – L5 SNR

  • s6 (numpy array of floats) – L6 SNR

  • s7 (numpy array of floats) – L7 SNR

  • s8 (numpy array of floats) – L8 SNR

  • snrE (numpy array of booleans) – whether SNR exists