定义初始条件的udf

#include "udf.h" DEFINE_INIT(domainInititial, d) {  cell_t c;  Thread *t;  real xc[ND_ND];  real x;  real y;  thread_loop_c(t, d)  {   begin_c_loop_all(c, t)   {    C_CENTROID(xc, c, t);    x = xc[0];

>>阅读原文<<