103 |
do i = imin,imax |
do i = imin,imax |
104 |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
105 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
106 |
cph gtNm1(i,j,k,bi,bj) = gtNm1(i,j,k,bi,bj) + |
if(theta(i,j,k,bi,bj).lt.-2.0) |
107 |
cph & fac*tmpfld3d(i,j,k,bi,bj) |
& theta(i,j,k,bi,bj)= -2.0 |
108 |
enddo |
enddo |
109 |
enddo |
enddo |
110 |
enddo |
enddo |
128 |
do i = imin,imax |
do i = imin,imax |
129 |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
130 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
|
cph gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) + |
|
|
cph & fac*tmpfld3d(i,j,k,bi,bj) |
|
131 |
enddo |
enddo |
132 |
enddo |
enddo |
133 |
enddo |
enddo |
151 |
do i = imin,imax |
do i = imin,imax |
152 |
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
153 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
|
cph gtr1Nm1(i,j,k,bi,bj) = gtr1Nm1(i,j,k,bi,bj) + |
|
|
cph & fac*tmpfld3d(i,j,k,bi,bj) |
|
154 |
enddo |
enddo |
155 |
enddo |
enddo |
156 |
enddo |
enddo |