inadequate RHOHV values in the radar file

[MohitDahliya]

• Py-ART version: 1.16.1
• Python version: 3.11.5
• Operating System: Wndows

Description

I want to add the Hydrometeor Classification Field into my radar object and I have which does that by considering the Z, ZDR, RHOHV & KDP fileds. The problem is that the RHOHV filed in my radar file is incorrect. maybe when they made the file at the source, they made a mistake. Due to wrong RHOHV values, hydrometeor classification field is not coming out correct. How can I overcome this issue?

What I Did

Initialize lists to store heights and Z values for each profile

all_heights = []
all_HCA_values = []
all_times = [] # Added to store time information

Iterate over each file and apply the QVP code

for file in files:
file_path = os.path.join(directory_path, file)

# Read radar data
radar = pyart.io.read(file_path)

# Extract desired angle
desired_angle = radar.fixed_angle['data']

# Interpolate the sounding to the radar heights
radar_T, radar_z = interpolate_sounding_to_radar(sounding, radar)

wh0 = np.where(np.isclose(np.abs(sounding['temp']),0.0,atol=1.0))
expected_ML = np.array(sounding['height'])[wh0[0]][0]/1000.

rband = 'C'                
scan_type = radar.scan_type 
use_temp = True             
verbose = False             
minRH = 0.5                 
                        
rheights = radar_z/1000.    

return_scores = False   

dz = np.ma.masked_array(
radar.fields['Z']['data'],
mask=np.logical_or(
    radar.fields['Z']['data'] < 5,
    radar.fields['Z']['data'] > 65
  )
)
dr = np.ma.masked_array(
radar.fields['ZDR']['data'],
mask=np.logical_or(
    radar.fields['ZDR']['data'] < -4,
    radar.fields['ZDR']['data'] > 6
  )
)
kd = np.ma.masked_array(radar.fields['KDP']['data'])
rh = np.ma.masked_array(
radar.fields['RHOHV']['data'],
mask=(radar.fields['RHOHV']['data'] > 0.8)
)

hcasummer = csu_fhc.csu_fhc_summer(dz=dz, zdr=dr, rho=rh, kdp=kd, use_temp=True,T=radar_T, band=rband,
                            verbose=False,use_trap=False,method='hybrid')

radar = add_field_to_radar_object(hcasummer, radar,field_name='HCA',dz_field='Z')
                        

# Create quasi-vertical profile
profile = pyart.retrieve.quasi_vertical_profile(radar, desired_angle=desired_angle, fields='HCA')

all_heights.append(profile['height'])
all_HCA_values.append(profile['HCA'])

# Extract time information (assuming 'start_time' is present in radar.metadata)
time_info = radar.metadata.get('start_time', None)

if time_info is not None:
    # Convert the time to a datetime object if needed
    # Adjust the format to include milliseconds
    time_obj = pd.to_datetime(time_info, format='%Y-%m-%d %H:%M:%S.%f')
    
    # Append time, heights, and Z values to the lists
    all_times.append(time_obj)

Convert the lists to numpy arrays for plotting

heights_array1 = np.vstack(all_heights)
HCA_values_array = np.vstack(all_HCA_values)

Choose a specific height from heights_array

selected_heights1 = heights_array1[0,: ]

Convert the time list to a NumPy array and then to matplotlib date format

times_array = date2num(np.array(all_times))

[mgrover1]

@MohitDahliya - without the data file mentioned here, I am not sure what we can do to help. You can use gatefilter to quality control your data, to see if that improves the analysis
https://arm-doe.github.io/pyart/notebooks/masking_data_with_gatefilters.html

But beyond that, we cannot make any specific recommendations

[MohitDahliya]

I have another data file for another day in which all the variables are correct. Is there some way I can compare with that file to improve my RHOHV values.?

[MohitDahliya]

If you want, I can provide you my data file...

[mgrover1]

@MohitDahliya - unfortunately this is a data quality issue that one cannot be resolved from using Py-ART. Py-ART cannot re-compute these fields.

[MohitDahliya]

No problem, I understand.

…

On Mon, Jun 17, 2024 at 7:20 AM Max Grover ***@***.***> wrote: @MohitDahliya <https://github.com/MohitDahliya> - unfortunately this is a data quality issue that one cannot be resolved from using Py-ART. Py-ART cannot re-compute these fields. — Reply to this email directly, view it on GitHub <#1593 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AQACDYWLB4OSP6O3TJZRVQ3ZH3WD7AVCNFSM6AAAAABIZ4S74GVHI2DSMVQWIX3LMV43SRDJONRXK43TNFXW4Q3PNVWWK3TUHM4TOOJVGM4DK> . You are receiving this because you were mentioned.Message ID: ***@***.***>

[scollis]

<!-- /* Font Definitions */ @font-face {font-family:"Cambria Math"; panose-1:2 4 5 3 5 4 6 3 2 4;} @font-face {font-family:Aptos; panose-1:2 11 0 4 2 2 2 2 2 4;} /* Style Definitions */ p.MsoNormal, li.MsoNormal, div.MsoNormal {margin:0in; font-size:12.0pt; font-family:"Aptos",sans-serif;} a:link, span.MsoHyperlink {mso-style-priority:99; color:blue; text-decoration:underline;} span.EmailStyle19 {mso-style-type:personal-reply; font-family:"Aptos",sans-serif; color:windowtext;} .MsoChpDefault {mso-style-type:export-only; font-size:10.0pt; mso-ligatures:none;} @page WordSection1 {size:8.5in 11.0in; margin:1.0in 1.0in 1.0in 1.0in;} div.WordSection1 {page:WordSection1;} -->If you want to discuss this issue more generally I suggest the Open Radar Science discoursehttps://openradar.discourse.group/ From: MohitDahliya ***@***.***>Date: Tuesday, June 18, 2024 at 02:08To: ARM-DOE/pyart ***@***.***>Cc: Subscribed ***@***.***>Subject: Re: [ARM-DOE/pyart] inadequate RHOHV values in the radar file (Discussion #1593)No problem, I understand. On Mon, Jun 17, 2024 at 7:20 AM Max Grover ***@***.***> wrote: > @MohitDahliya <https://github.com/MohitDahliya> - unfortunately this is a > data quality issue that one cannot be resolved from using Py-ART. Py-ART > cannot re-compute these fields. > > — > Reply to this email directly, view it on GitHub > <#1593 (comment)>, > or unsubscribe > <https://github.com/notifications/unsubscribe-auth/AQACDYWLB4OSP6O3TJZRVQ3ZH3WD7AVCNFSM6AAAAABIZ4S74GVHI2DSMVQWIX3LMV43SRDJONRXK43TNFXW4Q3PNVWWK3TUHM4TOOJVGM4DK> > . > You are receiving this because you were mentioned.Message ID: > ***@***.***> > —Reply to this email directly, view it on GitHub, or unsubscribe.You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>