求随位置不同源项产生速率udf

xintianli

想要模拟底板不同区域CO产生源项，由于利用建模软件DesignModeler划分多个区域效果并不好，于是想用udf实现源项随位置变化而变化。实际运行过程中，udf无报错，但是调节时发现CO源项并没有随位置变化。求大佬帮忙看看#include "udf.h"
// CO产生
DEFINE_SOURCE(CO_consumption, cell,thread,dS,i)  
{
        real r[ND_ND];
        real x,v;
        C_CENTROID(r,cell,thread);
        x=r[0];
        if(20<=x<=150)
        {
        v=0.00001;
        dS[i] = 0;
        }
        else
        {
                v=0;
                dS[i] = 0;
                }
                return v;
}