anchor.diff:           IfW_WriteOutput = y_AD%IfW_Outputs%WriteOutput
   !<<<         
   
  +      END IF
  +      
  +      
  +      IF (p_FAST%CompHydro) THEN !jmj doesn't think we need this here (because it's taken care of in ED_HD_InputOutputSolve)
  +         CALL HD_InputSolve( p_FAST, HD_Input(1), ED_Output(1), MeshMapData, ErrStat, ErrMsg )
  +            CALL CheckError( ErrStat, 'Message from HD_InputSolve: '//NewLine//ErrMsg  )
  +            
  +         CALL HydroDyn_CalcOutput( this_time, HD_Input(1), p_HD, x_HD_this, xd_HD_this, z_HD_this, OtherSt_HD, y_HD, ErrStat, ErrMsg )
  +            CALL CheckError( ErrStat, 'Message from HD_CalcOutput: '//NewLine//ErrMsg  )
  +            
         END IF
   
         IF ( p_FAST%CompServo ) THEN
  
positive.code:       !............................................................................................................................
        ! Set exit error code if there was an error;
        !............................................................................................................................
        IF (Error) CALL ProgAbort( ' Simulation error level: '//TRIM(GetErrStr(ErrLev)), &   !This assumes PRESENT(ErrID) is .TRUE.
                                    TrapErrors=.FALSE., TimeWait=3._ReKi )  ! wait 3 seconds (in case they double-clicked and got an error)
        
        !............................................................................................................................
        !  Write simulation times and stop
  
negative.code:       INTEGER(IntKi)                       :: ErrStat2                                    ! Error status
        CHARACTER(LEN(ErrMsg))               :: ErrMsg2                                     ! Error message
  
        
        
        !...............................................................................................................................
  