9 |
I tracer, |
I tracer, |
10 |
O vT, |
O vT, |
11 |
I myThid ) |
I myThid ) |
12 |
|
C /==========================================================\ |
13 |
|
C | SUBROUTINE GAD_FLUXLIMIT_ADV_Y | |
14 |
|
C | o Compute Meridional advective Flux of Tracer using | |
15 |
|
C | Flux Limiter Scheme | |
16 |
|
C |==========================================================| |
17 |
IMPLICIT NONE |
IMPLICIT NONE |
18 |
|
|
19 |
C == GLobal variables == |
C == GLobal variables == |
36 |
|
|
37 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
38 |
vT(i,1-Oly)=0. |
vT(i,1-Oly)=0. |
39 |
|
vT(i,2-Oly)=0. |
40 |
|
vT(i,sNy+Oly)=0. |
41 |
ENDDO |
ENDDO |
42 |
DO j=1-Oly+1,sNy+Oly |
DO j=1-Oly+2,sNy+Oly-1 |
43 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
44 |
Rjp=(tracer(i,j+1)-tracer(i,j))*maskS(i,j+1,k,bi,bj) |
Rjp=(tracer(i,j+1)-tracer(i,j))*maskS(i,j+1,k,bi,bj) |
45 |
Rj=(tracer(i,j)-tracer(i,j-1))*maskS(i,j,k,bi,bj) |
Rj=(tracer(i,j)-tracer(i,j-1))*maskS(i,j,k,bi,bj) |