code: def calc_2dsplinecoeffs_c(array2d):
      
      #Set up result arrays
      out= copy.copy(array2d)
      out= numpy.require(out,dtype=numpy.float64,requirements=['C','W'])
  
      #Set up the C code
      ndarrayFlags= ('C_CONTIGUOUS','WRITEABLE')
      interppotential_calc_2dsplinecoeffs= _lib.samples_to_coefficients
      interppotential_calc_2dsplinecoeffs.argtypes= [ndpointer(dtype=numpy.float64,flags=ndarrayFlags),
                                                     ctypes.c_int,
                                                     ctypes.c_int]
  
      #Run the C code
      interppotential_calc_2dsplinecoeffs(out,out.shape[1],out.shape[0])
  
      return out
docs: NAME:
         calc_2dsplinecoeffs_c
      PURPOSE:
         Use C to calculate spline coefficients for a 2D array
      INPUT:
         array2d
      OUTPUT:
         new array with spline coeffs
      HISTORY:
         2013-01-24 - Written - Bovy (IAS)