PyObject * PyArray_FromDims(int n_dimensions,
                            int dimensions[n_dimensions],
                            int type_num);


int dims[2]; dims[0] = nx; dims[2] = ny;
PyArrayObject *a;  int dims[2];
dims[0] = 10;  dims[1] = 21;
a = (PyArrayObject *) PyArray_FromDims(2, dims, PyArray_DOUBLE);