!:qinit:! Field variable initialization
    SUBROUTINE qinit(Info)
      TYPE (NodeInfo) :: Info  ! Data associated with this grid
      ! Internal declarations
      INTEGER i,iQ
      REAL (KIND=xPrec) :: x
      ! iQ=1 inside Hermes, at right; iQ=2 vacuum approximation
      x=Info%Xlower(1)+Info%dX(1)/2.0
      DO i=1,Info%mX(1)
        IF (x >= xBreach) THEN
           iQ = 1
        ELSE
           iQ = 2
        ENDIF
        Info%q(i,:,1,1,1) = rhoQ(iQ)
        Info%q(i,:,1,1,2) = rhoQ(iQ)*uQ(iQ)
        Info%q(i,:,1,1,3) = rhoQ(iQ)*vQ(iQ)
        Info%q(i,:,1,1,4) = pQ(iQ)/gamma1 + 0.5d0*rhoQ(iQ)*(uQ(iQ)**2 + vQ(iQ)**2)
        x=x+Info%dX(1)
      END DO
    END SUBROUTINE qinit