148 |
cgg exactly balance the volume flux, since we are dealing with |
cgg exactly balance the volume flux, since we are dealing with |
149 |
cgg the baroclinic velocity structure.. |
cgg the baroclinic velocity structure.. |
150 |
utop = tmpfldyz(j,k,bi,bj)* |
utop = tmpfldyz(j,k,bi,bj)* |
151 |
& maskW(i+ip1,j,k,bi,bj) * delZ(k) + utop |
& maskW(i+ip1,j,k,bi,bj) * delR(k) + utop |
152 |
cgg Add the barotropic velocity component. |
cgg Add the barotropic velocity component. |
153 |
if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then |
if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then |
154 |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
156 |
enddo |
enddo |
157 |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
158 |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
159 |
& - utop / delZ(1) |
& - utop / delR(1) |
160 |
enddo |
enddo |
161 |
enddo |
enddo |
162 |
enddo |
enddo |
183 |
cgg exactly balance the volume flux, since we are dealing with |
cgg exactly balance the volume flux, since we are dealing with |
184 |
cgg the baroclinic velocity structure.. |
cgg the baroclinic velocity structure.. |
185 |
utop = tmpfldyz(j,k,bi,bj)* |
utop = tmpfldyz(j,k,bi,bj)* |
186 |
& maskS(i,j,k,bi,bj) * delZ(k) + utop |
& maskS(i,j,k,bi,bj) * delR(k) + utop |
187 |
cgg Add the barotropic velocity component. |
cgg Add the barotropic velocity component. |
188 |
if (maskS(i,j,k,bi,bj) .ne. 0.) then |
if (maskS(i,j,k,bi,bj) .ne. 0.) then |
189 |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
191 |
enddo |
enddo |
192 |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
193 |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
194 |
& - utop / delZ(1) |
& - utop / delR(1) |
195 |
enddo |
enddo |
196 |
enddo |
enddo |
197 |
enddo |
enddo |
265 |
cgg exactly balance the volume flux, since we are dealing with |
cgg exactly balance the volume flux, since we are dealing with |
266 |
cgg the baroclinic velocity structure.. |
cgg the baroclinic velocity structure.. |
267 |
utop = tmpfldyz(j,k,bi,bj)* |
utop = tmpfldyz(j,k,bi,bj)* |
268 |
& maskW(i+ip1,j,k,bi,bj) * delZ(k) + utop |
& maskW(i+ip1,j,k,bi,bj) * delR(k) + utop |
269 |
cgg Add the barotropic velocity component. |
cgg Add the barotropic velocity component. |
270 |
if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then |
if (maskW(i+ip1,j,k,bi,bj) .ne. 0.) then |
271 |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
273 |
enddo |
enddo |
274 |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
275 |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
276 |
& - utop / delZ(1) |
& - utop / delR(1) |
277 |
enddo |
enddo |
278 |
enddo |
enddo |
279 |
enddo |
enddo |
300 |
cgg exactly balance the volume flux, since we are dealing with |
cgg exactly balance the volume flux, since we are dealing with |
301 |
cgg the baroclinic velocity structure.. |
cgg the baroclinic velocity structure.. |
302 |
utop = tmpfldyz(j,k,bi,bj)* |
utop = tmpfldyz(j,k,bi,bj)* |
303 |
& maskS(i,j,k,bi,bj) * delZ(k) + utop |
& maskS(i,j,k,bi,bj) * delR(k) + utop |
304 |
cgg Add the barotropic velocity component. |
cgg Add the barotropic velocity component. |
305 |
if (maskS(i,j,k,bi,bj) .ne. 0.) then |
if (maskS(i,j,k,bi,bj) .ne. 0.) then |
306 |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
tmpfldyz(j,k,bi,bj) = tmpfldyz(j,k,bi,bj)+ ubaro |
308 |
enddo |
enddo |
309 |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
cgg Compute the baroclinic velocity at level 1. Should balance flux. |
310 |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
tmpfldyz(j,1,bi,bj) = tmpfldyz(j,1,bi,bj) |
311 |
& - utop / delZ(1) |
& - utop / delR(1) |
312 |
enddo |
enddo |
313 |
enddo |
enddo |
314 |
enddo |
enddo |