rotation matrix for ICRF <-> Mean Ecliptic
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(ecliptic_frame), | intent(inout) | :: | me | |||
class(ephemeris_class), | intent(inout) | :: | eph | |||
logical, | intent(in) | :: | to_icrf | |||
real(kind=wp), | intent(out), | dimension(3,3) | :: | c | ||
real(kind=wp), | intent(out), | optional, | dimension(3,3) | :: | cdot | |
logical, | intent(out) | :: | status_ok |
subroutine get_c_cdot_ecliptic(me,eph,to_icrf,c,cdot,status_ok) !! rotation matrix for ICRF <-> Mean Ecliptic use obliquity_module implicit none class(ecliptic_frame),intent(inout) :: me class(ephemeris_class),intent(inout) :: eph logical,intent(in) :: to_icrf real(wp),dimension(3,3),intent(out) :: c real(wp),dimension(3,3),intent(out),optional :: cdot logical,intent(out) :: status_ok if (to_icrf) then c = mean_ecliptic_to_equatorial_rotmat() else c = equatorial_to_mean_ecliptic_rotmat() end if if (present(cdot)) cdot = zero status_ok = .true. end subroutine get_c_cdot_ecliptic