59 |
c & -recip_hFacW(i,j,k,bi,bj)*( Zip-Zij ) ) |
c & -recip_hFacW(i,j,k,bi,bj)*( Zip-Zij ) ) |
60 |
del2u(i,j) = |
del2u(i,j) = |
61 |
& + ( Dij-Dmj )*recip_DXC(i,j,bi,bj) |
& + ( Dij-Dmj )*recip_DXC(i,j,bi,bj) |
62 |
& -recip_hFacW(i,j,k,bi,bj)*( Zip-Zij )*recip_DYU(i,j,bi,bj) |
& -recip_hFacW(i,j,k,bi,bj)*( Zip-Zij )*recip_DYG(i,j,bi,bj) |
63 |
|
|
64 |
c del2v(i,j) = recip_rAs(i,j,bi,bj)*( |
c del2v(i,j) = recip_rAs(i,j,bi,bj)*( |
65 |
c & recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij ) |
c & recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij ) |
68 |
c & recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij ) |
c & recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij ) |
69 |
c & + ( Dij-Dim ) ) |
c & + ( Dij-Dim ) ) |
70 |
del2v(i,j) = |
del2v(i,j) = |
71 |
& recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij )*recip_DXV(i,j,bi,bj) |
& recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij )*recip_DXG(i,j,bi,bj) |
72 |
& + ( Dij-Dim )*recip_DYC(i,j,bi,bj) |
& + ( Dij-Dim )*recip_DYC(i,j,bi,bj) |
73 |
|
|
74 |
ENDDO |
ENDDO |