def curvature_mean(curvatures_per_phase):
    """ Given the curvatures along pp for each phase, calculate mean curvature
    """
    # Mean curvature per phase (1 value per phase)
    mean_curvature_per_phase = []
    for curvatures in curvatures_per_phase:
        mean_curvature_per_phase.append(float(curvatures.mean()))
    
    # stats
    mean_curvature_phase_min = min(mean_curvature_per_phase)
    mean_curvature_phase_max = max(mean_curvature_per_phase)
    mean_curvature_phase_change = mean_curvature_phase_max - mean_curvature_phase_min
    
    return (mean_curvature_per_phase, mean_curvature_phase_min, mean_curvature_phase_max, 
            mean_curvature_phase_change )