run('met3Toff47G_2Hz.m');                %simulated data from jmrui-NMR-SCOPE-B
DATAc = DATA(:,1)+1i.*(DATA(:,2));
DATAr = reshape(DATAc,[SIZE(1) SIZE(3)]);

sz = size(DATAr);
metab_list = [{'Cho'},{'Cr'},{'NAA'}];


for mt=1:sz(2)

mystruct.sim = DATAr(:,mt);
mystruct.fids = DATAr(:,mt);
mystruct.seq = {'PRESS'};
mystruct.name = metab_list(mt);
mystruct.spectralwidth = 1.0000e+03;
mystruct.Bo = 3.0000;
mystruct.te = 135;
mystruct.n = 1024;
mystruct.linewidth=2.0000;
mystruct.t= -linspace(1.0000e-03, 1024*1.0000e-03, 1024);
mystruct.sz = [1024 1];
mystruct.centerFreq = 0;


RF=io_writelcmraw(mystruct,[cell2mat(metab_list(mt)),'.RAW'],cell2mat(metab_list(mt)));

end
