piecewise_grace_maps.py

Reads in GRACE/GRACE-FO spatial files and fits a piecewise regression model at each grid point for breakpoint analysis

Calling Sequence

usage: piecewise_grace_maps.py [-h] [--output-directory OUTPUT_DIRECTORY]
                               [--file-prefix FILE_PREFIX] [--lmax LMAX]
                               [--mmax MMAX] [--start START] [--end END]
                               [--breakpoint BREAKPOINT]
                               [--missing MISSING [MISSING ...]]
                               [--radius RADIUS] [--destripe]
                               [--units {1,2,3,4,5}]
                               [--spacing dlon [dlat ...]]
                               [--interval {1,2,3}]
                               [--bounds lon_min lon_max lat_min lat_max]
                               [--format {ascii,netCDF4,HDF5}]
                               [--redistribute-removed]
                               [--cycles CYCLES [CYCLES ...]] [--log]
                               [--verbose] [--mode MODE]

Named Arguments

--output-directory, -O

Output directory for spatial files

--file-prefix, -P

Prefix string for input and output files

--lmax, -l

Maximum spherical harmonic degree

--mmax, -m

Maximum spherical harmonic order

--start, -S

Starting GRACE/GRACE-FO month for time series regression

--end, -E

Ending GRACE/GRACE-FO month for time series regression

--breakpoint, -B

Breakpoint GRACE/GRACE-FO month for piecewise regression

--missing, -N

Missing GRACE/GRACE-FO months

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

Spatial resolution of output data

--interval

Output grid interval

1: (0:360, 90:-90)
2: (degree spacing/2)
3: non-global grid (set with defined bounds)

--bounds

Bounding box for non-global grid

--format, -F

Possible choices: ascii, netCDF4, HDF5

Input/output data format

--redistribute-removed

Redistribute removed mass fields over the ocean

--cycles

Regression fit cyclical terms

--log

Output log file for each job

--verbose, -V

Verbose output of run

--mode, -M

Permissions mode of output files