? MJDGGES requires two square real matrices of the same

In dr1.m you should have (lines 209 to 217):

    if ~isreal(jacobia_)
        if max(max(abs(imag(jacobia_)))) < 1e-15
            jacobia_ = real(jacobia_);
        else
            info(1) = 6;
            info(2) = sum(sum(imag(jacobia_).^2));
            return
        end
    end

If the jacobian matrix is complex dr1 returns info(1)=6. Can you post (or send) your mod file if you do not fall in the trap ?

Stéphane.