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

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

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

revision 1.38 by jmc, Mon Mar 24 00:32:53 2008 UTC revision 1.39 by jmc, Tue Apr 22 22:20:31 2008 UTC
# Line 103  C     i,j,k         :: Loop counters Line 103  C     i,j,k         :: Loop counters
103        PARAMETER( iMin = 1 , iMax = sNx )        PARAMETER( iMin = 1 , iMax = sNx )
104        PARAMETER( jMin = 1 , jMax = sNy )        PARAMETER( jMin = 1 , jMax = sNy )
105  CEOP  CEOP
106    #ifdef ALLOW_DIAGNOSTICS
107          LOGICAL diagDiss, diagAdvec
108          LOGICAL  DIAGNOSTICS_IS_ON
109          EXTERNAL DIAGNOSTICS_IS_ON
110    #endif /* ALLOW_DIAGNOSTICS */
111    
112  C Catch barotropic mode  C--   Catch barotropic mode
113        IF ( Nr .LT. 2 ) RETURN        IF ( Nr .LT. 2 ) RETURN
114    
115    #ifdef ALLOW_DIAGNOSTICS
116          IF ( useDiagnostics ) THEN
117            diagDiss  = DIAGNOSTICS_IS_ON( 'Wm_Diss ', myThid )
118            diagAdvec = DIAGNOSTICS_IS_ON( 'Wm_Advec', myThid )
119          ELSE
120            diagDiss  = .FALSE.
121            diagAdvec = .FALSE.
122          ENDIF
123    #endif /* ALLOW_DIAGNOSTICS */
124    
125  C--   Initialise gW to zero  C--   Initialise gW to zero
126        DO k=1,Nr        DO k=1,Nr
127          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
# Line 488  C--          prepare for next level (k+1 Line 503  C--          prepare for next level (k+1
503    
504  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
505    
506    #ifdef ALLOW_DIAGNOSTICS
507          IF ( diagDiss )  CALL DIAGNOSTICS_FILL( gwDiss, 'Wm_Diss ',
508         &                                        k, 1, 2, bi,bj, myThid )
509          IF ( diagAdvec ) CALL DIAGNOSTICS_FILL( gW, 'Wm_Advec',
510         &                                        k,Nr, 1, bi,bj, myThid )
511    #endif /* ALLOW_DIAGNOSTICS */
512    
513  C--   Dissipation term inside the Adams-Bashforth:  C--   Dissipation term inside the Adams-Bashforth:
514          IF ( momViscosity .AND. momDissip_In_AB) THEN          IF ( momViscosity .AND. momDissip_In_AB) THEN
515            DO j=jMin,jMax            DO j=jMin,jMax
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.39
  ViewVC Help
Powered by ViewVC 1.1.22