subroutine orbsmoo
  use iso_c_binding, only: fp => c_double
  !
  ! smooth radial profile outputs of monte carlo calculation
  ! --the profiles just normalized in outptb
  !
  ! (dmc mar 2004 -- smoothing moved to generated code (nbo_load) ...
  !  only "hole patch" smooth is left here).
  !
  ! subroutines used:  comput/bsmooh.for, subroutine r8bsmoop
  !
  !   orbsmoo smooths the profiles for the **current specie** that were
  !   renormalized by a prior call to outptb
  !
  use nbi_com
  !
  implicit none
  real(fp) :: z1(mj)
  !
  !-----------
  !
  z1=1.0_fp
  !
  if(.not.nlfbmfpp(lsbeam)) then
     !
     call r8bsmoop(omegabs(1:mj,lsbeam),z1)
     call r8bsmoop(ebavps(1:mj,lsbeam),z1)
     call r8bsmoop(ebapp2s(1:mj,lsbeam),z1)
     call r8bsmoop(ebapl2s(1:mj,lsbeam),z1)
     call r8bsmoop(xzbmav(1:mj,lsbeam),z1)
     call r8bsmoop(wttmf(1:mj,lsbeam),z1)
     !
  end if
  !
  return
end subroutine orbsmoo
