| 查看: 3693 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
为什么Fortran调试无错误运行不出结果 已有2人参与
|
||
|
program main parameter n=600 double precision u(0:n),u1(0:n),x0,x1,x,t0,t1,h,dett,lambda,y(0:n) integer i,j x0=0.0d0 x1=6.0d0 t0=0.0d0 lambda=0.3d0 t1=0.5d0 h=0.01d0 dett=0.003d0 x=x0 !do i=0,n,1 !if(i<=250)then !y(i)=1.0 !else ! y(i)=0 !end if ! end do do i=0,n,1 if(x>=1 .and. x<=3) then u(i)=1.0d0 else u(i)=0.0d0 end if x=x+h end do open(1,file="yingfeng.dat"  t=t0+dett x=x0 do while(t<=t1) x=x0+h do j=1,n,1 u1(j)=u(j)-3*lambda*(u(j)-u(j-1)) if(abs(t-0.5)<=0.0001)then write(1,*)x,u1(j),y(j) end if x=x+h end do do j=1,n-1 u(j)=u1(j) u(0)=1 u(1)=0 end do t=t+dett end do stop end 求大神指导，为啥这个程序运行后yingfeng.dat文件里面啥数据都没有 |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有252人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Intel fortran与visual fortran
已经有11人回复
- FORTRAN
已经有0人回复
- Fortran
已经有7人回复
- 已经安装studio，为什么还是不能安装Fortran呢？
已经有0人回复
- 开始学习FORTRAN，结果杀毒软件总是认为生成的exe程序是木马，要求删掉，为什么呢？
已经有7人回复
- fortran
已经有4人回复
- Fortran77 代码调试中提示有Syntax error，但却不知道为什么出错。望高人指点！
已经有4人回复
- fortran循环求助
已经有3人回复
- Fortran 求解复数矩阵SVD时，用到cgesvd，结果与matlab不同，为什么？
已经有15人回复
- fortran里的intrisic function: precision()x结果表示什么意思
已经有5人回复
- fortran 程序修改后改回来，计算结果不同？？？为什么会这样
已经有2人回复
- data 给fortran 的二维数组如何命名？为什么老出错，解决送金币
已经有1人回复
- fortran输出结果为0
已经有3人回复
- fortran 程序有write和没有write两种情况计算的结果为什么会不一样?
已经有14人回复
- fortran输出求助
已经有6人回复
- 帮忙把公式转化为FORTRAN子程序
已经有6人回复
- fortran中计算数据写入文本为 NAN的问题
已经有4人回复
- 刚开始学Fortran 帮忙看一下这个程序 为什么算不出结果
已经有3人回复
- fortran断点
已经有2人回复
- fortran 运行问题
已经有2人回复
- 【求助】请问有什么软件可以将C语言代码转换为fortran吗？
已经有7人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
夜月丨噩梦: 金币+30, ★★★很有帮助 2016-01-05 11:42:26
夜月丨噩梦: 金币+30, ★★★很有帮助 2016-01-05 11:42:26
|
你的程序中最重要的部分如下： double precision ...t0,t1,...,dett... t0=0.0d0 t1=0.5d0 dett=0.003d0 ... t=t0+dett do while(t<=t1) ... if(abs(t-0.5)<=0.0001) then write(1,*)x,u1(j),y(j) endif ... 可以看到，这个程序能否向yingfeng.dat中写入数据的关键是：t的数值相距0.5的差值是否小到0.0001以内，和其它问题完全无关。 如果分析我给出的这些语句的话，就会发现问题所在： t从0.003开始累加，要和0.5比？？？能达到相差小于0.0001的程度么？ 相当于500/3，得到的余数是多少？ 条件达不到，当然不会输出。程序在逻辑上是没错的。 希望楼主能再审视一下相应的条件，让程序输出你要的结果。 另外，希望楼主能注意以下几个方面： 1 t1和t0都是双精度，但是t却没有定义啊！！！按照默认的I-N规则，t是单精度型，可能会引起误差，从而引起不必要的麻烦哦。 2 open之后没有对应的close语句，这点很不好！！很不好！！很不好！！楼主目前写的是小程序不觉得，程序一大，很容易搞混的，到时候想写这个文件却写到那个文件去了，或者想写却发现文件没打开之类的。而且这种错误会很难调试。 3 stop语句一般只用在这样两个方面：调试程序 或 出错后的强行停止。不建议用stop来结束程序，让程序自己运行完毕后停止岂不更好？ 4 前面有program xxx的话，最后建议用end program xxx结束。这样构造一个封闭的block，今后也会很容易读代码，便于维护管理。 5 楼主的程序使用了默认的I-N规则，建议使用Implicit None，取消默认规则，使每个变量都要定义类型（和C一样强类型）。使用默认规则的话，一旦有什么问题都是非常隐蔽极难调试的。 最后，楼主如果有程序方面的问题，欢迎移步“程序语言”版块提问。 |
夜月丨噩梦
5楼2016-01-02 19:32:34
bennett_14
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 7204.8
- 红花: 4
- 沙发: 1
- 帖子: 932
- 在线: 421.5小时
- 虫号: 1073877
- 注册: 2010-08-13
- 性别: GG
- 专业: 计算机科学
2楼2016-01-02 16:09:59
bennett_14
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 7204.8
- 红花: 4
- 沙发: 1
- 帖子: 932
- 在线: 421.5小时
- 虫号: 1073877
- 注册: 2010-08-13
- 性别: GG
- 专业: 计算机科学
3楼2016-01-02 16:35:26
4楼2016-01-02 17:51:51
50