oceansdb icon indicating copy to clipboard operation
oceansdb copied to clipboard

how to extract annual data?

Open uliw opened this issue 3 years ago • 0 comments

I tried

db = oceansdb.WOA()
db.tscale="annual"
o = db["dissolved_oxygen"].extract(
  var="mean", depth=depth, lat=latitude, lon=longitude 
)

which complaints about the missing doy. Adding a value for `doy' then results in this error

File "oxygen.py", line 21, in <module>
    o = db["dissolved_oxygen"].extract(
  File "/home/uliw/anaconda3/lib/python3.8/site-packages/oceansdb/woa.py", line 386, in extract
    output = self.interpolate(doy, depth, lat, lon, var)
  File "/home/uliw/anaconda3/lib/python3.8/site-packages/oceansdb/woa.py", line 269, in interpolate
    f = interp1d(dims['time'], subset[v], axis=0)
  File "/home/uliw/anaconda3/lib/python3.8/site-packages/scipy/interpolate/interpolate.py", line 554, in __init__
    raise ValueError("x and y arrays must have at "
ValueError: x and y arrays must have at least 2 entries

uliw avatar Mar 04 '21 15:03 uliw