版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

lijunjiexm

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 22.8
红花: 1
帖子: 182
在线: 87.1小时
虫号: 1198843
注册: 2011-02-03
专业: 常微分方程与动力系统

[求助] 请教一个程序小问题，谢谢

程序如下：

PROGRAM PROCESS3

C *** Program to draw velocity profiles.
C
   IMPLICIT    NONE
   INTEGER*4 ITIME, I, II, NUMNP, NTIME, COUNTER, IDUM, IDUM1
   INTEGER*4 WAIT
   PARAMETER (NUMNP = 17, NTIME = 10)
   REAL*4    STORE1, STORE2, VX(21), VY(21), DUM1
   REAL*4    XNODE(NUMNP), YNODE(NUMNP)
C
   OPEN (10, FILE = 'ns2dev.ve2')
   OPEN (11, FILE = 'mesh.dat')
   OPEN (12, FILE = 'vprof.dat')
c
c read nodes
c
   READ(11,*,END=300) IDUM, DUM1
   DO I = 1 , NUMNP
      READ(11, *, END = 300) IDUM, XNODE(I), YNODE(I)
   ENDDO
c
c read velocity field
c
   WAIT = 4
   READ(10,*,END=300) IDUM, DUM1
   DO ITIME = 1 , NTIME
      READ(10,*,END=300) IDUM, DUM1
      IF (ITIME .GT. WAIT) COUNTER = COUNTER + 1
      DO I = 1 , NUMNP
      READ(10, *, END=300) IDUM, STORE1, STORE2
      IF (ITIME .GT. WAIT) THEN
         IF ((I .GE. 7) .AND. (I .LE. 17)) THEN
            VX(I - 7 + 1) = VX(I - 7 + 1) + STORE1
            VY(I - 7 + 1) = VY(I - 7 + 1) + STORE2
         ENDIF
      ENDIF
      ENDDO
   ENDDO
300  CONTINUE
   DO II = 1 , 11
      VX(II) = VX(II) / COUNTER
      VY(II) = VY(II) / COUNTER
      WRITE (12, 1001) YNODE(II), VX(II), VY(II),' y vx vy'
1001 FORMAT (3(E14.5),A)
   ENDDO
   STOP
   END

为什么counter的数值在debug过程一直是零？ Parameter的用法有问题吗？
附：ns2dev.ve2文件：
2047       5000  numnp , itmax
         0  0.000000000000000E+000  it , time
1    0.00000E+00 0.00000E+00 node , vx , vy
2    0.00000E+00 0.00000E+00 node , vx , vy
3    0.00000E+00 0.00000E+00 node , vx , vy
4    0.00000E+00 0.00000E+00 node , vx , vy
5    0.00000E+00 0.00000E+00 node , vx , vy
6    0.00000E+00 0.00000E+00 node , vx , vy
7    0.00000E+00 0.00000E+00 node , vx , vy
8    0.00000E+00 0.00000E+00 node , vx , vy
9    0.00000E+00 0.00000E+00 node , vx , vy
10    0.00000E+00 0.00000E+00 node , vx , vy
11    0.00000E+00 0.00000E+00 node , vx , vy
12    0.00000E+00 0.00000E+00 node , vx , vy
13    0.00000E+00 0.00000E+00 node , vx , vy
14    0.00000E+00 0.00000E+00 node , vx , vy
15    0.00000E+00 0.00000E+00 node , vx , vy
16    0.00000E+00 0.00000E+00 node , vx , vy
17    0.00000E+00 0.00000E+00 node , vx , vy

mesh.dat
17       11  numnp , numel
1    0.00000E+00 0.00000E+00 1 node, x , y , ref. num.
2    0.00000E+00 0.35803E-03 1 node, x , y , ref. num.
3    0.00000E+00 0.70174E-03 1 node, x , y , ref. num.
4    0.00000E+00 0.10317E-02 1 node, x , y , ref. num.
5    0.00000E+00 0.13485E-02 1 node, x , y , ref. num.
6    0.00000E+00 0.16526E-02 1 node, x , y , ref. num.
7    0.00000E+00 0.19445E-02 1 node, x , y , ref. num.
8    0.00000E+00 0.22247E-02 1 node, x , y , ref. num.
9    0.00000E+00 0.24938E-02 1 node, x , y , ref. num.
10    0.00000E+00 0.27521E-02 1 node, x , y , ref. num.
11    0.00000E+00 0.30000E-02 1 node, x , y , ref. num.

[ Last edited by lijunjiexm on 2011-11-13 at 11:03 ]

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

1楼2011-11-13 10:44:41

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

zcius

至尊木虫 (知名作家)

应助: 5 (幼儿园)
金币: 16742.4
散金: 810
红花: 8
帖子: 8863
在线: 406.9小时
虫号: 1123409
注册: 2010-10-15
专业: 计算机科学的基础理论

【答案】应助回帖

★
余泽成(金币+1): 谢谢参与应助！ 2011-11-13 21:39:53

COUNTER 没有赋初值能行吗？
做除法前不做判断？

赞一下(1人)

回复此楼

2楼2011-11-13 11:05:09

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

lijunjiexm

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 22.8
红花: 1
帖子: 182
在线: 87.1小时
虫号: 1198843
注册: 2011-02-03
专业: 常微分方程与动力系统

ns2dev.ve2文件若改为：
2047       5000  numnp , itmax
         0  0.000000000000000E+000  it , time
  1    0.10296E+12  -0.13594E+10 node , vx , vy
2    0.14944E+11  -0.13925E+10 node , vx , vy
3    0.70706E+11  -0.29397E+09 node , vx , vy
4    0.77626E+11 0.58172E+09 node , vx , vy
5    0.62009E+11 0.70982E+09 node , vx , vy
6    0.00000E+00 0.00000E+00 node , vx , vy
7    -0.55675E+11 0.00000E+00 node , vx , vy
8    0.17093E+11 0.00000E+00 node , vx , vy
9    -0.27307E+11 0.00000E+00 node , vx , vy
  10    0.16832E+12 0.00000E+00 node , vx , vy
11    0.14461E+11 0.00000E+00 node , vx , vy
  12    0.72585E+11 0.00000E+00 node , vx , vy
  13    -0.26770E+11 0.00000E+00 node , vx , vy
  14    0.28571E+11 0.00000E+00 node , vx , vy
  15    0.48992E+11 0.00000E+00 node , vx , vy
  16    0.36732E+11 0.00000E+00 node , vx , vy
  17    0.00000E+00 0.00000E+00 node , vx , vy

counter最好结果居然不为零，且值为NTIME-wait=6

赞一下

回复此楼

3楼2011-11-13 11:10:42

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

lijunjiexm

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 22.8
红花: 1
帖子: 182
在线: 87.1小时
虫号: 1198843
注册: 2011-02-03
专业: 常微分方程与动力系统

引用回帖:

2楼: Originally posted by zcius at 2011-11-13 11:05:09:
COUNTER 没有赋初值能行吗？
做除法前不做判断？

初值不是默认是零吗？你看把把文件ns2dev.ve2改成上面的形式，counter就有数值了，这个很奇怪，我想不通，帮忙看一下哈，谢谢

赞一下

回复此楼

4楼2011-11-13 11:15:28

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

maomao1210

金虫 (正式写手)

程序强帖: 5
应助: 2 (幼儿园)
金币: 1431.3
散金: 242
红花: 16
沙发: 1
帖子: 991
在线: 441.5小时
虫号: 253215
注册: 2006-05-20
性别: MM
专业: 考古理论

★ ★
余泽成(金币+2): 谢谢参与应助！ 2011-11-13 21:40:12

引用回帖:

1楼: Originally posted by lijunjiexm at 2011-11-13 10:44:41:
程序如下：

PROGRAM PROCESS3

C *** Program to draw velocity profiles.
C
   IMPLICIT    NONE
   INTEGER*4 ITIME, I, II, NUMNP, NTIME, COUNTER, IDUM, IDUM1
   INTEGER*4 WAIT
...

不知道该表达的意思是什么，有太多的错误。。。。
变量定义类型有很多不匹配
读文件的时候，有太多的end =300, 其中读完文件mesh.dat，程序就跳到了300行了，而文件mesh.dat的数据显然没有用到， xnode,ynode.
貌似楼主是修改他人的程序，但貌似没理清楚逻辑顺序，而且变量混乱的很。

赞一下(1人)

回复此楼

5楼2011-11-13 14:42:04

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

相关版块跳转我要订阅楼主 lijunjiexm 的主题更新

返回列表