GKSguamam_NFLRscan.m
41e7a23b16aenew_coulomb

GKSguamam_NFLRscan.m
View Options

			clear all; close all;
			dir = '/home/bjfrei/Documents/MoliSolver/COSOlver/iCa/';
			f0 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_3_Pmaxe_-1_Jmaxe_-1_Pmaxi_3_Jmaxi_3_pamaxx_10_NFLR_0_kperp_0.5000.h5';
			f3 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_3_Pmaxe_-1_Jmaxe_-1_Pmaxi_3_Jmaxi_3_pamaxx_10_NFLR_3_kperp_0.5000.h5';
			f5 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_3_Pmaxe_-1_Jmaxe_-1_Pmaxi_3_Jmaxi_3_pamaxx_10_NFLR_5_kperp_0.5000.h5';
			f7 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_3_Pmaxe_-1_Jmaxe_-1_Pmaxi_3_Jmaxi_3_pamaxx_10_NFLR_7_kperp_0.5000.h5';
			f10 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_3_Pmaxe_-1_Jmaxe_-1_Pmaxi_3_Jmaxi_3_pamaxx_10_NFLR_10_kperp_0.5000.h5';

			Pmax = 3;
			Jmax = 3;
			bari = @(pp,jj) (Jmax +1)*pp + 1 + jj ;
			NFLR = [0,1,2,3,4,6,7,8];
			Cpj= zeros(length(NFLR),1);
			p =1;
			j =0;
			for iNFLR=1:length(NFLR)
			fi = ['self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_3_Pmaxe_-1_Jmaxe_-1_Pmaxi_3_Jmaxi_3_pamaxx_10_NFLR_',num2str(NFLR(iNFLR)),'_kperp_0.5000.h5']
			CiNFLR = h5read(fullfile(dir,fi),'/Caapj/Ciipj');
			Cpj(iNFLR) = CiNFLR(1,1)
			end

			figure; hold all;
			h1 = plot(NFLR(1:end-1),abs(Cpj(1:end-1) - Cpj(end)));
			set(gca,'YScale','log')