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

Diff of /MITgcm/pkg/mom_vecinv/mom_vi_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:33 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_CORIOLIS(
7         I        bi,bj,K,
8         I        uFld,vFld,omega3,r_hFacZ,
9         O        uCoriolisTerm,vCoriolisTerm,
10         I        myThid)
11          IMPLICIT NONE
12    C     /==========================================================\
13    C     | S/R VORTICITY_X_V                                        |
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 vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
27          _RL omega3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28          _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29          _RL uCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
30          _RL vCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
31          INTEGER myThid
32    
33    C     == Local variables ==
34          INTEGER I,J
35          _RL vBarXY,uBarXY,vort3u,vort3v
36    
37           DO J=1-Oly,sNy+Oly-1
38            DO I=2-Olx,sNx+Olx
39    c        vBarXY=0.25*(vFld( i ,j)+vFld( i ,j+1)
40    c    &               +vFld(i-1,j)+vFld(i-1,j+1))
41    c        vort3u=0.5*(omega3(i,j)+omega3(I,J+1))
42             vBarXY=0.25*(
43         &       vFld( i , j )*dxG( i , j ,bi,bj)*hFacS( i , j ,k,bi,bj)
44         &      +vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacS( i ,j+1,k,bi,bj)
45         &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*hFacS(i-1, j ,k,bi,bj)
46         &      +vFld(i-1,j+1)*dxG(i-1,j+1,bi,bj)*hFacS(i-1,j+1,k,bi,bj))
47             vort3u=0.5*(omega3(i,j)*r_hFacZ(i,j)
48         &              +omega3(i,j+1)*r_hFacZ(i,j+1))
49             uCoriolisTerm(i,j)=
50    C high order vorticity advection term
51    c    &   +vort3u*vBarXY*recip_dxc(i,j,bi,bj)
52    C linear Coriolis term (enstrophy conserving)
53         &   +0.5*(fCoriG(i,j,bi,bj)*r_hFacZ(i,j)
54         &        +fCoriG(i,j+1,bi,bj)*r_hFacZ(i,j+1))*vBarXY
55         *        *recip_dxc(i,j,bi,bj)
56    C full nonlinear Coriolis term
57    c    &   +0.5*(omega3(i,j)+omega3(i,j+1))*vBarXY
58    C correct energy conserving form of Coriolis term
59    c    &   +0.25*( fCori( i ,j,bi,bj)*(vFld( i ,j)+vFld( i ,j+1)) +
60    c    &           fCori(i-1,j,bi,bj)*(vFld(i-1,j)+vFld(i-1,j+1)) )
61    C original form of Coriolis term (copied from calc_mom_rhs)
62    c    &   +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY
63         &   *_maskW(I,J,K,bi,bj)
64            ENDDO
65           ENDDO
66    
67           DO J=2-Oly,sNy+Oly
68            DO I=1-Olx,sNx+Olx-1
69    c        uBarXY=0.25*( uFld(i, j )+uFld(i+1, j )
70    c    &                +uFld(i,j-1)+uFld(i+1,j-1))
71    c        vort3v=0.5*(omega3(i,j)+omega3(I+1,J))
72             uBarXY=0.25*(
73         &       uFld( i , j )*dyG( i , j ,bi,bj)*hFacW( i , j ,k,bi,bj)
74         &      +uFld(i+1, j )*dyG(i+1, j ,bi,bj)*hFacW(i+1, j ,k,bi,bj)
75         &      +uFld( i ,j-1)*dyG( i ,j-1,bi,bj)*hFacW( i ,j-1,k,bi,bj)
76         &      +uFld(i+1,j-1)*dyG(i+1,j-1,bi,bj)*hFacW(i+1,j-1,k,bi,bj))
77             vort3v=0.5*(omega3(i,j)*r_hFacZ(i,j)
78         &              +omega3(i+1,j)*r_hFacZ(i+1,j))
79             vCoriolisTerm(i,j)=
80    C high order vorticity advection term
81    c    &   -vort3v*uBarXY*recip_dyc(i,j,bi,bj)
82    C linear Coriolis term (enstrophy conserving)
83         &   -0.5 *(fCoriG(i,j,bi,bj)*r_hFacZ(i,j)
84         &         +fCoriG(i+1,j,bi,bj)*r_hFacZ(i+1,j))*uBarXY
85         &         *recip_dyc(i,j,bi,bj)
86    C full nonlinear Coriolis term
87    c    &   -0.5*(omega3(i,j)+omega3(i+1,j))*uBarXY
88    C correct energy conserving form of Coriolis term
89    c    &   -0.25*( fCori(i,j  ,bi,bj)*(uFld(i, j )+uFld(i+1,j)) +
90    c    &           fCori(i,j-1,bi,bj)*(uFld(i,j-1)+uFld(i+1,j)) )
91    C original form of Coriolis term (copied from calc_mom_rhs)
92    c    &   -0.5*(fCori(i,j,bi,bj)+fCori(i,j-1,bi,bj))*uBarXY
93         &   *_maskS(I,J,K,bi,bj)
94            ENDDO
95           ENDDO
96    
97    
98          RETURN
99          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22