为什么不能计算受力问题

贝丹

DEFINE_EXECUTE_AT_END(force)
{
Domain *d;                 
Thread *t ;
face_t f;
real   x[ND_ND];
real   f_glob[ND_ND],m_glob[ND_ND],x_cg[ND_ND];
real   total_f[ND_ND];
real   time = RP_Get_Real("flow-time");
int i;
#if !RP_NODE   
FILE *file= fopen("force.txt", "a+");
#endif
if(time>100.0)
{
#if !RP_HOST
  d=Get_Domain(1);
  t=Lookup_Thread(d,6);
  if (!Data_Valid_P ())   
   return;
   NV_S(f_glob, =, 0.0);
   NV_S(m_glob, =, 0.0);
   NV_S(x_cg, =, 0.0);
   NV_S(total_f,   =, 0.0);

  for(i=0;i<=ND_ND;i++)
   x_cg = 0.0;
  Compute_Force_And_Moment (d,t,x_cg,f_glob,m_glob,TRUE);

    #if RP_3D
  {
        total_f[0]=f_glob[0];
        total_f[1]=f_glob[1]-9.8*body_mass;
                total_f[2]=f_glob[2];
    Message ("================================================================);
    Message (" ================================================================);
    Message (" ================================================================);
    Message ("  ===============================================================);
  }
   #endif
#endif /*if!RP_HOST*/
  node_to_host_real(total_f, ND_ND);
  node_to_host_real(f_glob, ND_ND);
  node_to_host_real(m_glob, ND_ND);
  node_to_host_real(x_cg, ND_ND);
   
#if RP_HOST
   if (file != NULL)
       {
        #if(RP_3D)
       {
            fprintf (file, "=============================================");
                fclose (file);
          }
                 #endif
          }
       else
          Message ("\nWarning: cannot write %s file");
   #endif  /*#if RP_HOST×/
}
}

求助，这段程序用了老是会出错，且force.txt 在time 100时没有数据输入，且会卡住不能计算，用的是计算三维物体受力问题