| 查看: 925 | 回复: 6 | ||
gelitao1234金虫 (正式写手)
|
[求助]
程序运行的结果很奇怪，不知道原因是什么
|
|
程序如下，程序运行结果是TTT的值是NaN，数组TEMP的结果也是这样的，原因是什么？ C C....<< CACULATION OF TEMPERATURE>> SUBROUTINE TEMPER PARAMETER(MX=181,MY=151,MZ=161) IMPLICIT DOUBLE PRECISION(A-H,O-Z) COMMON /PARA3/ NX,NY,NZ COMMON /PARA4/ DT,TIME COMMON /TEMP/ TEMP(0:MX+2,0:MY+2,0:MZ+2) &, TEMPH(0:MZ+2) &, SGST(0:MX+2,0:MY+2,0:MZ+2) COMMON /FLOWS/ UU(0:MX+2,0:MY+2,0:MZ+2) &, VV(0:MX+2,0:MY+2,0:MZ+2) &, WW(0:MX+2,0:MY+2,0:MZ+2) COMMON /DELT1/ DXU(0:MX+1),DYU(0:MY+1),DZU(0:MZ+1) COMMON /DELT2/ DXV(0:MX+1),DYV(0:MY+1),DZV(0:MZ+1) COMMON /DELT3/ DXW(0:MX+1),DYW(0:MY+1),DZW(0:MZ+1) COMMON /DELT3/ DXT(0:MX+1),DYT(0:MY+1),DZT(0:MZ+1) COMMON /DELT3/ DXP(0:MX+1),DYP(0:MY+1),DZP(0:MZ+1) C....<INITIAL CONDITION > DO K=0,NZ+2 DO J=0,NY+2 DO I=0,NX+2 TEMP(I,J,K)=1.2D0 END DO END DO END DO DO J=0,NY+2 DO I=0,NX+2 TEMP(I,J,2)=6.D0 TEMP(I,J,1)=6.D0 END DO END DO C....<< TEMPERATURE B.C.>> DO K=1,NZ+1 DO J=1,NY+1 TEMP(1,J,K)=TEMP(2,J,K) TEMP(NX+1,J,K)=TEMP(NX,J,K) END DO END DO DO J=1,NY+1 DO I=1,NX+1 TEMP(I,J,1)=TEMP(I,J,2) TEMP(I,J,NZ+1)=TEMP(I,J,NZ) END DO END DO DO K=1,NZ+1 DO I=1,NX+1 TEMP(I,1,K)=TEMP(I,2,K) TEMP(I,NY+1,K)=TEMP(I,NY,K) END DO END DO C C....<TEMPERATURE CACULATION > DO K=2,NZ DO J=2,NY DO I=2,NX C....<<SPACING LENGTH>> H1=DXP(I-1) H2=DXP(I) H3=DYP(J-1) H4=DYP(J) H5=DZP(K-1) H6=DZP(K) DX=DXU(I-1) DY=DYV(J-1) DZ=DZW(K-1) C....<<CACULATION>> YIX=DT*UU(I,J,K)*(TEMP(I+1,J,K)-TEMP(I-1,J,K)) YIY=DT*VV(I,J,K)*(TEMP(I,J+1,K)-TEMP(I,J-1,K)) YIZ=DT*WW(I,J,K)*(TEMP(I,J,K+1)-TEMP(I,J,K-1)) ERX=TEMP(I+1,J,K)/(DX*H2)-TEMP(I,J,K)*(1.D0/(DX*H1)+1.D0/(DX*H2)) ERY=TEMP(I,J+1,K)/(DY*H4)-TEMP(I,J,K)*(1.D0/(DY*H3)+1.D0/(DY*H4)) ERY=TEMP(I,J,K+1)/(DZ*H6)-TEMP(I,J,K)*(1.D0/(DZ*H5)+1.D0/(DZ*H6)) TEMP(I,J,K)=TEMP(I,J,K)+YIX/(DXT(I)+DXT(I+1)) & +YIY/(DYT(J)+DYT(J+1)) & +YIZ/(DZT(K)+DZT(K+1)) & -DT*SGST(I,J,K)*(ERX+TEMP(I-1,J,K)/(DX*H1)) & -DT*SGST(I,J,K)*(ERY+TEMP(I,J-1,K)/(DY*H3)) & -DT*SGST(I,J,K)*(RYZ+TEMP(I,J,K-1)/(DZ*H5)) END DO END DO END DO C....<<HORIZENTAL AVERAGING >> DO K=0,NZ+1 TEMPH(K)= 0 END DO DO K=1,NZ+1 TEMPH(K-1)= TEMPH(K)/DFLOAT(NX+1)*(NY+1) DO J=1,NY+1 DO I=1,NX+1 TEMPH(K)=TEMPH(K)+TEMP(I,J,K) END DO END DO END DO C....<<AVERAGING & SUBTRACTION>> TTT=0.D0 DPLAY=0.D0 DO K=1,NZ+1 DO J=1,NY+1 DO I=1,NX+1 TTT=TTT+TEMP(I,J,K) END DO END DO END DO DO K=1,NZ+1 DO J=1,NY+1 DO I=1,NX+1 TEMP(I,J,K)=TEMP(I,J,K)-TTT/DFLOAT((NX+1)*(NY+1)*(NZ+1)) END DO END DO END DO DPLAY=TTT/DFLOAT((NX+1)*(NY+1)*(NZ+1)) WRITE(*,601)DPLAY 601 FORMAT(' TTT',F8.1) RETURN END |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
论文撤稿了
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
河北省自然科学基金
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有8人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 并行计算没有输出？怎么回事
已经有16人回复
- PID速度调节，实测时得到的图形很奇怪，是什么原因
已经有9人回复
- 关于VASP一些琐事的求助（K点，吸附优化等）
已经有12人回复
- 【讨论】A2/O工艺，采用间歇进水、间歇曝气为何系统还能正常运行
已经有16人回复
- 【求助】运行ms的时候跳出“由于另一个程序正在运行中，此操作无法完成。请选择切换…
已经有28人回复
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
【答案】应助回帖
感谢参与，应助指数 +1
|
ERY=TEMP(I,J,K+1)/(DZ*H6)-TEMP(I,J,K)*(1.D0/(DZ*H5)+1.D0/(DZ*H6)) TEMP(I,J,K)=TEMP(I,J,K)+YIX/(DXT(I)+DXT(I+1)) & +YIY/(DYT(J)+DYT(J+1)) & +YIZ/(DZT(K)+DZT(K+1)) & -DT*SGST(I,J,K)*(ERX+TEMP(I-1,J,K)/(DX*H1)) & -DT*SGST(I,J,K)*(ERY+TEMP(I,J-1,K)/(DY*H3)) & -DT*SGST(I,J,K)*(RYZ+TEMP(I,J,K-1)/(DZ*H5)) 这两段里，首先ERY可能是ERZ（重复了）？其次，RYZ是ERZ？ 你的RYZ完全没有赋值啊。 |
2楼2013-06-16 20:18:20
gelitao1234
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1787.6
- 散金: 366
- 红花: 4
- 帖子: 577
- 在线: 1339.1小时
- 虫号: 1206537
- 注册: 2011-02-19
- 性别: GG
- 专业: 边界层大气物理学和大气湍
3楼2013-06-16 23:56:43
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
4楼2013-06-17 05:54:43
gelitao1234
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1787.6
- 散金: 366
- 红花: 4
- 帖子: 577
- 在线: 1339.1小时
- 虫号: 1206537
- 注册: 2011-02-19
- 性别: GG
- 专业: 边界层大气物理学和大气湍
5楼2013-06-17 10:30:57
gelitao1234
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1787.6
- 散金: 366
- 红花: 4
- 帖子: 577
- 在线: 1339.1小时
- 虫号: 1206537
- 注册: 2011-02-19
- 性别: GG
- 专业: 边界层大气物理学和大气湍
6楼2013-06-17 10:31:15
gelitao1234
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1787.6
- 散金: 366
- 红花: 4
- 帖子: 577
- 在线: 1339.1小时
- 虫号: 1206537
- 注册: 2011-02-19
- 性别: GG
- 专业: 边界层大气物理学和大气湍
7楼2013-06-17 10:31:37
10