Hi all,

I was trying to use meshgrid with three arrays and got some odd results. Here's a simple example:

xt = np.array([1,2,3,4])
yt = np.array([6,7,8])
zt = np.array([12,13])
xxx,yyy,zzz = np.meshgrid(xt,yt,zt)

So I would expect that xxx[0,0,:] = array([1,2,3,4])

instead I get xxx[0,0,:] = array([1,1]) and xxx[0,:,0] = array([1,2,3,4])

also yyy[:,0,0] = array([6,7,8]), whereas I would expect yyy[0,:,0] = array([6,7,8])

So what's going on? This seems like a bug to me.

Any suggestions for getting what I wanted -- i.e. xxx.shape = (2,3,4), with values as appropriate?

Thanks,

Jon

Jonathan D. Slavin (he/him)

Astrophysicist - High Energy Astrophysics Division

Center for Astrophysics | Harvard & Smithsonian

Office: (617) 496-7981 | Cell: (781) 363-0035
60 Garden Street | MS 04 | Cambridge, MA 02138