Any nice way to access XBT and CTD casts from WOD?

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.

3 Likes

Ha ha! I navigated to this thread curious to see what the answers might be.

1 Like