anchor.diff:     endif
   
   
  +   nP = p%nWings * (  (p%nSpan+1)*(m%nNW-1+2) +(FWnSpan+1)*(m%nFW+1) )
  +   nFWEff = min(m%nFW, p%nFWFree)
  +   ! --- Display some status to screen
  +   if (mod(n,10)==0) print'(A,F10.3,A,I0,A,I0,A,I0,A,I0,A,I0,A,I0,A,I0,A,I0,A,I0)','FVW status - t:',t,'  Step:',n,'  nNW:',m%nNW-1,'/',p%nNWMax-1,'  nFW:',nFWEff, '+',m%nFW-nFWEff,'=',m%nFW,'/',p%nFWMax,'  nP:',nP
      if (DEV_VERSION) print'(A,F10.3,A,F10.3,A,I0,A,I0,A,I0,A,L1)','Update states, t:',t,'  t_u:', utimes(1),'            Step:',n,' nNW:',m%nNW,' nFW:',m%nFW,' ComputeWake: ',m%ComputeWakeInduced
   
   
  
positive.code:    ErrStat = ErrID_None
     ErrMsg  = "" 
  
     dt = real(p%DTaero,ReKi)
     ! Compute "right hand side"
     CALL FVW_CalcContStateDeriv( t, u, p, x, xd, z, OtherState, m, dxdt, ErrStat2, ErrMsg2); if (Failed()) return
  
  
negative.code:    character(*),                  intent(  out) :: ErrMsg     !< Error message if ErrStat /= ErrID_None
     ! local variables
     type(FVW_ContinuousStateType) :: dxdt ! time derivatives of continuous states
     real(ReKi)     :: dt
     integer(IntKi)       :: ErrStat2      ! temporary error status of the operation
     character(ErrMsgLen) :: ErrMsg2       ! temporary error message
  