309 |
|
|
310 |
C-- Compute vertical flux fVerT(kUp) at interface k (between k-1 & k): |
C-- Compute vertical flux fVerT(kUp) at interface k (between k-1 & k): |
311 |
C- Advective flux in R |
C- Advective flux in R |
312 |
|
#ifdef ALLOW_AIM |
313 |
|
C- a hack to prevent Water-Vapor vert.transport into the stratospheric level Nr |
314 |
|
IF (calcAdvection .AND. .NOT.implicitAdvection .AND. K.GE.2 .AND. |
315 |
|
& (.NOT.useAIM .OR.tracerIdentity.NE.GAD_SALINITY .OR.K.LT.Nr) |
316 |
|
& ) THEN |
317 |
|
#else |
318 |
IF (calcAdvection .AND. .NOT.implicitAdvection .AND. K.GE.2) THEN |
IF (calcAdvection .AND. .NOT.implicitAdvection .AND. K.GE.2) THEN |
319 |
|
#endif |
320 |
C- Compute vertical advective flux in the interior: |
C- Compute vertical advective flux in the interior: |
321 |
IF (advectionScheme.EQ.ENUM_CENTERED_2ND) THEN |
IF (advectionScheme.EQ.ENUM_CENTERED_2ND) THEN |
322 |
CALL GAD_C2_ADV_R(bi,bj,k,rTrans,tracer,af,myThid) |
CALL GAD_C2_ADV_R(bi,bj,k,rTrans,tracer,af,myThid) |