%% HeLaZ data
filename = 'results_00.h5';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Load the data
moment = 'Ni00';

kx       = h5read(filename,['/data/var2d/' moment '/coordkx']);
ky       = h5read(filename,['/data/var2d/' moment '/coordky']);
timeNi   = h5read(filename,'/data/var2d/time');
Nipj     = zeros(numel(timeNi),numel(kx),numel(ky));
for it = 1:numel(timeNi)
    tmp          = h5read(filename,['/data/var2d/', moment,'/', num2str(it,'%06d')]);
    Nipj(it,:,:) = tmp.real + 1i * tmp.imaginary; 
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Plot growth rate vs ky
K_RICCI = 1; %% add a sqrt(1+tau) to the kperps
ikx     = 1; %% Fix the kx value
plot_gamma_vs_k;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Plot Ni00 evolution
ikx     = 1; %% Fix the kx value
plot_Ni00_t_evolution;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%