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

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

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

revision 1.3 by heimbach, Wed Sep 5 17:46:03 2001 UTC revision 1.4 by jmc, Sun Aug 3 02:49:35 2003 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6        SUBROUTINE MOM_VI_U_CORIOLIS(        SUBROUTINE MOM_VI_U_CORIOLIS(
7       I        bi,bj,K,       I        bi,bj,k,
8       I        vFld,omega3,r_hFacZ,       I        vFld,omega3,hFacZ,r_hFacZ,
9       O        uCoriolisTerm,       O        uCoriolisTerm,
10       I        myThid)       I        myThid)
11        IMPLICIT NONE        IMPLICIT NONE
12  C     /==========================================================\  C     *==========================================================*
13  C     | S/R MOM_VI_U_CORIOLIS                                    |  C     | S/R MOM_VI_U_CORIOLIS
14  C     |==========================================================|  C     | o Calculate meridional flux of vorticity at U point
15  C     \==========================================================/  C     *==========================================================*
16    
17  C     == Global variables ==  C     == Global variables ==
18  #include "SIZE.h"  #include "SIZE.h"
# Line 24  C     == Routine arguments == Line 24  C     == Routine arguments ==
24        INTEGER bi,bj,K        INTEGER bi,bj,K
25        _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
26        _RL omega3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL omega3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
27          _RS hFacZ (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28        _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29        _RL uCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
30        INTEGER myThid        INTEGER myThid
31    
32  C     == Local variables ==  C     == Local variables ==
33          LOGICAL upwindVort3, use_original_hFac
34        INTEGER I,J        INTEGER I,J
35        _RL vBarXY,vort3u        _RL vBarXY,vort3u
36        LOGICAL upwindVort3        _RS epsil
37          PARAMETER ( upwindVort3      =.FALSE. )
38        upwindVort3=.FALSE.        PARAMETER ( use_original_hFac=.FALSE. )
39    
40         DO J=1-Oly,sNy+Oly-1        epsil = 1. _d -9
41          DO I=2-Olx,sNx+Olx  
42          DO J=1-Oly,sNy+Oly-1
43           DO I=2-Olx,sNx+Olx
44            IF ( use_original_hFac ) THEN
45           vBarXY=0.25*(           vBarXY=0.25*(
46       &       vFld( i , j )*dxG( i , j ,bi,bj)*hFacS( i , j ,k,bi,bj)       &       vFld( i , j )*dxG( i , j ,bi,bj)*hFacS( i , j ,k,bi,bj)
47       &      +vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacS( i ,j+1,k,bi,bj)       &      +vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacS( i ,j+1,k,bi,bj)
48       &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*hFacS(i-1, j ,k,bi,bj)       &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*hFacS(i-1, j ,k,bi,bj)
49       &      +vFld(i-1,j+1)*dxG(i-1,j+1,bi,bj)*hFacS(i-1,j+1,k,bi,bj))       &      +vFld(i-1,j+1)*dxG(i-1,j+1,bi,bj)*hFacS(i-1,j+1,k,bi,bj))
 c        vBarXY=0.25*(vFld( i ,j)+vFld( i ,j+1)  
 c    &               +vFld(i-1,j)+vFld(i-1,j+1))  
50           IF (upwindVort3) THEN           IF (upwindVort3) THEN
51            IF (vBarXY.GT.0.) THEN            IF (vBarXY.GT.0.) THEN
52             vort3u=omega3(I,J)*r_hFacZ(i,j)             vort3u=omega3(I,J)*r_hFacZ(i,j)
# Line 51  c    &               +vFld(i-1,j)+vFld(i Line 54  c    &               +vFld(i-1,j)+vFld(i
54             vort3u=omega3(I,J+1)*r_hFacZ(i,j+1)             vort3u=omega3(I,J+1)*r_hFacZ(i,j+1)
55            ENDIF            ENDIF
56           ELSE           ELSE
57            vort3u=0.5*(omega3(i,j)*r_hFacZ(i,j)             vort3u=0.5*(omega3(i,j)*r_hFacZ(i,j)
58       &               +omega3(i,j+1)*r_hFacZ(i,j+1))       &                +omega3(i,j+1)*r_hFacZ(i,j+1))
59           ENDIF           ENDIF
60           uCoriolisTerm(i,j)=          ELSE
61       &   +vort3u*vBarXY*recip_dxc(i,j,bi,bj)  c--      test a different formulation (relatively to hFac)
62       &   *_maskW(I,J,K,bi,bj)           vBarXY=0.5*(
63         &       vFld( i , j )*dxG( i , j ,bi,bj)*hFacZ(i,j)
64         &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*hFacZ(i,j)
65         &      +vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacZ(i,j+1)
66         &      +vFld(i-1,j+1)*dxG(i-1,j+1,bi,bj)*hFacZ(i,j+1)
67         &              )/MAX( epsil, hFacZ(i,j)+hFacZ(i,j+1) )
68             IF (upwindVort3) THEN
69              IF (vBarXY.GT.0.) THEN
70               vort3u=omega3(i,j)
71              ELSE
72               vort3u=omega3(i,j+1)
73              ENDIF
74             ELSE
75               vort3u=0.5*(omega3(i,j)+omega3(i,j+1))
76             ENDIF
77            ENDIF
78    
79            uCoriolisTerm(i,j)=
80         &   +vort3u*vBarXY*recip_dxC(i,j,bi,bj)*_maskW(i,j,k,bi,bj)
81  cph *note* put these comments after end of continued line  cph *note* put these comments after end of continued line
82  cph        to ensure TAMC compatibility  cph        to ensure TAMC compatibility
83  C high order vorticity advection term  C high order vorticity advection term
# Line 70  c    &   +0.5*( fCori(I  ,J,bi,bj)*vBarY Line 91  c    &   +0.5*( fCori(I  ,J,bi,bj)*vBarY
91  c    &          fCori(I-1,J,bi,bj)*vBarY(I-1,J,K,bi,bj)  )  c    &          fCori(I-1,J,bi,bj)*vBarY(I-1,J,K,bi,bj)  )
92  C original form of Coriolis term (copied from calc_mom_rhs)  C original form of Coriolis term (copied from calc_mom_rhs)
93  c    &   +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY  c    &   +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY
         ENDDO  
94         ENDDO         ENDDO
95          ENDDO
96    
97        RETURN        RETURN
98        END        END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22