half_2_full_cc_2D.m
59fa66532c42master

half_2_full_cc_2D.m
View Options

			function [ FF_f ] = half_2_full_cc_2D( FF_h )
			%half_2_full_cc_2D Retrieve full domain frequency from half domain
			% Meant for complex conjugate symmetric fields

			[Nkr,Nkz] = size(FF_h);

			if Nkr > Nkz
			FF_f = zeros(max(size(FF_h)));
			FF_f(1:Nkr,1:Nkz) = FF_h;

			for ikr = 1:Nkr
			for ikz = Nkr/2+2:Nkr
			FF_f(ikr,ikz) = FF_f(Nkr-ikr+1,Nkr-ikz+1);
			end
			end
			else
			FF_f = zeros(max(size(FF_h)));
			FF_f(1:Nkr,1:Nkz) = FF_h;

			for ikz = 1:Nkz
			for ikr = Nkz/2+2:Nkz
			FF_f(ikr,ikz) = FF_f(Nkz-ikr+1,Nkz-ikz+1);
			end
			end
			end
			end