Draw_B_Ellip_bifreq.m
6b37315661933D

Draw_B_Ellip_bifreq.m
View Options

			% Computational grid
			z=linspace(0.,1.2,2000);
			r=0;
			[Z,R] = meshgrid(z,r);


			% Default current in the coils
			I = [69.5248 -4.8];
			%final filename
			name='bifreq';

			%allocate memory for speed
			res=zeros([size(Z),3]);
			for i=1:size(R,2)
			res(:,i,:) = B_Ellip_bifreq({'aphi','bz','br'},'normal',I,R(:,i),Z(:,i));
			end
			Aphi=res(:,:,1);
			Bz=res(:,:,2);
			Br=res(:,:,3);


			% Plot the magnetic field lines
			f=figure;
			plot(Bz,z)
			ylabel('Z [m]')
			xlabel('B_z [T]')

			grid on
			set(gca,'fontsize',14)
			f.PaperUnits='centimeters';
			f.PaperSize=[12,8];
			f.Position=[1 1 0.3 1.5].*f.Position;
			namefig='phiBprofile_bifreq';
			exportgraphics(f,sprintf('%s.png',name),'Resolution',300)
			print(f,name,'-dpdf','-fillpage')
			set(f, 'Color', 'w');
			exportgraphics(f,sprintf('%s.eps',name))
			savefig(f,'phiBprofile_bifreq')
			hold off


			%% Save magnetic field to disk
			save=true;
			overwrite=true;
			if save
			savemagtoh5([name,'.h5'],r,z,Aphi,Br,Bz,overwrite);
			end