Explore HadGEM3 CMIP6 data¶

In [1]:

import fsspec
import xarray as xr
import hvplot.xarray
import hvplot.pandas

In [2]:

fs = fsspec.filesystem('s3', anon=True, client_kwargs={'endpoint_url':'https://mghp.osn.xsede.org'})

In [3]:

ds = xr.open_dataset(fs.open('s3://rsignellbucket1/testing/test.nc'))

In [4]:

ds['so']

Plot surface salinity values¶

In [5]:

ds.isel(lev=0).to_dataframe().hvplot.points(x='longitude', y='latitude', color='so', cmap='turbo',  geo=True, tiles='OSM')

/home/conda/users/222098f819851fd4ec9d4c9445d950452119627537a69aa8cd965d7dbb1ef7e2-20230602-144731-342097-199-pangeo/lib/python3.10/site-packages/geoviews/operation/__init__.py:14: HoloviewsDeprecationWarning: 'ResamplingOperation' is deprecated and will be removed in version 1.17, use 'ResampleOperation2D' instead.
  from holoviews.operation.datashader import (

Out[5]:

Vertical Section plot with Matplotlib¶

In [6]:

ds['so'].plot()

Out[6]:

<matplotlib.collections.QuadMesh at 0x7f2fdd7274c0>

Vertical Section plot with hvplot !¶

In [7]:

ds["so"] = ds["so"].reindex(trajectory=ds["so"]["trajectory"])

In [8]:

ds['so'].hvplot.quadmesh(cmap='viridis', frame_width=400, frame_height=400)

Out[8]:

In [ ]: