gnssrefl.quickPhase module

gnssrefl.quickPhase.count_result_arcs(result_path)

Count non-comment lines in a result file.

gnssrefl.quickPhase.main()

gnssrefl.quickPhase.parse_arguments()

gnssrefl.quickPhase.process_phase_day(year, doy, args)

Process a single day for sequential mode.

gnssrefl.quickPhase.process_phase_day_worker(worker_args)

Worker for parallel phase processing. Returns arc count for progress bar.

gnssrefl.quickPhase.process_phase_sequential(year, year_end, doy, doy_end, args)

Sequential processing function

gnssrefl.quickPhase.quickphase(station: str, year: int, doy: int, year_end: int | None = None, doy_end: int | None = None, snr: int = 66, fr: str | None = None, e1: float | None = None, e2: float | None = None, plt: bool = False, screenstats: bool = False, gzip: bool | None = None, extension: str = '', par: int | None = None, midnite: bool = True, ampl: float | None = None, savearcs: bool = False, savearcs_format: str | None = None, dec: int | None = None, dbhz: bool | None = None, recompute_lsp: bool = False)

quickphase computes phase, which are subquently used in vwc. The command line call is phase (which maybe we should change).

Examples

phase p038 2021 4

analyzes data for year 2021 and day of year 4

phase p038 2021 1 -doy_end 365

analyzes data for the whole year

phase p038 2021 1 -doy_end 365 -par 5

analyzes data for the whole year using 5 parallel processes

Parameters:

• station (str) – 4 character ID of the station.

• year (int) – full Year to evaluate.

• doy (int) – day of year to evaluate.

• year_end (int, optional) – year to end analysis. Using this option will create a range from year-year_end. Default is None.

• doy_end (int, optional) – Day of year to end analysis. Using this option will create a range of doy-doy_end. If also using year_end, then this will be the day to end analysis in the year_end requested. Default is None.

• snr (int, optional) –

  SNR format. This tells the code what elevation angles are in the SNR file value options:

  66 (default) : data with elevation angles less than 30 degrees

  99 : data with elevation angles between 5 and 30 degrees

  88 : data with all elevation angles

  50 : data with elevation angles less than 10 degrees

• fr (str, optional) – GNSS frequency. Currently only supports L2C. Default is 20 (l2c)

• e1 (float, optional) – Elevation angle lower limit in degrees for the LSP. default is from json (typically 5)

• e2 (float, optional) – Elevation angle upper limit in degrees for the LSP. default is from json (typically 25)

• plt (bool, optional) – Whether to plot results. Default is False

• screenstats (bool, optional) – Whether to print stats to the screen. Default is False

• gzip (bool, optional) – gzip the SNR file after use. Default is True

• par (int, optional) – Number of parallel processes to spawn (up to 10). Default is 1 (single process).

• midnite (bool, optional) – Allow midnight crossings. When True, loads +/- 2 hours from adjacent days. Default is True.

Returns:

• Saves a file for each day in the doy-doy_end range ($REFL_CODE/<year>/phase/<station>/<doy>.txt)

• columns in files – year doy hour phase nv azimuth sat ampl emin emax delT aprioriRH freq estRH pk2noise LSPAmp