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