00001 subroutine interpolation_metric(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   implicit none
00007   real(long), pointer :: fnc(:,:,:)
00008   real(long), pointer :: fncca(:,:,:)
00009   real(long) :: xc, yc, zc, cfn
00010   integer :: ix, iy, iz
00011 
00012   do iz = 1, nz
00013     zc = z(iz)
00014     do iy = 1, ny
00015       yc = y(iy)
00016       do ix = 1, nx
00017         xc = x(ix)
00018         call interpo_gr2cgr_4th(fnc,cfn,xc,yc,zc)
00019         fncca(ix,iy,iz) = cfn
00020       end do
00021     end do
00022   end do
00023 end subroutine interpolation_metric