pyrad.proc.process_phidp_kdp_Maesaka#

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

Estimates PhiDP and KDP using the method by Maesaka. This method only retrieves data in rain (i.e. below the melting layer)

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 string. Dataset keyword
    The input data types, must contain,
    "dBZ" or "dBZc", and,
    "PhiDP" or "PhiDPc" or "uPhiDP", and
    "TEMP" or "H_ISO0" (Optional)
rmin : float. Dataset keyword
    The minimum range where to look for valid data [m]. Default 1000.
rmax : float. Dataset keyword
    The maximum range where to look for valid data [m]. Default 50000.
rcell : float. Dataset keyword
    The length of a continuous cell to consider it valid precip [m].
    Default 1000.
Zmin : float. Dataset keyword
    The minimum reflectivity [dBZ]. Default 20
Zmax : float. Dataset keyword
    The maximum reflectivity [dBZ]. Default 40
fzl : float. Dataset keyword
    The freezing level height [m]. 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.
ml_thickness : float. Dataset keyword
    The melting layer thickness in meters. Default 700.
beamwidth : float. Dataset keyword
    the antenna beamwidth [deg]. If None that of the keys
    radar_beam_width_h or radar_beam_width_v in attribute
    instrument_parameters of the radar object will be used. If the key
    or the attribute are not present the beamwidth will be set to None

radar_list (list of Radar objects) – Optional. list of radar objects

Returns:

new_dataset (dict) – dictionary containing the output field “PhiDPc” and “KDPc”
ind_rad (int) – radar index