5 |
CStartOfInterFace |
CStartOfInterFace |
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,maskC, |
9 |
I K13,K23,KappaZT,KapGM, |
I K13,K23,KappaZT,KapGM, |
10 |
U af,df,fZon,fMer,fVerT, |
U af,df,fZon,fMer,fVerT, |
11 |
I myThid ) |
I myThid ) |
52 |
C fVerT - Flux of temperature (T) in the vertical |
C fVerT - Flux of temperature (T) in the vertical |
53 |
C direction at the upper(U) and lower(D) faces of a cell. |
C direction at the upper(U) and lower(D) faces of a cell. |
54 |
C maskUp - Land mask used to denote base of the domain. |
C maskUp - Land mask used to denote base of the domain. |
55 |
|
C maskC - Land mask for theta cells (used in TOP_LAYER only) |
56 |
C xA - Tracer cell face area normal to X |
C xA - Tracer cell face area normal to X |
57 |
C yA - Tracer cell face area normal to X |
C yA - Tracer cell face area normal to X |
58 |
C uTrans - Zonal volume transport through cell face |
C uTrans - Zonal volume transport through cell face |
72 |
_RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
73 |
_RL wTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL wTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
74 |
_RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
75 |
|
_RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
76 |
_RL K13 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL K13 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
77 |
_RL K23 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL K23 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
78 |
_RL KappaZT(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL KappaZT(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
260 |
DO j=jMin,jMax |
DO j=jMin,jMax |
261 |
DO i=iMin,iMax |
DO i=iMin,iMax |
262 |
gT(i,j,k,bi,bj)=gT(i,j,k,bi,bj) |
gT(i,j,k,bi,bj)=gT(i,j,k,bi,bj) |
263 |
& +maskUp(i,j)*( |
& +maskC(i,j)*( |
264 |
& -lambdaThetaClimRelax*(theta(i,j,k,bi,bj)-SST(i,j,bi,bj)) |
& -lambdaThetaClimRelax*(theta(i,j,k,bi,bj)-SST(i,j,bi,bj)) |
265 |
& -Qnet(i,j,bi,bj) ) |
& -Qnet(i,j,bi,bj) ) |
266 |
ENDDO |
ENDDO |