anchor.diff:       temp_ambient = model_S_state(npts_model_S,itemp_model)
   
   
  -     ! ambient X
  +     ! Ambient X
        xn_ambient(:) = model_S_state(npts_model_S,ispec_model:ispec_model-1+nspec)
   
   
  -     ! get the rest of the ambient thermodynamic state
  +     ! Given the inputs of small_dens and small_temp, figure out small_pres.
  +     ! Use the ambient gas composition (we don't need to worry about saving
  +     ! our result in eint_ambient since it will be reset in the next call).
  +
  +     call eos_given_RTX(eint_ambient,small_pres,small_dens, &
  +                        small_temp,xn_ambient)
  +
  +     ! Get the rest of the ambient thermodynamic state
        call eos_given_RTX(eint_ambient,pres_ambient,dens_ambient, &
                           temp_ambient,xn_ambient)
   
  -
  -     ! compute the radius of the primary
  +     ! Compute the radius of the primary
        radius_P_initial = -1.0d0
        do i = 1, npts_model_P
           if (model_P_state(i,idens_model) <= 1.1*dens_ambient) then
  
positive.code:      use model_parser_module
       use bl_constants_module
       use fundamental_constants_module
       use eos_module
       use com, only: mass_p, mass_s, com_loc_p, com_loc_s
  
  
negative.code:      end do
  
  
       ! Set namelist defaults
       denerr = 1.d20
       dengrad = 1.d20
       max_denerr_lev = 10
  