grace_raster_grids.py

  • Reads in GRACE/GRACE-FO spherical harmonic coefficients and exports projected spatial fields

  • Correct spherical harmonics with the specified GIA model group

  • Filters and smooths data with specified processing algorithms [33, 64]

  • Converts data to specified units and performs a Clenshaw harmonic summation to convert to the spatial domain [26, 70]

Source code

Calling Sequence

usage: grace_raster_grids.py [-h] [--directory DIRECTORY]
                             [--output-directory OUTPUT_DIRECTORY]
                             [--file-prefix FILE_PREFIX] --center PROC
                             [--release DREL] [--product DSET] [--lmin LMIN]
                             [--lmax LMAX] [--mmax MMAX] [--start START]
                             [--end END] [--missing MISSING [MISSING ...]]
                             [--love {0,1,2,3,4}] [--reference {CF,CM,CE}]
                             [--radius RADIUS] [--destripe]
                             [--units {1,2,3,4,5}]
                             [--spacing SPACING [SPACING ...]]
                             [--bounds xmin xmax ymin ymax]
                             [--projection PROJECTION] [--gia GIA]
                             [--gia-file GIA_FILE] [--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}]
                             [--format {netCDF4,HDF5}] [--mean-file MEAN_FILE]
                             [--mean-format {ascii,netCDF4,HDF5,gfc}]
                             [--remove-file REMOVE_FILE [REMOVE_FILE ...]]
                             [--remove-format {ascii,netCDF4,HDF5,index-ascii,index-netCDF4,index-HDF5} [{ascii,netCDF4,HDF5,index-ascii,index-netCDF4,index-HDF5} ...]]
                             [--redistribute-removed] [--mask MASK] [--log]
                             [--verbose] [--mode MODE]

Named Arguments

--directory, -D

Working data directory

--output-directory, -O

Output directory for raster files

--file-prefix, -P

Prefix string for input and output files

--center, -c

GRACE/GRACE-FO Processing Center

--release, -r

GRACE/GRACE-FO Data Release

--product, -p

GRACE/GRACE-FO Level-2 data product

--lmin

Minimum spherical harmonic degree

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

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

Use decorrelation (destriping) filter

--units, -U

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

Output units

  • 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

Output grid spacing

--bounds

Output grid extents

--projection

Spatial projection as EPSG code or PROJ4 string

  • 4326: latitude and longitude coordinates on WGS84 reference ellipsoid

--gia, -G

Possible choices: IJ05-R2, W12a, SM09, ICE6G, Wu10, AW13-ICE6G, AW13-IJ05, Caron, ICE6G-D, ascii, netCDF4, HDF5

GIA model type to read

  • 'IJ05-R2': Ivins R2 GIA Models [29]

  • 'W12a': Whitehouse GIA Models [82]

  • 'SM09': Simpson/Milne GIA Models [57]

  • 'ICE6G': ICE-6G GIA Models [48]

  • 'Wu10': Wu (2010) GIA Correction [83]

  • 'AW13-ICE6G': Geruo A ICE-6G GIA Models [1]

  • 'AW13-IJ05': Geruo A IJ05-R2 GIA Models [1]

  • 'Caron': Caron JPL GIA Assimilation [7]

  • 'ICE6G-D': ICE-6G Version-D GIA Models [49]

  • 'ascii': reformatted GIA in ascii format

  • 'netCDF4': reformatted GIA in netCDF4 format

  • 'HDF5': reformatted GIA in HDF5 format

--gia-file

GIA file to read

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

--format, -F

Possible choices: netCDF4, HDF5

Input/output data format

--mean-file

GRACE/GRACE-FO mean file to remove from the harmonic data

--mean-format

Possible choices: ascii, netCDF4, HDF5, gfc

Input data format for GRACE/GRACE-FO mean file

--remove-file

Monthly files to be removed from the GRACE/GRACE-FO data

--remove-format

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

Input data format for files to be removed

--redistribute-removed

Redistribute removed mass fields over the ocean

--mask

Land-sea mask for redistributing land water flux

--log

Output log file for each job

--verbose, -V

Verbose output of run

--mode, -M

Permissions mode of output files