anchor.diff:  	         ierr = 0
   			 
   			 original_diffusion_dt_limit = s% diffusion_dt_limit			          
  +			 s% other_wind => VW_superwind
   	      end subroutine extras_controls
               
   	      integer function extras_startup(s, id, restart, ierr)
  
positive.code: 			 rot_full_off = s% job% extras_rpar(3) !1.2
  			 rot_full_on = s% job% extras_rpar(4) !1.8
  			 
  	         if (s% star_mass < rot_full_off) then
  				 frac2 = 0
  				 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  				 write(*,*) 'no rotation'
  				 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  	         else if (s% star_mass >= rot_full_off .and. s% star_mass <= rot_full_on) then
  				 frac2 = (s% star_mass - rot_full_off) / &
  	                    (rot_full_on - rot_full_off)
  	             frac2 = 0.5d0*(1 - cos(pi*frac2))
  				 s% job% set_near_zams_omega_div_omega_crit_steps = 10
  				 s% job% new_omega_div_omega_crit = s% job% extras_rpar(5) * frac2
  				 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  				 write(*,*) 'new omega_div_omega_crit, fraction', s% job% new_omega_div_omega_crit, frac2
  				 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  			 else
  				 frac2 = 1.0
  				 s% job% set_near_zams_omega_div_omega_crit_steps = 10
  				 s% job% new_omega_div_omega_crit = s% job% extras_rpar(5) * frac2 !nominally 0.4
  				 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  				 write(*,*) 'new omega_div_omega_crit, fraction', s% job% new_omega_div_omega_crit, frac2
  				 write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  	         end if
  
  			 !set VARCONTROL: for massive stars, turn up varcontrol gradually to help them evolve
  			 vct30 = 1e-4
  
negative.code:       use star_lib
        use star_def
        use const_def
  	  use crlibm_lib
        ! use run_star_support
        
        implicit none
  