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

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

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

revision 1.1 by adcroft, Wed May 9 13:59:24 2001 UTC revision 1.1.2.1 by adcroft, Wed May 9 13:59:24 2001 UTC
# Line 0  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "CPP_OPTIONS.h"
5    
6          SUBROUTINE MOM_VI_U_CORIOLIS_C4(
7         I        bi,bj,K,
8         I        vFld,omega3,r_hFacZ,
9         O        uCoriolisTerm,
10         I        myThid)
11          IMPLICIT NONE
12    C     /==========================================================\
13    C     | S/R MOM_VI_U_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 vFld(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 uCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29          INTEGER myThid
30    
31    C     == Local variables ==
32          INTEGER I,J
33          _RL vBarXY,vort3u,Rjp,Rjm
34          LOGICAL upwindVort3
35          LOGICAL fourthVort3
36    
37          upwindVort3=.FALSE.
38          fourthVort3=.TRUE.
39    
40           DO J=1-Oly,sNy+Oly-1
41            DO I=2-Olx,sNx+Olx
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    c        vBarXY=0.25*(vFld( i ,j)+vFld( i ,j+1)
48    c    &               +vFld(i-1,j)+vFld(i-1,j+1))
49             IF (upwindVort3) THEN
50              IF (vBarXY.GT.0.) THEN
51               vort3u=omega3(I,J)*r_hFacZ(i,j)
52              ELSE
53               vort3u=omega3(I,J+1)*r_hFacZ(i,j+1)
54              ENDIF
55             ELSEIF (fourthVort3) THEN
56              Rjp=omega3(i,j+2)*r_hFacZ(i,j+2)
57         &       -omega3(i,j+1)*r_hFacZ(i,j+1)
58              Rjm=omega3(i,j)*r_hFacZ(i,j)
59         &       -omega3(i,j-1)*r_hFacZ(i,j-1)
60              vort3u=0.5*(omega3(i,j)*r_hFacZ(i,j)
61         &               +omega3(i,j+1)*r_hFacZ(i,j+1)
62         &               -1./12.*(Rjp-Rjm)
63         &               )
64             ELSE
65              vort3u=0.5*(omega3(i,j)*r_hFacZ(i,j)
66         &               +omega3(i,j+1)*r_hFacZ(i,j+1))
67             ENDIF
68             uCoriolisTerm(i,j)=
69    C high order vorticity advection term
70         &   +vort3u*vBarXY*recip_dxc(i,j,bi,bj)
71    C linear Coriolis term
72    c    &   +0.5*(fCoriG(I,J,bi,bj)+fCoriG(I,J+1,bi,bj))*vBarXY
73    C full nonlinear Coriolis term
74    c    &   +0.5*(omega3(I,J)+omega3(I,J+1))*vBarXY
75    C correct energy conserving form of Coriolis term
76    c    &   +0.5*( fCori(I  ,J,bi,bj)*vBarY(I  ,J,K,bi,bj) +
77    c    &          fCori(I-1,J,bi,bj)*vBarY(I-1,J,K,bi,bj)  )
78    C original form of Coriolis term (copied from calc_mom_rhs)
79    c    &   +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY
80         &   *_maskW(I,J,K,bi,bj)
81            ENDDO
82           ENDDO
83    
84          RETURN
85          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.2.1
  ViewVC Help
Powered by ViewVC 1.1.22