pkg/cal/cal_toseconds.F

C $Header: /u/gcmpack/MITgcm/pkg/cal/cal_toseconds.F,v 1.2 2002/12/19 13:12:42 cheisey Exp $

#include "CAL_CPPOPTIONS.h"

      subroutine cal_ToSeconds(
     I                          date,
     O                          timeint,
     I                          mythid
     &                        )

c     ==================================================================
c     SUBROUTINE cal_ToSeconds
c     ==================================================================
c
c     o Given a time interval as a date array return the number of
c       seconds in that time interval.
c
c       If one wanted to use calendar dates in this routine, then
c       the date should be after the calendar's refdate and timeint
c       would be the number of seconds that have elapsed since the
c       refdate. Of course this can also be done by first calling
c       sub cal_TimePassed and then calling this routine with the
c       resulting time interval array.
c              
c     started: Christian Eckert eckert@mit.edu  30-Jun-1999
c
c     changed: Christian Eckert eckert@mit.edu  29-Dec-1999
c
c              - restructured the original version in order to have a
c                better interface to the MITgcmUV.
c
c              Christian Eckert eckert@mit.edu  03-Feb-2000
c
c              - Introduced new routine and function names, cal_<NAME>,
c                for verion 0.1.3.
c
c              21-Sep-2003: fixed check_sign logic to work with
c              negative intervals (menemenlis@jpl.nasa.gov)
c
c     ==================================================================
c     SUBROUTINE cal_ToSeconds
c     ==================================================================

      implicit none

c     == global variables ==

#include "cal.h"

c     == routine arguments ==

      integer date(4)
      _RL     timeint
      integer mythid

c     == local variables ==

      integer fac
      integer ndays
      integer nsecs
      integer hhmmss
      integer ierr
      _RL check_sign

c     == end of interface ==
c     print *,'cal_toseconds: date',date
c     print *,'cal_toseconds: timeint',timeint

      check_sign = sign(1,dble(date(1))*dble(date(2)))
      if (((date(4) .eq. -1) .and.
     &    (date(3) .eq.  0) .and.
     &    (check_sign .ge. 0)) .or.
     &    usingModelCalendar) then
        if ((date(1) .lt. 0) .or.
     &      (date(2) .lt. 0)) then
          ndays  = -date(1)
          hhmmss = -date(2)
          fac    = -1
        else
          ndays  = date(1)
          hhmmss = date(2)
          fac    = 1
        endif
        nsecs   = ndays*secondsperday +
     &            (hhmmss/10000)*secondsperhour +
     &            mod(hhmmss/100,100)*secondsperminute +
     &            mod(hhmmss,100)
        timeint = fac*nsecs
      else

        ierr = 1001
        call cal_PrintError( ierr, mythid )
        stop ' stopped in cal_ToSeconds.'

      endif

      return
      end
1	dimitri	1.3	C $Header: /u/gcmpack/MITgcm/pkg/cal/cal_toseconds.F,v 1.2 2002/12/19 13:12:42 cheisey Exp $
2	heimbach	1.1
3			#include "CAL_CPPOPTIONS.h"
4
5			subroutine cal_ToSeconds(
6			I date,
7			O timeint,
8			I mythid
9			& )
10
11			c ==================================================================
12			c SUBROUTINE cal_ToSeconds
13			c ==================================================================
14			c
15			c o Given a time interval as a date array return the number of
16			c seconds in that time interval.
17			c
18			c If one wanted to use calendar dates in this routine, then
19			c the date should be after the calendar's refdate and timeint
20			c would be the number of seconds that have elapsed since the
21			c refdate. Of course this can also be done by first calling
22	dimitri	1.3	c sub cal_TimePassed and then calling this routine with the
23	heimbach	1.1	c resulting time interval array.
24			c
25			c started: Christian Eckert eckert@mit.edu 30-Jun-1999
26			c
27			c changed: Christian Eckert eckert@mit.edu 29-Dec-1999
28			c
29			c - restructured the original version in order to have a
30			c better interface to the MITgcmUV.
31			c
32			c Christian Eckert eckert@mit.edu 03-Feb-2000
33			c
34			c - Introduced new routine and function names, cal_<NAME>,
35			c for verion 0.1.3.
36			c
37	dimitri	1.3	c 21-Sep-2003: fixed check_sign logic to work with
38			c negative intervals (menemenlis@jpl.nasa.gov)
39			c
40	heimbach	1.1	c ==================================================================
41			c SUBROUTINE cal_ToSeconds
42			c ==================================================================
43
44			implicit none
45
46			c == global variables ==
47
48			#include "cal.h"
49
50			c == routine arguments ==
51
52			integer date(4)
53			_RL timeint
54			integer mythid
55
56			c == local variables ==
57
58			integer fac
59			integer ndays
60			integer nsecs
61			integer hhmmss
62			integer ierr
63	dimitri	1.3	_RL check_sign
64	heimbach	1.1
65			c == end of interface ==
66	cheisey	1.2	c print *,'cal_toseconds: date',date
67			c print *,'cal_toseconds: timeint',timeint
68	heimbach	1.1
69	dimitri	1.3	check_sign = sign(1,dble(date(1))*dble(date(2)))
70	cheisey	1.2	if (((date(4) .eq. -1) .and.
71	heimbach	1.1	& (date(3) .eq. 0) .and.
72	dimitri	1.3	& (check_sign .ge. 0)) .or.
73	cheisey	1.2	& usingModelCalendar) then
74	heimbach	1.1	if ((date(1) .lt. 0) .or.
75			& (date(2) .lt. 0)) then
76			ndays = -date(1)
77			hhmmss = -date(2)
78			fac = -1
79			else
80			ndays = date(1)
81			hhmmss = date(2)
82			fac = 1
83			endif
84			nsecs = ndays*secondsperday +
85			& (hhmmss/10000)*secondsperhour +
86			& mod(hhmmss/100,100)*secondsperminute +
87			& mod(hhmmss,100)
88			timeint = fac*nsecs
89			else
90
91			ierr = 1001
92			call cal_PrintError( ierr, mythid )
93			stop ' stopped in cal_ToSeconds.'
94
95			endif
96
97			return
98			end