anchor.diff:        INTEGER(IntKi)                                     :: J                    !< Generic counter
         INTEGER(IntKi)                                     :: K                    !< Generic counter
         TYPE(FFT_DataType)                                 :: FFT_Data             !< Temporary array for the FFT module we're using
  -
  +      CHARACTER(*), PARAMETER                            :: RoutineName = 'NewmanApp_InitCalc'
  + 
   
            ! Wave information and QTF temporary
         COMPLEX(SiKi)                                      :: QTF_Value            !< Temporary complex number for QTF
  
positive.code:             IF ( InitInp%WaveDirMin < MINVAL(MnDriftData%Data4D%WvDir1) ) THEN
                 CALL SetErrStat( ErrID_Fatal,' Minimum wave direction required of '//TRIM(Num2LStr(InitInp%WaveDirMin))//' is not'//&
                       'found in the WAMIT data file '//TRIM(MnDriftData%Filename)//' for the first wave direction.', &
                       ErrStat, ErrMsg, 'MnDrift_InitCalc')
              ENDIF
              IF ( InitInp%WaveDirMax < MAXVAL(MnDriftData%Data4D%WvDir1) ) THEN
                 CALL SetErrStat( ErrID_Fatal,' Maximum wave direction required of '//TRIM(Num2LStr(InitInp%WaveDirMax))//' is not'//&
                       'found in the WAMIT data file '//TRIM(MnDriftData%Filename)//' for the first wave direction.', &
                       ErrStat, ErrMsg, 'MnDrift_InitCalc')
              ENDIF
  
  
  
negative.code: 
              ! If we need to create the 3D data set, then
           CALL Copy_InitData4Dto3D( NewmanAppData%Data4D, NewmanAppData%Data3D, ErrStatTmp, ErrMsgTmp )
           CALL SetErrStat(ErrStatTmp,ErrMsgTmp,ErrStat,ErrMsg,RoutineName)
           IF (ErrStat >= AbortErrLev) RETURN
  
           NewmanAppData%DataIs3D = .TRUE.       ! Set flag to indicate we now have the 3D data.
  