00001 module radial_green_fn_helmholtz
00002   use phys_constant, only : long
00003   implicit none
00004   real(long), pointer  ::  bsjy(:,:,:,:), sbsjy(:,:,:), sbsjyp(:,:,:)
00005 
00006   contains
00007   subroutine allocate_radial_green_fn_helmholtz
00008     use grid_parameter, only : nrg, nlg
00009     use make_array_3d
00010     use make_array_4d
00011     implicit none 
00012     call alloc_array4d(bsjy,1,nrg,0,nlg,0,nlg,0,nrg)
00013     call alloc_array3d(sbsjy,0,nlg,0,nlg,0,nrg)
00014     call alloc_array3d(sbsjyp,0,nlg,0,nlg,0,nrg)
00015   end subroutine allocate_radial_green_fn_helmholtz
00016 end module radial_green_fn_helmholtz