56 |
IF (tempStepping) THEN |
IF (tempStepping) THEN |
57 |
DO j=jMin,jMax |
DO j=jMin,jMax |
58 |
DO i=iMin,iMax |
DO i=iMin,iMax |
59 |
theta(i,j,k,bi,bj)=theta(i,j,k,bi,bj) |
gTNm1(i,j,k,bi,bj)=theta(i,j,k,bi,bj) |
60 |
& +deltaTtracer*(ab15*gT(i,j,k,bi,bj)+ab05*gTNm1(i,j,k,bi,bj)) |
& +deltaTtracer*(ab15*gT(i,j,k,bi,bj)+ab05*gTNm1(i,j,k,bi,bj)) |
|
gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj) |
|
61 |
ENDDO |
ENDDO |
62 |
ENDDO |
ENDDO |
63 |
ENDIF |
ENDIF |
66 |
IF (saltStepping) THEN |
IF (saltStepping) THEN |
67 |
DO j=jMin,jMax |
DO j=jMin,jMax |
68 |
DO i=iMin,iMax |
DO i=iMin,iMax |
69 |
salt(i,j,k,bi,bj)=salt(i,j,k,bi,bj) |
gSNm1(i,j,k,bi,bj)=salt(i,j,k,bi,bj) |
70 |
& +deltaTtracer*(ab15*gS(i,j,k,bi,bj)+ab05*gSNm1(i,j,k,bi,bj)) |
& +deltaTtracer*(ab15*gS(i,j,k,bi,bj)+ab05*gSNm1(i,j,k,bi,bj)) |
|
gSNm1(i,j,k,bi,bj)=gS(i,j,k,bi,bj) |
|
71 |
ENDDO |
ENDDO |
72 |
ENDDO |
ENDDO |
73 |
ENDIF |
ENDIF |