/[MITgcm]/MITgcm_contrib/bbl/code/mypackage_calc_rhs.F
ViewVC logotype

Diff of /MITgcm_contrib/bbl/code/mypackage_calc_rhs.F

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

revision 1.1 by dimitri, Thu Nov 18 04:00:05 2010 UTC revision 1.2 by dimitri, Tue Nov 23 07:06:25 2010 UTC
# Line 58  C     fMer      :: Meridional flux Line 58  C     fMer      :: Meridional flux
58        _RL     bbl_rho ( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly )        _RL     bbl_rho ( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly )
59        _RL     fZon    ( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly )        _RL     fZon    ( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly )
60        _RL     fMer    ( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly )        _RL     fMer    ( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly )
61          CHARACTER*(MAX_LEN_MBUF) msgBuf
62  CEOP  CEOP
63    
64  C     Initialize tendency terms and  C     Initialize tendency terms and
# Line 197  C     for botommost wet grid box. Line 198  C     for botommost wet grid box.
198         ENDDO         ENDDO
199        ENDDO        ENDDO
200    
201    #ifdef ALLOW_DEBUG
202          IF ( debugLevel .GE. debLevB ) THEN
203    C     Check salinity conservation
204           bbl_tend=0
205           DO j=1,sNy
206            DO i=1,sNx
207             kBot = kLowC(i,j,bi,bj)
208             IF ( kBot .GT. 0 ) THEN
209              bbl_tend = bbl_tend + bbl_TendSalt(i,j,bi,bj) *
210         &         hFacC(i,j,kBot,bi,bj) * drF(kBot) *rA(i,j,bi,bj)
211             ENDIF
212            ENDDO
213           ENDDO
214           _GLOBAL_SUM_RL( bbl_tend, myThid )
215           WRITE(msgBuf,'(A,E10.2)') 'total salt tendency = ', bbl_tend
216           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
217         &      SQUEEZE_RIGHT, myThid )
218          ENDIF
219    #endif /* ALLOW_DEBUG */
220    
221    
222        CALL EXCH_XY_RL( bbl_theta, myThid )        CALL EXCH_XY_RL( bbl_theta, myThid )
223        CALL EXCH_XY_RL( bbl_salt , myThid )        CALL EXCH_XY_RL( bbl_salt , myThid )
224    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22