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: 	 	            s% diffusion_dt_limit = original_diffusion_dt_limit
  	 	         endif
  	 	     endif
  		  
  			 !turn off DIFFUSION: diffusion isn't super important post-MS   
               min_center_h1_for_diff = 1d-4
               if (s% do_element_diffusion .and. s% center_h1 < min_center_h1_for_diff) then
  				 if (s% do_element_diffusion) then !only print the first time
  				 	write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  					write(*,*) 'turning off diffusion'
  					write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'  
  				 end if
                   s% do_element_diffusion = .false.
               endif
  			  
  	      end function extras_check_model
  
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
  