DrawFRCs

Contents

Draw FRCs
Draw 1/4 of the FRC - inside
Outside the FRC (no mirror)
Draw the whole FRC

Draw FRCs

See also FRCFieldLines, DrawFRCLines, NewFig

%--------------------------------------------------------------------------
%	 Copyright (c) 2019 Princeton Satellite Systems, Inc.
%  All rights reserved.
%--------------------------------------------------------------------------

Draw 1/4 of the FRC - inside

kappa = 5;
r = linspace(0,1);
z = linspace(0,kappa);
[R,Z] = meshgrid(r,z);
psi = FRCFieldLines( R, Z, 1, 1, kappa );
NewFig('FRCFieldLines')
contour(R,Z,psi,linspace(0,max(psi(:)),15));
ylabel('z')
zlabel('r')
title(['Solov''ev Model'])
grid on

Outside the FRC (no mirror)

kappa = 2.5;
r = linspace(0,2);
z = linspace(0,2*kappa);
[R,Z] = meshgrid(r,z);
[~,psiOut] = FRCFieldLines( R, Z, 1, 1, kappa );
NewFig('FRCFieldLines - Outside')
contour(Z',R',psiOut',100);

Draw the whole FRC

kappa = 2.5;
rs = 1;
DrawFRCLines( rs, kappa );


%--------------------------------------
% $Date$
% $Id: e6e6cb9d56dfdb5e776f42717c386bcd8ae263f5 $