function [ sensor_data, sensor] = COMP_SPLdata( kgrid, medium, source )
%UNTITLED13 Summary of this function goes here
% ========================================================================
% MICROPHONES2
% =========================================================================
% MICROPHONES POSITIONS
%sensor.mask = [1, 1, 1, kgrid.Nx, kgrid.Ny,kgrid.Nz].';
sensor.mask = zeros(kgrid.Nx,kgrid.Ny,kgrid.Nz);
sensor.mask(1:kgrid.Nx, 1:kgrid.Ny, round(kgrid.Nz/2))=1;
sensor.record = {'p','p_max','p_rms'};
% define the frequency response of the sensor elements
%bandwidth = 23;  
%sensor.frequency_response = [center_freq, bandwidth ];
% =========================================================================
% COMPUTING 2
% =========================================================================
display_mask = source.p_mask;
input_args = {'DisplayMask', display_mask, 'PMLInside', false, 'PlotPML', false,'PlotSim', false};
sensor_data = kspaceFirstOrder3D(kgrid, medium, source, sensor, input_args{:});

end

