anchor.diff:  			 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
  +			 if (s% job% extras_rpar(5) > 0.0) then
  +	         	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
  +			else
  +				write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  +				write(*,*) 'no rotation'
  +				write(*,*) '+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'
  +			end if
  +				
   
   			 !set VARCONTROL: for massive stars, turn up varcontrol gradually to help them evolve
   			 vct30 = 1e-4
  
positive.code: 			 endif
  	      end function extras_finish_step
        
  	      subroutine extras_after_evolve(s, id, id_extra, ierr)
  	         type (star_info), pointer :: s
  	         integer, intent(in) :: id, id_extra
  
negative.code: 	         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)
  