anchor.diff:      type(gridedgevar_t), intent(inout) :: U_l, U_r
       real (kind=dp_t),    intent(in   ) :: dt
   
  -    integer :: i
  +    type(gridvar_t) :: Q
  +    type(gridedgevar_t) :: Q_l, Q_r
  +    type(gridvar_t) :: ldelta, temp
  +
  +    real (kind=dp_t) :: e
  +    real (kind=dp_t) :: test
  +
  +    integer :: i, n
   
       ! piecewise linear slopes -- this is a 1-d version of the Colella
       ! 2nd order unsplit Godunov scheme (Colella 1990).  See also
  
positive.code:     ! compute the monotonized central differences
      !-------------------------------------------------------------------------
  
  
  
  
      !-------------------------------------------------------------------------
      ! compute left and right primitive variable states
      !-------------------------------------------------------------------------
  
  
  
  
      !-------------------------------------------------------------------------
      ! transform the states into conserved variables
      !-------------------------------------------------------------------------
  
  
  
  
      !-------------------------------------------------------------------------
  
negative.code: 
  
  
  
    end subroutine make_interface_states_plm
  
  