107 |
do i = imin,imax |
do i = imin,imax |
108 |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
109 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
110 |
gtNm1(i,j,k,bi,bj) = gtNm1(i,j,k,bi,bj) + |
cph gtNm1(i,j,k,bi,bj) = gtNm1(i,j,k,bi,bj) + |
111 |
& fac*tmpfld3d(i,j,k,bi,bj) |
cph & fac*tmpfld3d(i,j,k,bi,bj) |
112 |
enddo |
enddo |
113 |
enddo |
enddo |
114 |
enddo |
enddo |
132 |
do i = imin,imax |
do i = imin,imax |
133 |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
134 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
135 |
gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) + |
cph gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) + |
136 |
& fac*tmpfld3d(i,j,k,bi,bj) |
cph & fac*tmpfld3d(i,j,k,bi,bj) |
137 |
enddo |
enddo |
138 |
enddo |
enddo |
139 |
enddo |
enddo |
157 |
do i = imin,imax |
do i = imin,imax |
158 |
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
159 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
160 |
gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) + |
cph gtr1Nm1(i,j,k,bi,bj) = gtr1Nm1(i,j,k,bi,bj) + |
161 |
& fac*tmpfld3d(i,j,k,bi,bj) |
cph & fac*tmpfld3d(i,j,k,bi,bj) |
162 |
enddo |
enddo |
163 |
enddo |
enddo |
164 |
enddo |
enddo |