[求助] 求教一个UDF中源项设置的问题

DEFINE_ADJUST(adj_uds0,mixturedomain)
{
        Thread *thread;
        cell_t cell;
        Domain *gasdomain;

        gasdomain = DOMAIN_SUB_DOMAIN(mixturedomain,1);
        thread_loop_c(thread, gasdomain)
        {
                begin_c_loop(cell, thread)
                {
                        C_UDSI(cell, thread, 0)=C_VOF(cell, thread);                               
                }
                end_c_loop (cell, thread)
        }
}
先是通过这个获得一个全流场的气含率信息。
然后：
DEFINE_SOURCE(lateralforce_cfd1,cell,thread,dS,eqn)
{
        Thread *mix_thread, *thread_gas, *thread_liq;
        float kl,  vof_gas, y_vel_liq, y_vel_gas, vel_slip;
        float x_vel_liq, x_vel_gas, y_vel_slip, x_vel_slip, klg;
        float FTD=0, CWi=0.1, source=0,sign=0;       
        float delta_gas_a=0, CTD=0.7, Cvm=0.5;
        float x[ND_ND];
        float Re, CD, Eo, FCL, CL, FEo;
        float Eo2, FTD2, FTD3, vf_ratio, VF_TRANS=0.14;
        FILE *fp_gas;
        float current_time=1;
            
        mix_thread= THREAD_SUPER_THREAD(thread);
        thread_gas= THREAD_SUB_THREAD(mix_thread, 1);
           thread_liq= THREAD_SUB_THREAD(mix_thread, 0);
            x_vel_liq=C_U(cell, thread_liq);
        x_vel_gas=C_U(cell, thread_gas);
        y_vel_liq=C_V(cell, thread_liq);
        y_vel_gas=C_V(cell, thread_gas);
        y_vel_slip=fabs(y_vel_gas-y_vel_liq);
        x_vel_slip=fabs(x_vel_gas-x_vel_liq);
        vel_slip=sqrt(y_vel_slip*y_vel_slip+x_vel_slip*x_vel_slip);       
        vof_gas=C_VOF(cell, thread_gas);
       

            kl=C_K(cell,thread_liq);
        klg=0.5*vof_gas*Cvm*vel_slip*vel_slip;
        delta_gas_a=C_UDSI_G(cell, thread_gas, 0)[1];
        FTD=-CTD*roul*delta_gas_a*(kl+klg);       

   
           if (thread==thread_liq)
           sign=-1.0;       
        else if (thread==thread_gas)
           sign=1.0;       
        else Message("error in lift force");

        source=(FTD)*sign;
        dS[eqn]=0;
        fp_gas=fopen("gas_holdup.dat","a");
        fprintf(fp_gas,"%f\t",current_time);
    return source;
}