Help with RHI plotting

[Arunas Kuciauskas]

Hi:

I'm new here. I have a volumetric radar file that I would like to plot RHIs.

From viewing the Pyart RHI examples in: 

https://arm-doe.github.io/pyart/source/auto_examples/index.html#plotting-examples,

the data formata are only stored in MDV or Sigmet.  The file that I have is in h5 and can be read using odimh5.

I have attached the code and the resultant output which is not what I expected.

The data file is in the 2nd attachment

I want  plot along the 135 degree azimuth angle.

I appreciate any help.

Thanks,

Arunas

[Zach Sherman]

Hi Arunas,

Your file is a ppi format file, from the plot I created and from what the scan_type says.

Zach S.

[Zach Sherman]

Use for getting an rhi for the specific azimuth using your ppi radar:
rhi = pyart.util.cross_section_ppi(radar, [145.0])

fig = plt.figure(figsize=(6,6))

display = pyart.graph.RadarMapDisplay(rhi)

display.plot_rhi('reflectivity', vmin=-20, vmax=20,

                 fig=fig)

[Zach Sherman]

Just use 135 instead as I realized I used 145...

[Arunas Kuciauskas]

Zach:

Thanks very much for your help and insight. I will apply your code and let you know

of my results.

Arunas

[Arunas Kuciauskas]

Is there a way of setting the y axis limits (distance above radar)?  I looked at 

pyart.graph.RadarMapDisplay.set_limits(ylim=40)

but that was unsuccessful.

Arunas

[Zach Sherman]

Yeah, to set the limits just use:

plt.xlim(0, 20)

plt.ylim(0, 50)

for setting the limits.

Zach S.

[Arunas Kuciauskas]

Hi: 

I'm not sure if I am posting this properly...

I have a new request. I was able to plot RHI's from the PPI file that I provided you, thanks!

My RHI plot is in the first  attachment.

I would like to plot an RHI similar to the Py-ART RHI figure in:http://arm-doe.github.io/pyart/

(see 2nd attachment).

The corresponding code is in the 3rd attachment.

Arunas

[Zach Sherman]

Arunas,

Do you mean same colormap?The tick labels should mostly be the same, the first image they are hard to see, did you save the first image with plt.savefig? I don't see it in the first code.

Zach S.

[Arunas Kuciauskas]

No, No, No...The plotting format, where I want the vertical cross section to have radar echoes extending to the top of the atmosphere throughout the entire scope, from behind the radar to ahead of the radar

[Zach Sherman]

Scott or Bobby, you might have better input on that, my guess would be, being a PPI your not getting those higher elevation angles, like the RHI in image 2

Zach S.

[Kai Mühlbauer]

Hi Zach, Arunas,

Zach, correct, this particular Cross-section plot is only possible with real RHI data. From volume PPI the pseudo RHI has the so called cone of silence where no measurement was taken.

Think about your highest elevation in the PPI is say 30 degrees, then you are blind the complete area from 30 degree on the one side moving up to and over 90 degree (vertically pointing) and again down to 30 degrees.

HTH,
Kai

Am 15. Februar 2020 00:33:20 MEZ schrieb Zach Sherman <sherma...@gmail.com>:

Scott or Bobby, you might have better input on that, my guess would be, being a PPI your not getting those higher elevation angles, like the RHI in image 2

Zach S.

--
Diese Nachricht wurde von meinem Android-Gerät mit K-9 Mail gesendet.

[Arunas Kuciauskas]

Kai and Zach:

Thanks for your analysis of the dataset. My goal was to capture the height of some of these fires that turned into PyroCb’s (fire thunderstorms) that can reach the lower stratosphere. There’s a lot of current research going on during some of the most intensive fire days along the SE coast of Australia, during 28-31 December and the 4th of January. Looks like my radar dataset will not be able to capture these echo tops. Oh well, but thank you very much,

Arunas 

Sent from my iPhone

On Feb 15, 2020, at 4:28 AM, Kai Mühlbauer <kai.mue...@uni-bonn.de> wrote:

Hi Zach, Arunas,

--
You received this message because you are subscribed to the Google Groups "Py-ART Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pyart-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/pyart-users/AEE2F203-8BE7-4655-994E-E0D025944009%40uni-bonn.de.

[Zach Sherman]

Thanks Kai! Much appreciated!

Zach S.