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

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

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

revision 1.1 by adcroft, Fri Mar 30 21:13:34 2001 UTC revision 1.2 by adcroft, Tue May 29 14:01:39 2001 UTC
# Line 0  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "CPP_OPTIONS.h"
5    
6          SUBROUTINE MOM_VI_V_CORIOLIS(
7         I        bi,bj,K,
8         I        uFld,omega3,r_hFacZ,
9         O        vCoriolisTerm,
10         I        myThid)
11          IMPLICIT NONE
12    C     /==========================================================\
13    C     | S/R MOM_VI_V_CORIOLIS                                    |
14    C     |==========================================================|
15    C     \==========================================================/
16    
17    C     == Global variables ==
18    #include "SIZE.h"
19    #include "EEPARAMS.h"
20    #include "GRID.h"
21    #include "PARAMS.h"
22    
23    C     == Routine arguments ==
24          INTEGER bi,bj,K
25          _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
26          _RL omega3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
27          _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28          _RL vCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29          INTEGER myThid
30    
31    C     == Local variables ==
32          INTEGER I,J
33          _RL uBarXY,vort3v
34          LOGICAL upwindVort3
35    
36          upwindVort3=.FALSE.
37    
38           DO J=2-Oly,sNy+Oly
39            DO I=1-Olx,sNx+Olx-1
40             uBarXY=0.25*(
41         &       uFld( i , j )*dyG( i , j ,bi,bj)*hFacW( i , j ,k,bi,bj)
42         &      +uFld(i+1, j )*dyG(i+1, j ,bi,bj)*hFacW(i+1, j ,k,bi,bj)
43         &      +uFld( i ,j-1)*dyG( i ,j-1,bi,bj)*hFacW( i ,j-1,k,bi,bj)
44         &      +uFld(i+1,j-1)*dyG(i+1,j-1,bi,bj)*hFacW(i+1,j-1,k,bi,bj))
45    c        uBarXY=0.25*( uFld(i, j )+uFld(i+1, j )
46    c    &                +uFld(i,j-1)+uFld(i+1,j-1))
47             IF (upwindVort3) THEN
48              IF (uBarXY.GT.0.) THEN
49               vort3v=omega3(i,j)*r_hFacZ(i,j)
50              ELSE
51               vort3v=omega3(i+1,j)*r_hFacZ(i+1,j)
52              ENDIF
53             ELSE
54              vort3v=0.5*(omega3(i,j)*r_hFacZ(i,j)
55         &               +omega3(i+1,j)*r_hFacZ(i+1,j))
56             ENDIF
57             vCoriolisTerm(i,j)=
58    C high order vorticity advection term
59         &   -vort3v*uBarXY*recip_dyc(i,j,bi,bj)
60    C linear Coriolis term
61    c    &   -0.5 *(fCoriG(I,J,bi,bj)+fCoriG(I+1,J,bi,bj))*uBarXY
62    C full nonlinear Coriolis term
63    c    &   -0.5*(omega3(I,J)+omega3(I+1,J))*uBarXY
64    C correct energy conserving form of Coriolis term
65    c    &   -0.5 *( fCori(I,J  ,bi,bj)*uBarX(I,J  ,K,bi,bj) +
66    c    &           fCori(I,J-1,bi,bj)*uBarX(I,J-1,K,bi,bj)  )
67    C original form of Coriolis term (copied from calc_mom_rhs)
68    c    &   -0.5*(fCori(i,j,bi,bj)+fCori(i,j-1,bi,bj))*uBarXY
69         &   *_maskS(I,J,K,bi,bj)
70            ENDDO
71           ENDDO
72    
73          RETURN
74          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22