Rainfall accumulation

I thought I posted here a few weeks ago, but couldn’t find a post from me when I went to look back for it. Apologies if I’ve missed it.

I have a problem I’m trying to help a colleague with and am comimg up short. He has two weeks’ worth of radar data from a site in Hawaii. From this data, he needs to know the total rainfall at a certain lat/lon within that two week period. I know Xarray is involved, but I’m still stuck on how to approach this.


Thanks for asking this question @davidmoranccm - I would start by checking out the Xradar documentation

This package will enable you to read in the data (what format is the dataset?)

Once you read in the data, you will either need to use one of the existing rainfall fields, or estimate the rainfall rate from reflectivity using some relationship.

If you are interested in using Py-ART for this, an example of calculating Quantitative Precipitation Estimation (QPE) from radar data was covered at the recent ERAD workshop


Thanks! I will check these out!

OK, I feel stupid for asking, but I started installing Pyart and followed the instructions here: The Python ARM Radar Toolkit - Py-ART — Py-ART 1.14.1.post1 documentation. I grabbed a script off the Pyart website and when I ran it, I get:

Traceback (most recent call last):
File “C:\Users\David\Desktop\plot_nexrad_data_aws.py”, line 18, in
import pyart
ModuleNotFoundError: No module named ‘pyart’

Where could I have gone wrong?


There are no stupid questions! Did you activate your conda environment? What are you using to manage your environments?

I don’t think so; how do I do that? I’m using Anaconda.

I really like the description/videos from the NCAR Earth System Data Science site.


I would follow those videos/commands. Once you activate your environment

conda activate my-environment

You will need to install Py-ART

conda install -c conda-forge arm_pyart

Then, you can run

python -m my_script.py

The conda cheatsheet might help here too.