gemb_smb_harmonics.py

Reads GEMB surface mass balance variables and converts to spherical harmonic coefficients [29]
Estimates the grid values as point masses for calculating the gravitational spherical harmonic coefficients

Calling Sequence

usage: gemb_smb_harmonics.py [-h] [--product {accum_SMB,dFAC}]
                             [--mask MASK [MASK ...]] [--area AREA]
                             [--lmax LMAX] [--mmax MMAX] [--love {0,1,2,3,4}]
                             [--reference {CF,CM,CE}]
                             [--format {ascii,netCDF4,HDF5}] [--verbose]
                             [--mode MODE]
                             infile

Positional Arguments

infile: GEMB SMB file to run

Named Arguments

--product, -P

Possible choices: accum_SMB, dFAC

GEMB product to calculate

--mask

netCDF4 masks file for reducing to regions

--area

netCDF4 area file for calculating mass

--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 [8]
1: Gegout (2005) values from PREM [6]
2: Wang et al. (2012) values from PREM [30]
3: Wang et al. (2012) values from PREM with hard sediment [30]
4: Wang et al. (2012) values from PREM with soft sediment [30]

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

--format, -F

Possible choices: ascii, netCDF4, HDF5

Input and output data format

--verbose, -V

Verbose output of processing run

--mode, -M

Permission mode of directories and files