/[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.12 by jmc, Sat Jan 3 00:51:42 2004 UTC
# Line 435  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K Line 435  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K
435  #endif  #endif
436    
437  C--   Vertical shear terms (-w*du/dr & -w*dv/dr)  C--   Vertical shear terms (-w*du/dr & -w*dv/dr)
438         CALL MOM_VI_U_VERTSHEAR(bi,bj,K,uVel,wVel,uCf,myThid)         IF ( .NOT. momImplVertAdv ) THEN
439         DO j=jMin,jMax          CALL MOM_VI_U_VERTSHEAR(bi,bj,K,uVel,wVel,uCf,myThid)
440          DO i=iMin,iMax          DO j=jMin,jMax
441           gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+uCf(i,j)           DO i=iMin,iMax
442              gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+uCf(i,j)
443             ENDDO
444          ENDDO          ENDDO
445         ENDDO          CALL MOM_VI_V_VERTSHEAR(bi,bj,K,vVel,wVel,vCf,myThid)
446         CALL MOM_VI_V_VERTSHEAR(bi,bj,K,vVel,wVel,vCf,myThid)          DO j=jMin,jMax
447         DO j=jMin,jMax           DO i=iMin,iMax
448          DO i=iMin,iMax            gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+vCf(i,j)
449           gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+vCf(i,j)           ENDDO
450          ENDDO          ENDDO
451         ENDDO         ENDIF
452    
453  C--   Bernoulli term  C--   Bernoulli term
454         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.12
  ViewVC Help
Powered by ViewVC 1.1.22