gnssrefl.vwc module

gnssrefl.vwc.main()
gnssrefl.vwc.parse_arguments()
gnssrefl.vwc.vwc(station: str, year: int, year_end: int = None, fr: int = 20, plt: bool = True, screenstats: bool = False, min_req_pts_track: int = 50, polyorder: int = -99, minvalperday: int = 10, snow_filter: bool = False, circles: bool = False, subdir: str = None, tmin: str = None, tmax: str = None, warning_value: float = 5.5, auto_removal: bool = False)

Computes VWC from GNSS-IR phase estimates. It concatenates previously computed phase results, makes plots for the four geographic quadrants, computes daily average phase files before converting to volumetric water content (VWC).

Examples

vwc p038 2017

one year for station p038

vwc p038 2015 -year_end 2017

three years for station p038

vwc p038 2015 -year_end 2017 -warning_value 6

warns you about tracks greater than 6 degrees rms

vwc p038 2015 -year_end 2017 -warning_value 6 -auto_removal

makes new list of tracks based on your new warning value

Parameters:

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

  • year (int) – full Year

  • year_end (int, optional) – last year for analysis

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

  • plt (bool, optional) – Whether to produce plots to the screen. Default is True

  • min_req_pts_track (int, optional) – how many points needed to keep a satellite track default is 50

  • polyorder (int) – polynomial order used for leveling. Usually the code picks it but this allows to users to override. Default is -99 which means let the code decide

  • minvalperday (integer) – how many phase measurements are needed for each daily measurement default is 10

  • snow_filter (boolean) – whether you want to attempt to remove points contaminated by snow default is False

  • circles (boolean) – whether you want circles in the final plot (lines are default)

  • subdir (str) – subdirectory in $REFL_CODE/Files for plots and text file outputs

  • tmin (str) – minimum soil texture value, e.g. 0.05

  • tmax (str) – maximum soil texture value, e.g. 0.45

  • warning_value (float) – screen warning about bad tracks (phase rms, in degrees). default is 5.5

  • auto_removal (boolean, optional) – whether to automatically remove tracks that hit your bad track threshold default value is false

Returns:

  • Daily phase results in a file at $REFL_CODE/<year>/phase/<station>_phase.txt – with columns: Year DOY Ph Phsig NormA MM DD

  • VWC results in a file at $$REFL_CODE/<year>/phase/<station>_vwc.txt – with columns: FracYr Year DOY VWC Month Day