-
Notifications
You must be signed in to change notification settings - Fork 237
/
mom_calc_absvort3.F
50 lines (44 loc) · 1.24 KB
/
mom_calc_absvort3.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
#include "MOM_COMMON_OPTIONS.h"
SUBROUTINE MOM_CALC_ABSVORT3(
I bi,bj,k,
I vort3,
O omega3,
I myThid)
IMPLICIT NONE
C /==========================================================\
C | S/R MOM_CALC_ABSVORT3 |
C |==========================================================|
C \==========================================================/
C == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
C == Routine arguments ==
C myThid - Instance number for this innvocation of CALC_MOM_RHS
INTEGER bi,bj,k
_RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
_RL omega3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
INTEGER myThid
C == Local variables ==
INTEGER i,j
_RL nonLinFac,useCoriolisFac
IF (momAdvection) THEN
nonLinFac=1.
ELSE
nonLinFac=0.
ENDIF
IF (useCoriolis) THEN
useCoriolisFac=1.
ELSE
useCoriolisFac=0.
ENDIF
DO j=1-Oly,sNy+Oly
DO i=1-Olx,sNx+Olx
omega3(I,J)=
& fCoriG(i,j,bi,bj)*useCoriolisFac
& +vort3(i,j)*nonLinFac
ENDDO
ENDDO
RETURN
END