GRACE gravity field solutions¶

[1]:

# Imports
from gfeatpy.gravity import GravityField, EquivalentWaterHeight

from gfeatpy import plotting
import matplotlib.pyplot as plt

data_root = "../../../data/"

[2]:

# Load static gravity field
static = GravityField(90).load(f"{data_root}/gravity/static/ITSG-Grace2018s.gfc")

Non-resonant case¶

[3]:

# Load monthly static gravity field
monthly = GravityField(90).load(f"{data_root}/gravity/monthly/ITSG-Grace2018_n96_2007-03.gfc")
# Remove static field
monthly.coefficients = monthly.coefficients - static.coefficients

[4]:

# Plot EWH
plotting.synthesis(monthly, 360, 180, EquivalentWaterHeight(0))
plt.show()

../../_images/notebooks_gravity-field-data_monthly_5_0.png

[5]:

# Plot EWH + Gaussian spatial smoothing filter r=750 km
plotting.synthesis(monthly, 360, 180, EquivalentWaterHeight(750e3))
plt.show()

../../_images/notebooks_gravity-field-data_monthly_6_0.png

Resonant case¶

[6]:

# Load monthly static gravity field
monthly = GravityField(90).load(f"{data_root}/gravity/monthly/ITSG-Grace2018_n96_2004-09.gfc")
# Remove static field
monthly.coefficients = monthly.coefficients - static.coefficients

[7]:

# Plot EWH
plotting.synthesis(monthly, 360, 180, EquivalentWaterHeight(0))
plt.show()

../../_images/notebooks_gravity-field-data_monthly_9_0.png

[8]:

# Plot EWH + Gaussian spatial smoothing filter r=750 km
plotting.synthesis(monthly, 360, 180, EquivalentWaterHeight(750e3))
plt.show()

../../_images/notebooks_gravity-field-data_monthly_10_0.png