nodal_adj {TideHarmonics}

Calculate nodal corrections.

Description

Calculates nodal corrections from astronomical longitudes (lambdas). Mainly for internal use.

Usage

  nodal_adj(lambp, lambN, lambph, indegree = TRUE, outdegree = TRUE) 

Arguments

lambp

Lambdas for the lunar perigee (p).

lambN

Lambdas for the lunar nodal (N).

lambph

Lambdas for the sun's perihelion (ph).

indegree

If FALSE, the lambdas are in radians. Degrees are assumed by default.

outdegree

If FALSE, the returned phase corrections are in radians. Degrees are produced by default.

Value

A list with two elements, where both are matrices with 409 rows, which correspond to the 409 harmonic constituents available in the package. The first element fn gives the amplitude corrections. The second element un gives the phase corrections.

See Also

ftide, lambdas

Examples

days <- seq(as.Date("2012-12-30"), as.Date("2013-01-08"), 1)
lamb <- lambdas(days)
nodal_adj(lamb[3,], lamb[4,], lamb[5,])