This page was generated from examples/ssebop.ipynb. Interactive online version:

SSEBop

import matplotlib.pyplot as plt
from pynhd import NLDI

import pygeohydro as gh

The daily actual evapotranspiration can be retrieved from SEEBop database. Note that since this service does not offer a web service and data are available as raster files on the server, so this function is not as fast as other functions and download speed might be the bottleneck.

You can get the actual ET for location using ssebopeta_byloc and for a region using ssebopeta_bygeom. Let’s get a watershed geometry using NLDI and then get the actual ET.

geometry = NLDI().get_basins("01031500").geometry[0]

dates = ("2005-10-01", "2005-10-05")
coords = (geometry.centroid.x, geometry.centroid.y)
eta_p = gh.ssebopeta_byloc(coords, dates=dates)
eta_g = gh.ssebopeta_bygeom(geometry, dates=dates)

ax = eta_g.isel(time=4).plot(size=5)
ax.figure.savefig("../_static/eta.png", bbox_inches="tight", dpi=100)