Draw_B_Ellip_gt118.m
a494effb5a96master

Draw_B_Ellip_gt118.m
View Options

			% Computational grid
			z=linspace(0.225,0.525,100);
			r=linspace(0.00,0.10,100);
			[Z,R] = meshgrid(z,r);


			% Default current in the coils
			I = [9.85 -15.80 33.59 33.21];
			%final filename
			name='gt118_magnet';

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


			% Plot the magnetic field lines
			f=figure;
			axis equal
			[~,cont1]=contour(Z,R,R.*Aphi,15,'r:');
			hold on
			contourf(Z,R,sqrt(Bz.^2+Br.^2))
			c=colorbar;
			c.Label.String='\|B\| [T]';
			xlabel('Z [m]')
			ylabel('R [m]')



			f.PaperUnits='centimeters';
			f.PaperSize=[12,8];
			xlabel('z [m]')
			ylabel('r [m]')
			print(f,'phiBprofile_gt118','-dpdf','-fillpage')
			savefig(f,'phiBprofile_gt118')
			hold off


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