anchor.diff:                 J = J + 1
               END DO
            END IF
  -         
  -      !END IF
  -      
  +
         IF ( p%OutSwtch == 1 .OR. p%OutSwtch == 3 ) THEN
            CALL HDOut_OpenOutput( HydroDyn_ProgDesc, InitInp%OutRootName, p, InitOut, ErrStat, ErrMsg )
            IF (ErrStat >= AbortErrLev ) RETURN
  
positive.code:    
     CALL HDOUT_ChkOutLst( InitInp%OutList(1:p%NumOuts), y, p, ErrStat, ErrMsg )
     IF ( ErrStat /= 0 ) RETURN
     
     
    !IF ( ALLOCATED( p%OutParam ) .AND. p%NumOuts > 0 ) THEN           ! Output has been requested so let's open an output file            
    !    
    !    ALLOCATE( y%WriteOutput( p%NumOuts ),  STAT = ErrStat )
    !    IF ( ErrStat /= ErrID_None ) THEN
    !       ErrMsg  = ' Error allocating space for WriteOutput array.'
    !       ErrStat = ErrID_Fatal
    !       RETURN
    !    END IF
    !    y%WriteOutput = 0.0_ReKi
    !    
    !      ALLOCATE ( InitOut%WriteOutputHdr(p%NumOuts), STAT = ErrStat )
    !    IF ( ErrStat /= ErrID_None ) THEN
    !       ErrMsg  = ' Error allocating space for WriteOutputHdr array.'
    !       ErrStat = ErrID_Fatal
    !       RETURN
    !    END IF
    !    
    !    ALLOCATE ( InitOut%WriteOutputUnt(p%NumOuts), STAT = ErrStat )
    !    IF ( ErrStat /= ErrID_None ) THEN
    !       ErrMsg  = ' Error allocating space for WriteOutputHdr array.'
    !       ErrStat = ErrID_Fatal
    !       RETURN
    !    END IF   
    !
    !    DO I = 1,p%NumOuts
    !       
    !       InitOut%WriteOutputHdr(I) = TRIM( p%OutParam(I)%Name  )
    !       InitOut%WriteOutputUnt(I) = TRIM( p%OutParam(I)%Units )      
    !    
    !    END DO   
    !    
    ! END IF   ! there are any requested outputs   
  
        
       
           ! Aggregate the sub-module initialization outputs for the glue code
        !IF ( p%OutSwtch == 2 .OR. p%OutSwtch == 3 ) THEN
           
           !hasWAMITOuts   = .FALSE.
           hasWAMIT2Outs  = .FALSE.
           hasWaves2Outs  = .FALSE.
           hasMorisonOuts = .FALSE.
  
negative.code:             
           J = p%NumOuts + 1
           
           IF ( hasWaves2Outs ) THEN
              DO I=1, p%Waves2%NumOuts
                 InitOut%WriteOutputHdr(J) = InitOut%Waves2%WriteOutputHdr(I)
  