anchor.diff:  
     ! binary properties
     double precision, save :: period
  -  double precision, save :: a_A, a_B
  +  double precision, save :: a_P_initial, a_S_initial
     
  -  double precision, save :: x_cen_A, y_cen_A, z_cen_A
  -  double precision, save :: x_cen_B, y_cen_B, z_cen_B
  +  double precision, dimension(3), save :: center_P_initial, center_S_initial
     
   end module probdata_module
  
positive.code:   integer         , save ::  max_raderr_lev,   max_radgrad_lev
    
    ! model info
    character (len=80), save :: model_A_name
    character (len=80), save :: model_B_name
    
    double precision, allocatable, save :: model_A_r(:), model_A_state(:,:)
    double precision, allocatable, save :: model_B_r(:), model_B_state(:,:)
    
    integer, save :: npts_model_A, npts_model_B
    
    double precision, save :: mass_A, mass_B
    double precision, save :: radius_A, radius_B
    
    double precision, save :: dens_ambient, temp_ambient
    double precision, save :: xn_ambient(nspec)
  
negative.code:   logical, save :: inertial
  
    ! damping
    logical, save    :: damping
    double precision :: damping_alpha
  
    ! grid info
  