bounce_ripple_jd

 Draws the bounce and ripple loss domains

0001 function bounce_ripple_jd(bounce_mode,pnmax,xmhubounce2,ir_display,Zmripple,fontsize,ax)
0002 %
0003 % Draws the bounce and ripple loss domains
0004 %
0005 if nargin < 7,
0006     ax = gca;
0007 end
0008 if nargin < 6,
0009     fontsize = 20+14i;
0010 end
0011 %
0012 nmhu_ripple = 100;
0013 %
0014 if bounce_mode >= 1,%Draw bounce domain boundaries
0015     %
0016     pparmax = pnmax*sqrt(xmhubounce2(ir_display));
0017     pperpmax = pnmax*sqrt(1 - xmhubounce2(ir_display));
0018     %
0019     graph1D_jd([0,pparmax],[0,pperpmax],0,0,'','','',NaN,NaN,NaN,'-','none','k',2,fontsize,ax);
0020     graph1D_jd([0,-pparmax],[0,pperpmax],0,0,'','','',NaN,NaN,NaN,'-','none','k',2,fontsize,ax);
0021     %
0022     if isempty(Zmripple) == 0,
0023         %if isnan(Zmripple) == 0,%Draw the magnetic ripple loss cone
0024             %
0025             xpndetrap = Zmripple.xpndetrap;
0026             xmhusupertrap = cos((90 - Zmripple.xthetaloss)*pi/180);
0027             %
0028             mhu_ripple = linspace(-xmhusupertrap(ir_display),xmhusupertrap(ir_display),nmhu_ripple);           
0029             pparsupertrap = xpndetrap(ir_display)*mhu_ripple;%For contour plot of the supertrapped domain
0030             pperpsupertrap = xpndetrap(ir_display)*sqrt(1-mhu_ripple.^2);%For contour plot of the supertrapped domain
0031             %
0032             pparmin = xpndetrap(ir_display)*xmhusupertrap(ir_display);
0033             pperpmin = xpndetrap(ir_display)*sqrt(1-xmhusupertrap(ir_display)^2);
0034             pparmax = pnmax*xmhusupertrap(ir_display);
0035             pperpmax = pnmax*sqrt(1-xmhusupertrap(ir_display)^2);
0036             %
0037             graph1D_jd(pparsupertrap,pperpsupertrap,0,0,'','','',NaN,NaN,NaN,'--','none','k',2,fontsize,ax);
0038             graph1D_jd([pparmin,pparmax],[pperpmin,pperpmax],0,0,'','','',NaN,NaN,NaN,'--','none','k',2,fontsize,ax);
0039             graph1D_jd([-pparmin,-pparmax],[pperpmin,pperpmax],0,0,'','','',NaN,NaN,NaN,'--','none','k',2,fontsize,ax);
0040         %end
0041     end
0042 end
0043 %
0044 %--------------------------------------------------------------------------
0045 %