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

[Nkx,Nky] = size(FF_h);

if Nkx > Nky
    FF_f = zeros(max(size(FF_h)));
    FF_f(1:Nkx,1:Nky) = FF_h;

    for ikx = 1:Nkx
        for iky = Nkx/2+2:Nkx
            FF_f(ikx,iky) = FF_f(Nkx-ikx+1,Nkx-iky+1);
        end
    end
else
    FF_f = zeros(max(size(FF_h)));
    FF_f(1:Nkx,1:Nky) = FF_h;

    for iky = 1:Nky
        for ikx = Nky/2+2:Nky
            FF_f(ikx,iky) = FF_f(Nky-ikx+1,Nky-iky+1);
        end
    end
end
end