Function IntegratePhase 

pub unsafe extern "C" fn IntegratePhase(
    psi: *mut *mut f32,
    phi: *mut *mut f32,
    flows: *mut *mut c_short,
    nrow: c_long,
    ncol: c_long,
) -> c_int

This function takes row and column flow information and integrates wrapped phase to create an unwrapped phase field. The unwrapped phase field will be the same size as the wrapped field. The array rowflow should have size N-1xM and colflow size NxM-1 where the phase fields are NxM. Output is saved to a file.