/home/uryu/Cocal_svn_checkinout/Cocal/code/Analysis/Subroutine/interpolation_matter.f90

Go to the documentation of this file.

subroutine interpolation_matter(fnc,fncca)
  use phys_constant, only : long
  use grid_parameter_cartesian, only : nx, ny, nz
  use coordinate_grav_xyz, only : x, y, z
  use interface_modules_cartesian, ignore_me => interpolation_matter
  implicit none
  real(long), pointer :: fnc(:,:,:)
  real(long), pointer :: fncca(:,:,:)
  real(long) :: xc, yc, zc, cfn
  integer :: ix, iy, iz
!
  do iz = 1, nz
    zc = z(iz)
    do iy = 1, ny
      yc = y(iy)
      do ix = 1, nx
        xc = x(ix)
        call interpo_fl2cgr_4th(fnc,cfn,xc,yc,zc)
        fncca(ix,iy,iz) = cfn
      end do
    end do
  end do
end subroutine interpolation_matter

---