GHCN Daily Precip Data

[jdossgollin]

• Add the GHCN data to the Snakefile
• Should download all stations and the index file
• Individual users should write their own scripts to extract what they need
• There is no way to incrementally re-run this so will need to think about how to update as needed (every year?)

[jdossgollin]

Steps:

1. Download https://www.ncei.noaa.gov/data/global-historical-climatology-network-daily/archive/
2. Unzip and move to correct location

[jdossgollin]

If we were ambitious, we could turn into a netcdf file (indexed by station, time) -- perhaps 1 netcdf4 file per station?

for each station, have the following data variables:

• precip_in (indexed by time)
• quality_flag (indexed by time)
• station ID (not indexed)
• lon (not indexed)
• lat (not indexed)

[jdossgollin]

Then provide an example of how to use xarray and dask (xr.open_mfdataset) to read in precip data from multiple stations at once

[jdossgollin]

Also note https://www.ncei.noaa.gov/data/global-historical-climatology-network-daily/doc/ghcnd-stations.txt