forked from altMITgcm/MITgcm
-
Notifications
You must be signed in to change notification settings - Fork 0
/
mom_vi_u_grad_ke.F
39 lines (32 loc) · 1.05 KB
/
mom_vi_u_grad_ke.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
C $Header: /u/gcmpack/MITgcm/pkg/mom_vecinv/mom_vi_u_grad_ke.F,v 1.4 2015/09/10 18:08:51 jmc Exp $
C $Name: $
#include "MOM_VECINV_OPTIONS.h"
SUBROUTINE MOM_VI_U_GRAD_KE(
I bi,bj,k,
I KE,
O dKEdx,
I myThid)
IMPLICIT NONE
C /==========================================================\
C | S/R MOM_U_GRAD_KE |
C |==========================================================|
C \==========================================================/
C == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "GRID.h"
C == Routine arguments ==
INTEGER bi,bj,k
_RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
_RL dKEdx(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
INTEGER myThid
C == Local variables ==
INTEGER i,j
DO j=1-OLy,sNy+OLy
DO i=2-OLx,sNx+OLx
dKEdx(i,j)=-_recip_dxC(i,j,bi,bj)*(KE(i,j)-KE(i-1,j))
& *_maskW(i,j,k,bi,bj)*recip_deepFacC(k)
ENDDO
ENDDO
RETURN
END