anchor.diff:      call eos(eos_input_rt, eos_state, .false.)
   
       small_pres = eos_state % p
  -!    small_ener = eos_state % e
  +    small_ener = eos_state % e
   
     end subroutine set_small
   
  
positive.code:        call init_1d(model_S_r, model_S_state, npts_model, dx, radius_S_initial, mass_S, &
                      central_density_S, stellar_temp, stellar_comp, ambient_state)
  
         if (ioproc == 1) then
            print *, "Generated initial model for secondary WD of mass", mass_S, &
                     ", central density", central_density_S, ", and radius", radius_S_initial
         endif
  
negative.code:     call init_1d(model_P_r, model_P_state, npts_model, dx, radius_P_initial, mass_P, &
                   central_density_P, stellar_temp, stellar_comp, ambient_state)
  
      if (ioproc == 1 .and. init == 1) then
          print *, "Generated initial model for primary WD of mass", mass_P, &
                   ", central density", central_density_P, ", and radius", radius_P_initial
      endif
  