[原创]（请教）非稳态计算，使用那个宏将求出的流场随时间输出？

lijuanwm

[watermark]非稳态计算，使用那个宏将求出的流场值随时间输出？
是DEFINE_ADJUST还是DEFINE_EXECUTE_AT_END?
DEFINE_EXECUTE_AT_END的变量是domain;而我想在某个面上循环求出值。
是否先在domain中找到thread,然后在thread中的face循环。（二维轴对称计算）
具体如下;
#include "stdio.h"
#include "udf.h"
#include "metric.h"
#include "math.h"

DEFINE_EXECUTE_AT_END(Load)
{
  Domain *d;
  Thread *x_tip;
  Thread *t0;

  cell_t c0;
  face_t f0;
   
  real A[ND_ND];
  real Load;
  FILE *fp = fopen("pintleLoad", "w");
  d=Get_Domain(1);
  x_tip=Lookup_Thread(d,7);
  Load=0.;

Message(" eeeee hhh \n");
  
thread_loop_f(x_tip, d)/* loops over all face threads in a domain*/
{  
begin_f_loop(f0, x_tip)    /* loops over faces in a face thread  */
{   

     c0 = F_C0(f0,x_tip);
     t0 = THREAD_T0(x_tip);

     F_AREA(A,f0,x_tip);
  
   
     Load +=F_P(f0,x_tip)* A[1];
   
}                        
end_f_loop(f0, x_tip)
     Message(" fffff  \n");
     Message("Load=%g (N) \n ",Load);
  if (fp != NULL)
    {
      fprintf (fp, "%e", Load);
      fclose (fp);
    }
  else
    Message ("\nWarning: cannot write %s file", Load);
  
}
}   
请各位帮忙[/watermark][br][br][以下内容由 lijuanwm 在 2006年05月26日 03:32pm 时添加] [br]
请帮忙[br][br][以下内容由 lijuanwm 在 2006年05月26日 03:32pm 时添加] [br]
请帮忙