00001 subroutine interpolation_matter(fnc,fncca)
00002 use phys_constant, only : long
00003 use grid_parameter_cartesian, only : nx, ny, nz
00004 use coordinate_grav_xyz, only : x, y, z
00005 use interface_modules_cartesian
00006
00007 implicit none
00008 real(long), pointer :: fnc(:,:,:)
00009 real(long), pointer :: fncca(:,:,:)
00010 real(long) :: xc, yc, zc, cfn
00011 integer :: ix, iy, iz
00012
00013 do iz = 1, nz
00014 zc = z(iz)
00015 do iy = 1, ny
00016 yc = y(iy)
00017 do ix = 1, nx
00018 xc = x(ix)
00019 call interpo_fl2cgr_4th(fnc,cfn,xc,yc,zc)
00020 fncca(ix,iy,iz) = cfn
00021 end do
00022 end do
00023 end do
00024 end subroutine interpolation_matter