mod_time::mod Interface Reference

Public Member Functions
type(time) function	modulo_time (A, B)

Detailed Description

Definition at line 69 of file mod_time.f90.

Member Function/Subroutine Documentation

modulo_time()

type(time) function mod_time::mod::modulo_time	(	type(time), intent(in)	A,
		type(time), intent(in)	B
	)

Definition at line 177 of file mod_time.f90.

! COMMENTS: THIS PROGRAM IS ROBUST BUT SLOW.....
     IMPLICIT NONE
     TYPE(TIME), INTENT(IN) ::A,B
     real(DP)               :: DPA,DPB,DIV,CNT
     TYPE(TIME) ::T, ZT

     zt%mjd=0
     zt%MuSOD=0

     dpa = days(a)
     dpb = days(b)
     div = dpa/dpb

     cnt = 1.0_dp

     IF(abs(div) .GT. 1e5_dp) cnt = log10(abs(div))

     IF( abs(div) .GT. 1e16_dp) THEN
        ! THE DIVISOR EXCEEDS THE NUMERICAL ACCURACY - RETURN A BAD RESULT!
        modulo_time%MJD = -huge(zt%mjd)
        modulo_time%MUSOD = -huge(zt%musod)
        RETURN
     END IF

     

     IF (b .EQ. zt) THEN
     ! SECOND ARGUMENT IS ZERO - ESCAPE
        modulo_time%MJD = huge(zt%mjd)
        modulo_time%MUSOD = huge(zt%musod)

     ELSE IF(abs(a) .LT. abs(b))THEN
     ! THE FIRST IS SMALLER THAN THE SECOND
        modulo_time = a

     ELSE IF(abs(a) .EQ. abs(b))THEN
    ! THE FIRST IS SMALLER THAN THE SECOND
        modulo_time = zt
        
     ELSE IF (a .GT. zt .and. b .GT. zt) THEN
        ! BOTH ARE POSITIVE

        t = b * int((div-cnt),itime)
        DO WHILE (a .GE. t+b)
           t = t +b
        END DO
        
        modulo_time = a - t
     ELSE IF(a .LT. zt .and. b .LT. zt) THEN
        ! BOTH ARE NEGATIVE

        t = b * int((div-cnt),itime)
        DO WHILE(a .LE. t+b)
           t = t +b
        END DO
        modulo_time = a - t
     ELSE IF (a .LT. zt .and. b .GT. zt) THEN
        ! A IS NEGATIVE AND B IS POSITIVE

        t = b *  int((div+cnt),itime)
        DO WHILE(a .LE. t-b)
           t = t -b
        END DO
        modulo_time = a - t

     ELSE IF (a .GT. zt .and. b .LT. zt) THEN
        ! A IS POSITIVE AND B IS NEGATIVE

        t = b *  int((div+cnt),itime)
        DO WHILE(a .GE. t-b)
           t = t-b
        END DO
        modulo_time = a - t
     ELSE

        ! THIS SHOULD NEVER HAPPEN!
        modulo_time%MJD = -huge(zt%mjd)
        modulo_time%MUSOD = -huge(zt%musod)

     END IF

---

The documentation for this interface was generated from the following file:

• FVCOMg01/mod_time.f90