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