anchor.diff:  	 	         endif
   	 	     endif
   		  
  -			 !turn off diffusion post MS   
  +			 !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
  
positive.code: 	            call unpack_extra_info(s)
  	         end if
  			 
  			 !set Zbase for Type 2 Opacities automatically to the Z for the star
  			 s% Zbase = 1.0 - (s% job% initial_h1 + s% job% initial_h2 + &
  			 s% job% initial_he3 + s% job% initial_he4)
  			 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  
negative.code: 	         extras_finish_step = keep_going
  	         call store_extra_info(s)
  
  			 !set NET: evolve 100 steps then change to approx21_extras net
               if (s% model_number == 100) then
  				
  				write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  