% combinemattomat_fixeddt.m
% A MATLAB script for combining the files saved by running "calcsumcurr_manyareagsynmediumtau_parts_fixeddt.py 20 0.025 0.000042 10000 10000 2 myseed 200"
% Expects that the variable "myseed" has been initialized
% Tuomo MakiMarttunen, 20142016
synloctype = 2;
nsegs = 20;
dt = 0.025;
tstop = 10000;
Nsynlocs = 10000;
singleSimT = 200;
syngmaxes = [nan 0.000042 nan];
syngmax = syngmaxes(synloctype);
Nsims = floor(1.0*tstop/singleSimT+0.9999);
Nparts = 13;
ina = [];
ik = [];
ica = [];
ih = [];
il = [];
VtimesA = [];
imemb = [];
Vsoma = [];
icap = [];
times = [];
dt_int = 0.1;
tic
for isim=1:Nsims
disp(['isim=' num2str(isim)]);
disp(['Loading isim=' num2str(isim) ', myseed=' num2str(myseed) ', toc=' num2str(toc)]);
A = load(['currsums_parts_' num2str(Nsynlocs) 'areagsynsmediumtau_fixeddt_type' num2str(synloctype) '_amp' num2str(syngmax) '_tstop' num2str(tstop) '.0_nseg' num2str(nsegs) '_dt' num2str(dt) '_seed' num2str(myseed) '_sim' num2str(isim1) 'x200.0.mat']);
times = [times;A.times];
ina = [ina, A.ina];
ik = [ik, A.ik];
ica = [ica, A.ica];
ih = [ih, A.ih];
il = [il, A.il];
VtimesA = [VtimesA, A.VtimesA];
imemb = [imemb, A.imemb];
icap = [icap, A.icap];
Vsoma = [Vsoma; A.Vsoma];
if isim==1
ts_syn = A.ts_syn;
part_syn = A.part_syn;
end
end
clear A
clear isim
clear Nsims
clear dt_int
clear ts
save(['currsums_parts_' num2str(Nsynlocs) 'areagsynsmediumtau_fixeddt_type' num2str(synloctype) '_amp' num2str(syngmax) '_tstop' num2str(tstop) '.0_nseg' num2str(nsegs) '_dt' num2str(dt) '_seed' num2str(myseed) '_comb200.0.mat']);
