1 |
C $Header: /u/gcmpack/models/MITgcmUV/pkg/mom_vecinv/Attic/mom_vi_calc_absvort3.F,v 1.1.2.2 2001/04/09 18:36:45 adcroft Exp $ |
2 |
C $Name: pre38-close $ |
3 |
|
4 |
#include "CPP_OPTIONS.h" |
5 |
|
6 |
SUBROUTINE MOM_VI_CALC_ABSVORT3( |
7 |
I bi,bj,k, |
8 |
I vort3, |
9 |
O omega3, |
10 |
I myThid) |
11 |
IMPLICIT NONE |
12 |
C /==========================================================\ |
13 |
C | S/R MOM_CALC_ABSVORT3 | |
14 |
C |==========================================================| |
15 |
C \==========================================================/ |
16 |
|
17 |
C == Global variables == |
18 |
#include "SIZE.h" |
19 |
#include "EEPARAMS.h" |
20 |
#include "PARAMS.h" |
21 |
#include "GRID.h" |
22 |
C == Routine arguments == |
23 |
C myThid - Instance number for this innvocation of CALC_MOM_RHS |
24 |
INTEGER bi,bj,k |
25 |
_RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
26 |
_RL omega3(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
27 |
INTEGER myThid |
28 |
|
29 |
C == Local variables == |
30 |
INTEGER i,j |
31 |
_RL nonLinFac,useCoriolisFac |
32 |
|
33 |
nonLinFac=0. |
34 |
IF (momAdvection) nonLinFac=0. |
35 |
useCoriolisFac=0. |
36 |
IF (useCoriolis) useCoriolisFac=1. |
37 |
|
38 |
DO J=1-Oly,sNy+Oly |
39 |
DO I=1-Olx,sNx+Olx |
40 |
omega3(I,J)= |
41 |
& fCoriG(I,J,bi,bj)*useCoriolisFac |
42 |
& +vort3(i,j)*nonLinFac |
43 |
ENDDO |
44 |
ENDDO |
45 |
|
46 |
RETURN |
47 |
END |