graph_ohm_jd

 This function plots the characteristics of ohmic data

0001 function graph_ohm_jd(equil,ohm,axs,font)
0002 %
0003 % This function plots the characteristics of ohmic data
0004 %
0005 if nargin < 4
0006     font = struct;
0007 end
0008 if nargin < 3 || ~ishandle(axs(1))
0009     figure
0010     ax = gca;
0011 else
0012     ax = axs(1);
0013 end
0014 %
0015 vloop = ohm.epsi*(2*pi*equil.Rp);
0016 %
0017 ylim = [min([vloop(:);0]),max([vloop(:)+eps;0])]*1.2;
0018 %
0019 graph1D_jd(ohm.rho,vloop,0,0,'r/a','V_{loop} (V)','',NaN,[0,1],ylim,'-','none','r',2,font,ax);
0020 %
0021 if nargin < 3 || ~ishandle(axs(2))
0022     figure
0023     ax = gca;
0024 else
0025     ax = axs(2);
0026 end
0027 %
0028 E_c = Ecritical_jd(equil,ohm.rho);
0029 EoEc = abs(ohm.epsi)./E_c;
0030 %
0031 graph1D_jd(ohm.rho,EoEc,0,0,'r/a','|E|/E_c','',NaN,[0,1],NaN,'-','none','r',2,font,ax);
0032