ft_defaults
fnam='/home/pellizze/tzagarak/msitransfer/Sall_f20_5_pttest_025_normean_alltrial_41515_voxels_allroi_clustdics_cue_pwranalytic_negc_finf_labmatrix.txt';% begin: various preparative things no calls to fieldtrip functions
lab=h_prepare_plilab(fnam,[1 3 33 57 59 61 65]);

sub=[13 15 16 17 18 19 23 26 28 29];
outf=cell(size(sub,2),1,codistributor('1d',1));
szsub=size(sub,2);
     
% end  
     
     
     
for a=drange(1:szsub)
    ft_defaults
[outf(a,1),~]=h_multitime_phase(1,sub(1,a));% this part works fine on a non-parallel call - contains a  lot of fieldtrip code
end
outff=gather(outf);
tst=h_save_phase_par(outff,'/home/pellizze/tzagarak/analysis/results');