ticktime.f90

CLM集合卡曼滤波数据同化算法

       SUBROUTINE ticktime (dtime, jday, iyear, mcsec)

!=======================================================================
!      Source file: ticktime.f90
! Original version: Yongjiu Dai, September 15, 1999
!
!      Notes: Greenwich time
!
!=======================================================================

      IMPLICIT NONE

      real, INTENT(in) :: dtime            ! model time step (second)
      integer, INTENT(inout) :: iyear,jday,mcsec ! year and julian day, respectively

!     local
      integer maxday                       ! days of one year (365 or 366)
!-----------------------------------------------------------------------

      mcsec = mcsec + nint(dtime)
      if(mcsec>86400)then
         mcsec = mcsec - 86400
         if((mod(iyear,4) == 0 .AND. mod(iyear,100) /=0) .OR. mod(iyear,400) == 0)then
             maxday = 366
         else
             maxday = 365
         endif
         jday = jday + 1
         if(jday > maxday) then
            iyear = iyear + 1
            jday = 1
         endif
      endif

      END SUBROUTINE ticktime