/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vecinv.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vecinv.F

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

revision 1.11 by edhill, Tue Nov 4 19:51:54 2003 UTC revision 1.13 by dimitri, Sun Jan 25 00:31:52 2004 UTC
# Line 426  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K Line 426  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K
426         ENDDO         ENDDO
427    
428  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
429    #ifndef HRCUBE
430         IF (taveFreq.GT.0.) THEN         IF (taveFreq.GT.0.) THEN
431           CALL TIMEAVE_CUMUL_1K1T(uZetatave,vCf,deltaTClock,           CALL TIMEAVE_CUMUL_1K1T(uZetatave,vCf,deltaTClock,
432       &                           Nr, k, bi, bj, myThid)       &                           Nr, k, bi, bj, myThid)
433           CALL TIMEAVE_CUMUL_1K1T(vZetatave,uCf,deltaTClock,           CALL TIMEAVE_CUMUL_1K1T(vZetatave,uCf,deltaTClock,
434       &                           Nr, k, bi, bj, myThid)       &                           Nr, k, bi, bj, myThid)
435         ENDIF         ENDIF
436  #endif  #endif /* ALLOW_TIMEAVE */
437    #endif /* ndef HRCUBE */
438    
439  C--   Vertical shear terms (-w*du/dr & -w*dv/dr)  C--   Vertical shear terms (-w*du/dr & -w*dv/dr)
440         CALL MOM_VI_U_VERTSHEAR(bi,bj,K,uVel,wVel,uCf,myThid)         IF ( .NOT. momImplVertAdv ) THEN
441         DO j=jMin,jMax          CALL MOM_VI_U_VERTSHEAR(bi,bj,K,uVel,wVel,uCf,myThid)
442          DO i=iMin,iMax          DO j=jMin,jMax
443           gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+uCf(i,j)           DO i=iMin,iMax
444              gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+uCf(i,j)
445             ENDDO
446          ENDDO          ENDDO
447         ENDDO          CALL MOM_VI_V_VERTSHEAR(bi,bj,K,vVel,wVel,vCf,myThid)
448         CALL MOM_VI_V_VERTSHEAR(bi,bj,K,vVel,wVel,vCf,myThid)          DO j=jMin,jMax
449         DO j=jMin,jMax           DO i=iMin,iMax
450          DO i=iMin,iMax            gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+vCf(i,j)
451           gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+vCf(i,j)           ENDDO
452          ENDDO          ENDDO
453         ENDDO         ENDIF
454    
455  C--   Bernoulli term  C--   Bernoulli term
456         CALL MOM_VI_U_GRAD_KE(bi,bj,K,KE,uCf,myThid)         CALL MOM_VI_U_GRAD_KE(bi,bj,K,KE,uCf,myThid)
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22