racmo_smb_harmonics.py

Reads monthly RACMO variables and converts to spherical harmonic coefficients [29]
Estimates the grid values as point masses for calculating the gravitational spherical harmonic coefficients

Calling Sequence

model_file: full path to input RACMO netCDF4 file

Command Line Options

usage: racmo_smb_harmonics.py [-h] [--product PRODUCT]
                              [--mask MASK [MASK ...]] [--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: RACMO SMB file to run

Named Arguments

--product, -P

Possible choices: precip, rainfall, refreeze, runoff, smb, sndiv, snowfall, snowmelt, subl

RACMO SMB product to calculate

'precip': Precipitation
'rainfall': Rainfall
'refreeze': Meltwater Refreeze
'runoff': Meltwater Runoff
'smb': Surface Mass Balance
'sndiv': Snow Drift Erosion
'snowfall': Snowfall
'snowmelt': Snowmelt
'subl': Sublimation

--mask

netCDF4 masks file for reducing to regions

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