function [X,Z,C_slice]=open_single_dat_file()

Nx=201;
Nz=101;

n_lines=Nx*Nz;

filename='poisson_slice.dat';
fid=fopen(filename);

X=0*(1:n_lines);
Z=0*(1:n_lines);
C_slice=0*(1:n_lines);

for i=1:n_lines
    c1=fgetl(fid);
    vec_temp=sscanf(c1,'%f');
    x_temp=vec_temp(1);
    y_temp=vec_temp(2);
    C_temp=vec_temp(3);
    
    X(i)=x_temp;
    Z(i)=y_temp;
    C_slice(i)=C_temp;
end

X=reshape(X,Nx,Nz);
Z=reshape(Z,Nx,Nz);
C_slice=reshape(C_slice,Nx,Nz);

fclose(fid);
    
end