write_linear_sweet_spot Subroutine

  subroutine write_linear_sweet_spot(dim, report_unit, nfacs_dim, facs_dim, &
       previous_dim_block, nmesh, ncut)
    implicit none
    character, intent(in) :: dim
    integer, intent(in) :: report_unit, nfacs_dim, previous_dim_block, nmesh, ncut
    integer, dimension(:), intent(in) :: facs_dim
    integer :: i, npe
    ! We start at 2 to avoid the duplicate factor coming from the previous dimension
    do i = 2, nfacs_dim
       npe = facs_dim(i) * previous_dim_block
       if (nmesh / npe > ncut) write (report_unit, fmt="('  npe = ',i8,'  (',a,')')") npe, dim
    end do
  end subroutine write_linear_sweet_spot