Changeset 5965 for branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/TOOLS/OBSTOOLS/src/date_utils.F90

Timestamp:

2015-12-01T16:35:30+01:00 (8 years ago)

Author:

timgraham

Message:

Upgraded branch to r5518 of trunk (v3.6 stable revision)

File:

• branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/TOOLS/OBSTOOLS/src/date_utils.F90 (modified) (2 diffs)

branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/TOOLS/OBSTOOLS/src/date_utils.F90

@@ -34 +34 @@
 
    END SUBROUTINE add_date
+
+
+   SUBROUTINE add_days_to_date(initial_date,days,final_date)
+
+      ! Add a number of days to initial_date and return it in final_date
+
+      IMPLICIT NONE
+
+
+      !! Arguments
+      INTEGER,INTENT(in) :: initial_date ! Initial date (YYYYMMDD)
+      INTEGER,INTENT(in) :: days         ! Number of days to add
+      INTEGER,INTENT(out) :: final_date  ! Final date (YYYYMMDD)
+
+      !! Local variables
+
+      INTEGER :: isec,imin,ihours,iyear,imon,iday ! temporary results
+      REAL(dp):: juld
+
+      ! Account for lack of hours in date format (initial_date*100)
+      CALL split_date(initial_date*100,iyear,imon,iday,ihours)
+
+      CALL greg2jul(0,0,ihours,iday,imon,iyear,juld)
+
+      juld=juld+REAL(days)
+
+      CALL jul2greg(isec,imin,ihours,iday,imon,iyear,juld)
+
+      final_date=(iyear*1000000+imon*10000+iday*100+ihours)/100
+
+   END SUBROUTINE add_days_to_date
+
 
    SUBROUTINE split_date(iyyyymmddhh,iyyyy,imm,idd,ihh)
@@ -379 +411 @@
    END FUNCTION diffdate
 
+   INTEGER FUNCTION difftime(itime1,itime2)
+
+      ! Compute difference in minutes between times
+      ! Assumes HHMM or HMM or MM or M format for dates
+      !
+      ! ORDER MATTERS - itime1 is ealier time
+      ! Result is an integer number of minutes
+
+      IMPLICIT NONE
+      INTEGER, INTENT(IN) :: itime1,itime2    ! Times to be diffed.
+      INTEGER :: imin1, imin2, ihr1, ihr2
+
+      ihr1 = (itime1/100)
+      ihr2 = (itime2/100)
+      
+      imin1 = (itime1 - ihr1*100) + (60 * ihr1)
+      imin2 = (itime2 - ihr2*100) + (60 * ihr2)
+      
+      ! Assume that itime2 is later, so wrap around midnight if necessary.
+      IF (imin2 < imin1) THEN
+         imin2 = imin2 + 24*60
+      END IF
+      
+      difftime = imin2 - imin1
+      
+   END FUNCTION difftime
+
+
+   INTEGER FUNCTION add_mins_to_time(itime1, imin_add)
+
+      ! Add number of minutes onto given time
+      ! Assumes time in HHMM or HMM or MM or M format
+      !
+      ! Result is in HHMM format
+
+      IMPLICIT NONE
+      INTEGER, INTENT(IN) :: itime1,imin_add
+      INTEGER :: imin1, ihr1, imin2, ihr2
+
+      ihr1 = (itime1/100)
+      
+      ! itime1 in minutes from previous midnight
+      imin1 = (itime1 - ihr1*100) + (60 * ihr1)
+      
+      imin1 = imin1 + imin_add
+      
+      ! Add 1day if time went nagative
+      IF (imin1 < 0) THEN
+         imin1 = imin1 + 24*60
+      END IF
+      
+      ! Turn number of minutes back into HHMM
+      ihr2 = imin1/60
+      imin2 = imin1 - ihr2*60
+
+      DO
+         IF (ihr2<0) THEN
+            ihr2 = ihr2 + 24
+         ELSE IF (ihr2>=24) THEN
+            ihr2 = ihr2 - 24
+         END IF
+         IF ((ihr2>=0).OR.(ihr2<24)) EXIT
+      END DO
+
+      add_mins_to_time = ihr2*100 + imin2
+   
+   END FUNCTION add_mins_to_time
+
 
 END MODULE date_utils