    function [seti,uTotRX,uIncRX,uScaRX] = matchIncFieldTrans(filename,frequencyHz)

    % Read Fresnel data
    %filename = 'inexpdata/fresnel_opus_1/twodielTM_8f.exp';
    [uTotRX, uIncRX, frequencies, rTX, nTX, rRX, nRX] = readRAWData(filename);

    [frequencyId,frequencyHz] = freqID(frequencies,frequencyHz); % Choose data to specific frequency, e.g. 5 GHz via 5*1E9
    seti.k = freqToWaveNumber(frequencyHz); % Compute wave number

    % Choose fields at specific frequency and conjugate the fields to adapt
    % from time dependence exp(iwt) to exp(-iwt):
    [uTotRX, uIncRX, uScaRX] = uSca(uTotRX,uIncRX,frequencyId);

    % Transfer experimental set-up into structural array seti
    seti.radSrc = rTX;
    seti.incNb = nTX;
    seti.radMeas = rRX;
    seti.measNb = nRX;

    % Experimental Setup
    seti.dim = 2;               % two dimensional problem
    % Set transmitters positions:
    seti.incType = 'pointSource';
    seti.incPntsType = 'circle';
    seti = setIncPnts(seti);    % seti.incPnts contains the coordinates
    % Set receivers positions:
    seti.measType = 'nearField';
    seti.measPntsType = 'circle';
    seti = setMeasPnts(seti);   % seti.measPnts contains the coordinates

    % Set grid
    seti.rCD = 0.2;         % size of computational domain [-rCD,rCD)^dim
    seti.nCD = 256;         % number of discretization points for each dimension of CD
    seti = setGrid(seti);   % grid is stored in seti.gridROI and discretization points in each dimension in seti.nROI

    end