proc_rundke_splitting

 script proc_rundke_splitting

0001 % script proc_rundke_splitting
0002 %
0003 clear all
0004 %
0005 Dr0 = 1;
0006 %
0007 load(['DKE_RESULTS_TScyl_LH_karney_3D_splitting_',num2str(Dr0),'.mat'])
0008 %
0009 % total current calculation
0010 %
0011 ifig = 0;
0012 %
0013 %color = {'b','g','r','m'};
0014 %style = {'-','-','-','-'};
0015 %marker = {'none','none','none','none'};
0016 %lsize = [2,2,2,2];
0017 %
0018 %red = 0.85;
0019 %bspace = 0.55;
0020 %
0021 color = {'k','k','k','k'};
0022 style = {'-','--','-.','none'};
0023 marker = {'none','none','none','o'};
0024 lsize = [0.5,2,2,2];
0025 %
0026 red = 0.9;
0027 bspace = 0.7;
0028 %
0029 xlim = [0,tfinal];
0030 ylim = [1e-1,1e2];
0031 xlab = 't/\tau_c';
0032 ylab = 'I(t) (kA)';
0033 tit = '';%str;
0034 %
0035 for icoll = 1:ncoll,
0036     for iopt = 1:nopt,
0037         %
0038         ifig = ifig + 1;
0039         %
0040         str = ['D_r = ',num2str(Dr0),' m^2/s; coll. mode = ',num2str(coll_mode_list(icoll)),'; op. split = ',num2str(optsplit_list(iopt))];
0041         figure(ifig),clf
0042         set(gcf,'name',str)
0043         %
0044         for idtn = 1:ndtn - 1, 
0045             graph1D_jd(dtn_list(idtn):dtn_list(idtn):tfinal,abs(I{idtn,iopt,icoll}),1,1,'','','',NaN,xlim,ylim,style{idtn},marker{idtn},color{idtn},lsize(idtn),20,gca);
0046             %
0047             leg{idtn} = ['\Deltat/\tau_c = ',num2str(dtn_list(idtn))];
0048         end
0049         leg{length(dtn_list)} = ['\Deltat/\tau_c = ',num2str(dtn_list(end))];
0050         leg = [leg,'Location','SouthEast'];
0051         %
0052         graph1D_jd(dtn_list(idtn):dtn_list(idtn):tfinal,abs(I{idtn,iopt,icoll}),1,1,xlab,ylab,tit,leg,xlim,ylim,style{ndtn},marker{ndtn},color{ndtn},lsize(ndtn),20,gca,red,0.7,bspace);
0053         %
0054         set(gca,'ytick',10.^[-1:2])
0055         set(gca,'xtick',10.^[0:5])
0056         set(gca,'xminorgrid','off')
0057         set(gca,'yminorgrid','off')
0058         %
0059         print_jd(p_opt,['Fig_splitting_',num2str(Dr0),'_coll_',num2str(coll_mode_list(icoll)),'_split_',num2str(optsplit_list(iopt))],'./figures',ifig)
0060         %
0061     end
0062 end
0063 %