grace_mean_harmonics.py

  • Calculates the temporal mean of the GRACE/GRACE-FO spherical harmonics for a specified date range

  • Used to estimate the static gravitational field over a given date rage

Source code

Calling Sequence

usage: grace_mean_harmonics.py [-h] [--directory DIRECTORY] --center PROC
                               [--release DREL] [--product DSET] [--lmax LMAX]
                               [--mmax MMAX] [--start START] [--end END]
                               [--missing MISSING [MISSING ...]]
                               [--atm-correction] [--pole-tide]
                               [--geocenter DEG1]
                               [--geocenter-file GEOCENTER_FILE]
                               [--interpolate-geocenter]
                               [--slr-c20 {CSR,GFZ,GSFC}]
                               [--slr-21 {CSR,GFZ,GSFC}] [--slr-22 {CSR,GSFC}]
                               [--slr-c30 {CSR,GFZ,GSFC,LARES}]
                               [--slr-c40 {CSR,GSFC,LARES}]
                               [--slr-c50 {CSR,GSFC,LARES}]
                               [--mean-file MEAN_FILE]
                               [--mean-format {ascii,netCDF4,HDF5,gfc}]
                               [--log] [--verbose] [--mode MODE]

Named Arguments

--directory, -D

Working data directory

--center, -c

GRACE/GRACE-FO Processing Center

--release, -r

GRACE/GRACE-FO Data Release

--product, -p

GRACE/GRACE-FO Level-2 data product

--lmax, -l

Maximum spherical harmonic degree

--mmax, -m

Maximum spherical harmonic order

--start, -S

Starting GRACE/GRACE-FO month

--end, -E

Ending GRACE/GRACE-FO month

--missing, -N

Missing GRACE/GRACE-FO months

--atm-correction

Apply atmospheric jump correction coefficients

--pole-tide

Correct for pole tide drift

--geocenter

Possible choices: Tellus, SLR, SLF, UCI, Swenson, GFZ

Update Degree 1 coefficients with SLR or derived values

  • None

  • 'Tellus': GRACE/GRACE-FO TN-13 coefficients from PO.DAAC

  • 'SLR': satellite laser ranging coefficients from CSR

  • 'UCI': Sutterley and Velicogna coefficients, Remote Sensing (2019)

  • 'Swenson': GRACE-derived coefficients from Sean Swenson

  • 'GFZ': GRACE/SLR derived coefficients from GFZ GravIS

--geocenter-file

Specific geocenter file if not default

--interpolate-geocenter

Least-squares model missing Degree 1 coefficients

--slr-c20

Replace C20 coefficients with SLR values

  • None: use original values

  • 'CSR': use values from CSR (TN-07, TN-09, TN-11)

  • 'GFZ': use values from GFZ

  • 'GSFC': use values from GSFC (TN-14)

--slr-21

Possible choices: CSR, GFZ, GSFC

Replace C21 and S21 coefficients with SLR values

--slr-22

Replace C22 and S22 coefficients with SLR values

  • None: use original values

  • 'CSR': use values from CSR

  • 'GSFC': use values from GSFC

--slr-c30

Replace C30 coefficients with SLR values

  • None: use original values

  • 'CSR': use values from CSR (5×5 with 6,1)

  • 'GFZ': use values from GFZ GravIS

  • 'GSFC': use values from GSFC (TN-14)

  • 'LARES': use filtered values from CSR

--slr-c40

Replace C40 coefficients with SLR values

  • None: use original values

  • 'CSR': use values from CSR (5×5 with 6,1)

  • 'GSFC': use values from GSFC

  • 'LARES': use filtered values from CSR

--slr-c50

Replace C50 coefficients with SLR values

  • None: use original values

  • 'CSR': use values from CSR (5×5 with 6,1)

  • 'GSFC': use values from GSFC

  • 'LARES': use filtered values from CSR

--mean-file

Output GRACE/GRACE-FO mean file

--mean-format

Possible choices: ascii, netCDF4, HDF5, gfc

Output data format for GRACE/GRACE-FO mean file

--log

Output log file for each job

--verbose, -V

Verbose output of run

--mode, -M

Permissions mode of output files