Extracting data

[davidmoranccm]

I’m a meteorologist who’s been programming in Python for a little over a year, but only recently started looking at PyART. I have a problem I’m not sure if I can use with PyART and wanted to ask.

I have two weeks’ worth of radar data from a NEXRAD site. I need to know the total rainfall at a certain lat/lon over the two week period. Can this be done with PyART?

Thanks!
David

[zssherman]

Hi @davidmoranccm, you could always try to do a a loop with pyarts nexrad reader, and use one of the estimate rain rate functions in pyarts retrieval module. For extracting at a specific lat and lon, I think the column sect might work after the rain rate retrieval, if not, once things get ported to xarray the .sel will help.

[davidmoranccm]

Thanks, Zach! I’ll play around with it later. David David Moran, CCM AMS Certified Consulting Meteorologist E-Mail: ***@***.*** Phone: (405) 693-5632

…

On Nov 4, 2022, at 6:17 PM, Zach Sherman ***@***.***> wrote:  Hi @davidmoranccm, you could always try to do a a loop with pyarts nexrad reader, and use one of the estimate rain rate functions in pyarts retrieval module. For extracting at a specific lat and lon, I think the column sect might work, if not, once things get ported to xarray the .sel will help. — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you were mentioned.

[wxwolff]

Just a thought if you are familiar with Xarray… For each volume scan, use the radar.to_grid method and then grid.to_xarray. Create an empty lists, gr_list, and append each grid. Use xarray.concat(gr_list) to create an Xarray dataset From there, you can use Xarray’s powerful tools to calculate aggregate statistics. D

_____________________________________ David B. Wolff, Code 612 GPM Deputy Project Scientist for GV GPM Ground Validation System Manager NASA Wallops Flight Facility 34200 Fulton Street Building N159, Room E220 Wallops Island, VA 23337 Ph: 757-824-1515; Cell: 443-262-1146 From: davidmoranccm ***@***.***> Reply-To: ARM-DOE/pyart ***@***.***> Date: Friday, November 4, 2022 at 7:22 PM To: ARM-DOE/pyart ***@***.***> Cc: Subscribed ***@***.***> Subject: [EXTERNAL] Re: [ARM-DOE/pyart] Extracting data (Discussion #1316) Thanks, Zach! I’ll play around with it later. David David Moran, CCM AMS Certified Consulting Meteorologist E-Mail: ***@***.*** Phone: (405) 693-5632

On Nov 4, 2022, at 6:17 PM, Zach Sherman ***@***.***> wrote: Hi @davidmoranccm, you could always try to do a a loop with pyarts nexrad reader, and use one of the estimate rain rate functions in pyarts retrieval module. For extracting at a specific lat and lon, I think the column sect might work, if not, once things get ported to xarray the .sel will help. — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you were mentioned.

— Reply to this email directly, view it on GitHub<https://gcc02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FARM-DOE%2Fpyart%2Fdiscussions%2F1316%23discussioncomment-4061271&data=05%7C01%7CDavid.B.Wolff%40nasa.gov%7Cedc795c9450449fa0a5308dabebb62ee%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C638032009270185695%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=8KOdRh2MybRfogH4isnEZUhNCb%2BF7LVVsvA%2BFoHy9jw%3D&reserved=0>, or unsubscribe<https://gcc02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FADJRUYSEZACSOETD4HN3A6LWGWLBTANCNFSM6AAAAAARXT2QHU&data=05%7C01%7CDavid.B.Wolff%40nasa.gov%7Cedc795c9450449fa0a5308dabebb62ee%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C638032009270185695%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=itu9o8NrkTP3ZD%2FnU8EMoJbVZDzTK2wW6%2FYXzF9GGXM%3D&reserved=0>. You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

[scollis]

Also for general QPE questions highly suggest joining https://openradar.discourse.group/