| 查看: 1549 | 回复: 16 | ||
[求助]
求助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，出来的也是奇奇怪怪的数字，请教各位大神应该怎么写才对呢？另外，读入语句怎么写？ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有52人回复
-
时间戳又变了8-15
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】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人回复
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 64
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
3楼2012-10-09 19:16:09
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 64
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
2楼2012-10-09 19:14:04
4楼2012-10-10 14:37:09
5楼2012-10-10 14:37:39
6楼2012-10-10 14:45:08
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 64
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
7楼2012-10-10 14:54:00
|
目前是写在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); } |
8楼2012-10-10 15:31:12
9楼2012-10-10 15:35:42
10楼2012-10-10 16:02:26
