00001 close all; 00002 clear all; 00003 00004 load Kp.dat; 00005 Bf=Kp; 00006 Bf1=Kp; 00007 Bf3=Kp 00008 Kp_baseline=1; 00009 Kp_24_max=squeeze(zeros(length(Kp(:,2)),1)); 00010 for i=1:1:length(Kp(:,2)) 00011 if (i>(23+24)) 00012 Kp_24_max(i)=max(Kp(i-(23+24):i,2)); 00013 else 00014 Kp_24_max(i)=max(Kp(1:i,2)); 00015 end 00016 end 00017 00018 B_unnorm=10.^((-1)*Kp(:,2)/Kp_baseline); 00019 Bf(:,2)=B_unnorm(:)* length(B_unnorm)/sum(B_unnorm); 00020 00021 B_unnorm=10.^((-1)*Kp_24_max/Kp_baseline); 00022 Bf2(:,2)=B_unnorm(:)* length(B_unnorm)/sum(B_unnorm); 00023 00024 B_unnorm=10.^((-1)*Kp_24_max/Kp_baseline); 00025 Bf3(:,2)=B_unnorm(:)* length(B_unnorm)/sum(B_unnorm); 00026 00027 00028 save 'bf.dat' Bf -ASCII 00029 save 'bf2.dat' Bf1 -ASCII 00030 save 'bf_strong_boundary_drops.dat' Bf3 -ASCII 00031 figure 00032 plot (Bf(:,1),Bf(:,2),':r') 00033 hold on 00034 plot (Bf2(:,1),Bf2(:,2),'black') 00035 hold on 00036 semilogy (Bf3(:,1),Bf3(:,2), 'LineWidth', 3 ) 00037 hold on 00038