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
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
ind_rad (int) – radar index