On making CFAD over a specific area/domain

[gene_cgha]

I’m currently working on generating a CFAD (Contoured Frequency by Altitude Diagram) using Py-ART over a specific geographic domain, similar to this visualization of Chen et al. (2019) in Figure 6. I’ve already written an initial Python script that attempts to process and visualize the CFAD, but the output doesn’t appear to be correct or meaningful.

May I kindly ask if someone could provide a working example or correct code that shows how to properly ‘subset’ radar data (CFRadial) over a defined lat-lon domain to generate a CFAD?

I’m happy to share my code snippet and sample data for this case.

Code: Jupyter Notebook

Sample Data: CFRadial Data

[mgrover1]

Hello!

Any reason for subsetting over a lat/lon area as opposed to subsetting using the native radial coordiantes ie range, azimuth?

One suggestion would be to feed gatefilters in using the desired region, similar to this example

Fun With Gatefilters - ARM Open Science in the Forest Summer School

Then feed that into the CFAD example in Py-ART

Creating a CFAD diagram — Py-ART 2.0.1 documentation