grid modelがしんどいわけ　その３

staggered gridを使うと物理領域内のgridの数が変数によって異なる

do i = 0, nx
  dp_dx(i) = ( p(i+1) - p(i) ) / dx
end do

do i = 1, nx
  du_dx(i) = ( u(i) - u(i-1) ) / dx
end do