用Fortran写的代码如下:
program Moon_Angles
implicit none
real :: T, az, alt
T = get_time()
call Moon_Position(T, az, alt)
print *, az, alt
end program Moon_Angles
subroutine Moon_Position(T, az, alt)
implicit none
real, intent(in) :: T real, intent(out) :: az, alt
real :: a, b, c, d, e, f, g, h