anchor.diff:    ! Expose only the needed variables and functions
     public create_body, set_state 
     public get_rotation
  -  public get_angrate, get_angrate_dot
  +  public get_angrate, get_angrate_dot, get_angrate_inv
     public print_body
   
     !*******************************************************************!
  
positive.code: 
    end function get_angrate_dot_cosines
  
    ! ********************************************************
    ! routine that prints the state and properties of the body
    ! ********************************************************
    subroutine print_body(alpha)
  
      use dispmodule
  
negative.code:     ! Assuming the body axis to be at the centre of mass of the body
      ! Inertial or body frame?
      alpha%KE = 0.5_dp*(mass * alpha%v *  alpha%v &
           & + alpha%omega*alpha%J* alpha%omega) ! + coupling term
  
      ! potential energy due to position
      ! include strain energy later
  