graph_transp_jd

 This function plots the characteristics of transpfaste data

0001 function graph_transp_jd(equil,transpfaste,axs,font)
0002 %
0003 % This function plots the characteristics of transpfaste data
0004 %
0005 if nargin < 4,
0006     font = struct;
0007 end
0008 if nargin < 3 || any(~ishandle(axs)),
0009     figure
0010     axs(1) = gca;
0011     figure
0012     axs(2) = gca;
0013 end
0014 %
0015 xrho = equil.ptx(:,1)/equil.ptx(end,1);
0016 %
0017 xDr = transpfaste.Dr0*(1 + transpfaste.pDr*xrho.^2);
0018 xVr = transpfaste.Vr0*(1 + transpfaste.pVr*xrho.^2);
0019 %
0020 ylim1 = [min([xDr(:);0]),max([xDr(:)+eps;0])]*1.2;
0021 ylim2 = [min([xVr(:);0]),max([xVr(:)+eps;0])]*1.2;
0022 %
0023 graph1D_jd(xrho,xDr,0,0,'r/a','D_r (m^2/s)','',NaN,[0,1],ylim1,'-','none','r',2,font,axs(1));
0024 graph1D_jd(xrho,xVr,0,0,'r/a','V_r (m/s)','',NaN,[0,1],ylim2,'-','none','r',2,font,axs(2));