1 |
C $Header$ |
C $Header$ |
2 |
|
C $Name$ |
3 |
|
|
4 |
#include "GMREDI_OPTIONS.h" |
#include "GMREDI_OPTIONS.h" |
5 |
|
|
41 |
C I, J - Loop counters |
C I, J - Loop counters |
42 |
INTEGER I, J |
INTEGER I, J |
43 |
|
|
44 |
IF (use_GMRedi) THEN |
IF (useGMRedi) THEN |
45 |
|
|
46 |
C-- Area integrated meridional flux |
C-- Area integrated meridional flux |
47 |
DO j=jMin,jMax |
DO j=jMin,jMax |
50 |
& -yA(i,j) |
& -yA(i,j) |
51 |
& *(GM_background_K |
& *(GM_background_K |
52 |
#ifdef GM_VISBECK_VARIABLE_K |
#ifdef GM_VISBECK_VARIABLE_K |
53 |
& +0.5*(VisbeckK(i,j,myThid)+VisbeckK(i,j-1,myThid)) |
& +0.5*(VisbeckK(i,j,bi,bj)+VisbeckK(i,j-1,bi,bj)) |
54 |
#endif |
#endif |
55 |
& ) |
& ) |
56 |
#ifdef GM_NON_UNITY_DIAGONAL |
#ifdef GM_NON_UNITY_DIAGONAL |
57 |
& *Kvy(i,j,k,myThid) |
& *Kvy(i,j,k,bi,bj) |
58 |
#endif |
#endif |
59 |
& *_recip_dyC(i,j,bi,bj) |
& *_recip_dyC(i,j,bi,bj) |
60 |
& *(Tracer(i,j,k,bi,bj)-Tracer(i,j-1,k,bi,bj)) |
& *(Tracer(i,j,k,bi,bj)-Tracer(i,j-1,k,bi,bj)) |