| 查看: 561 | 回复: 5 | |||
[交流]
【求助】子程序传值出错
|
|||
|
我在主程序中调用了一个子程序test，然后再子程序中又调用了一个子程序test1.在子程序中的write语句都能得到ihbnum=1，可是主程序write语句的ihbnum去等于0，请大家帮忙看看，先谢了 program main -- call test（a1,a2,----ihbnum1） write (*,*) ihbnum1 enddo 子程序为： subroutine test & (io,ih,ia,xxx,yyy,zzz,nattot,cell,atomname,resname,ihbnum) implicit none include 'constants.h' integer io,ih,ia,nattot,ihbnum double precision xxx(nattot),yyy(nattot),zzz(nattot) double precision cell(*) character*4 atomname(nattot) character*3 resname(nattot) integer iat ia = 0 do iat=1, nattot if((resname(iat).eq.'DC5'.or.resname(iat).eq.'DC ').and.& (atomname(iat).eq.' N1 '.or.atomname(iat).eq.' N4 '.or.& atomname(iat).eq.' N3 '.or.atomname(iat).eq.' O2 ')then call test1 & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) if(ihbnum.eq.1) then ia=iat write(21,*) 'test',ihbnum,ia endif enddo return end subroutine test1 & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) ------- ihbnum=0 call bondcalculation(io,ih,xxx,yyy,zzz,nattot,cell,dx0,dy0,dz0,r0) call bondcalculation(io,iat,xxx,yyy,zzz,nattot,cell,dx1,dy1,dz1,r1) if (r1.lt.roamax) then theta = dacos(dx0*dx1+dy0*dy1+dz0*dz1) if(theta.lt.thetamax*pi/180.0) then call bondcalculation & (ih,iat,xxx,yyy,zzz,nattot,cell,dx2,dy2,dz2,r2) if (r2.lt.rhamax) then ihbnum = 1 write(21,*) 'find',ihbnum,iat endif endif endif return end |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
【求助】急求！！！
已经有5人回复
- 【求助】子程序的参数传递是传址还是传数
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
求助紫外可见分光光度计
+1/99
-
太原科技大学材料科学与工程学院招收金属材料类博士一名
+1/82
-
长江师范学院 材料与化工专业 金磊老师课题组招收调剂学生
+1/44
-
江苏科技大学能源材料化学课题组郭兴梅教授招收博士研究生1-2名
+1/40
-
烟台大学精准材料高等研究院26年材料方向研究生招生
+1/39
-
沈阳农业大学生物与医药专硕招收调剂生！！！！！！！！！！
+1/37
-
招收2026级环境科学方向（专业代码083001）硕士研究生1-2名
+1/37
-
武汉工程大学资源生物技术课题组生物/环境/化工/矿物专业硕士、博士研究生招生及调剂
+1/17
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/16
-
南京林业大学-国家级青年人才团队 招2026级博士、硕士（合成化学、植物化学方向）
+1/12
-
天津理工大学国家杰青王铁团队招聘学术骨干、博士后
+1/12
-
南京医科大学-谭磊课题组-诚招博士后
+1/10
-
诚挚招生（硕士+博士）
+1/9
-
三峡大学 计算力学课题组招收调剂
+1/7
-
广东医科大学 李雪萌教授团队 招收申请考核制博士（临床检验诊断）
+1/7
-
大湾区大学张效勇课题组诚招理论和计算化学方向2026级博士研究生和博士后
+1/7
-
半导体光电子器件与集成
+1/4
-
Top-88悉尼科技大学数据科学/AI 招收2027年入学 校奖 博士生1到2名(国际和本地学生）
+1/2
-
湘潭大学教师招聘
+1/1
-
北京高校副校长团队招收机械类，环境类学硕和专硕
+1/1
2楼2010-11-23 06:06:52
|
再问一下，对于子程序，这两个write语句输出的ihbnum的值是不是write（20）ihbnum=1，而write（21）的ihbnum=0？ subroutine test1 & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) ------- ihbnum=0 call bondcalculation(io,ih,xxx,yyy,zzz,nattot,cell,dx0,dy0,dz0,r0) call bondcalculation(io,iat,xxx,yyy,zzz,nattot,cell,dx1,dy1,dz1,r1) if (r1.lt.roamax) then theta = dacos(dx0*dx1+dy0*dy1+dz0*dz1) if(theta.lt.thetamax*pi/180.0) then call bondcalculation & (ih,iat,xxx,yyy,zzz,nattot,cell,dx2,dy2,dz2,r2) if (r2.lt.rhamax) then ihbnum = 1 write(20,*) 'find',ihbnum endif endif endif write(21,*) 'find',ihbnum return end |
3楼2010-11-23 06:32:40
4楼2010-11-23 06:42:52
6楼2010-11-23 19:20:11
简单回复
我非我亦我5楼
2010-11-23 08:02
回复
astringent(金币+1):谢谢参与
nono2009:专业区请勿灌水。谢谢！ 2010-11-25 07:35:16
astringent(金币+3): 2010-11-29 15:56:24
祝福祝福
nono2009:专业区请勿灌水。谢谢！ 2010-11-25 07:35:16
astringent(金币+3): 2010-11-29 15:56:24
