#!/usr/bin/env python
# coding: utf-8

# # 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')


# ## Vertical Section plot with Matplotlib

# In[6]:


ds['so'].plot()


# ## 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)


# In[ ]: