39 |
& )*_maskW(i,j,k,bi,bj) |
& )*_maskW(i,j,k,bi,bj) |
40 |
ENDDO |
ENDDO |
41 |
ENDDO |
ENDDO |
42 |
|
|
43 |
C Step forward meridional velocity (store in Gv) |
C Step forward meridional velocity (store in Gv) |
44 |
DO j=jMin,jMax |
DO j=jMin,jMax |
45 |
DO i=iMin,iMax |
DO i=iMin,iMax |
51 |
& )*_maskS(i,j,k,bi,bj) |
& )*_maskS(i,j,k,bi,bj) |
52 |
ENDDO |
ENDDO |
53 |
ENDDO |
ENDDO |
54 |
|
|
55 |
C Step forward temperature |
C Step forward temperature |
56 |
|
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) |
theta(i,j,k,bi,bj)=theta(i,j,k,bi,bj) |
61 |
gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj) |
gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj) |
62 |
ENDDO |
ENDDO |
63 |
ENDDO |
ENDDO |
64 |
|
ENDIF |
65 |
|
|
66 |
|
C Step forward salt |
67 |
|
IF (saltStepping) THEN |
68 |
|
DO j=jMin,jMax |
69 |
|
DO i=iMin,iMax |
70 |
|
salt(i,j,k,bi,bj)=salt(i,j,k,bi,bj) |
71 |
|
& +deltaTtracer*(ab15*gS(i,j,k,bi,bj)+ab05*gSNm1(i,j,k,bi,bj)) |
72 |
|
gSNm1(i,j,k,bi,bj)=gS(i,j,k,bi,bj) |
73 |
|
ENDDO |
74 |
|
ENDDO |
75 |
|
ENDIF |
76 |
|
|
77 |
RETURN |
RETURN |
78 |
END |
END |