wqp.get_results with multiple siteid(s)

Created by: jsta

I am trying to pull WQP data from multiple sites as follows:

from dataretrieval import nwis,wqp

sites, md = wqp.get_results(siteid=["USGS-04024315", "USGS-01594440"], characteristicName = "pH")

but the result is empty:

sites.shape

(0, 63)

Maybe a list is not the correct data structure to pass to siteid? Or, perhaps this functionality not exist?

This is possible with the dataRetrieval R package which creates the following urlCall with this query:

"https://www.waterqualitydata.us/data/Result/search?siteid=USGS-04024315%3BUSGS-01594440&characteristicName=pH&zip=yes&mimeType=tsv"