anchor.diff:  	         end if
   
   	         s% overshoot_f_above_burn_h = frac * s% overshoot_f_above_burn_h
  -	         write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  +			 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
   			 write(*,*) 'core convective overshoot fraction: ', frac
   			 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
   
  
positive.code: 			 !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
  			   	 extras_check_model = terminate
  			 end if
  			 
  
negative.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
  