openTauContourFile Subroutine

    subroutine openTauContourFile
        implicit none

        ! Only root writes out the trajectory
        if (iProcIndex == root) then
            iunitCycLog = get_free_unit()
            call get_unique_filename('tauContour', .true., .true., 0, trajFile)

            open(iunitCycLog, file=trajFile, status='new')
        end if

    end subroutine openTauContourFile