/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/flgrad_midpoint.f90

Go to the documentation of this file.

subroutine flgrad_midpoint(fnc,dfdx,dfdy,dfdz)
  use phys_constant, only : long
  use grid_parameter, only : nrf, ntf, npf
  implicit none
  real(long), pointer :: fnc(:,:,:)
  real(long), pointer :: dfdx(:,:,:), dfdy(:,:,:), dfdz(:,:,:)
  real(long) :: dfdx0, dfdy0, dfdz0
  integer :: irf, itf, ipf
!
! --- Compute the gradient of a function.
! --- The gradient is evaluated at mid points.
!
! --- r, theta, phi derivatives.
!
  do ipf = 1, npf
    do itf = 1, ntf
      do irf = 1, nrf
        call flgrad_midpoint_type0(fnc,dfdx0,dfdy0,dfdz0,irf,itf,ipf)
        dfdx(irf,itf,ipf) = dfdx0
        dfdy(irf,itf,ipf) = dfdy0
        dfdz(irf,itf,ipf) = dfdz0
      end do
    end do
  end do
!
end subroutine flgrad_midpoint

---