function [ nspecSQ, nspecMQ, swSQ, swMQ, Spectrum, n_spec ] = read_ga_data2( file_name )
fp = fopen(file_name);
n_spec = fread(fp,1,'int');
nspecMQ = fread(fp,1,'int');
nspecSQ = fread(fp,1,'int');
swMQ = fread(fp,1,'double');
swSQ = fread(fp,1,'double');
%input = zeros(nspecSQ*nspecMQ*n_spec,1);
input = fread(fp,nspecSQ*nspecMQ*n_spec,'double');
Spec=reshape(input,nspecSQ,nspecMQ,n_spec);
for i=1:n_spec
    Spectrum(:,:,i)=Spec(:,:,i)';
end
end