I noticed that too, and I think that's a bug, in fact.  At least if nx != ny, which is never the case currently, as you noticed.  I think the phase factor between gamma and kappa should be defined in terms of k, not ik.  But given that we create the PSR with nx=ngrid and ny=ngrid, it's not a problem currently.