这是个表示物性密度随高度（y方向）变化的udf，大家看看哪里有问题 谢谢

milaoshu

这是个表示物性密度随高度（y方向）变化的udf，大家看看哪里有问题 谢谢
#include "udf.h"

DEFINE_PROPERTY(cell_density,c,t)
{
    real mu_den;
    real x[ND_ND];
    real y;
    begin_c_loop(c,t)
      {
       C_CENTROID(x,c,t);
       y =x[1];
       mu_den=100*y;
      
      }
    end_c_loop(c,t)
    return mu_den;
  }