function [ fig ] = plot_metric( data, options )

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
NPLOT = options.SHOW_FLUXSURF + options.SHOW_METRICS;
if NPLOT > 0
    fig = figure;
    if options.SHOW_METRICS
    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
    if options.SHOW_FLUXSURF
        R = [geo_arrays(1,:,12),geo_arrays(1,1,12)];
        Z = [geo_arrays(1,:,13),geo_arrays(1,1,13)];
    plot(R,Z); 
    axis equal
    end
end
end