58 |
C == Local variables == |
C == Local variables == |
59 |
C I, J - Loop counters |
C I, J - Loop counters |
60 |
INTEGER I, J |
INTEGER I, J |
61 |
_RL tmpFac |
_RL tmpFac, tmpFac1 |
62 |
|
|
63 |
|
|
64 |
|
#ifdef ALLOW_SALT_PLUME |
65 |
|
IF ( useSALT_PLUME ) THEN |
66 |
|
tmpFac1 = 1. _d 0 |
67 |
|
ELSE |
68 |
|
tmpFac1 = 0. _d 0 |
69 |
|
ENDIF |
70 |
|
#endif /* ALLOW_SALT_PLUME */ |
71 |
|
|
72 |
#ifdef ALLOW_GMREDI |
#ifdef ALLOW_GMREDI |
73 |
|
|
83 |
& *KPPghat(i,j,km1,bi,bj) |
& *KPPghat(i,j,km1,bi,bj) |
84 |
& *( surfaceForcingS(i,j,bi,bj) |
& *( surfaceForcingS(i,j,bi,bj) |
85 |
#ifdef ALLOW_SALT_PLUME |
#ifdef ALLOW_SALT_PLUME |
86 |
& + saltPlumeFlux(i,j,bi,bj)*recip_rhoConst |
& + tmpFac1*saltPlumeFlux(i,j,bi,bj)*recip_rhoConst |
87 |
& *(1.-KPPplumefrac(i,j,bi,bj)) |
& *(1.-KPPplumefrac(i,j,bi,bj)) |
88 |
#endif /* ALLOW_SALT_PLUME */ |
#endif /* ALLOW_SALT_PLUME */ |
89 |
& ) |
& ) |
99 |
& *KPPghat(i,j,km1,bi,bj) |
& *KPPghat(i,j,km1,bi,bj) |
100 |
& *( surfaceForcingS(i,j,bi,bj) |
& *( surfaceForcingS(i,j,bi,bj) |
101 |
#ifdef ALLOW_SALT_PLUME |
#ifdef ALLOW_SALT_PLUME |
102 |
& + saltPlumeFlux(i,j,bi,bj)*recip_rhoConst |
& + tmpFac1*saltPlumeFlux(i,j,bi,bj)*recip_rhoConst |
103 |
& *(1.-KPPplumefrac(i,j,bi,bj)) |
& *(1.-KPPplumefrac(i,j,bi,bj)) |
104 |
#endif /* ALLOW_SALT_PLUME */ |
#endif /* ALLOW_SALT_PLUME */ |
105 |
& ) |
& ) |