anchor.diff:       end do
   
   
  -     ! set namelist defaults
  +     ! Set namelist defaults
        denerr = 1.d20
        dengrad = 1.d20
        max_denerr_lev = 10
  
positive.code: 
       dens_ambient = model_S_state(npts_model_S,idens_model)
  
       ! check that the cutoff temperature match
       if (model_S_state(npts_model_S,itemp_model) /= &
           model_P_state(npts_model_P,itemp_model)) then
          call bl_error("ERROR: Cutoff temperatures for primary and secondary models do not agree.")
  
negative.code:      temp_ambient = model_S_state(npts_model_S,itemp_model)
  
  
       ! Ambient X
       xn_ambient(:) = model_S_state(npts_model_S,ispec_model:ispec_model-1+nspec)
  
  
       ! 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
       radius_P_initial = -1.0d0
       do i = 1, npts_model_P
          if (model_P_state(i,idens_model) <= 1.1*dens_ambient) then
  