anchor.diff:  
         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 HD_CalculateLoads( ZTime,  HD_AllMarkers,  HydroDyn_data, HD_AllLoads, InitOutData_HD, 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()
  
positive.code:       HydroDyn_InitData%OutRootName = p_FAST%OutFileRoot
        HD_ConfigMarkers%Substructure%Position = (/0._ReKi, 0._ReKi, 0._ReKi/)   
  
        CALL HD_Init(HydroDyn_InitData, HD_ConfigMarkers, HD_AllMarkers, HydroDyn_data, ErrStat)
        CALL CheckError( ErrStat, 'Error initializing HydroDyn.' )
  
  
  
negative.code:    CALL SrvD_DestroyInitInput(  InitInData_SrvD, ErrStat, ErrMsg )
     CALL SrvD_DestroyInitOutput( InitOutData_SrvD, ErrStat, ErrMsg )
  
  
     !CALL HD_DestroyInitOutput( InitOutData_HD, ErrStat, ErrMsg )
     !
     !...............................................................................................................................
     ! loose coupling
     !...............................................................................................................................
  