anchor.diff:        real(dp), intent(out) :: w ! wind in units of Msun/year (value is >= 0)
         integer, intent(out) :: ierr
         real(dp) :: logP, P, reimers_w, pre_superwind_w, superwind_w, vexp, agb_w, center_h1, center_he4
  -      ierr = 0
  -
  +	  integer :: h1, he4
  +	  type (star_info), pointer :: s
  +      ierr = 0
  +      call star_ptr(id, s, ierr)
  +      if (ierr /= 0) return
  +	  
   !     first Reimers for the MS+RGB
         reimers_w = 4d-13*(Lsurf*Rsurf/Msurf)/(Lsun*Rsun/Msun)
         reimers_w = reimers_w * s% Reimers_wind_eta
  
positive.code:       
        subroutine Reimers_then_VW(id, Lsurf, Msurf, Rsurf, Tsurf, w, ierr)
        use star_def
        integer, intent(in) :: id
        real(dp), intent(in) :: Lsurf, Msurf, Rsurf, Tsurf ! surface values (cgs)
  !     NOTE: surface is outermost cell. not necessarily at photosphere.
  
negative.code:       end if
        end if
        
  !      define STOPPING CRITERION: stopping criterion for C burning, massive stars.
        if ((s% center_h1 < 1d-4) .and. (s% center_he4 < 1d-4) .and. (s% center_c12 < 1d-4)) then
           termination_code_str(t_xtra1) = 'central C12 mass fraction below 1e-4'
           s% termination_code = t_xtra1
  