140 |
INTEGER nipass,ipass |
INTEGER nipass,ipass |
141 |
INTEGER myTile, nCFace |
INTEGER myTile, nCFace |
142 |
LOGICAL N_edge, S_edge, E_edge, W_edge |
LOGICAL N_edge, S_edge, E_edge, W_edge |
143 |
|
#ifdef ALLOW_DIAGNOSTICS |
144 |
|
CHARACTER*8 diagName |
145 |
|
CHARACTER*4 GAD_DIAG_SUFX, diagSufx |
146 |
|
EXTERNAL GAD_DIAG_SUFX |
147 |
|
#endif |
148 |
CEOP |
CEOP |
149 |
|
|
150 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
168 |
endif |
endif |
169 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
170 |
|
|
171 |
|
#ifdef ALLOW_DIAGNOSTICS |
172 |
|
C-- Set diagnostic suffix for the current tracer |
173 |
|
IF ( useDiagnostics ) THEN |
174 |
|
diagSufx = GAD_DIAG_SUFX( tracerIdentity, myThid ) |
175 |
|
ENDIF |
176 |
|
#endif |
177 |
|
|
178 |
C-- Set up work arrays with valid (i.e. not NaN) values |
C-- Set up work arrays with valid (i.e. not NaN) values |
179 |
C These inital values do not alter the numerical results. They |
C These inital values do not alter the numerical results. They |
180 |
C just ensure that all memory references are to valid floating |
C just ensure that all memory references are to valid floating |
602 |
ENDDO |
ENDDO |
603 |
ENDIF |
ENDIF |
604 |
|
|
605 |
|
#ifdef ALLOW_DIAGNOSTICS |
606 |
|
IF ( useDiagnostics ) THEN |
607 |
|
diagName = 'ADVx'//diagSufx |
608 |
|
CALL DIAGNOSTICS_FILL(afx,diagName, k,1, 2,bi,bj, myThid) |
609 |
|
diagName = 'ADVy'//diagSufx |
610 |
|
CALL DIAGNOSTICS_FILL(afy,diagName, k,1, 2,bi,bj, myThid) |
611 |
|
ENDIF |
612 |
|
#endif |
613 |
|
|
614 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
615 |
IF ( debugLevel .GE. debLevB |
IF ( debugLevel .GE. debLevB |
616 |
& .AND. tracerIdentity.EQ.GAD_TEMPERATURE |
& .AND. tracerIdentity.EQ.GAD_TEMPERATURE |
728 |
ENDDO |
ENDDO |
729 |
ENDDO |
ENDDO |
730 |
|
|
731 |
|
#ifdef ALLOW_DIAGNOSTICS |
732 |
|
IF ( useDiagnostics ) THEN |
733 |
|
diagName = 'ADVr'//diagSufx |
734 |
|
CALL DIAGNOSTICS_FILL( fVerT(1-Olx,1-Oly,kUp), |
735 |
|
& diagName, k,1, 2,bi,bj, myThid) |
736 |
|
ENDIF |
737 |
|
#endif |
738 |
|
|
739 |
C-- End of K loop for vertical flux |
C-- End of K loop for vertical flux |
740 |
ENDDO |
ENDDO |
741 |
C-- end of if not.implicitAdvection block |
C-- end of if not.implicitAdvection block |