[Python]对fluent输出的无序的压力分布系数进行排序

问题描述：利用fluent计算翼型，提取翼型压力分布系数时，得到的压力分布排列顺序混乱无规则，本文以计算超临界翼型RAE2822为例提取压力系数分布。

按下图输出翼型的压力分布：

则可得到翼型的压力系数分布如下图：

但由于与实验对比时，实验数据用散点图、验证结果通常使用线性图，由于输出数据排列无序，若直接画线性图，图像无序混乱。

利用如下pyhton代码，按照压力系数变化的梯度对其进行排列。

def cp_sort(cp):
    indicates_x=np.argsort(-cp.iloc[:,1])
    cp=cp.T[indicates_x].T
    cp = cp.reset_index(drop = True)#重新索引
    [m,n]=cp.shape
    num=m-1
    x=cp.iloc[:,1]
    cpp=cp.iloc[:,3]
    x_up=[]
    cp_up=[]
    x_down=[]
    cp_down=[]
    x_up.append(x[0])
    cp_up.append(cpp[0])
    for i in range(1,num) :
        if (cpp[i]-cpp[i-1])/abs(cpp[i]) < -0.25 :
            x_up.append(x[i])
            cp_up.append(cpp[i])
        elif abs(x[i]!=x[i-1]):
            x_down.append(x[i])
            cp_down.append(cpp[i])
    x_down=x_down[::-1]
    cp_down=cp_down[::-1]
    x_down[-1]=x[0]
    cp_down[-1]=cpp[0]
    X=x_up+x_down
    CP=cp_up+cp_down
    X=np.array(X)
    CP=np.array(CP)
    data=np.column_stack((X,CP))
    return data

 对其重新排序后，可得到压力分布图像如下：

不同翼型得到的压力分布有所不同，所以根据其梯度变化 来进行排序的方法鲁棒性可能不高，仅提供自己进行排序的一种方法，具体操作可根据具体压力分布来进行调整。