anchor.diff:         ! velocity
          sum_xm = 0.0_dp_t
          sum_xp = 0.0_dp_t
  -       do n = 1, nprim
  +       do n = 1, nwaves
             sum_xm = sum_xm + beta_xm(n)*rvec(n,2)
             sum_xp = sum_xp + beta_xp(n)*rvec(n,2)
          enddo
  
positive.code: module interface_states_plm_module
  
    use datatypes_module
    use  grid_module
    use   variables_module
    use params_module
    use eos_module
  
    implicit none
  
  
negative.code:     type(gridedgevar_t) :: Q_l, Q_r
      type(gridvar_t) :: ldelta, temp
  
      real (kind=dp_t) :: rvec(nwaves,nprim), lvec(nwaves,nprim), eval(nwaves)
      real (kind=dp_t) :: dQ(nprim)
  
      real (kind=dp_t) :: r, ux, p, cs
      real (kind=dp_t) :: ldr, ldu, ldp
      real (kind=dp_t) :: r_xm, r_xp, u_xm, u_xp, p_xm, p_xp
      real (kind=dp_t) :: beta_xm(nwaves), beta_xp(nwaves)
      real (kind=dp_t) :: sum_xm, sum_xp
      real (kind=dp_t) :: sum
  
  