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: 				 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
  
negative.code: 				 s% termination_code = t_xtra2
  			   	 extras_check_model = terminate
  			 end if
  			 
  			 !check DIFFUSION: to determine whether or not diffusion should happen 
  	 	     if(s% do_element_diffusion) then !only check if diffusion is on
  	             if(abs(s% mass_conv_core - s% star_mass) < 1d-2) then ! => fully convective
  	 	            s% diffusion_dt_limit = huge_dt_limit
  