This subroutine returns the values for nutation in longitude and nutation in obliquity for a given TDB Julian date.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in) | :: | t |
tdb time in julian centuries since j2000.0 |
||
real(kind=wp), | intent(out) | :: | dpsi |
nutation in longitude in arcseconds |
||
real(kind=wp), | intent(out) | :: | deps |
nutation in obliquity in arcseconds |
subroutine nod (t,dpsi,deps) real(wp),intent(in) :: t !! tdb time in julian centuries since j2000.0 real(wp),intent(out) :: dpsi !! nutation in longitude in arcseconds real(wp),intent(out) :: deps !! nutation in obliquity in arcseconds real(wp) :: t1,dp,de ! get method/accuracy mode call getmod ( mode ) t1 = t * 36525.0_wp ! evaluate nutation series ! resulting nutation in longitude and obliquity in arcseconds ! call subroutine to evaluate nutation series if ( mod ( mode, 2 ) == 0 ) then ! high accuracy mode -- iers subroutine call nu2000a ( jd_j2000, t1, dp, de ) else ! low accuracy mode -- modification of iers subroutine call nu2000k ( jd_j2000, t1, dp, de ) end if dpsi = dp * seccon deps = de * seccon end subroutine nod