anchor.diff:     !call print_x_NW_FW(p, m, z, x,'Conv')
   
   
  -   ! --- 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_')
  +   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
  
positive.code: end subroutine FVW_InitConstraint
  ! ==============================================================================
  subroutine FVW_Init_Y( p, u, y, ErrStat, ErrMsg )
  !TODO: move the r_wind to Miscvars
     type(FVW_ParameterType),         intent(in   )  :: p              !< Parameters
     type(FVW_InputType),             intent(inout)  :: u              !< An initial guess for the input; input mesh must be defined
     type(FVW_OutputType),            intent(  out)  :: y              !< Constraints
  
negative.code:    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 )
  