Correct nfchoa_order() and added references
Added the possibility of Gauss grids for spherical arrays