I was wondering if anyone knows about a good way to access data from the World Ocean Database (CTD and XBT casts for example) or if you know of any tutorials for doing this optimally on a cluster.
From Thomas Moore in the OHW slack channel:
fs = s3fs.S3FileSystem(anon=True)
fs.ls('s3://noaa-wod-pds/')
f = fs.open('s3://noaa-wod-f2 = fs.open('s3://noaa-wod-pds/2022/wod_pfl_2022.nc')
wod_pfl_2022 = xr.open_dataset(f2, engine='h5netcdf')
plt.rcParams["figure.figsize"] = (20,10)
plt.scatter(x=wod_pfl_2022.lon, y=wod_pfl_2022.lat, linewidths=0, s=0.5)
Not sure if you find everything you need there.
2 Likes
Ha ha! I navigated to this thread curious to see what the answers might be.
1 Like