/[MITgcm]/MITgcm/pkg/timeave/timeave_cumul_2v.F
ViewVC logotype

Annotation of /MITgcm/pkg/timeave/timeave_cumul_2v.F

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

Revision 1.7 - (hide annotations) (download)
Fri Aug 19 22:50:26 2005 UTC (18 years, 9 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint62u, checkpoint62t, checkpoint57s_post, checkpoint58b_post, checkpoint57y_post, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint58d_post, checkpoint62s, checkpoint58a_post, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint57z_post, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint58t_post, checkpoint58m_post, checkpoint57t_post, checkpoint63g, checkpoint57v_post, checkpoint64, checkpoint65, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint58w_post, checkpoint57y_pre, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, mitgcm_mapl_00, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint58r_post, checkpoint58n_post, checkpoint65o, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpint57u_post, checkpoint57q_post, checkpoint58k_post, checkpoint62b, checkpoint58v_post, checkpoint64y, checkpoint64x, checkpoint58l_post, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint59j, checkpoint58h_post, checkpoint58j_post, checkpoint61q, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61z, checkpoint61x, checkpoint61y, HEAD
Changes since 1.6: +12 -12 lines 
Move deltaT to deltaTloc to avoid conflict with PARAMS.h
1 heimbach 1.7 C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_cumul_2v.F,v 1.6 2004/07/26 16:57:29 jmc Exp $
2 jmc 1.1 C $Name: $
3 edhill 1.5 #include "TIMEAVE_OPTIONS.h"
4 jmc 1.1
5     CStartofinterface
6     SUBROUTINE TIMEAVE_CUMUL_2V(
7     O fldtave,
8 heimbach 1.7 I fld1, fld2, Ksize, dir, deltaTloc,
9 jmc 1.1 I bi, bj, myThid )
10 jmc 1.3 C /==========================================================*
11     C | SUBROUTINE TIMEAVE_CUMUL_2V
12     C | o Sum over time a product of two arrays depending on the
13     C | relative position of the 2 fields.
14     C | (tracer point, u, v, w ...)
15     C \==========================================================*
16 jmc 1.1 IMPLICIT NONE
17    
18     C == Global variables ===
19     #include "SIZE.h"
20     #include "EEPARAMS.h"
21 jmc 1.4 #include "GRID.h"
22 jmc 1.1
23     C == Routine arguments ==
24     C myThid - Thread number for this instance of the routine.
25     C fldtave - time averaged Field
26     C fld1,fld2 - Input Field
27 jmc 1.2 C dir - type of grid for 2nd array relatively to the 1rst array
28     C 0: same grid ; 1: dX/2 shift ; 2: dY/2 shift ; 3: dr/2 shift
29 jmc 1.3 C (2 digits => also shift the 1rst array)
30 jmc 1.1 C Ksize - 3rd dimension of local arrays (Input and Output fields)
31     INTEGER Ksize, dir
32     _RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
33     _RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
34     _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
35 heimbach 1.7 _RL deltaTloc
36 jmc 1.1 INTEGER bi, bj, myThid
37    
38     CEndofinterface
39    
40 jmc 1.3 #ifdef ALLOW_TIMEAVE
41 jmc 1.1 C == Local variables ==
42     C i,j,k,bi,bj - Loop counters
43     INTEGER i, j, k
44 jmc 1.6 INTEGER km1
45 jmc 1.1
46     IF ( dir.eq.0 ) THEN
47 jmc 1.4 c- both fields at the same location :
48 jmc 1.1
49     C DO bj = myByLo(myThid), myByHi(myThid)
50     C DO bi = myBxLo(myThid), myBxHi(myThid)
51     DO k=1,Ksize
52     DO j=1,sNy
53     DO i=1,sNx
54     fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj)
55 heimbach 1.7 & + fld1(i,j,k,bi,bj)*fld2(i,j,k,bi,bj)*deltaTloc
56 jmc 1.1 ENDDO
57     ENDDO
58 jmc 1.3 ENDDO
59 jmc 1.1 C ENDDO
60 jmc 1.3 C ENDDO
61 jmc 1.1
62     ELSEIF ( dir.eq.1 ) THEN
63 jmc 1.4 c- 2nd field shifted by -dX/2 (e.g.: 1=T, 2=U) :
64 jmc 1.1
65     DO k=1,Ksize
66     DO j=1,sNy
67     DO i=1,sNx
68     fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj)
69     & + .5 * ( fld1(i-1,j,k,bi,bj) + fld1(i,j,k,bi,bj) )
70     & * fld2(i,j,k,bi,bj)
71 heimbach 1.7 & * deltaTloc
72 jmc 1.1 ENDDO
73     ENDDO
74 jmc 1.3 ENDDO
75 jmc 1.1
76     ELSEIF ( dir.eq.2 ) THEN
77 jmc 1.4 c- 2nd field shifted by -dY/2 (e.g.: 1=T, 2=V) :
78 jmc 1.1
79     DO k=1,Ksize
80     DO j=1,sNy
81     DO i=1,sNx
82     fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj)
83     & + .5 * ( fld1(i,j-1,k,bi,bj) + fld1(i,j,k,bi,bj) )
84 jmc 1.2 & * fld2(i,j,k,bi,bj)
85 heimbach 1.7 & * deltaTloc
86 jmc 1.2 ENDDO
87     ENDDO
88 jmc 1.3 ENDDO
89 jmc 1.2
90     ELSEIF ( dir.eq.3 ) THEN
91 jmc 1.4 c- 2nd field shifted by -dR/2 (e.g.: 1=T, 2=W) :
92 jmc 1.2
93 jmc 1.6 DO k=1,Ksize
94     km1 = MAX(k-1,1)
95 jmc 1.2 DO j=1,sNy
96     DO i=1,sNx
97     fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj)
98 jmc 1.6 & + .5 * ( fld1(i,j,km1,bi,bj) + fld1(i,j,k,bi,bj) )
99 jmc 1.1 & * fld2(i,j,k,bi,bj)
100 heimbach 1.7 & * deltaTloc
101 jmc 1.1 ENDDO
102     ENDDO
103 jmc 1.3 ENDDO
104    
105     ELSEIF ( dir.eq.12 ) THEN
106 jmc 1.4 c- 1rst & 2nd fields shifted by -dY/2 & -dX/2
107     c (e.g.: 1=U, 2=V, product at the corner) :
108 jmc 1.3
109     DO k=1,Ksize
110     DO j=1,sNy
111     DO i=1,sNx
112     fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
113     & + .25 _d 0*( fld1(i,j-1,k,bi,bj) + fld1(i,j,k,bi,bj) )
114     & *( fld2(i-1,j,k,bi,bj) + fld2(i,j,k,bi,bj) )
115 heimbach 1.7 & * deltaTloc
116 jmc 1.4 ENDDO
117     ENDDO
118     ENDDO
119    
120     ELSEIF ( dir.eq.13 ) THEN
121     c- 1rst & 2nd fields shifted by -dR/2 & -dX/2 (e.g.: 1=U, 2=W):
122    
123 jmc 1.6 DO k=1,Ksize
124     km1 = MAX(k-1,1)
125 jmc 1.4 DO j=1,sNy
126     DO i=1,sNx
127     fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
128 jmc 1.6 & + .25 _d 0*( fld1(i,j,km1,bi,bj) + fld1(i,j,k,bi,bj) )
129 jmc 1.4 & *( fld2(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)
130     & +fld2( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
131     & )*recip_rAw(i,j,bi,bj)
132 heimbach 1.7 & * deltaTloc
133 jmc 1.4 ENDDO
134     ENDDO
135     ENDDO
136    
137     ELSEIF ( dir.eq.23 ) THEN
138     c- 1rst & 2nd fields shifted by -dR/2 & -dY/2 (e.g.: 1=V, 2=W):
139    
140 jmc 1.6 DO k=1,Ksize
141     km1 = MAX(k-1,1)
142 jmc 1.4 DO j=1,sNy
143     DO i=1,sNx
144     fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
145 jmc 1.6 & + .25 _d 0*( fld1(i,j,km1,bi,bj) + fld1(i,j,k,bi,bj) )
146 jmc 1.4 & *( fld2(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj)
147     & +fld2(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
148     & )*recip_rAs(i,j,bi,bj)
149 heimbach 1.7 & * deltaTloc
150 jmc 1.4 ENDDO
151     ENDDO
152     ENDDO
153    
154     ELSEIF ( dir.eq.-13 ) THEN
155     c- gradient of the 1rst field * 2nd fields, shifted by -dR/2 & -dX/2 resp.
156     c- (e.g.: used for advective form of vertical advection: w.du/dr)
157    
158     DO k=2,Ksize
159     DO j=1,sNy
160     DO i=1,sNx
161     fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
162     & + .5 _d 0*( fld1(i,j,k-1,bi,bj) - fld1(i,j,k,bi,bj) )
163     & *( fld2(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)
164     & +fld2( i ,j,k,bi,bj)*rA( i ,j,bi,bj)
165     & )*recip_rAw(i,j,bi,bj)
166 heimbach 1.7 & * deltaTloc
167 jmc 1.4 ENDDO
168     ENDDO
169     ENDDO
170    
171     ELSEIF ( dir.eq.-23 ) THEN
172     c- gradient of the 1rst field * 2nd fields, shifted by -dR/2 & -dY/2 resp.
173     c- (e.g.: used for advective form of vertical advection: w.dv/dr)
174    
175     DO k=2,Ksize
176     DO j=1,sNy
177     DO i=1,sNx
178     fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
179     & + .5 _d 0*( fld1(i,j,k-1,bi,bj) - fld1(i,j,k,bi,bj) )
180     & *( fld2(i,j-1,k,bi,bj)*rA(i,j-1,bi,bj)
181     & +fld2(i, j ,k,bi,bj)*rA(i, j ,bi,bj)
182     & )*recip_rAs(i,j,bi,bj)
183 heimbach 1.7 & * deltaTloc
184 jmc 1.3 ENDDO
185     ENDDO
186     ENDDO
187 jmc 1.1
188     ENDIF
189 jmc 1.3
190     #endif /* ALLOW_TIMEAVE */
191 jmc 1.1
192     RETURN
193     END
  ViewVC Help
Powered by ViewVC 1.1.22