6 |
SUBROUTINE CALC_GT( |
SUBROUTINE CALC_GT( |
7 |
I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown, |
I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown, |
8 |
I xA,yA,uTrans,vTrans,wTrans,maskup, |
I xA,yA,uTrans,vTrans,wTrans,maskup, |
9 |
I K13,K23,K33,KapGM, |
I K13,K23,KappaZT,KapGM, |
10 |
U af,df,fZon,fMer,fVerT, |
U af,df,fZon,fMer,fVerT, |
11 |
I myThid ) |
I myThid ) |
12 |
C /==========================================================\ |
C /==========================================================\ |
72 |
_RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
73 |
_RL K13 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL K13 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
74 |
_RL K23 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL K23 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
75 |
_RL K33 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL KappaZT(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
76 |
_RL KapGM (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL KapGM (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
77 |
_RL af (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL af (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
78 |
_RL df (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL df (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
198 |
DO j=jMin,jMax |
DO j=jMin,jMax |
199 |
DO i=iMin,iMax |
DO i=iMin,iMax |
200 |
df(i,j) = _zA(i,j,bi,bj)*( |
df(i,j) = _zA(i,j,bi,bj)*( |
|
& -(diffKzT+KapGM(i,j)*K33(i,j,k))*rdzC(k) |
|
|
& *(theta(i,j,kM1,bi,bj)-theta(i,j,k,bi,bj)) |
|
201 |
& -KapGM(i,j)*K13(i,j,k)*dTdx(i,j) |
& -KapGM(i,j)*K13(i,j,k)*dTdx(i,j) |
202 |
& -KapGM(i,j)*K23(i,j,k)*dTdy(i,j) |
& -KapGM(i,j)*K23(i,j,k)*dTdy(i,j) |
203 |
& ) |
& ) |
204 |
ENDDO |
ENDDO |
205 |
ENDDO |
ENDDO |
206 |
|
IF (.NOT.implicitDiffusion) THEN |
207 |
|
DO j=jMin,jMax |
208 |
|
DO i=iMin,iMax |
209 |
|
df(i,j) = df(i,j) + _zA(i,j,bi,bj)*( |
210 |
|
& -KappaZT(i,j,k)*rdzC(k) |
211 |
|
& *(theta(i,j,kM1,bi,bj)-theta(i,j,k,bi,bj)) |
212 |
|
& ) |
213 |
|
ENDDO |
214 |
|
ENDDO |
215 |
|
ENDIF |
216 |
C Net vertical flux |
C Net vertical flux |
217 |
DO j=jMin,jMax |
DO j=jMin,jMax |
218 |
DO i=iMin,iMax |
DO i=iMin,iMax |