!:GlobalVariables:! Problem specific global variables
  ! Parameters
  INTEGER, PARAMETER :: OK=0
  REAL (KIND=qPrec), PARAMETER :: zero=0., half=0.5d0, one=1.d0, two=2.d0
  ! Variables
  LOGICAL efix                                     ! Flag: efix=T to apply entropy fix
  REAL (KIND=qPrec) :: gamma,gamma1                ! Adiabatic coefficient gamma, gamma-1
  REAL (KIND=xPrec) :: xBreach                     ! Breach location
  REAL (KIND=qPrec) :: vHermes,mHermes             ! Hermes velocity, mass
  REAL (KIND=qPrec), DIMENSION(2) :: pQ,rhoQ,uQ,vQ ! Initial values

  !=======
  CONTAINS
  !=======