FLUENT UDF怎么获取某个面液相的面积

[梅尔比恩]

在做液滴撞击的模拟，因为要用到动态接触角所以需要获取三相点位置以及速度。目前是想用两步迭代的位移差来计算速度，但是现在不知道怎么去算这个三相点的位置。代码是这么写的

    Thread* thread = Lookup_Thread(domain, 9); 
    Thread** pt = THREAD_SUB_THREADS(thread); 
    cell_t cell; 
    face_t f; 
    real x[ND_ND], area[ND_ND]; 
    double max_x = 0, f_Hoff_inverse, x_hoff, temp, Ca, volume; 
    int n; 
    sum = 0.0; 
    begin_c_loop_all(cell, pt[1]) 
    { 
        if (C_VOF(cell, pt[1]) != 0) 
        { 
            C_CENTROID(x, cell, pt[1]); 
            if (x[0] > max_x) 
            max_x = x[0]; 
        } 
        sum += C_VOF(cell, pt[1]); 
    } 
    end_c_loop_all(cell, pt[1])

我想问一下怎么把SUM变成面积啊