def obtainValues():
url = 'http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.11&request=ProcessCoverages&query=for c in (river_discharge_forecast_opt2) return encode (c[ansi("2008-04-08T00:00"),Lat(41.32),Long(-89.0)],"csv")'
print(" hardcoded value for url "+ url)
r = requests.get(url,
proxies={'http':None}
)
r.raise_for_status()
x = np.array(eval(r.text.replace('{','[').replace('}',']')))
listOfMedian=[]
for i in range(30):
listOfEnsambles=x[i]
## print(listOfEnsambles)
listOfEnsambles_masked = np.ma.masked_where(listOfEnsambles == 0, listOfEnsambles)
## print(listOfEnsambles_masked)
medianValue=np.nanmedian(listOfEnsambles_masked)
## print(medianValue)
listOfMedian.append(medianValue)
return(listOfMedian)
评论列表
文章目录