anchor.diff:      use meth_params_module, only: rot_period
       use initial_model_module, only: init_1d
       use prob_params_module, only: center
  +    use meth_params_module, only: rot_axis
   
       implicit none
   
       double precision :: dx
       type (eos_t) :: ambient_state
  +
  +    double precision :: loc_P(3), loc_S(3)
   
       call get_ambient(ambient_state)
   
  
positive.code: 
  
  
    ! This routine checks to see if the primary mass is actually larger
    ! than the secondary mass, and switches them if not.
  
  
negative.code:   double precision, dimension(3) :: center_P_initial, center_S_initial
  
    integer :: star_axis
    integer :: initial_motion_dir
  
    ! Density of ambient medium
    double precision :: ambient_density
  