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

登录注册

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

24小时热门版块排行榜

返回列表

霹雳旋风

荣誉版主 (著名写手)

应助: 3 (幼儿园)
贵宾: 0.95
金币: 6428.7
散金: 100
红花: 4
帖子: 1124
在线: 175.7小时
虫号: 62568
注册: 2005-03-17
性别: GG
专业: 爆炸与冲击动力学

[求助] 求高人帮忙看一下我这个Fortran程序的错误出在哪里？

程序代码如下：

program ganeiliwucha
            implicit none
External FS
Integer I,L,N
Dimension X(2),Y(2),C(2)
Double Precision X, Y, EPS,C
open (2000, file='test.dat', status='old' )
Data X/0.00,0.01/
   EPS=1.0E-06
CALL DSNSE(2,EPS, X, Y, FS, L)
IF(L.GT.0)then
   C(1)=cos(X(2))/sin(X(1)+X(2))*tan(20.0/180*3.14)-1.0 ! AC杆的内力误差
C(2)=cos(X(1))/sin(X(1)+X(2))*sin(20.0/180*3.14)-1.0 ! BC杆的内力误差
do I=1,2
write(2000,*) X(I)
end do
end IF
close(2000)
end

SUBROUTINE DSNSE(N,EPS,X,Y,FS,L)
DIMENSION X(N),Y(N)
DOUBLE PRECISION X,Y,F,D,S
L=500
5 CALL FS(X,N,F,Y)
IF (F.GE.EPS) THEN
   L=L-1
   IF (L.EQ.0) RETURN
   D=0.0
   DO 20 J=1,N
20    D=D+Y(J)*Y(J)
   IF (D+1.0.EQ.1.0) THEN
      L=-1
      RETURN
   END IF
   S=F/D
   DO 30 I=1,N
30    X(I)=X(I)-S*Y(I)
   GOTO 5
END IF
RETURN
END

SUBROUTINE FS(X, N, F, Y)
DIMENSION X(N), Y(N)
Double precision X,Y,F,F1,F2,DF1,DF2,A,B
   A=20.0/180*3.14
   K=0.5
   F1=X(1)+X(2)-cos(A)
   F2=sin(X(1))*cos(A)*(1+K*cotan(A))+SIN(X(2))*(1.0-K/sin(A))-cos(A)
   F=F1*F1+F2*F2
   DF1=1.0
   DF2=cos(X(1))*cos(A)*(1+K*cotan(A))
   Y(1)=2.0*(F1*DF1+F2*DF2)
   DF1=1.0
   DF2=COS(X(2))*(1.0-K/sin(A))
   Y(2)=2.0*(F1*DF1+F2*DF2)
RETURN
END

  程序编译通过，但是出现计算数据无法输出的问题，请高手帮忙看一下。不胜感激。

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有7人回复
申请2026年博士已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

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

FORTRAN属性声明出错已经有8人回复
有会看fortran程序的，能帮忙看看这个程序吗？已经有10人回复
求助计算机高人，编写个Fortran程序已经有13人回复
请教一个fortran小程序编译出错的问题，谢谢已经有9人回复
写了一个fortran90的小程序，编译通不过，请大侠帮忙已经有59人回复
【求助】读文本数据出现段错误[fortran] 已经有14人回复
【求助】如何写FORTRAN程序实现求平均最近邻距离已经有34人回复
【求助】运行fortran程序出现错误 : error LNK2019 怎么回事啊已经有8人回复
【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙已经有8人回复
【求助】fortran程序未出错，但运行不出结果，中途中断！已经有6人回复
【求助】此FORTRAN90程序没有错误，但是运行不出来。已经有4人回复
【求助】用fortran求解大型线性方程组时出现的错误【已解决】已经有11人回复
【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙已经有11人回复

1楼2012-02-18 21:24:39

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

回帖置顶 ( 共有2个 )

霹雳旋风

荣誉版主 (著名写手)

应助: 3 (幼儿园)
贵宾: 0.95
金币: 6428.7
散金: 100
红花: 4
帖子: 1124
在线: 175.7小时
虫号: 62568
注册: 2005-03-17
性别: GG
专业: 爆炸与冲击动力学

引用回帖:

楼: Originally posted by ygl820 at 2012-02-21 10:31:22:
运行后，test.dat生成了没有。

运行后，生成了test.dat

赞一下

回复此楼

8楼2012-02-22 13:05:55

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

霹雳旋风

荣誉版主 (著名写手)

应助: 3 (幼儿园)
贵宾: 0.95
金币: 6428.7
散金: 100
红花: 4
帖子: 1124
在线: 175.7小时
虫号: 62568
注册: 2005-03-17
性别: GG
专业: 爆炸与冲击动力学

引用回帖:

楼: Originally posted by snoopyzhao at 2012-02-23 13:54:41:
这个程序用 gfortran 没有办法编译，显示 cotan 未定义……

你这个 test.dat 是空的很正常啊，说明 L .LE. 0 啊，呵呵……

谢谢楼上，
程序中的cotan（）的含义是求角度的余切值呀。请高人指点哦。

赞一下

回复此楼

10楼2012-02-23 19:37:37

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

普通回帖

夕阳西下

金虫 (著名写手)

应助: 138 (高中生)
贵宾: 0.041
金币: 9247.1
散金: 73
红花: 28
沙发: 1
帖子: 1268
在线: 1164.4小时
虫号: 220750
注册: 2006-03-20
专业: 凝聚态物性 II ：电子结构
管辖: 物理

【答案】应助回帖

感谢参与，应助指数 +1
jjdg(金币+1): 感谢应助 2012-02-20 01:39:45

open (2000, file='test.dat', status='old' )
write(2000,*) X(I)
write 2000这样写是不是有问题呢？

赞一下

回复此楼

2楼2012-02-19 09:17:19

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

haizhfly

木虫 (正式写手)

应助: 16 (小学生)
金币: 3134.8
帖子: 316
在线: 67.4小时
虫号: 925614
注册: 2009-12-11
专业: 天气学

【答案】应助回帖

感谢参与，应助指数 +1
jjdg(金币+1): 感谢应助 2012-02-20 01:39:55

open (2000, file='test.dat', status='old' )
write(2000,*) X(I)

status='old' 要求输出文本 test.dat是已经存在的文件，删去这一文件状态控制语句试一下看，个人觉得最大问题在这里

赞一下

回复此楼

3楼2012-02-19 11:21:31

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

霹雳旋风

荣誉版主 (著名写手)

应助: 3 (幼儿园)
贵宾: 0.95
金币: 6428.7
散金: 100
红花: 4
帖子: 1124
在线: 175.7小时
虫号: 62568
注册: 2005-03-17
性别: GG
专业: 爆炸与冲击动力学

引用回帖:

楼: Originally posted by haizhfly at 2012-02-19 11:21:31:
open (2000, file='test.dat', status='old' )
write(2000,*) X(I)

status='old' 要求输出文本 test.dat是已经存在的文件，删去这一文件状态控制语句试一下看，个人觉得最大问题在这里

这个我已经调试过。文件能编译完，也能运行，就是没有输出结果。

赞一下

回复此楼

4楼2012-02-19 19:48:38

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