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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】大家帮我看一下这个程序的错误
91/1返回列表
查看: 586  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

zyj8119

木虫 (著名写手)

[交流] 【求助】大家帮我看一下这个程序的错误 已有1人参与

program main
        parameter (natom0=18331,nN=5000)
        character atomname*2,nname(3)*3
        integer i,j
        real coor(natom0,3),xn,yn,zn,n(natom0,3)
        real dl
        open(10,file='1.pdb',status='old')
        if((i.le.natom0).and.(i.ge.1))then
        read(10,*)i,atomname,coor(i,1),coor(i,2),coor(i,3)
      end if
      close (10)
        do 25 i=1,nN
         do 26 j=1,3
         n(i,j)=0
26    continue
25    continue

   

      do 30 i=1,nN
50          xn=64.185*RAN2(INUM)
          yn=149.764*RAN2(INUM)
          zn=42.79*RAN2(INUM)
         do 40 j=1,natom0
         dl=sqrt((xn-coor(j,1))**2+(yn-coor(j,2))**2+(zn-coor(j,3))**2)
      if (dl.lt.0.15) goto 50
40        continue
            
        if(i.ge.1)then
         do 60 k=1,i
         dl=sqrt((xn-n(k,1))**2+(yn-n(k,2))**2+(zn-n(k,3))**2)
      if(dl.lt.0.15) goto 50
60    continue
      end if
      
             n(i,1)=xn
        n(i,2)=yn
        n(i,3)=zn

30        continue

        do 110 i=1,nN
         do 120 j=1,3
        nname(1)='N'
        nname(2)='N'
        nname(3)='N'
120   continue
110   continue
     
      open(60,file='2.pdb',access='append')
        do 100 j=1,3
      write(60,600)i,nname(j),n(i,j)
        write(60,600)i,atomname,coor(i,j)
100        continue
600   format(1X,I3,A2,3F8.3)

      close(60)
        end


        FUNCTION RAN2(INUM)
        INTEGER idum,IM1,IM2,IMM1,IA1,IA2,IQ1,IQ2,IR1,IR2,NTAB,NDIV
        REAL ran2,AM,eps,RNMX
        PARAMETER(IM1=214783563,IM2=214783399,AM=1./IM1,IMM1=IM1-1,
     *IA1=40014,IA2=40692,IQ1=53668,IQ2=52774,IR1=12211,IR2=3791,
     *NTAB=32,NDIV=1+IMM1/NTAB,EPS=1.2E-7,RNMX=1.-EPS)
        INTEGER idum2,j,k,iv(NTAB),iy
         SAVE iv,iy,idum2
         DATA idum2/123456789/,iv/NTAB*0/,iy/0/
        if(idum.le.0)then
        idum=max(-idum,1)
        idum2=idum
        do 11 j=NTAB+8,1,-1
         k=idum/iq1
        idum=IA1*(idum-k*IQ1)-k*IR1
        if(idum.lt.0)idum=idum+IM1
        if(j.le.NTAB)iv(j)=idum
11     continue
       iy=iv(1)
        end if
        k=idum/IQ1
        idum=IA1*(idum-k*IQ1)-k*IR1
        if(idum.lt.0)idum=idum+IM1
        k=idum2/IQ2
        idum2=IA2*(idum2-k*IQ2)-k*IR2
        if(idum2.lt.0)idum2=idum2+IM2
        j=1+iy/NDIV
        iy=iv(j)-idum2
        iv(j)=idum
        if(iy.lt.1)iy=iy+IMM1
        ran2=min(AM*iy,RNMX)
        return
        END
说的是数组越界，在iy=iv(j)-idum2这里，使用的是COMPAQ VISUAL FORTRAN的DEBUGGING 调试的，输入文件的格式是
1        SI        38.349        21.095        16.446
2        O        30.974        20.423        14.478
3        O        34.674        13.2        5.204
4        SI        33.557        9.684        4.419
5        O        32.251        18.364        8.928
6        O        33.171        11.24        4.572
7        O        32.016        12.941        0.731
8        SI        36.701        18.416        16.43
9        O        27.251        20.294        2.714
10        O        30.31        15.536        13.347
11        O        32.387        5.65        15.653
12        SI        37.202        20.618        21.333
13        O        28.179        16.802        6.719
14        O        30.904        19.77        10.618
15        O        30.236        13.134        20.145
回复此楼

» 猜你喜欢
好好学习，天天向上。
1楼 2010-03-27 21:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
能编译通过，但是运行不出来
赞一下 回复此楼
好好学习，天天向上。
2楼2010-03-27 21:21:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
请把错误一并贴上来，帮助分析。
赞一下(1人) 回复此楼
3楼2010-03-29 14:32:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanjj

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
请问RAN2(INUM)中参数INUM是干什么的？
赞一下(1人) 回复此楼
在下才疏学浅，请多多指教了。
4楼2010-04-01 22:25:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
引用回帖:
Originally posted by wanjj at 2010-04-01 22:25:46:
请问RAN2(INUM)中参数INUM是干什么的？

调用子程序的参数。。
赞一下 回复此楼
好好学习，天天向上。
5楼2010-04-02 00:36:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

figgers

铜虫 (初入文坛)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+1):谢谢参与应助！ 2010-04-03 15:19
wangen994(金币+2):感谢参与应助 2010-04-03 15:20
wangen994(金币+2):活动期间内额外奖励 2010-04-03 15:21
这么几行程序，就单步一步步看就是了，肯定是j超过32了或者小于1了
btw: 这个程序的风格太糟糕了. 三十年前的风格啊...
如果整个工程就这么些程序，完全可以用90格式写的清清楚楚，跟C一样的
赞一下(4人) 回复此楼
6楼2010-04-03 14:44:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
引用回帖:
Originally posted by figgers at 2010-04-03 14:44:10:
这么几行程序，就单步一步步看就是了，肯定是j超过32了或者小于1了
btw: 这个程序的风格太糟糕了. 三十年前的风格啊...
如果整个工程就这么些程序，完全可以用90格式写的清清楚楚，跟C一样的

兄弟，我是新手，以后还请多多指教，我的QQ是675096471.
赞一下 回复此楼
好好学习，天天向上。
7楼2010-04-03 23:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
jjdg:引用帖子或直接PM可以更快获得答复！ 2010-04-03 23:47
j超过32了或者小于1了
是什么意思？？？
赞一下(1人) 回复此楼
好好学习，天天向上。
8楼2010-04-03 23:29:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
引用回帖:
Originally posted by figgers at 2010-04-03 14:44:10:
这么几行程序，就单步一步步看就是了，肯定是j超过32了或者小于1了
btw: 这个程序的风格太糟糕了. 三十年前的风格啊...
如果整个工程就这么些程序，完全可以用90格式写的清清楚楚，跟C一样的

肯定是j超过32了或者小于1了
?什么意思？
赞一下 回复此楼
好好学习，天天向上。
9楼2010-04-04 09:18:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zyj8119 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-16 22:25 by lindoudou
[考博] 现在不知道怎么办，感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[基金申请] 重磅！青年科学基金项目（C类）资助增幅预计超过50% +6 水和泥不是水泥 2026-05-13 8/400 2026-05-16 19:50 by Equinoxhua
[基金申请] 精华III:  评审感受-评审感受-评审感受 +15 ferrarichen 2026-05-11 19/950 2026-05-16 19:44 by Equinoxhua
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 19:10 by j6b2pdz07o
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-16 18:45 by j6b2pdz07o
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 3/150 2026-05-16 17:30 by j6b2pdz07o
[有机交流] 求有机合成大神指点三硫酸乙烯酯（CAS：2793408-99-6）的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 16:35 by x28q7dxf75
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 考博自荐 +3 科研狗111 2026-05-13 4/200 2026-05-16 11:45 by 科研狗111
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 11:36 by h3oerqvkv9
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 cjf4bx70cj 2026-05-14 6/300 2026-05-16 11:16 by h3oerqvkv9
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 04:29 by k37jurhrau
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-16 01:49 by k37jurhrau
[基金申请] 青C资助名额大幅增加！ +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
关闭
请填处理意见
关闭