anchor.diff:           
      ELSEIF ( p_FAST%CompAero == Module_AD ) THEN
         
  +!debug_print2_unit = 92   
         CALL AD_InputSolve_NoIfW( p_FAST, AD%Input(1), ED%Output(1), MeshMapData, ErrStat2, ErrMsg2 )   
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )        
  +!debug_print2_unit = 90      
   
            ! because we're not calling InflowWind_CalcOutput, this probably can be skipped; 
            ! TODO: alternatively, we could call InflowWind_CalcOutput, too.
  
positive.code:       
     ELSE IF ( p_FAST%CompAero == Module_AD ) THEN 
                          
        CALL AD_InputSolve_NoIfW( p_FAST, AD%Input(1), ED%Output(1), MeshMapData, ErrStat2, ErrMsg2 )
           CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName ) 
           
     END IF
     
  
negative.code:       
        InitInData_BD%RootName     = p_FAST%OutFileRoot     
        InitInData_BD%gravity      = (/ 0.0_ReKi, 0.0_ReKi, -InitOutData_ED%Gravity /)       ! "Gravitational acceleration" m/s^2
        InitInData_BD%GlbPosHub    = ED%Output(1)%HubPtMotion%Position(:,1)                  ! "Initial Position Vector of the hub (center of rotation)" 
        
           ! now initialize IceD for additional legs (if necessary)
        dt_BD = p_FAST%dt_module( MODULE_BD )
                          
  