combine_harmonics.py

  • Converts a file from the spherical harmonic domain into the spatial domain [73]

Source code

Calling Sequence

usage: combine_harmonics.py [-h] [--lmax LMAX] [--mmax MMAX]
                            [--love {0,1,2,3,4}] [--reference {CF,CM,CE}]
                            [--radius RADIUS] [--destripe]
                            [--units {0,1,2,3,4,5}]
                            [--spacing dlon [dlat ...]] [--interval {1,2,3}]
                            [--bounds lon_min lon_max lat_min lat_max]
                            [--redistribute-mass] [--mask MASK] [--mean MEAN]
                            [--format {ascii,netCDF4,HDF5,index-ascii,index-netCDF4,index-HDF5}]
                            [--date] [--verbose] [--mode MODE]
                            [infile] [outfile]

Positional Arguments

infile

Input harmonic file

outfile

Output spatial file

Named Arguments

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

--destripe, -d

Verbose output of run

--units, -U

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

Output units

  • 0: no unit conversion

  • 1: cm of water thickness

  • 2: mm of geoid height

  • 3: mm of elastic crustal deformation

  • 4: μGal gravitational perturbation

  • 5: mbar equivalent surface pressure

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

  • 1: (0:360, 90:-90)

  • 2: (degree spacing/2)

  • 3: non-global grid (set with defined bounds)

--bounds, -B

Bounding box for non-global grid

--redistribute-mass

Redistribute total mass over the ocean

--mask

Land-sea mask for redistributing over the ocean

--mean

Mean file to remove from the harmonic data

--format, -F

Possible choices: ascii, netCDF4, HDF5, index-ascii, index-netCDF4, index-HDF5

Input and output data format

--date, -D

Input and output files have date information

--verbose, -V

Verbose output of run

--mode, -M

Permissions mode of output files