想定义一个随空间变化的扩散系数

shadow501

有一关于电场计算的方程，在矩形块内与外场要定义两个不同扩散系数，我试着写了这样一个但是调用不起来啊，求高手来帮我分析分析啊
矩形域为0<x<0.003,-0.001<y<0
DEFINE_DIFFUSIVITY(my_diff,c,t,p)
{
    real xc[ND_ND];
       
    cell_t c;
    real x1=1;
    real x2=2.7;
    begin_c_loop(c,t)
   {
   C_CENTROID(xc,c,t);
     if((xc[0]>0)&&(xc[0]<0.0030))
     {
       if((xc[1]>0)&&(xc[1]<-0.0010))
       {
         return x2;
       }
       else return x1;
     }
     else return x1;
   }
   end_c_loop(c,t)
}