Skip to content

Walling_documentation

Walling9 edited this page Dec 7, 2022 · 3 revisions

Codes

GRACEpy

Description:

This script contains multiple functions listed below with their descriptions:

Modules used

GRACEconstants


sc2cs

Description:

The function converts the rectangular (L+1)x(2L+1) matrix field, containing spherical harmonics coefficients in /S|C\ storage format into a square (L+1)x(L+1) matrix in |C\S| format. (inverse of cs2sc)

Inputs

field     the rectangular (L+1)x(2L+1) matrix FIELD, containing spherical harmonics coefficients in /S|C\ storage format

Outputs

cs     square (L+1)x(L+1) matrix in |C\S| format

cs2sc

Description:

The function converts the storage format into a square (L+1)x(L+1) matrix in |C\S| format to rectangular (L+1)x(2L+1) matrix field, containing spherical harmonics coefficients in /S|C\ . (inverse of sc2cs)

Inputs

field     the square (L+1)x(L+1) matrix FIELD , containing spherical harmonics coefficients in |C\S| storage format

Outputs

cs     rectangular (L+1)x(2L+1) matrix in /S|C\format

upwcon

Description:

This function returns the upward continuation terms, which are the field values at a height above the earth's radius.

Inputs

degree  the spherical harmonic degree in a scalar or vector format, integer
height  the height above the mean Earth radius in [m] in a scalar or vector format

Outputs

uc   the upward continuation terms

lovenr

Description:

returns the interpolated Love numbers of the Earth up to degree 200. Love numbers are a representation of the rigidity of the Earth.

Inputs

lmax     the maximum degree

Outputs

kn     the Love number k

lovenrPREM

Description:

for a degree lmax, it computes the Love and Shida number (l) of the Earth for different reference frames (provided by Olivier Francis from the PREM Earth model for selected degrees)

Inputs

lmax     the maximum degree
frame    a string referring to the frame of reference

Outputs

kn    the Love number k
hn    the Love number h
ln    the Love number l

eigengrav

Description:

This script contains multiple functions listed below with their descriptions:

Description:

The function yields the eigenvalues (isotropic spectral transfer) for switching between the following gravity related quantities:
• gravitational potential [m2/s 2 ]
• geoid height [m]
• gravity disturbances [mGal]
• gravity anomaly [mGal]
• 2nd radial derivative of the field [E]
• slope [arcsec]
• equivalent water thickness [m]
• surface mass density
• coefficient dependent

Inputs

lmax     the maximum degree
fstr     string for the relevant gravity related quantity with the following parameters
   none        least squares
   geoid      weighted least squares
   dg or gravity   gravity anomaly, potential
   tr          gravity disturbance
   trr          (d^2/dr^2)
   slope       block mean values
   water       equivalent water thickness
   smd       surface mass density
   height       vertical displacements

Outputs

tf      the transfer function

Modules used

GRACEconstants, GRACEpy