anchor.diff:           IF ( EqualRealNos( ZTime, OutTime ) )  THEN 
                              
                  ! Generate glue-code output file
  -            CALL WrOutputLine( ZTime, p_FAST, y_FAST, EDOutput=y_ED%WriteOutput, SrvDOutput=y_SrvD%WriteOutput, IfWOutput=IfW_WriteOutput, ErrStat=ErrStat, ErrMsg=ErrMsg )
  +            CALL WrOutputLine( ZTime, p_FAST, y_FAST, IfW_WriteOutput, y_ED%WriteOutput, y_SrvD%WriteOutput, InitOutData_HD%WriteOutput, ErrStat, ErrMsg )
               CALL CheckError( ErrStat, ErrMsg )
               
                  ! Generate AeroDyn's element data if desired:
  
positive.code: 
        IF ( p_FAST%CompHydro ) THEN
           CALL HD_InputSolve( p_ED, x_ED, OtherSt_ED, u_ED, y_ED, ErrStat, ErrMsg )
           CALL HD_CalculateLoads( ZTime,  HD_AllMarkers,  HydroDyn_data, HD_AllLoads,  ErrStat )
              CALL CheckError( ErrStat, 'Error calculating hydrodynamic loads in HydroDyn.'  )
        END IF
        
           ! User Tower Loading
        IF ( p_FAST%CompUserTwrLd ) THEN !bjj: array below won't work... routine needs to be converted to UsrTwr_CalcOutput()
  
negative.code: REAL(ReKi)                     :: IfW_WriteOutput(3)                         ! Temporary hack for getting wind speeds from InflowWind
     
     ! Data for the HydroDyn module:
  TYPE(HD_InitOutputType)         :: InitOutData_HD                            ! Initialization output for HydroDyn
     
  
     ! Other/Misc variables
  REAL(DbKi)                     :: TiLstPrn                                   ! The time of the last print
  