verify_parameters Subroutine

  subroutine verify_parameters
    use geometry, only: bishop, local_eq
    use file_utils, only: error_unit
    implicit none
    integer :: ierr

    ierr = error_unit()

    !Disable scan in s_hat when incompatible bishop used
    if (n_shat > 1) then
       select case(bishop)
       case(2,3,4,5,8,9)
          !These are valid values
       case default
          if( (bishop /= 0) .and. local_eq) then
             !This is also ok
          else
             !This is not ok
             n_shat = 1
             write(ierr,'("ERROR : Cannot scan in shat for current combination of bishop and equilbrium type")')
          end if
       end select
    end if

    !Disable scan in beta_prime when incompatible bishop used
    if (n_beta > 1) then
       select case(bishop)
       case(2,3,4,5,6,7,8,9)
          !These are valid values
       case default
          !This is not ok
          n_beta = 1
          write(ierr,'("ERROR : Cannot scan in beta_prime for current combination of bishop and equilbrium type")')
       end select
    end if
  end subroutine verify_parameters