set_vgrid Subroutine

  subroutine set_vgrid
    use species, only: nspec
    use egrid, only: setvgrid, setvgrid_genquad
    allocate (energy(negrid,nspec), w(negrid,nspec))
    if (genquad) then
       call setvgrid_genquad (negrid, energy, w)
    else
       call setvgrid (vcut, negrid, energy, w, nesub)
    end if
  end subroutine set_vgrid