anchor.diff:      vector_gas, &
       flux_olg, &
       flux_gas, &
  -    gaussrootsweights
  -
  +    gaussrootsweights, &
  +    reconstruct_velocity_quadrature
       !use refinter, only: &
       !andean_mountain_data, &
       !smooth_andean_mountain, &
  
positive.code:     call zero_vector(tracereq)
  
      !===============================================================
      !Calculate temperature tendency
      !===============================================================
  
  
negative.code:     do i=1,nodes
          node(i)%phi_exa=node(i)%phi_new2
      enddo
      !print*, minval(node(1:nodes)%phi_new2), maxval(node(1:nodes)%phi_new2)
  
  end subroutine highorder_adv_vars
  
  