327 |
DO j = 1,sNy |
DO j = 1,sNy |
328 |
DO i = 1,sNx |
DO i = 1,sNx |
329 |
tmpMk(i,j,k,bi,bj) = tmp1k(i,j,bi,bj) |
tmpMk(i,j,k,bi,bj) = tmp1k(i,j,bi,bj) |
330 |
& + (rF(k)-rC(k))*hFacC(i,j,k,bi,bj) |
& + (rF(k+1)-rC(k))*hFacC(i,j,k,bi,bj)*rkSign |
331 |
C above: more general (setInterFDr/setCenterDr) than line below |
C above: more general (setInterFDr/setCenterDr) than line below |
332 |
c & + drF(k)*hFacC(i,j,k,bi,bj)*0.5 _d 0 |
c & + drF(k)*hFacC(i,j,k,bi,bj)*0.5 _d 0 |
333 |
tmp1k(i,j,bi,bj) = tmp1k(i,j,bi,bj) |
tmp1k(i,j,bi,bj) = tmp1k(i,j,bi,bj) |