!:setprob:! Read problem global data
    SUBROUTINE setprob
      INTEGER i
      ! Read gas parameters from file
      OPEN(UNIT=7,FILE='setprob.data',STATUS='old',FORM='formatted')
      READ(7,*)efix                           ! set to T to apply entropy fix in Riemann solver
      READ(7,*) gamma                         ! adiabatic constant
      gamma1 = gamma - one
      READ(7,*) xBreach                       ! position of Hermes breach
      READ(7,*) mHermes
      DO i=1,2
        READ(7,*)pQ(i),rhoQ(i),uQ(i),vQ(i)    ! initial values (primitive variables)
      END DO
      CLOSE(7)
      ! Initialize Hermes velocity, reaction force of escaping gas
      vHermes=zero
    END SUBROUTINE setprob