20 |
C dir - type of grid for 2nd array (0: tracer, 1: zonal, 2 merid) |
C dir - type of grid for 2nd array (0: tracer, 1: zonal, 2 merid) |
21 |
C fld - Input Field |
C fld - Input Field |
22 |
C fldtave - time averaged Field |
C fldtave - time averaged Field |
23 |
_RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy) |
_RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
24 |
_RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy) |
_RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
25 |
_RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz,nSx,nSy) |
_RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
26 |
INTEGER dir |
INTEGER dir |
27 |
REAL deltaT |
REAL deltaT |
28 |
INTEGER bi, bj, K, myThid |
INTEGER bi, bj, K, myThid |
37 |
|
|
38 |
C DO bj = myByLo(myThid), myByHi(myThid) |
C DO bj = myByLo(myThid), myByHi(myThid) |
39 |
C DO bi = myBxLo(myThid), myBxHi(myThid) |
C DO bi = myBxLo(myThid), myBxHi(myThid) |
40 |
C DO k=1,Nz |
C DO k=1,Nr |
41 |
DO j=1,sNy |
DO j=1,sNy |
42 |
DO i=1,sNx |
DO i=1,sNx |
43 |
fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj) |
fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj) |
52 |
|
|
53 |
C DO bj = myByLo(myThid), myByHi(myThid) |
C DO bj = myByLo(myThid), myByHi(myThid) |
54 |
C DO bi = myBxLo(myThid), myBxHi(myThid) |
C DO bi = myBxLo(myThid), myBxHi(myThid) |
55 |
C DO k=1,Nz |
C DO k=1,Nr |
56 |
DO j=1,sNy |
DO j=1,sNy |
57 |
DO i=1,sNx |
DO i=1,sNx |
58 |
fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj) |
fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj) |
69 |
|
|
70 |
C DO bj = myByLo(myThid), myByHi(myThid) |
C DO bj = myByLo(myThid), myByHi(myThid) |
71 |
C DO bi = myBxLo(myThid), myBxHi(myThid) |
C DO bi = myBxLo(myThid), myBxHi(myThid) |
72 |
C DO k=1,Nz |
C DO k=1,Nr |
73 |
DO j=1,sNy |
DO j=1,sNy |
74 |
DO i=1,sNx |
DO i=1,sNx |
75 |
fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj) |
fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj) |