/[MITgcm]/MITgcm/pkg/mom_common/mom_w_coriolis_nh.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_w_coriolis_nh.F

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

revision 1.1 by jmc, Mon Mar 20 14:40:31 2006 UTC revision 1.2 by jmc, Mon Mar 12 23:54:26 2007 UTC
# Line 7  CBOP Line 7  CBOP
7  C !ROUTINE: MOM_W_CORIOLIS_NH  C !ROUTINE: MOM_W_CORIOLIS_NH
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE MOM_W_CORIOLIS_NH(        SUBROUTINE MOM_W_CORIOLIS_NH(
11       I               bi,bj,k,       I               bi,bj,k,
12       I               uFld, vFld,       I               uFld, vFld,
13       U               wCoriolisTerm,       U               wCoriolisTerm,
14       I               myThid)       I               myThid )
15    
16  C !DESCRIPTION:  C !DESCRIPTION:
17  C Calculates the Coriolis term in the vertical momentum equation:  C Calculates the Coriolis term in the vertical momentum equation:
# Line 52  C Energy conserving discretization of 2* Line 52  C Energy conserving discretization of 2*
52        IF ( k.GT.1 .AND. k.LE.Nr ) THEN        IF ( k.GT.1 .AND. k.LE.Nr ) THEN
53          DO j=1-Oly,sNy+Oly-1          DO j=1-Oly,sNy+Oly-1
54           DO i=1-Olx,sNx+Olx-1           DO i=1-Olx,sNx+Olx-1
55             wCoriolisTerm(i,j) =             wCoriolisTerm(i,j) =
56       &      -gravitySign*fCoriCos(i,j,bi,bj)*       &      -gravitySign*fCoriCos(i,j,bi,bj)*
57       &       ( angleCosC(i,j,bi,bj)*(       &       ( angleCosC(i,j,bi,bj)*(
58       &                (uFld(i,j,k-1,bi,bj)+uFld(i+1,j,k-1,bi,bj))       &                (uFld(i,j,k-1,bi,bj)+uFld(i+1,j,k-1,bi,bj))
# Line 62  C Energy conserving discretization of 2* Line 62  C Energy conserving discretization of 2*
62       &                (vFld(i,j,k-1,bi,bj)+vFld(i,j+1,k-1,bi,bj))       &                (vFld(i,j,k-1,bi,bj)+vFld(i,j+1,k-1,bi,bj))
63       &              + (vFld(i,j, k ,bi,bj)+vFld(i,j+1, k ,bi,bj))       &              + (vFld(i,j, k ,bi,bj)+vFld(i,j+1, k ,bi,bj))
64       &                              )*0.25 _d 0       &                              )*0.25 _d 0
65       &       )*horiVertRatio       &       )*wUnit2rVel(k)
66           ENDDO           ENDDO
67          ENDDO          ENDDO
68        ELSE        ELSE
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22