anchor.diff:  			 !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
  +				 s% termination_code = t_xtra1
   			   	 extras_check_model = terminate
   			 end if
   			 
  
positive.code: 				 if (((envelope_mass_fraction < 0.15) .and. (L_He/L_tot < 0.05)) .or. &
  					 ((s% center_h1 < 1d-4) .and. (s% center_he4 < 1d-4) .and. (s% Teff > 10**4.5) .and. (s% L_phot > 4))) then
  					 if (s% category_factors(3) > 0) then !only print the first time
                           write(*,*) '++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  						 write(*,*) 'now at post AGB phase, turning off all burning except for H'
  						 write(*,*) '++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  					 end if
  
negative.code: 				 end if
  				 s% varcontrol_target = 1d-3
  			 end if
  			 !set MAX MDOT: cap the mass loss at 1e-3 M/yr
  			 s% star_mdot = max(-1e-3, s% star_mdot)
  			 
  			 !suppress LATE BURNING:turn off burning post-AGB
  			 envelope_mass_fraction = (s% star_mass - s% c_core_mass)/(s% star_mass)
  			 L_He = s% power_he_burn*Lsun
  			 L_tot = s% photosphere_L*Lsun
  