Plots maps of the real and imaginary geocentric pole tide admittance functions from Desai et al. (2002)
io.ocean_pole_tide.py
: Read ocean pole load tide map from IERSutilities.py
: download and management utilities for filesimport matplotlib.pyplot as plt
import cartopy.crs as ccrs
import pyTMD.io
import pyTMD.utilities
# read ocean pole tide map from Desai (2002)
ocean_pole_tide_file = pyTMD.utilities.get_data_path(['data','opoleloadcoefcmcor.txt.gz'])
iur,iun,iue,ilon,ilat = pyTMD.io.ocean_pole_tide(ocean_pole_tide_file)
fig,(ax1,ax2) = plt.subplots(ncols=2,sharex=True,sharey=True,figsize=(10,4),
subplot_kw=dict(projection=ccrs.PlateCarree()))
ax1.imshow(iur.real.T,extent=(ilon[0],ilon[-1],ilat[0],ilat[-1]),origin='lower')
ax2.imshow(iur.imag.T,extent=(ilon[0],ilon[-1],ilat[0],ilat[-1]),origin='lower')
ax1.set_title('Radial Ocean Pole Tide (real component)')
ax2.set_title('Radial Ocean Pole Tide (imaginary component)')
ax1.coastlines(); ax2.coastlines()
fig.subplots_adjust(left=0.01, right=0.99, bottom=0.10, top=0.95, wspace=0.05)
plt.show()