80 & LBi, UBi, LBj, UBj, &
81 & IminS, ImaxS, JminS, JmaxS, &
87 & pm, pmon_p, pmon_r, &
88 & pn, pnom_p, pnom_r, &
97 integer,
intent(in) :: ng, tile
98 integer,
intent(in) :: LBi, UBi, LBj, UBj
99 integer,
intent(in) :: IminS, ImaxS, JminS, JmaxS
100 integer,
intent(in) :: nrhs, nnew
104 real(r8),
intent(in) :: pmask(LBi:,LBj:)
106 real(r8),
intent(in) :: Hz(LBi:,LBj:,:)
107 real(r8),
intent(in) :: pm(LBi:,LBj:)
108 real(r8),
intent(in) :: pmon_p(LBi:,LBj:)
109 real(r8),
intent(in) :: pmon_r(LBi:,LBj:)
110 real(r8),
intent(in) :: pn(LBi:,LBj:)
111 real(r8),
intent(in) :: pnom_p(LBi:,LBj:)
112 real(r8),
intent(in) :: pnom_r(LBi:,LBj:)
114 real(r8),
intent(inout) :: tl_u(LBi:,LBj:,:,:)
115 real(r8),
intent(inout) :: tl_v(LBi:,LBj:,:,:)
118 real(r8),
intent(in) :: pmask(LBi:UBi,LBj:UBj)
120 real(r8),
intent(in) :: Hz(LBi:UBi,LBj:UBj,N(ng))
121 real(r8),
intent(in) :: pm(LBi:UBi,LBj:UBj)
122 real(r8),
intent(in) :: pmon_p(LBi:UBi,LBj:UBj)
123 real(r8),
intent(in) :: pmon_r(LBi:UBi,LBj:UBj)
124 real(r8),
intent(in) :: pn(LBi:UBi,LBj:UBj)
125 real(r8),
intent(in) :: pnom_p(LBi:UBi,LBj:UBj)
126 real(r8),
intent(in) :: pnom_r(LBi:UBi,LBj:UBj)
128 real(r8),
intent(inout) :: tl_u(LBi:UBi,LBj:UBj,N(ng),2)
129 real(r8),
intent(inout) :: tl_v(LBi:UBi,LBj:UBj,N(ng),2)
138 real(r8),
dimension(IminS:ImaxS,JminS:JmaxS) :: UFe
139 real(r8),
dimension(IminS:ImaxS,JminS:JmaxS) :: UFx
140 real(r8),
dimension(IminS:ImaxS,JminS:JmaxS) :: VFe
141 real(r8),
dimension(IminS:ImaxS,JminS:JmaxS) :: VFx
143# include "set_bounds.h"
158 k_loop :
DO k=1,n(ng)
167 & (tl_u(i+1,j,k,nrhs)- &
174 & 0.25_r8*(hz(i,j ,k)+hz(i-1,j ,k)+ &
175 & hz(i,j-1,k)+hz(i-1,j-1,k))* &
176 & (tl_u(i,j ,k,nrhs)- &
177 & tl_u(i,j-1,k,nrhs))
179 ufe(i,j)=ufe(i,j)*pmask(i,j)
186 & 0.25_r8*(hz(i,j ,k)+hz(i-1,j ,k)+ &
187 & hz(i,j-1,k)+hz(i-1,j-1,k))* &
188 & (tl_v(i ,j,k,nrhs)- &
189 & tl_v(i-1,j,k,nrhs))
191 vfx(i,j)=vfx(i,j)*pmask(i,j)
199 & (tl_v(i,j+1,k,nrhs)- &
210 tl_u(i,j,k,nnew)=tl_u(i,j,k,nnew)+ &
211 & cff*(pm(i-1,j)+pm(i,j))* &
212 & (pn(i-1,j)+pn(i,j))* &
213 & (ufx(i,j)-ufx(i-1,j)+ &
214 & ufe(i,j+1)-ufe(i,j))
220 tl_v(i,j,k,nnew)=tl_v(i,j,k,nnew)+ &
221 & cff*(pm(i,j)+pm(i,j-1))* &
222 & (pn(i,j)+pn(i,j-1))* &
223 & (vfx(i+1,j)-vfx(i,j)+ &
224 & vfe(i,j)-vfe(i,j-1))
subroutine tl_uv3drelax_tile(ng, tile, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, nrhs, nnew, pmask, hz, pm, pmon_p, pmon_r, pn, pnom_p, pnom_r, tl_u, tl_v)