pyart.correct.selfconsistency_bias#

pyart.correct.selfconsistency_bias(radar, zdr_kdpzh_dict, min_rhohv=0.92, filter_rain=True, max_phidp=20.0, smooth_wind_len=5, doc=None, fzl=None, thickness=700.0, min_rcons=20, dphidp_min=2, dphidp_max=16, parametrization='None', refl_field=None, phidp_field=None, zdr_field=None, temp_field=None, iso0_field=None, hydro_field=None, rhohv_field=None, temp_ref='temperature', check_wet_radome=True, wet_radome_refl=25.0, wet_radome_len_min=4, wet_radome_len_max=8, wet_radome_ngates_min=180, valid_gates_only=False, keep_points=False, kdp_wind_len=12)[source]#

Estimates reflectivity bias at each ray using the self-consistency algorithm by Gourley

Parameters:

  • radar (Radar) – radar object

  • zdr_kdpzh_dict (dict) – dictionary containing a look up table relating ZDR with KDP/Zh for different elevations

  • min_rhohv (float) – minimum RhoHV value to consider the data valid

  • filter_rain (bool) – If True the hydrometeor classification is going to be used to filter out all gates that are not rain

  • max_phidp (float) – maximum PhiDP value to consider the data valid

  • smooth_wind_len (int) – length of the smoothing window

  • doc (float) – Number of gates at the end of each ray to to remove from the calculation.

  • fzl (float) – Freezing layer, gates above this point are not included in the correction.

  • thickness (float) – assumed melting layer thickness [m]

  • min_rcons (int) – minimum number of consecutive gates to consider a valid segment of PhiDP

  • dphidp_min (float) – minimum differential phase shift in a segment

  • dphidp_max (float) – maximum differential phase shift in a segment

  • parametrization (str) – The type of parametrization for the self-consistency curves. Can be ‘None’, ‘Gourley’, ‘Wolfensberger’, ‘Louf’, ‘Gorgucci’ or ‘Vaccarono’. ‘None’ will use tables contained in zdr_kdpzh_dict.

  • refl_field, phidp_field, zdr_field (str) – Field names within the radar object which represent the reflectivity, differential phase and differential reflectivity fields. A value of None will use the default field name as defined in the Py-ART configuration file.

  • temp_field, iso0_field, hydro_field, rhohv_field (str) – Field names within the radar object which represent the temperature, the height relative to the iso0, the hydrometeor classification and the co-polar correlation fields. A value of None will use the default field name as defined in the Py-ART configuration file. They are going to be used only if available.

  • kdpsim_field, phidpsim_field (str) – Field names which represent the estimated specific differential phase and differential phase. A value of None will use the default field name as defined in the Py-ART configuration file.

  • temp_ref (str) – the field use as reference for temperature. Can be either temperature, height_over_iso0 or fixed_fzl

  • check_wet_radome (Bool) – if True the average reflectivity of the closest gates to the radar is going to be check to find out whether there is rain over the radome. If there is rain no bias will be computed

  • wet_radome_refl (Float) – Average reflectivity of the gates close to the radar to consider the radome as wet

  • wet_radome_len_min, wet_radome_len_max (int) – Mim and max gate indices of the disk around the radome used to decide whether the radome is wet

  • wet_radome_ngates_min (int) – Minimum number of valid gates to consider that the radome is wet

  • valid_gates_only (Bool) – If True the reflectivity bias obtained for each valid ray is going to be assigned only to gates of the segment used. That will give more weight to longer segments when computing the total bias.

  • keep_points (Bool) – If True the ZDR, ZH and KDP of the gates used in the self- consistency algorithm are going to be stored for further analysis

  • kdp_wind_len (int) – The length of the window used to compute KDP with the single window least square method

Returns:

  • refl_bias_dict (dict) – the bias at each ray field and metadata

  • selfconsistency_dict (dict) – If keep_poinst set, a dictionary containing the measured valid values of ZDR, Zh and KDP. None otherwise