24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

renzy0315

铜虫 (初入文坛)

应助: 3 (幼儿园)
金币: 53.4
帖子: 37
在线: 10.5小时
虫号: 1949253
注册: 2012-08-21
专业: 内流流体力学

[求助] 求助udf的输入输出

我现在在程序里面获得了三个1维数组real x1[1000],x2[1000].pre[1000],想把他们写入一个txt或dat文件，然后在以后的迭代还要从程序里面把x,y,p从那个文件读回来，下面是我写的输出语句，可是输出来数字都是乱的：
FILE *fp;
fp=fopen("faces.txt","w"

;
for(i=0;i<=count-1;i++) //count=240
{
fprintf(fp,"%g ",x1);
fprintf(fp,"%g ",x2);
fprintf(fp,"%g ",pre);
}
fprintf(fp,"\n total number = %d\n",count);
fclose(fp);

我试了把%g改成%f,%d，出来的也是奇奇怪怪的数字，请教各位大神应该怎么写才对呢？另外，读入语句怎么写？

回复此楼

» 猜你喜欢

垃圾破二本职称评审标准已经有18人回复
职称评审没过，求安慰已经有53人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
EST投稿状态问题已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】FLUENT UDF 气液两相流汽蚀汽化潜热(问题详见内容) 已经有21人回复
边界条件UDF求助已经有4人回复
求助：UDF模型的编译问题已经有15人回复
【求助】编一个UDF 已经有3人回复
【求助】udf读取混合物密度失败、及平均流速读取？已经有8人回复
【求助】单机多核fluent并行udf 已经有13人回复
【求助】求UDF高手帮忙看看我的问题已经有18人回复
【求助】Warning: incorrect cg motion UDF value on zone 35 已经有5人回复
【求助】fluent计算中，udf编译通过，初始化出错已经有11人回复
【求助】物料属性可以用UDF修改吗已经有8人回复
【求助】udf出错已经有13人回复

1楼 2012-10-09 16:19:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

renzy0315

铜虫 (初入文坛)

应助: 3 (幼儿园)
金币: 53.4
帖子: 37
在线: 10.5小时
虫号: 1949253
注册: 2012-08-21
专业: 内流流体力学

引用回帖:

7楼: Originally posted by fangdian at 2012-10-10 14:54:00
你在每个fprintf后面加一个Message验证一下此时的值是不是你想要的。如：fprintf(fp,"%g ",x1);后面加：
Message("%g ",x1);
还有，你这段代码写在哪个宏中的？

目前是写在DEFINE_ON_DEMAND里面的，完整代码是：
#include "udf.h"
#include
/*get the pressure of each cell*/

face_t f;
Domain *domain;
int ID=5;
int i;
Thread *t;
FILE *fp;
int n;
real pre[1000];
real x1[1000];
real x2[1000];
real x3[1000];
real x[ND_ND];
DEFINE_ON_DEMAND(get_p)
{
int count=0;/*****20*/
  domain = Get_Domain(1);

  t=Lookup_Thread(domain,ID);

begin_f_loop(f,t)
{
pre[count]=F_P(f,t);
F_CENTROID(x,f,t); //get the position of f
x1[count]=x[0];
x2[count]=x[1];
if (count==120)
Message("\n%g  ",pre[count]);
else Message("%g  ",pre[count]);

count=count+1;
}
end_f_loop(f,t)
fp=fopen("faces.txt","w"

;
for(i=0;i<=count-1;i++)
{
fprintf(fp,"%g ",x1);
fprintf(fp,"%g ",x2);
fprintf(fp,"%g ",pre);
}
fprintf(fp,"\n total number = %d\n",count);
fclose(fp);
Message("count=%d\n",count);
}