Gravityfield2DisplacementTimeSeries

This program computes a time series of displacements of a list of stations (grid) due to the effect of time variable loading masses. The displacement $\M u$ of a station is calculated according to \[\label{eq:displacement} \M u(\M r) = \frac{1}{\gamma}\sum_{n=0}^\infty \left[\frac{h_n}{1+k_n}V_n(\M r)\,\M e_{up} + R\frac{l_n}{1+k_n}\left( \frac{\partial V_n(\M r)}{\partial \M e_{north}}\M e_{north} +\frac{\partial V_n(\M r)}{\partial \M e_{east}} \M e_{east}\right)\right], \]where $\gamma$ is the normal gravity, the load Love and Shida numbers $h_n,l_n$ are given by inputfileDeformationLoadLoveNumber and the load Love numbers $k_n$ are given by inputfilePotentialLoadLoveNumber. The $V_n$ are the spherical harmonics expansion of the full time variable gravitational potential (potential of the loading mass + deformation potential): \[ V(\M r) = \sum_{n=0}^\infty V_n(\M r). \]Deformations due to Earth tide and due to polar tides are computed using the IERS conventions. Eq. \eqref{eq:displacement} is not used in these cases.

The outputfileTimeSeries is an instrument file, MISCVALUES. The data columns contain the deformation of each station in $x,y,z$ in a global terrestrial reference frame or alternatively in a local elliposidal frame (north, east, up) if localReferenceFrame is set.