24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2042  |  回复: 1

liangsg521

新虫 (初入文坛)

[求助] udf能读入,但是一初始化计算就闪退!

用udf编写了如下程序,希望达到读取速度分量作为大涡模拟的入口条件,可是每次一初始化就闪退,求问各位大神问题出在哪里?

# include <udf.h>
DEFINE_PROFILE(vel_X,t,index)
{
float V_x;
FILE *fp;
char ch;
face_t f;
int i,j;
int count=1;

fp=fopen("E:\LHC\LES\velocity_U","r";
if(fp==NULL){
        printf("Can not open file!\n";       
}
begin_f_loop(f,t)
{
        int step=RP_Get_Integer("time-step";
        for(i=0;i<count;i++){
                for(j=0;j<step+1;j++){
                        fscanf(fp,"%lf%c",&V_x,&ch);
                }
        }
        F_PROFILE(f,t,index)=V_x;
        count++;
}
end_f_loop(f,t);
if(fclose(fp)){
        printf("Can not close the file!\n";
       
}
}
//V_y
DEFINE_PROFILE(vel_y,t,index)
{
float V_y;
FILE *fp;
char ch;
face_t f;
int i,j;
int count=1;

fp=fopen("E:\LHC\LES\velocity_V","r";
if(fp==NULL){
        printf("Can not open file!\n";       
}
begin_f_loop(f,t)
{
        int step=RP_Get_Integer("time-step";
        for(i=0;i<count;i++){
                for(j=0;j<step+1;j++){
                        fscanf(fp,"%lf%c",&V_y,&ch);
                }
        }
        F_PROFILE(f,t,index)=V_y;
        count++;
}
end_f_loop(f,t);
if(fclose(fp)){
        printf("Can not close the file!\n";
       
}
}
//V_z
DEFINE_PROFILE(vel_z,t,index)
{
float V_z;
FILE *fp;
char ch;
face_t f;
int i,j;
int count=1;

fp=fopen("E:\LHC\LES\velocity_W","r";
if(fp==NULL){
        printf("Can not open file!\n";       
}
begin_f_loop(f,t)
{
        int step=RP_Get_Integer("time-step";
        for(i=0;i<count;i++){
                for(j=0;j<step+1;j++){
                        fscanf(fp,"%lf%c",&V_z,&ch);
                }
        }
        F_PROFILE(f,t,index)=V_z;
        count++;
}
end_f_loop(f,t);
if(fclose(fp)){
        printf("Can not close the file!\n";
       
}
}

非常感谢!!!
回复此楼

» 猜你喜欢

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

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

小虫子被鸟吃

新虫 (小有名气)

要不试试,把fscanf放到loop外面试试?
2楼2015-05-27 16:13:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liangsg521 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本9一志愿2 0854低分专硕286求调剂 +10 芒种111 2026-04-04 10/500 2026-04-10 12:31 by luosha500
[考研] 生物学调剂 可调剂到生物与医药 +6 李政莹 2026-04-06 7/350 2026-04-10 10:09 by 314126402
[考研] 268分085602化学工程调剂 +24 月照花林。 2026-04-09 24/1200 2026-04-10 08:09 by Sammy2
[考研] 269电子信息求调剂,可转专业 +9 独酌wl 2026-04-06 9/450 2026-04-09 20:55 by laoshidan
[考研] 0703调剂 +18 拾玖壹 2026-04-04 20/1000 2026-04-09 20:53 by zhouxiaoyu
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-09 20:49 by lbsjt
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +23 熊哥xtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +5 momo(上岸版) 2026-04-06 8/400 2026-04-09 15:07 by only周
[考研] 288求调剂 +15 没有答案_ 2026-04-05 15/750 2026-04-09 10:22 by 5268321
[考研] 土木水利专硕276分求调剂 +6 我想上学!!6 2026-04-05 9/450 2026-04-08 17:45 by 宋小宝HQ
[考研] 298求调剂 +4 manman511 2026-04-05 4/200 2026-04-08 16:50 by tjzhao
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 331求调剂 +5 张元一 2026-04-07 6/300 2026-04-07 22:13 by hemengdong
[考研] 085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干 +3 1021075758 2026-04-06 4/200 2026-04-07 09:23 by 蓝云思雨
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 求调剂 +5 wos666 2026-04-03 5/250 2026-04-06 10:13 by 蓝云思雨
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
信息提示
请填处理意见