plot_metric.m
View Options

			function [ fig ] = plot_metric( data )

			names = {'Jacobian','gradxB','gradyB','gradzB','gradz_coeff',...
			'gxx','gxy','gyy','gyz','gzz','hatB','hatR','hatZ'};
			geo_arrays = zeros(2,data.Nz,numel(names));

			for i_ = 1:numel(names)
			namae = names{i_};
			geo_arrays(:,:,i_) = h5read(data.outfilenames{end},['/data/metric/',namae])';
			end
			fig = figure;
			subplot(311)
			for i = 1:5
			plot(data.z, geo_arrays(1,:,i),'DisplayName',names{i}); hold on;
			end
			xlim([min(data.z),max(data.z)]); legend('show'); title('GYACOMO geometry');

			subplot(312)
			for i = 6:10
			plot(data.z, geo_arrays(1,:,i),'DisplayName',names{i}); hold on;
			end
			xlim([min(data.z),max(data.z)]); legend('show');

			subplot(313)
			for i = 11:13
			plot(data.z, geo_arrays(1,:,i),'DisplayName',names{i}); hold on;
			end
			xlim([min(data.z),max(data.z)]); legend('show');
			end

gyacomo/matlab/plot/plot_metric.m0597cb3c91dcmaster