time_module_test Subroutine

public subroutine time_module_test()

Test routine for the Julian date routines.

Arguments

None

Calls

proc~~time_module_test~~CallsGraph proc~time_module_test time_module_test interface~julian_date julian_date proc~time_module_test->interface~julian_date proc~calendar_date_realsec calendar_date_realsec proc~time_module_test->proc~calendar_date_realsec proc~calendar_date_to_et calendar_date_to_et proc~time_module_test->proc~calendar_date_to_et proc~jd_to_et jd_to_et proc~time_module_test->proc~jd_to_et proc~julian_date_intsec julian_date_intsec interface~julian_date->proc~julian_date_intsec proc~julian_date_realsec julian_date_realsec interface~julian_date->proc~julian_date_realsec proc~is_leap_year is_leap_year proc~calendar_date_to_et->proc~is_leap_year proc~julian_date_intsec->proc~julian_date_realsec proc~julian_day julian_day proc~julian_date_realsec->proc~julian_day

Source Code

    subroutine time_module_test()

    implicit none

    real(wp) :: jd, sec
    integer :: year,month,day,hrs,min
    real(wp) :: et1, et2

    write(*,*) ''
    write(*,*) '---------------'
    write(*,*) ' time_module_test'
    write(*,*) '---------------'
    write(*,*) ''

    ! JD = 2451545.0
    jd = julian_date(2000,1,1,12,0,0)

    call calendar_date_realsec(jd,year,month,day,hrs,min,sec)

    write(*,*) 'jd    ', jd
    write(*,*) 'year  ', year
    write(*,*) 'month ', month
    write(*,*) 'day   ', day
    write(*,*) 'hrs   ', hrs
    write(*,*) 'min   ', min
    write(*,*) 'sec   ', sec

    if (year/=2000)  error stop 'error: incorrect year'
    if (month/=1)    error stop 'error: incorrect month'
    if (day/=1)      error stop 'error: incorrect day'
    if (hrs/=12)     error stop 'error: incorrect hrs'
    if (min/=0)      error stop 'error: incorrect min'
    if (sec/=0.0_wp) error stop 'error: incorrect sec'

    ! compare two ways to convert calendar date to ephemeris time:
    ! et1 should be more accurate since it doesn't go through jd
    write(*,*) ''
    et1 = calendar_date_to_et(2026,3,4,8,9,10.12345678_wp)
    et2 = jd_to_et(julian_date(2026,3,4,8,9,10.12345678_wp))
    write(*,*) 'et1   ', et1
    write(*,*) 'et2   ', et2
    write(*,*) 'et diff: ', abs(et1-et2)

    write(*,*) ''
    write(*,*) 'PASSED'

    end subroutine time_module_test