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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 程序运行的结果很奇怪，不知道原因是什么
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 904  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-06-16 15:23:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gelitao1234

金虫 (正式写手)
引用回帖:
2楼: Originally posted by pippi6 at 2013-06-16 20:18:20
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))  
     &a ...

万分感谢。我尝试一下看看。
赞一下 回复此楼
3楼2013-06-16 23:56:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与，应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gelitao1234: 金币+20 2013-06-17 10:31:06
引用回帖:
3楼: Originally posted by gelitao1234 at 2013-06-16 23:56:43
万分感谢。我尝试一下看看。...

以后最好在subroutine开始时使用
implicit none
可以帮助你查出许多没有定义、没有赋值、写错变量这样的初级错误
赞一下 回复此楼
4楼2013-06-17 05:54:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gelitao1234

金虫 (正式写手)
非常感谢你，你的建议对我帮助很大
赞一下 回复此楼
5楼2013-06-17 10:30:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +5 懒羊羊轻置玉臀 2026-04-02 5/250 2026-04-02 13:47 by 星空星月
[考研] 265求调剂 +4 梁梁校校 2026-04-01 4/200 2026-04-02 13:06 by yulian1987
[考研] 一志愿9初试366 本双非求调剂 +4 运气来得若有似� 2026-04-02 4/200 2026-04-02 09:56 by guanxin1001
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考研] 805600专硕材料与化工348分求调剂 +6 上学啦！ 2026-04-01 6/300 2026-04-02 07:49 by 2026材料调剂
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +5 what张 2026-04-01 5/250 2026-04-01 22:35 by guanxin1001
[考研] 275求调剂 +12 waltzh 2026-04-01 12/600 2026-04-01 21:44 by chyhaha
[考研] 求调剂推荐 材料 304 +21 荷包蛋hyj 2026-03-26 21/1050 2026-04-01 21:09 by lijunpoly
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 【求调剂】085601材料工程专硕 | 总分272 | +10 脚滑的守法公民 2026-03-27 10/500 2026-04-01 17:23 by pies112
[考研] 英一数二生物信息学287分，本科生物科学，求调剂 +3 碧水xyz 2026-03-29 4/200 2026-04-01 15:35 by plum
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +7 @taotao 2026-03-30 7/350 2026-04-01 14:30 by chenqifeng666
[考研] 085410人工智能 初试316分 求调剂 +3 残星拂曙 2026-03-31 3/150 2026-04-01 11:09 by 小熊raider
[考研] 考研材料工程351分调剂 +5 整个好的 2026-03-31 5/250 2026-04-01 09:36 by topgun2009
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 0856 335分 +9 cccchenso 2026-03-29 9/450 2026-03-31 16:37 by lishahe
[考研] 279求调剂 +12 j的立方 2026-03-29 12/600 2026-03-30 20:30 by dick_runner
[考研] 340求调剂 +6 Amber00 2026-03-26 6/300 2026-03-29 12:06 by 无际的草原
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
信息提示
关闭
请填处理意见
关闭