135 |
|
|
136 |
C-- Unless we have already calculated the advection terms we initialize |
C-- Unless we have already calculated the advection terms we initialize |
137 |
C the tendency to zero. |
C the tendency to zero. |
138 |
IF (calcAdvection) THEN |
C <== now done earlier at the beginning of thermodynamics. |
139 |
DO j=1-Oly,sNy+Oly |
c IF (calcAdvection) THEN |
140 |
DO i=1-Olx,sNx+Olx |
c DO j=1-Oly,sNy+Oly |
141 |
gTracer(i,j,k,bi,bj)=0. _d 0 |
c DO i=1-Olx,sNx+Olx |
142 |
ENDDO |
c gTracer(i,j,k,bi,bj)=0. _d 0 |
143 |
ENDDO |
c ENDDO |
144 |
ENDIF |
c ENDDO |
145 |
|
c ENDIF |
146 |
|
|
147 |
C-- Pre-calculate del^2 T if bi-harmonic coefficient is non-zero |
C-- Pre-calculate del^2 T if bi-harmonic coefficient is non-zero |
148 |
IF (diffK4 .NE. 0.) THEN |
IF (diffK4 .NE. 0.) THEN |