anchor.diff:                   .NOT. EqualRealNos(AIMAG(Data4D%DataSet( TmpCoord(1), TmpCoord(2), TmpCoord(3), TmpCoord(4), TmpCoord(5) )), &
                           InitInp%RhoXg * InitInp%WAMITULEN**K * RawData4D(I,9)))  THEN
                  CALL SetErrStat( ErrID_Fatal, ' Line '//TRIM(Num2Lstr(NumHeaderLines+I))//'of '//TRIM(Filename4D)// &
  -                        ' contains different values for the real and imaginary part (columns 7 and 8) than was '// &
  +                        ' contains different values for the real and imaginary part (columns 8 and 9) than was '// &
                           'given earlier in the file for the same values of wave frequency and wave direction.', &
                           ErrStat, ErrMsg, 'Read_DataFile4D' )
                  IF (ALLOCATED(RawData4D))        DEALLOCATE(RawData4D,STAT=ErrStatTmp)
  
positive.code:          ENDIF
  
              ! Find which force component this belongs to
           TmpCoord(5) = INT(RawData4D(I,5))
  
  
              !> The data from the WAMIT file is non-dimensional, so we need to dimensionalize it here.  This
  
negative.code:          ! load component direction must be between 1 and 6 (translational: 1,2,3; rotational: 4,5,6).
        Data3D%LoadComponents = .FALSE.
        DO I=1,K
           IF ( NINT(TmpRealArr(I)) < 1 .OR. NINT(TmpRealArr(K)) > 6 ) THEN
              CALL SetErrStat( ErrID_Fatal, ' Load components listed in column 4 of '//TRIM(Filename3D)// &
                    ' must be between 1 and 6.', ErrStat,ErrMsg,'Read_DataFile3D')
              IF (ALLOCATED(TmpRealArr))       DEALLOCATE(TmpRealArr,STAT=ErrStatTmp)
  