anchor.diff:    kinetic_energy = ZERO
     internal_energy = ZERO
     potential_energy = ZERO
  +  mach_number = ZERO
  +  max_mach = ZERO
  +  !Radius of the star. Model radius for secondary is 7.60156250e8, Model radius for primary is 9.69531250e8
  +  rad_star = 7.60156250e8
   
     imask(:,:,:) = .true.
   
  
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: 
    ! gravitational potential
    phi_comp = plotfile_var_index(pf, "phiGrav")
    
    !Mach Number
    mach_comp = plotfile_var_index(pf, "MachNumber")
  
  
  
  