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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

fish8900

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 701.9
散金: 103
帖子: 105
在线: 34小时
虫号: 1496643
注册: 2011-11-17
性别: GG
专业: 零件加工制造

[求助] fortran的debug调试

1、调试过程中设置断点后，程序运行到断点后停止，接下来让程序继续运行应该按哪个图标啊？
2、有一个程序会循环很多次，前面都不会出问题，到后面某一步会出现错误，对于这种情况应该怎么设置断点？例如
do time=1,2000

PW(1)=1
DO 203 w1=1,100
DO 202 w0=1,100
AXA1=0.0
bxb1=0.0
   DO I=NIF1,NIF2
   DO J=0,JF(I)
   S1(I,J)=S(I,J)
   end do
   end do

   DO I=NIB1,NIB2
   DO J=0,JB(I)
   BS1(I,J)=BS(I,J)
   end do
   end do

do w2=1,20

  DO I=NIF1,NIF2

      S(I,-1)=S(I,1)
   S(I+1,-1)=S(I+1,1)
      S(I-1,-1)=S(I-1,1)

      DO J=0,JF(I)
         if(t(i,j,0)<=tm)then
   s(i,j)=0.0
   else
            SX=0.5*(S(I+1,J)-S(I-1,J))/DI
            SY=0.5*(S(I,J+1)-S(I,J-1))/DJ
      if((abs(sx)>=1.0e+7).or.(abs(sy)>=1.0e+7))then
   s(i,j)=tl
   bs(i,j)=0.0
else

         SYY=(S(I,J+1)-2.0*S(I,J)+S(I,J-1))/(DJ*DJ)
   SXX=(S(I+1,J)-2.0*S(I,J)+S(I-1,J))/(DI*DI)
SXY=0.25*(S(I+1,J+1)-S(I-1,J+1)-S(I+1,J-1)+S(I-1,J-1))/(DI*DJ)
            SS1=1.0+SX*SX
SS2=1.0+SY*SY
            SS3=SQRT(1.0+SX**2+SY**2)**3
SS4=2.0*SX*SY*SXY
           E1=SS3*P9*G9
     E2=SS3*(P(I,J)-PW(w1))-SGM(I,J,0)*SS4
G1=SGM(I,J,0)*SS1/(DJ*DJ)
G2=SGM(I,J,0)*SS2/(DI*DI)

S(I,J)=(E2+G1*(S(I,J+1)+S(I,J-1))+G2*(S(I+1,J)+S(I-1,J)))/(E1+2.0*G1+2.0*G2)
         end if
            end if


   if(w2/=20)goto 200

   AXA1=AXA1+S(I,J)*DI*DJ*1.0E+3
   AXA(w1)=AXA1

200          end do
   end do

end do

do w2=1,20

  DO I=NIB1,nib2

      BS(I,-1)=BS(I,1)
      BS(I+1,-1)=BS(I+1,1)
      BS(I-1,-1)=BS(I-1,1)

      DO J=0,JB(I)
         if(t(i,j,mz)<=tm)then
   bs(i,j)=0.0
   else
            BSX=0.5*(BS(I+1,J)-BS(I-1,J))/DI
            BSY=0.5*(BS(I,J+1)-BS(I,J-1))/DJ
            if((abs(bsx)>=1.0e+7).or.(abs(bsy)>=1.0e+7))then
   s(i,j)=tl
   bs(i,j)=0.0
else
         BSXX=(BS(I+1,J)-2.0*BS(I,J)+BS(I-1,J))/(DI*DI)
BSYY=(BS(I,J+1)-2.0*BS(I,J)+BS(I,J-1))/(DJ*DJ)
         BSXY=0.25*(BS(I+1,J+1)-BS(I-1,J+1)-BS(I+1,J-1)+BS(I-1,J-1))/(DI*DJ)
            BSS1=1.0+BSX*BSX
BSS2=1.0+BSY*BSY
BSS3=SQRT(1.0+BSX**2+BSY**2)**3
BSS4=2.0*BSX*BSY*BSXY

        BE1=BSS3*P9*G9
BE2=BSS3*(-PW(w1)-P9*G9*TL)+SGM(I,J,mz)*BSS4
BG1=SGM(I,J,mz)*BSS1/(DJ*DJ)
BG2=SGM(I,J,mz)*BSS2/(DI*DI)

   BS(I,J)=(BE2-BG1*(BS(I,J+1)+BS(I,J-1))-BG2*(BS(I+1,J)+BS(I-1,J)))/(BE1-2.0*BG1-2.0*BG2)
            end if

            end if

   if(w2/=20)goto 201

   BXB1=BXB1+BS(I,J)*DI*DJ*1.0E+3
   BXB(w1)=BXB1

201       end do
   end do

end do

   DO I=NIF1,NIF2
      DO J=0,JF(I)

      S2=S(I,J)-S1(I,J)
      IF(ABS(S2/MAX(ABS(S1(I,J)),1.0E-5)).GE.1.0E-3) GO TO 202
   end do
   end do

   do I=NIB1,NIB2
      do J=0,JB(I)
      BS2=BS(I,J)-BS1(I,J)
      IF(ABS(BS2/MAX(ABS(BS1(I,J)),1.0E-5)).GE.1.0E-3) GO TO 202
end do
   end do

GOTO 204

202 continue

204  IF(w1<3) THEN
      PW(w1+1)=PW(w1)+10
      ELSE
      SD1=(AXA(w1-1)-BXB(w1-1))-(AXA(w1)-BXB(w1))
   SD2=(PW(w1-1)-PW(w1))

   if(sd1==0.0)then
sd1=1.0E-10
end if

   PW(w1+1)=PW(w1)-(AXA(w1)-BXB(w1))*SD2/SD1

END IF

   if(w1.GT.6) THEN
   !SDDX=AXA(w1)-BXB(w1)
   SDDX=BXB(w1)-AXA(W1)

   IF(ABS(SDDX).LE.1.0E-8)then
   GOTO 205
   ELSE
   GOTO 203
      end if
   end if

203 continue

当time=1000时会出现出现报错，但前面都是好好的，我想从1000开始观察sqrt()函数里面的数值的变化，断点应该怎么设置啊
3、出现无穷大的数值时，无穷大的数能不能和具体的某个数进行比较啊？就像上面那个程序里的sx如果出现为无穷大，sx能和1.0e+8进行比较？

回复此楼

» 收录本帖的淘帖专辑推荐

Fortran程序

» 猜你喜欢

288环境专硕,求调材料方向已经有39人回复
353求调剂已经有4人回复
一志愿211电子信息347求调剂已经有3人回复
323求调剂已经有4人回复
285求调剂已经有7人回复
一志愿211，化学学硕，310分，本科重点双非，求调剂已经有9人回复
22408 266求调剂已经有4人回复
313求调剂已经有3人回复
305求调剂已经有3人回复
求调剂已经有11人回复

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

1楼 2013-04-25 14:46:15

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

oven1986

至尊木虫 (知名作家)

钓鱼岛岛主兼琉球地下党总书记

应助: 29 (小学生)
贵宾: 1.068
金币: 28863.4
散金: 9415
红花: 154
沙发: 14
帖子: 9876
在线: 1252.2小时
虫号: 1920718
注册: 2012-08-02
性别: GG
专业: 分离过程

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2013-04-29 08:13:54

先回答一个吧如果要调试找出错误，可以用print,write,pause等语句写出来，观察是否能运行出来。这样慢慢就能逼近发生错误的地方。放置点肯定是完整的do..enddo语句之前之后。

赞一下

回复此楼

优秀是一种习惯！

2楼2013-04-25 16:55:57

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

fish8900

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 701.9
散金: 103
帖子: 105
在线: 34小时
虫号: 1496643
注册: 2011-11-17
性别: GG
专业: 零件加工制造

你说的方法我也试过哈，出现的错误是出现了NAN，我也不太懂为什么会出现这个，所以想请教下。还有就是实在无法避免的话，用什么判据能继续接着运行后面的程序

赞一下

回复此楼

3楼2013-04-25 21:25:01

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

zhangghost

银虫 (初入文坛)

应助: 1 (幼儿园)
金币: 414
红花: 1
帖子: 12
在线: 26.5小时
虫号: 628364
注册: 2008-10-17
专业: 岩土力学与岩土工程

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2013-04-29 08:13:40

你应该首先说明你的编译环境，CVF6.5还是IVF，当然这些也不重要

1. 断点处继续运行直接按F5，如果是CVF的话是感叹号右边的书页形符号，IVF点三角符号。

2.解决方法很多，给一个简单易懂的吧，添加一个if语句，在里面添加断点
  if(time=1000) then
   !在这里随便写个语句然后添加断点
  endif
3. 出现NAN基本上是你代码中有一些小问题，先解决掉。

赞一下

回复此楼

4楼2013-04-27 16:32:31

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

fish8900

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 701.9
散金: 103
帖子: 105
在线: 34小时
虫号: 1496643
注册: 2011-11-17
性别: GG
专业: 零件加工制造

引用回帖:

4楼: Originally posted by zhangghost at 2013-04-27 16:32:31
你应该首先说明你的编译环境，CVF6.5还是IVF，当然这些也不重要

1. 断点处继续运行直接按F5，如果是CVF的话是感叹号右边的书页形符号，IVF点三角符号。

2.解决方法很多，给一个简单易懂的吧，添加一个if语句， ...

出现NAN我找了很长时间的错误，都没能发现是哪出错了，所以很纠结啊！其实主要的还是第三个问题一直困扰着我

赞一下

回复此楼

5楼2013-04-27 20:55:31

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

gaolugang

木虫 (小有名气)

应助: 12 (小学生)
金币: 1680.4
帖子: 132
在线: 65.1小时
虫号: 250274
注册: 2006-05-13
性别: GG
专业: 电力电子学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fish8900: 金币+30, ★★★很有帮助 2013-05-08 15:02:47

1. 看你的IDE，一般是<F5>, 或者点<运行>那个按钮。
2. 用PAUSE 语句来实现。
do time=1,2000
...
...
if (time.eq.100) then
pause
...
...
enddo
3. NaN (Not a Number) 出现的情况，是因为你的程序中有0/0的运算。找到这个bug就可以解决。
至于楼主说的无穷大，我没有理解是什么意思。 NaN 不是无穷大，他的意思是Not a Number.

赞一下

回复此楼

6楼2013-05-03 01:42:48

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

spoweryang

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 295.2
红花: 2
帖子: 136
在线: 92.9小时
虫号: 2445307
注册: 2013-05-02
专业: 安全科学与工程

【答案】应助回帖

查一查你的变量定义类型，在计算中是否出现常量与变量运算不一致的情况

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

fish8900

7楼2013-05-03 04:55:03

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

fish8900

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 701.9
散金: 103
帖子: 105
在线: 34小时
虫号: 1496643
注册: 2011-11-17
性别: GG
专业: 零件加工制造

送红花一朵

引用回帖:

7楼: Originally posted by spoweryang at 2013-05-03 04:55:03
查一查你的变量定义类型，在计算中是否出现常量与变量运算不一致的情况

常量与变量运算不一致指的是什么意思啊？

赞一下

回复此楼

8楼2013-05-08 15:00:12

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

fish8900

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 701.9
散金: 103
帖子: 105
在线: 34小时
虫号: 1496643
注册: 2011-11-17
性别: GG
专业: 零件加工制造

引用回帖:

6楼: Originally posted by gaolugang at 2013-05-03 01:42:48
1. 看你的IDE，一般是<F5>, 或者点<运行>那个按钮。
2. 用PAUSE 语句来实现。
do time=1,2000
...
...
if (time.eq.100) then
pause
...
...
enddo
3. NaN (Not a Number) 出现的情况，是 ...

那应该具体如何找呢？能不能详细的说一下？

赞一下

回复此楼

9楼2013-05-08 15:02:34

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

gaolugang

木虫 (小有名气)

应助: 12 (小学生)
金币: 1680.4
帖子: 132
在线: 65.1小时
虫号: 250274
注册: 2006-05-13
性别: GG
专业: 电力电子学

【答案】应助回帖

引用回帖:

9楼: Originally posted by fish8900 at 2013-05-08 15:02:34
那应该具体如何找呢？能不能详细的说一下？...

我更倾向于吧变量的值打印出来，写到一个log文件里。然后观察在哪一步，哪个值出的问题。你也可以，一步一步的调试，找到出错的那条语句，出错的变量。

赞一下

回复此楼

10楼2013-05-09 01:57:37

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

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

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 284求调剂 +13	梵@@ 2026-04-06	13/650	2026-04-07 21:53 by l_paradox
[考研] 0857大类环境工程B区求调剂 +3	龚禹铭 2026-04-05	3/150	2026-04-06 10:22 by 蓝云思雨
[考研] 一志愿南昌大学，085600，344分求调剂 +9	调剂上岸玘 2026-04-05	10/500	2026-04-06 09:30 by dongzh2009
[考研] 070300化学学硕311分求调剂 +11	梁富贵险中求 2026-04-04	13/650	2026-04-06 07:24 by houyaoxu
[考研] 329求调剂 +17	miaodesi 2026-04-02	20/1000	2026-04-05 18:33 by 蓝云思雨
[考研] 282求调剂 +3	aaa车辆 2026-04-01	3/150	2026-04-05 17:03 by yulian1987
[考研] 288求调剂，一志愿华南理工大学071005 +6	ioodiiij 2026-04-04	6/300	2026-04-05 10:09 by guoweigw
[考研] 0835学硕299求调剂 08大类可接受 +5	useryy 2026-04-03	5/250	2026-04-04 20:07 by 蓝云思雨
[考研] 一志愿沪9，求生物学调剂，326分 +6	刘墨墨 2026-04-04	6/300	2026-04-04 19:44 by 唐沐儿
[考研] 280求调剂 +21	咕噜晓晓 2026-04-02	22/1100	2026-04-04 11:12 by 猪会飞
[考研] 一志愿双非085502，267分，过四级求调剂 +3	再忙也要吃饭啊 2026-04-03	3/150	2026-04-04 05:03 by gswylq
[考研] 一志愿南昌大学324求调剂 +13	hanamiko 2026-04-01	13/650	2026-04-03 18:30 by ls刘帅
[考研] 考研调剂 +8	不爱喝饮料 2026-04-03	8/400	2026-04-03 16:40 by Mistake-J
[考研] 求调剂22408 288分 +5	new382 2026-04-02	5/250	2026-04-03 09:13 by 醉在风里
[考研] 085600，320分求调剂 +6	大馋小子 2026-04-02	6/300	2026-04-02 21:54 by dongzh2009
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4	1050389037 2026-04-02	5/250	2026-04-02 21:40 by dongzh2009
[考研] 283求调剂 +3	jiouuu 2026-04-02	4/200	2026-04-02 14:08 by 哒哒哒呱呱呱
[考研] 348环境工程调剂 +3	吴彦祖24k 2026-04-01	3/150	2026-04-02 09:14 by nanaliuyun
[考研] 一志愿北交大材料工程，总分358 +4	cs0106 2026-04-01	4/200	2026-04-02 07:42 by 尚水阁主
[考研] 安全工程 285 求调剂 +3	Xinyu56 2026-04-01	4/200	2026-04-01 21:50 by 静静静静静静静�