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

Go to the documentation of this file.

subroutine outer_boundary_d_poty(sou_surf,Bfun,dBfundy)
  use phys_constant, only  : long
  use grid_parameter, only : nrg, ntg, npg, rgin
  implicit none
  real(long), pointer :: sou_surf(:,:), Bfun(:,:,:), dBfundy(:,:,:)
  real(long) :: dBdx,dBdy,dBdz, bf
  integer    :: itg, ipg
! Reset boundary condition for BH
  do ipg = 1, npg
    do itg = 1, ntg
      dBdy = 0.25d0*(dBfundy(nrg,itg,ipg)   + dBfundy(nrg,itg-1,ipg)  &
                   + dBfundy(nrg,itg,ipg-1) + dBfundy(nrg,itg-1,ipg-1) )

      bf = 0.25d0*(Bfun(0,itg,ipg)   + Bfun(0,itg-1,ipg)  &
                 + Bfun(0,itg,ipg-1) + Bfun(0,itg-1,ipg-1) )
!
!      sou_surf(itg,ipg) = 0.0d0 + 0.25d0*dBdy/bf
!      sou_surf(itg,ipg) = 0.0d0 + 0.25d0*dBdy
      sou_surf(itg,ipg) = 0.0d0 
    end do
  end do
!
end subroutine outer_boundary_d_poty

---