anchor.diff:    use rpar_indices
     
     implicit none
  +  !$acc routine seq
   
     real(kind=dp_t), intent(in   ) :: t9, dens, y(nspec)
     real(kind=dp_t), intent(inout) :: rpar(n_rpar_comps)
  
positive.code:   return
  
  end subroutine rhs
  
  
  
  
  subroutine make_rates(t9,dens,y,rpar)
  
negative.code:   use rpar_indices
    
    implicit none
    !$acc routine seq
  
    real(kind=dp_t), intent(IN   ) :: ymol(nspec), t9
    real(kind=dp_t), intent(INOUT) :: rpar(n_rpar_comps)
  