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