/[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.9 by adcroft, Mon Jun 1 22:27:14 1998 UTC revision 1.10 by cnh, Mon Jun 8 21:43:00 1998 UTC
# Line 84  CEndOfInterface Line 84  CEndOfInterface
84  C     == Local variables ==  C     == Local variables ==
85  C     I, J, K - Loop counters  C     I, J, K - Loop counters
86        INTEGER i,j        INTEGER i,j
87          LOGICAL TOP_LAYER
88        _RL afFacT, dfFacT        _RL afFacT, dfFacT
89        _RL dTdx(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL dTdx(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
90        _RL dTdy(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL dTdy(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
91    
92        afFacT = 1. _d 0        afFacT = 1. _d 0
93        dfFacT = 1. _d 0        dfFacT = 1. _d 0
94          TOP_LAYER = K .EQ. 1
95    
96  C---  Calculate advective and diffusive fluxes between cells.  C---  Calculate advective and diffusive fluxes between cells.
97    
# Line 216  C           boundary condition. Line 218  C           boundary condition.
218  C     Net vertical flux  C     Net vertical flux
219        DO j=jMin,jMax        DO j=jMin,jMax
220         DO i=iMin,iMax         DO i=iMin,iMax
221          fVerT(i,j,kUp) = (afFacT*af(i,j) + dfFacT*df(i,j))*maskUp(i,j)          fVerT(i,j,kUp) = ( afFacT*af(i,j)+  dfFacT*df(i,j) )*maskUp(i,j)
222         ENDDO         ENDDO
223        ENDDO        ENDDO
224          IF ( TOP_LAYER ) THEN
225           DO j=jMin,jMax
226            DO i=iMin,iMax
227             fVerT(i,j,kUp) = afFacT*af(i,j)*freeSurfFac
228            ENDDO
229           ENDDO
230          ENDIF
231    
232  C--   Tendency is minus divergence of the fluxes.  C--   Tendency is minus divergence of the fluxes.
233  C     Note. Tendency terms will only be correct for range  C     Note. Tendency terms will only be correct for range
# Line 228  C           they are not algorithmically Line 237  C           they are not algorithmically
237  C           are not used.  C           are not used.
238        DO j=jMin,jMax        DO j=jMin,jMax
239         DO i=iMin,iMax         DO i=iMin,iMax
240    C    &   -_rhFacC(i,j,k,bi,bj)*rdzF(k)*_rdxF(i,j,bi,bj)*_rdyF(i,j,bi,bj)
241    C    &   -_rhFacC(i,j,k,bi,bj)*rdzF(k)/_zA(i,j,bi,bj)
242    C #define _rVolT(i,j,k,bi,bj) _rhFacC(i,j,k,bi,bj)*rdzF(k)*_rdxF(i,j,bi,bj)*_rdyF(i,j,bi,bj)
243    #define _rVolT(i,j,k,bi,bj) _rhFacC(i,j,k,bi,bj)*rdzF(k)/_zA(i,j,bi,bj)
244          gT(i,j,k,bi,bj)=          gT(i,j,k,bi,bj)=
245       &   -_rhFacC(i,j,k,bi,bj)*rdzF(k)*_rdxF(i,j,bi,bj)*_rdyF(i,j,bi,bj)       &   -_rVolT(i,j,k,bi,bj)
246       &   *(       &   *(
247       &    +( fZon(i+1,j)-fZon(i,j) )       &    +( fZon(i+1,j)-fZon(i,j) )
248       &    +( fMer(i,j+1)-fMer(i,j) )       &    +( fMer(i,j+1)-fMer(i,j) )
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22