function [ Gumb_log_mom ] = Gumb_mom(X,Y,mm,momX,momY)
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
for j = 1:numel(Y)
    for i = 1:numel(X)
        Gumb_log_mom(j,i) = exp( -( exp( -mm*(Y(j)-momY(1))/momY(2) ) +...
            exp( -mm*(X(i)-momX(1))/momX(2) ) ).^(1/mm) );
    end
end

end