calc_sensitivity_kernel.py

  • Calculates spatial sensitivity kernels through a least-squares mascon procedure following Jacob et al. [31], Tiwari et al. [67]

Source code

Calling Sequence

usage: calc_sensitivity_kernel.py [-h] [--output-directory OUTPUT_DIRECTORY]
                                  [--lmin LMIN] [--lmax LMAX] [--mmax MMAX]
                                  [--love {0,1,2,3,4}]
                                  [--reference {CF,CM,CE}] [--radius RADIUS]
                                  [--format {ascii,netCDF4,HDF5}]
                                  [--mascon-file MASCON_FILE]
                                  [--redistribute-mascons]
                                  [--fit-method {1,2}]
                                  [--solver {inv,lstsq,gelsd,gelsy,gelss}]
                                  [--mask MASK] [--spatial]
                                  [--spacing dlon [dlat ...]]
                                  [--interval {1,2,3}]
                                  [--bounds lon_min lon_max lat_min lat_max]
                                  [--log] [--verbose] [--mode MODE]

Named Arguments

--output-directory, -O

Output directory for mascon files

--lmin

Minimum spherical harmonic degree

--lmax, -l

Maximum spherical harmonic degree

--mmax, -m

Maximum spherical harmonic order

--love, -n

Possible choices: 0, 1, 2, 3, 4

Treatment of the Load Love numbers

  • 0: Han and Wahr (1995) values from PREM [23]

  • 1: Gegout (2005) values from PREM [21]

  • 2: Wang et al. (2012) values from PREM [80]

  • 3: Wang et al. (2012) values from PREM with hard sediment [80]

  • 4: Wang et al. (2012) values from PREM with soft sediment [80]

--reference

Possible choices: CF, CM, CE

Reference frame for load Love numbers

  • 'CF': Center of Surface Figure

  • 'CM': Center of Mass of Earth System

  • 'CE': Center of Mass of Solid Earth

--radius, -R

Gaussian smoothing radius (km)

--format, -F

Possible choices: ascii, netCDF4, HDF5

Input data format for auxiliary files

--mascon-file

Index file of mascons spherical harmonics

--redistribute-mascons

Redistribute mascon mass over the ocean

--fit-method

Possible choices: 1, 2

Method for fitting sensitivity kernel to harmonics

  • 1: mass coefficients

  • 2: geoid coefficients

--solver, -s

Least squares solver for sensitivity kernels

  • 'inv': matrix inversion

  • 'lstsq': least squares solution

  • 'gelsy': complete orthogonal factorization solution

  • 'gelss': singular value decomposition (SVD) solution

  • 'gelsd': singular value decomposition (SVD) solution with a divide and conquer method

--mask

Land-sea mask for redistributing mascon mass

--spatial, -o

Output spatial grid file for each mascon

--spacing, -S

Spatial resolution of output data

--interval, -I

Possible choices: 1, 2, 3

Output grid interval (1: global, 2: centered global, 3: non-global)

--bounds, -B

Bounding box for non-global grid

--log

Output log file for each job

--verbose, -V

Verbose output of processing run

--mode, -M

Permissions mode of output files