anchor.diff:  #ifdef SOLVE_OPTION_1_BEFORE_2
   ! used for Option 1 before Option 2:
   
  -   IF ( p_FAST%CompSub /= Module_None .OR. p_FAST%CompHydro == Module_HD ) THEN
  +   IF ( p_FAST%CompSub == Module_SD .OR. p_FAST%CompHydro == Module_HD ) THEN
      ! Because SubDyn needs a better initial guess from ElastoDyn, we'll add an additional call to ED_CalcOutput to get them:
      ! (we'll do the same for HydroDyn, though I'm not sure it's as critical)
      
  
positive.code:    !! Check to see if we should output data this time step:
     !----------------------------------------------------------------------------------------
  
     CALL WriteOutputToFile(n_t_global, m_FAST%t_global, p_FAST, y_FAST, ED, BD, AD14, AD, IfW, OpFM, HD, SD, ExtPtfm, &
                            SrvD, MAPp, FEAM, MD, Orca, IceF, IceD, MeshMapData, ErrStat2, ErrMsg2)
        CALL SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
     !----------------------------------------------------------------------------------------
     !! Display simulation status every SttsTime-seconds (i.e., n_SttsTime steps):
     !----------------------------------------------------------------------------------------   
        
     IF (p_FAST%WrSttsTime) then
        IF ( MOD( n_t_global + 1, p_FAST%n_SttsTime ) == 0 ) THEN
        
           if (.not. Cmpl4SFun) then   
              CALL SimStatus( m_FAST%TiLstPrn, m_FAST%PrevClockTime, m_FAST%t_global, p_FAST%TMax, p_FAST%TDesc )
           end if
        
        ENDIF
     ENDIF
       
  END SUBROUTINE FAST_Solution
  !----------------------------------------------------------------------------------------------------------------------------------
  
negative.code:    CHARACTER(*), PARAMETER                 :: RoutineName = 'WrVTK_AllMeshes'
     
     
     if (allocated(ED%Output)) then
        NumBl = SIZE(ED%Output(1)%BladeRootMotion)            
     else
        NumBl = 0
     end if
     
     
  