Your work is really great. It helps me a lot. And I have some special needs:
It seems the ge_contourf can't handle the cases when the matrix Z contains 'NaN'. But the matlab version contourf just ignore it entry. Can you do that as well? Thanks. God bless you.