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

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

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

revision 1.122 by jmc, Sat Jul 30 23:39:48 2005 UTC revision 1.125 by jmc, Fri Sep 16 21:40:58 2005 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "PACKAGES_CONFIG.h"  #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    #undef DYNAMICS_GUGV_EXCH_CHECK
7    
8  CBOP  CBOP
9  C     !ROUTINE: DYNAMICS  C     !ROUTINE: DYNAMICS
# Line 219  C         (1 + dt * K * d_zz) salt[n] = Line 220  C         (1 + dt * K * d_zz) salt[n] =
220  C---  C---
221  CEOP  CEOP
222    
223    #ifdef ALLOW_DEBUG
224          IF ( debugLevel .GE. debLevB )
225         &   CALL DEBUG_ENTER( 'DYNAMICS', myThid )
226    #endif
227    
228  C-- Call to routine for calculation of  C-- Call to routine for calculation of
229  C   Eliassen-Palm-flux-forced U-tendency,  C   Eliassen-Palm-flux-forced U-tendency,
230  C   if desired:  C   if desired:
# Line 428  CADJ STORE gU(:,:,:,bi,bj) = comlev1_bib Line 434  CADJ STORE gU(:,:,:,bi,bj) = comlev1_bib
434  #endif    /* ALLOW_AUTODIFF_TAMC */  #endif    /* ALLOW_AUTODIFF_TAMC */
435            CALL IMPLDIFF(            CALL IMPLDIFF(
436       I         bi, bj, iMin, iMax, jMin, jMax,       I         bi, bj, iMin, iMax, jMin, jMax,
437       I         0, KappaRU,recip_HFacW,       I         -1, KappaRU,recip_HFacW,
438       U         gU,       U         gU,
439       I         myThid )       I         myThid )
440  #ifdef    ALLOW_AUTODIFF_TAMC  #ifdef    ALLOW_AUTODIFF_TAMC
# Line 437  CADJ STORE gV(:,:,:,bi,bj) = comlev1_bib Line 443  CADJ STORE gV(:,:,:,bi,bj) = comlev1_bib
443  #endif    /* ALLOW_AUTODIFF_TAMC */  #endif    /* ALLOW_AUTODIFF_TAMC */
444            CALL IMPLDIFF(            CALL IMPLDIFF(
445       I         bi, bj, iMin, iMax, jMin, jMax,       I         bi, bj, iMin, iMax, jMin, jMax,
446       I         0, KappaRV,recip_HFacS,       I         -2, KappaRV,recip_HFacS,
447       U         gV,       U         gV,
448       I         myThid )       I         myThid )
449          ENDIF          ENDIF
# Line 488  C---+----1----+----2----+----3----+----4 Line 494  C---+----1----+----2----+----3----+----4
494  C--   Step forward W field in N-H algorithm  C--   Step forward W field in N-H algorithm
495        IF ( momStepping .AND. nonHydrostatic ) THEN        IF ( momStepping .AND. nonHydrostatic ) THEN
496  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
497            IF ( debugLevel .GE. debLevB )           IF ( debugLevel .GE. debLevB )
498       &     CALL DEBUG_CALL('CALC_GW',myThid)       &     CALL DEBUG_CALL('CALC_GW', myThid )
499  #endif  #endif
500           CALL TIMER_START('CALC_GW          [DYNAMICS]',myThid)           CALL TIMER_START('CALC_GW          [DYNAMICS]',myThid)
501           CALL CALC_GW( myTime, myIter, myThid )           CALL CALC_GW( myTime, myIter, myThid )
# Line 544  Cml) Line 550  Cml)
550        ENDIF        ENDIF
551  #endif  #endif
552    
553    #ifdef DYNAMICS_GUGV_EXCH_CHECK
554    C- jmc: For safety checking only: This Exchange here should not change
555    C       the solution. If solution changes, it means something is wrong,
556    C       but it does not mean that it is less wrong with this exchange.
557          IF ( debugLevel .GT. debLevB ) THEN
558           CALL EXCH_UV_XYZ_RL(gU,gV,.TRUE.,myThid)
559          ENDIF
560    #endif
561    
562    #ifdef ALLOW_DEBUG
563          IF ( debugLevel .GE. debLevB )
564         &   CALL DEBUG_LEAVE( 'DYNAMICS', myThid )
565    #endif
566    
567        RETURN        RETURN
568        END        END
Legend:
Removed from v.1.122  
changed lines
  Added in v.1.125
  ViewVC Help
Powered by ViewVC 1.1.22