136 |
DO j=jMin,jMax |
DO j=jMin,jMax |
137 |
DO i=iMin,iMax |
DO i=iMin,iMax |
138 |
dPhiHydX(i,j) = dPhiHydX(i,j) |
dPhiHydX(i,j) = dPhiHydX(i,j) |
139 |
& +factorP*(tFld(i-1,j,k,bi,bj)+tFld(i,j,k,bi,bj)) |
& +factorP*(tFld(i-1,j,k,bi,bj)+tFld(i,j,k,bi,bj)-2.*tRef(k)) |
140 |
& *(rStarFacC(i,j,bi,bj)-rStarFacC(i-1,j,bi,bj)) |
& *(rStarFacC(i,j,bi,bj)-rStarFacC(i-1,j,bi,bj)) |
141 |
& *rC(k)*recip_dxC(i,j,bi,bj) |
& *rC(k)*recip_dxC(i,j,bi,bj) |
142 |
dPhiHydY(i,j) = dPhiHydY(i,j) |
dPhiHydY(i,j) = dPhiHydY(i,j) |
143 |
& +factorP*(tFld(i,j-1,k,bi,bj)+tFld(i,j,k,bi,bj)) |
& +factorP*(tFld(i,j-1,k,bi,bj)+tFld(i,j,k,bi,bj)-2.*tRef(k)) |
144 |
& *(rStarFacC(i,j,bi,bj)-rStarFacC(i,j-1,bi,bj)) |
& *(rStarFacC(i,j,bi,bj)-rStarFacC(i,j-1,bi,bj)) |
145 |
& *rC(k)*recip_dyC(i,j,bi,bj) |
& *rC(k)*recip_dyC(i,j,bi,bj) |
146 |
ENDDO |
ENDDO |