00001 subroutine outer_boundary_d_Bfun(sou_surf) 00002 use phys_constant, only : long 00003 use grid_parameter, only : nrg, ntg, npg, rgin 00004 implicit none 00005 real(long), pointer :: sou_surf(:,:) 00006 integer :: itg, ipg 00007 ! 00008 do ipg = 1, npg 00009 do itg = 1, ntg 00010 ! sou_surf(itg,ipg) = 1.0d0 00011 sou_surf(itg,ipg) = 0.0d0 00012 end do 00013 end do 00014 ! 00015 end subroutine outer_boundary_d_Bfun