/[MITgcm]/MITgcm/model/src/calc_gt.F
ViewVC logotype

Diff of /MITgcm/model/src/calc_gt.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by cnh, Sat May 30 02:10:16 1998 UTC revision 1.9 by adcroft, Mon Jun 1 22:27:14 1998 UTC
# Line 6  CStartOfInterFace Line 6  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,
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     /==========================================================\
# Line 72  C     myThid - Instance number for this Line 72  C     myThid - Instance number for this
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)
# Line 198  C           boundary condition. Line 198  C           boundary condition.
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
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22