anchor.diff:     nMax = nMax + (FWnSpan+1) * (p%nFWMax+1) * p%nWings   ! Far wake points
      call AllocAry( m%r_wind, 3, nMax, 'Requested wind points', ErrStat2, ErrMsg2 );call SetErrStat ( ErrStat2, ErrMsg2, ErrStat,ErrMsg,'FVW_InitMisc' )
      m%r_wind = 0.0_ReKi     ! set to zero so InflowWind can shortcut calculations
  +   m%OldTime = 0.0_DbKi-p%DTfvw    ! First time for wake interaction calculations
   end subroutine FVW_InitMiscVars
   ! ==============================================================================
   subroutine FVW_InitStates( x, p, m, ErrStat, ErrMsg )
  
positive.code:    call AllocAry( dxdt%r_FW , 3   ,  FWnSpan+1  ,p%nFWMax+1,  p%nWings, 'Wind on FW ', ErrStat, ErrMsg ); dxdt%r_FW= -999999_ReKi;
  
     ! Only calculate freewake after start time and if on a timestep when it should be calculated.
     if (t> p%FreeWakeStart .and. m%ComputeWakeInduced) then
        nFWEff = min(m%nFW, p%nFWFree)
        
        ! --- Compute Induced velocities on the Near wake and far wake based on the marker postions:
        ! (expensive N^2 call)
        ! In  : x%r_NW,    r%r_FW 
  
negative.code:    !call print_x_NW_FW(p, m, z, x,'Conv')
  
  
     if (m%ComputeWakeInduced) then
        ! --- t+dtaero
        ! Propagation/creation of new layer of panels
        call PropagateWake(p, m, z, x, ErrStat2, ErrMsg2)
        !call print_x_NW_FW(p, m, z, x,'Prop_')
     endif
  
     ! Inputs at t+dtaero
     call FVW_Input_ExtrapInterp(u(1:size(utimes)),utimes,uInterp,t+dtaero, ErrStat2, ErrMsg2); if(Failed()) return
  