anchor.diff:        
         temp   = (1.0-1.3*AxInd_disk)*Un_disk
         
  -      tau1   = 1.1*u%R_disk/temp                                                                      ! Eqn. 1.2 (note that we've eliminated possibility of temp being 0)
  +      tau1   = 1.1*u%R_disk/temp          ! Eqn. 1.2 (note that we've eliminated possibility of temp being 0)
  +      tau1   = min(tau1, 100.0_ReKi)      ! put a limit on this time constant so it isn't unrealistically long (particularly at initialization)
         
      end if
   
  
positive.code:       if ( u%AxInd_disk > max_AxInd ) then
           AxInd_disk = max_AxInd
           if (m%FirstWarn_tau1) then
              call setErrStat( ErrID_Severe, 'Rotor-averaged axial induction factor is greater than 0.7; limiting time-varying tau1.' &
                                //' This message will not be repeated though the condition may persist.', ErrStat, ErrMsg, RoutineName ) ! don't print this error more than one time
              m%FirstWarn_tau1 = .false.
           end if
        else
  
negative.code:    real(ReKi)                                   :: temp
  
     real(ReKi)                                   :: AxInd_disk
     real(ReKi), parameter                        :: max_AxInd = 0.5_ReKi
  
     real(ReKi)                                   :: Un_disk
     real(ReKi), parameter                        :: min_Un = 0.1_ReKi
     
     character(*), parameter                      :: RoutineName = 'ComputeTau1'
  
  