$LastChangedDate: 27-Juil-2011 $

# CL_mod_nutationMatrix

Nutation matrix and rotation vectors (IERS 1996)

## Calling Sequence

       [NUT[,OM[,OMP]]] = CL_mod_nutationMatrix(K,KP,KPP[,flag])

## Description

- Computes the nutation matrix and the rotation vectors (standard IERS 1996).

## Parameters

K :

  [EPS0 ; DPSI ; DEPS ] (3xN)

KP :

  [EPS0P ; DPSIP ; DEPSP] (3xN)

KPP :

  [EPS0PP; DPSIPP; DEPSPP] (3xN)

flag :

  (optional) "n" -> matrix only; "p" -> matrix and rotation vector; "s" -> matrix, rotation vector and its time derivative (default depends on number of outputs)

NUT :

  Nutation matrix (3x3xN)

OM :

  (optional) Rotation vectors (3xN)

OMP :

  (optional) first time derivative of OM (3xN)

## Authors

- CNES - DCT/SB

## Bibliography

- 1) IERS Conventions (1996), Dennis D. McCarthy

- 2) Explanatory Supplement to the Astronomical Almanac, Seidelman (1992)

## See also

- CL_mod_nutationAngles

- CL_mod_nutationArg

- CL_mod_meanObliquity

- CL_rot_angles2matrix

- CL_rot_angularVelocity

## Examples

    jj_tai = [19500:1:20500];
    ss_tai = 43200*ones(jj_tai);
    // Nutation angles (dPsi,dEps) :
    [NUT,F,NUTP,FP,NUTPP,FPP]=CL_mod_nutationAngles(jj_tai,ss_tai,"s");

    // Mean obliquity :
    [eps0,eps0p,eps0pp]=CL_mod_meanObliquity(jj_tai,ss_tai,"s");
    K=[eps0; NUT]; Kp=[eps0p; NUTP]; Kpp=[eps0pp; NUTPP];

    [NUT,OM,OMP] = CL_mod_nutationMatrix(K,Kp,Kpp,"s");