pyrad.proc.process_selfconsistency_kdp_phidp#

pyrad.proc.process_selfconsistency_kdp_phidp(procstatus, dscfg, radar_list=None)[source]#

Computes specific differential phase and differential phase in rain using the selfconsistency between Zdr, Zh and KDP

Parameters:
  • procstatus (int) – Processing status: 0 initializing, 1 processing volume, 2 post-processing

  • dscfg (dictionary of dictionaries) –

    data set configuration. Accepted Configuration Keywords:

    datatype : list of strings. Dataset keyword
        The input data types, must contain
        "dBZ" or "dBZc", and,
        "ZDR" or "ZDRc", and,
        "PhiDP" or "PhiDPc", and,
        "uRhoHV" or "RhoHV" or "RhoHVc", and,
        "TEMP" (Optional), and,
        "H_ISO0" (Optional), and,
        "hydro" (Optional, only used if filter_rain)
    parametrization : str
        The type of parametrization for the self-consistency curves. Can
        be 'None', 'Gourley', 'Wolfensberger', 'Louf', 'Gorgucci' or
        'Vaccarono'
        'None' will use tables from config files. Default 'None'.
    rsmooth : float. Dataset keyword
        length of the smoothing window [m]. Default 2000.
    min_rhohv : float. Dataset keyword
        minimum valid RhoHV. Default 0.92
    filter_rain : Bool. Dataset keyword
        If True the hydrometeor classification is used to filter out gates
        that are not rain. Default True
    max_phidp : float. Dataset keyword
        maximum valid PhiDP [deg]. Default 20.
    ml_thickness : float. Dataset keyword
        assumed melting layer thickness [m]. Default 700.
    fzl : float. Dataset keyword
        The default freezing level height. It will be used if no
        temperature field name is specified or the temperature field is
        not in the radar object. Default 2000.
    sounding : str. Dataset keyword
        The nearest radiosounding WMO code (5 int digits). It will be used to
        compute the freezing level, if no temperature field name is specified,
        if the temperature field isin the radar object or if no freezing_level
        is explicitely defined.
    frequency : float. Dataset keyword
        the radar frequency [Hz]. If None that of the key
        frequency in attribute instrument_parameters of the radar
        object will be used. If the key or the attribute are not present
        the selfconsistency will not be computed
    
  • radar_list (list of Radar objects) – Optional. list of radar objects

Returns:

  • new_dataset (dict) – dictionary containing the output fields KDP and PhiDP “”

  • ind_rad (int) – radar index