Support and Examples for Using USGS Stations

[SorooshMani-NOAA]

@Delta-function, @Rjialceky during our last observation data meeting we discussed the interest in using searvey for retrieving USGS data as well as the lack of example for getting station data based on the station ID. We'd like to work with you to add such example and also support the use of searvey in your workflow, so we were wondering what sort of input and output data do you have in mind for your use-case, so that we can work on providing an example of achieving that.

Searvey USGS module extracts instantaneous value data for stations from NWIS API. Is it sufficient for example to demonstrate how a given list of station IDs (in the form of list of strings input) can be used to retrieve an xarray.Dataset object filled with station data? Or do you have anything else in mind?

Thank you in advance for your feedback!

[Delta-function]

Soroosh, to illustrate my use-case, I uploaded my Matlab script that does what I need: gets station data based on a list of stations IDs, and saves the data in netcdf files. Please see this scenario here:
https://github.com/noaa-ocs-modeling/VDatum-Spatially-Varying-Uncertainty/tree/main/utilities
If you have Matlab, you can play this use-case by copying both files (*.m and *.txt), creating a dir USGSdata for the output files, and then running *.m.
Thanks!

[SorooshMani-NOAA]

Thank you @Delta-function, I'll play with the Matlab code and write something with searvey to get similar results.

[SorooshMani-NOAA]

@Delta-function please see https://github.com/SorooshMani-NOAA/searvey/blob/enhance/usgs_examples/examples/USGS_by_id.ipynb for the example of using searvey to download data and format and save it as you do in the Matlab code.