get_c_cdot_iau_earth Subroutine

private subroutine get_c_cdot_iau_earth(me, eph, to_icrf, c, cdot, status_ok)

rotation matrix for IAU_EARTH <-> ICRF

Type Bound

iau_earth_rotating_frame

Arguments

Type IntentOptional Attributes Name
class(iau_earth_rotating_frame), intent(inout) :: me
class(ephemeris_class), intent(inout) :: eph

for ephemeris computations (assumed to have already been initialized)
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

Calls

proc~~get_c_cdot_iau_earth~~CallsGraph proc~get_c_cdot_iau_earth transformation_module::iau_earth_rotating_frame%get_c_cdot_iau_earth proc~icrf_to_iau_earth iau_orientation_module::icrf_to_iau_earth proc~get_c_cdot_iau_earth->proc~icrf_to_iau_earth proc~iau_rotation_matrix iau_orientation_module::iau_rotation_matrix proc~icrf_to_iau_earth->proc~iau_rotation_matrix proc~rotation_matrix vector_module::rotation_matrix proc~iau_rotation_matrix->proc~rotation_matrix

Source Code

    subroutine get_c_cdot_iau_earth(me,eph,to_icrf,c,cdot,status_ok)

    !! rotation matrix for IAU_EARTH <-> ICRF

    implicit none

    class(iau_earth_rotating_frame),intent(inout) :: me
    class(ephemeris_class),intent(inout)          :: eph     !! for ephemeris computations (assumed to have already been initialized)
    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

    c = icrf_to_iau_earth(me%et)

    !... don't have the cdot code yet... need to refactor iau code ...
    ! see also: ftp://naif.jpl.nasa.gov/pub/naif/toolkit_docs/FORTRAN/spicelib/tisbod.html

    !if (present(cdot)) cdot =
    error stop 'not yet supported'

    status_ok = .true.

    end subroutine get_c_cdot_iau_earth