求助：平均场协同角的UDF函数

dujuan826

各位高手：
我的问题是：我要求平均场协同角，在单机8核双CPU服务器上运行，用的是解释型，出现如下错：
mpt_read:failed:errno=10054
mpt_read:error:read failed trying to read 4 bytes:Invalid argument
MPI Application rank 0 exited before MPI_Finalize() with cortex received a fatal signal .

不知如何解决，程序序如下：
#include "udf.h"
DEFINE_ON_DEMAND(JIAO)
{
   real j,m,volume,vol_tot=0,stavg=0;
   Domain *d;
   Thread *t;
   cell_t c;
   d=Get_Domain(1);
   thread_loop_c(t,d)
    {
      begin_c_loop(c,t)
      {
        j=C_T_G(c,t)[0]*C_U(c,t)+C_T_G(c,t)[1]*C_V(c,t)+C_T_G(c,t)[2]*C_W(c,t);
        m=sqrt(C_T_G(c,t)[0]*C_T_G(c,t)[0]+C_T_G(c,t)[1]*C_T_G(c,t)[1]+C_T_G(c,t)[2]*C_T_G(c,t)[2])*sqrt(C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t)+C_W(c,t)*C_W(c,t));
        volume=C_VOLUME(c,t);
        stavg+=acos(j/m)*volume;
        vol_tot+=volume;
     }

    end_c_loop(c,t)
    stavg/=vol_tot;
    printf("\n Stavg=%g \n",stavg);
}

}



一定重谢，急求！

蜗牛0708

您好~敢问您弄懂了么？我是UDF小白，想请教一下你