query_disc issue

qdtest.py

NSIDE=8
vec = [ 0.17101007,  0.03015369,  0.98480775]
qd=healpy.query_disc(NSIDE, vec, 6, deg=True, nest=False)
m = np.zeros(healpy.nside2npix(NSIDE))
m[qd]=1
healpy.mollview(m)
print(qd)
#array([ 3,  0, 11,  4])

#IDL
#HIDL> query_disc, 8, [ 0.17101007,  0.03015369,  0.98480775],6,listpix,/DEG,NESTED=0           
#HIDL> print,listpix
#           4

#HIDL> query_disc, 8, [ 0.17101007,  0.03015369,  0.98480775],6,listpix,/DEG,NESTED=0,/inclusive
#HIDL> print,listpix
#           0           3           4           5          11          12          13          23