reanalysis_pressure_harmonics.py

Reads atmospheric surface pressure fields from reanalysis and calculates sets of spherical harmonics using a thin-layer 2D geometry with realistic earth [2, 24]

Calling Sequence

usage: reanalysis_pressure_harmonics.py [-h] [--directory DIRECTORY]
                                        [--year YEAR [YEAR ...]]
                                        [--mean START END] [--redistribute]
                                        [--lmax LMAX] [--mmax MMAX]
                                        [--love {0,1,2,3,4}]
                                        [--reference {CF,CM,CE}]
                                        [--format {ascii,netCDF4,HDF5}]
                                        [--verbose] [--mode MODE]
                                        {ERA-Interim,ERA5,MERRA-2,NCEP-DOE-2,NCEP-CFSR,JRA-55} [{ERA-Interim,ERA5,MERRA-2,NCEP-DOE-2,NCEP-CFSR,JRA-55} ...]

Positional Arguments

model

Possible choices: ERA-Interim, ERA5, MERRA-2, NCEP-DOE-2, NCEP-CFSR, JRA-55

Reanalysis Model

Named Arguments

--directory, -D

Working data directory

--year, -Y

Years of model outputs to run

--mean

Start and end year range for mean

--redistribute

Redistribute pressure values over the ocean

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