call query_disc( nside, vector0, radius, listpix, nlist[, nest, inclusive] )
Returns the NESTED pixel index of all pixels north of the equatorial line in a map.