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

Go to the documentation of this file.

subroutine grgrad_midpoint_r3rd_nsbh(fnc,dfdx,dfdy,dfdz,cobj)
  use phys_constant, only : long
  use grid_parameter, only : nrg, ntg, npg
  use interface_grgrad_midpoint_r3rd_type0_ns
  use interface_grgrad_midpoint_r3rd_type0_bh
  implicit none
  real(long), pointer :: fnc(:,:,:)
  real(long), pointer :: dfdx(:,:,:), dfdy(:,:,:), dfdz(:,:,:)
  real(long) :: dfncdx, dfncdy, dfncdz
  integer    :: irg, itg, ipg
  character(len=2), intent(in) :: cobj
!
! --- Compute the gradient of a function.
! --- The gradient is evaluated at mid points.
! --- r, theta, phi derivatives.
!
  if(cobj=='ns') then
    do irg = 1, nrg
      do itg = 1, ntg
        do ipg = 1, npg
          call grgrad_midpoint_r3rd_type0_ns(fnc,dfncdx,dfncdy,dfncdz,irg,itg,ipg)
          dfdx(irg,itg,ipg) = dfncdx
          dfdy(irg,itg,ipg) = dfncdy
          dfdz(irg,itg,ipg) = dfncdz
        end do 
      end do 
    end do
  else
    do irg = 1, nrg
      do itg = 1, ntg
        do ipg = 1, npg
          call grgrad_midpoint_r3rd_type0_bh(fnc,dfncdx,dfncdy,dfncdz,irg,itg,ipg)
          dfdx(irg,itg,ipg) = dfncdx
          dfdy(irg,itg,ipg) = dfncdy
          dfdz(irg,itg,ipg) = dfncdz
        end do
      end do
    end do
  end if
!
end subroutine grgrad_midpoint_r3rd_nsbh

---