anchor.diff:  
     ! gravitational potential
     phi_comp = plotfile_var_index(pf, "phiGrav")
  +  
  +  !Mach Number
  +  mach_comp = plotfile_var_index(pf, "MachNumber")
   
   
   
  
positive.code: 
    real(kind=dp_t), pointer :: p(:,:,:,:)
  
    real(kind=dp_t) :: kinetic_energy, potential_energy, internal_energy
  
    integer :: dens_comp, rhoe_comp, u_comp, v_comp, w_comp, phi_comp
  
    real(kind=dp_t) :: f_rotate
  
  
negative.code: 100 format(1x, 5(g20.10))
  
    write (*, 100) pf%tm, kinetic_energy, internal_energy, -potential_energy, &
         kinetic_energy + internal_energy - potential_energy, max_mach
  
  
    call destroy(pf)
  