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