gridloop2 in C; header

gridloop2: as gridloop1, but array a is returned

static PyObject *gridloop2(PyObject *self, PyObject *args)
{
  PyArrayObject *a, *xcoor, *ycoor;
  int a_dims[2];
  PyObject *func1, *arglist, *result;
  int nx, ny, i, j;
  double *a_ij, *x_i, *y_j;

  /* arguments: xcoor, ycoor, func1 */
  if (!PyArg_ParseTuple(args, "O!O!O:gridloop2",
			&PyArray_Type, &xcoor,
			&PyArray_Type, &ycoor,
			&func1)) {
    return NULL; /* PyArg_ParseTuple has raised an exception */
  } 
  nx = xcoor->dimensions[0];  ny = ycoor->dimensions[0];