Rotation matrix from J2000 to Mean Ecliptic.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
procedure(mean_obliquity_func), | optional | :: | obliquity_func |
optional function to compute the mean obliquity. If not present, then mean_obliquity_of_ecliptic_iau1980 is used. |
rotation matrix
pure function equatorial_to_mean_ecliptic_rotmat(obliquity_func) result(rot) implicit none real(wp),dimension(3,3) :: rot !! rotation matrix procedure(mean_obliquity_func),optional :: obliquity_func !! optional function to compute !! the mean obliquity. If not !! present, then !! [[mean_obliquity_of_ecliptic_iau1980]] !! is used. rot = transpose(mean_ecliptic_to_equatorial_rotmat(obliquity_func)) end function equatorial_to_mean_ecliptic_rotmat