function [z]=read_data(name)

f=fopen(name);
line1=fgets(f);
line2=fscanf(f,'%f',4);
line3=fscanf(f,'%d',2);
nx=line3(1);
ny=line3(2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
z=fscanf(f,'%f',nx*ny);
fclose(f);

% FOR DEBUG
size_z = size(z)
