subroutine setprob
      implicit double precision (a-h,o-z)
      common /param/ gamma, gamma1
c
c     # Set gamma and gamma1 = gamma-1 for Euler equations
c     # Passed to the Riemann solver rp1.f in a common block
c
      open(unit=7,file='setprob.data',status='old',form='formatted')


      read(7,*) gamma
      gamma1 = gamma - 1.d0

      return
      end