/[MITgcm]/MITgcm/pkg/cal/cal_addtime.F
ViewVC logotype

Contents of /MITgcm/pkg/cal/cal_addtime.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.1 - (show annotations) (download)
Mon May 14 22:07:26 2001 UTC (22 years, 11 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47e_post, checkpoint44e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint47c_post, release1_p13_pre, checkpoint50c_post, checkpoint46f_post, checkpoint48e_post, checkpoint50c_pre, checkpoint44f_post, checkpoint46b_post, checkpoint43a-release1mods, ecco_c50_e32, ecco_c50_e33, ecco_c50_e30, ecco_c50_e31, release1_p13, checkpoint40pre3, checkpoint40pre1, checkpoint40pre7, checkpoint40pre6, checkpoint48i_post, checkpoint40pre9, checkpoint40pre8, checkpoint46l_pre, chkpt44d_post, checkpoint51, checkpoint50, release1_p8, release1_p9, checkpoint50d_post, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, checkpoint50b_pre, checkpoint44e_pre, checkpoint51f_post, release1_b1, ecco_c51_e34d, ecco_c51_e34e, ecco_c51_e34f, ecco_c51_e34g, ecco_c51_e34a, ecco_c51_e34b, ecco_c51_e34c, checkpoint48b_post, checkpoint43, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, release1_chkpt44d_post, checkpoint47a_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint47d_post, icebear5, icebear4, icebear3, icebear2, checkpoint46d_pre, checkpoint40pre2, checkpoint48d_post, release1-branch_tutorials, checkpoint48f_post, checkpoint45d_post, checkpoint46j_pre, chkpt44a_post, checkpoint44h_pre, checkpoint48h_post, checkpoint40pre4, ecco_c50_e29, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, ecco_c50_e28, chkpt44c_pre, checkpoint48a_post, checkpoint45a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, release1_p12, release1_p10, release1_p11, release1_p16, release1_p17, release1_p14, release1_p15, checkpoint47j_post, ecco_c50_e33a, branch-exfmods-tag, checkpoint44g_post, branchpoint-genmake2, checkpoint46e_pre, checkpoint48c_post, checkpoint45b_post, checkpoint46b_pre, release1-branch-end, release1_final_v1, checkpoint51b_post, checkpoint51c_post, checkpoint46c_pre, checkpoint46, checkpoint47b_post, checkpoint44b_post, ecco_c51_e34, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint50g_post, checkpoint45c_post, ecco_ice2, ecco_ice1, checkpoint44h_post, checkpoint46g_post, release1_p12_pre, checkpoint39, ecco_c44_e22, checkpoint50h_post, checkpoint50e_pre, checkpoint50i_post, ecco_c44_e25, checkpoint51i_pre, checkpoint40pre5, checkpoint47f_post, checkpoint50e_post, chkpt44a_pre, checkpoint46i_post, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, checkpoint46c_post, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, checkpoint50d_pre, checkpoint46e_post, release1_beta1, checkpoint51e_post, checkpoint44b_pre, checkpoint42, checkpoint40, checkpoint41, checkpoint47, checkpoint44, checkpoint45, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint51f_pre, chkpt44c_post, checkpoint48g_post, checkpoint47h_post, checkpoint44f_pre, checkpoint51g_post, checkpoint46d_post, checkpoint50b_post, release1-branch_branchpoint, checkpoint51a_post
Branch point for: c24_e25_ice, branch-exfmods-curt, release1_final, release1-branch, branch-genmake2, release1, ecco-branch, release1_50yr, icebear, release1_coupled 
Added calendar package.
Not currently supported by mitgcm, i.e. disabled by default.
1 C $Header: /u/gcmpack/development/heimbach/ecco_env/pkg/cal/cal_addtime.F,v 1.5 2001/02/02 16:57:22 heimbach Exp $
2
3 #include "CAL_CPPOPTIONS.h"
4
5 subroutine cal_AddTime(
6 I date,
7 I interval,
8 O added,
9 I mythid
10 & )
11
12 c ==================================================================
13 c SUBROUTINE cal_AddTime
14 c ==================================================================
15 c
16 c o Add a time interval either to a calendar date or to a time
17 c interval.
18 c
19 c started: Christian Eckert eckert@mit.edu 30-Jun-1999
20 c
21 c changed: Christian Eckert eckert@mit.edu 29-Dec-1999
22 c
23 c - restructured the original version in order to have a
24 c better interface to the MITgcmUV.
25 c
26 c Christian Eckert eckert@mit.edu 03-Feb-2000
27 c
28 c - Introduced new routine and function names, cal_<NAME>,
29 c for verion 0.1.3.
30 c
31 c ralf.giering@fastopt.de 31-May-2000
32 c datesecs was computed at wrong place (cph)
33 c
34 c ==================================================================
35 c SUBROUTINE cal_AddTime
36 c ==================================================================
37
38 implicit none
39
40 c == global variables ==
41
42 #include "cal.h"
43
44 c == routine arguments ==
45
46 integer date(4)
47 integer interval(4)
48 integer added(4)
49 integer mythid
50
51 c == local variables ==
52
53 integer intsecs
54 integer datesecs
55 integer nsecs
56 integer hhmmss
57 integer yi,mi,di,si,li,wi
58 integer ndays
59 integer date_1,date_2
60 integer intv_1,intv_2
61 integer fac
62 integer iday
63 integer switch
64 integer ndayssub
65 integer ierr
66
67 c == external ==
68
69 integer cal_IsLeap
70 external cal_IsLeap
71
72 c == end of interface ==
73
74 date_1 = 0
75 date_2 = 0
76
77 fac = 1
78
79 if (interval(4) .eq. -1) then
80
81 if (date(4) .eq. -1) then
82 if (date(1) .ge. 0) then
83 date_1 = date(1)
84 date_2 = date(2)
85 intv_1 = interval(1)
86 intv_2 = interval(2)
87 else
88 if (interval(1) .lt. 0) then
89 date_1 = -date(1)
90 date_2 = -date(2)
91 intv_1 = -interval(1)
92 intv_2 = -interval(2)
93 fac = -1
94 else
95 date_1 = interval(1)
96 date_2 = interval(2)
97 intv_1 = date(1)
98 intv_2 = date(2)
99 fac = 1
100 endif
101 endif
102 else
103 if (interval(1) .ge. 0) then
104 intv_1 = interval(1)
105 intv_2 = interval(2)
106 else
107 intv_1 = -interval(1)
108 intv_2 = -interval(2)
109 fac = -1
110 endif
111 endif
112
113 intsecs = fac*(intv_2/10000*secondsperhour +
114 & (mod(intv_2/100,100)*secondsperminute +
115 & mod(intv_2,100)))
116
117 if (date(4) .eq. -1) then
118 datesecs = date_2/10000*secondsperhour +
119 & mod(date_2/100,100)*secondsperminute +
120 & mod(date_2,100)
121 date_1 = date_1 + intv_1
122 nsecs = datesecs + intsecs
123 if ((date_1 .gt. 0) .and.
124 & (nsecs .lt. 0)) then
125 date_1 = date_1 - 1
126 nsecs = nsecs + secondsperday
127 endif
128 nsecs = fac*nsecs
129 yi = 0
130 mi = 0
131 di = fac*date_1
132 li = 0
133 wi = -1
134 else
135 call cal_ConvDate( date,yi,mi,di,si,li,wi,mythid )
136 if ((interval(1) .ge. 0) .and.
137 & (interval(2) .ge. 0)) then
138 nsecs = si + intsecs
139 ndays = nsecs/secondsperday
140 nsecs = mod(nsecs,secondsperday)
141 do iday = 1,interval(1)+ndays
142 di = di + 1
143 if (di .gt. ndaymonth(mi,li)) then
144 di = 1
145 mi = mi + 1
146 endif
147 switch = (mi-1)/nmonthyear
148 yi = yi + switch
149 mi = mod(mi-1,nmonthyear)+1
150 if (switch .eq. 1) li = cal_IsLeap( yi, mythid )
151 enddo
152 wi = mod(wi+interval(1)+ndays-1,7)+1
153 else
154 nsecs = si + intsecs
155 if (nsecs .ge. 0) then
156 ndayssub = intv_1
157 else
158 nsecs = nsecs + secondsperday
159 ndayssub = intv_1 + 1
160 endif
161 do iday = 1,ndayssub
162 di = di - 1
163 if (di .eq. 0) then
164 mi = mod(mi+10,nmonthyear)+1
165 switch = mi/nmonthyear
166 yi = yi - switch
167 if (switch .eq. 1) li = cal_IsLeap( yi, mythid )
168 di = ndaymonth(mi,li)
169 endif
170 enddo
171 wi = mod(wi+6-mod(ndayssub,7),7)+1
172 endif
173 endif
174
175 c Convert to calendar format.
176 added(1) = yi*10000 + mi*100 + di
177 hhmmss = nsecs/secondsperminute
178 added(2) = hhmmss/minutesperhour*10000 +
179 & (mod(fac*hhmmss,minutesperhour)*100 +
180 & mod(fac*nsecs,secondsperminute))*fac
181 added(3) = li
182 added(4) = wi
183 else
184
185 ierr = 601
186 call cal_PrintError( ierr, mythid)
187 stop ' stopped in cal_AddTime.'
188
189 endif
190
191 return
192 end
193
  ViewVC Help
Powered by ViewVC 1.1.22