24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1408  |  回复: 3

bobo-scut

新虫 (初入文坛)

[求助] Fortran 77语言的数组定义问题

我编写了一个abaqus(有限元)子程序,使用Fortran 77语言调用一个数组(180*180),因为Fortran 77 单个数组最多只能有512行,一行6个数,所以需要平均分为12段赋值给新数组输出,为什么输出的数组混乱了?第10、11、12段数组跑到前面来了,导致输出的数组顺序颠倒。。。求大神指教
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cdut

铁虫 (正式写手)

把代码贴出来啊
2楼2017-08-02 19:32:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bobo-scut

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by cdut at 2017-08-02 19:32:30
把代码贴出来啊

SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
     1 COORDS,JLTYP,SNAME)

      INCLUDE 'ABA_PARAM.INC'

      DIMENSION TIME(2), COORDS (3)
      CHARACTER*80 SNAME     

      real:: a(32400)
      F=0
          data a(1:2700)/0,0,0,0,0,0,
     1 0,0.00220443,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0.00497663,0.00290721,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0.00480093,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0.00209624,0.00517915,
     1 0,0,0,0,0,0.00745109,
     1 0,0,0,0,0,0,
     1 0,0,0,0.00202873,0,0,
     1 0,0.00220443,0,0,0,0.00220443,
     1 0.0145517,0,0,0.00516617,0,0,
。。。。。。。。。。。中间数据太长,省略了(共12段数据)。。。。。。

         data a(29701:32400)/0.0168773,0.0358914,0,0,0,0,
     1 0,0,0,0,0,0.058191,
     1 0.285805,0.192845,0.178793,0.279138,0.817874,1.42453,
     1 0.568069,0.17191,0.0252345,0,0,0,
     1 0,0,0.0415708,0.0328483,0.0314557,0.006924,
     1 0.0146728,0.00747879,0.00206941,0.00923662,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0.0280474,0.0282101,
     1 0.134518,0.150935,0.192466,0.425515,0.438037,0.286305,
     1 0.196144,0.0212048,0.00206941,0.109296,0.143416,0.215023,
     1 0.130637,0.019311,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0,
     1 0,0,0,0,0,0/
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

      if (NOEL.ge.907201 .and. NOEL.le.939600)then
          if (a(NOEL-907200).gt.0.0000000000001)then
              F=a(NOEL-907200)*1E6   
          end if        
      end if
      RETURN
      END
3楼2017-08-02 19:59:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cdut

铁虫 (正式写手)

引用回帖:
3楼: Originally posted by bobo-scut at 2017-08-02 19:59:58
SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
     1 COORDS,JLTYP,SNAME)

      INCLUDE 'ABA_PARAM.INC'

      DIMENSION TIME(2), COORDS (3)
      CHARACTER*80 SNAME     

      re ...

你的问题我也不知道是怎么回事。
给你一个小建议:把数据保存到单独的数据文件中;然后在程序中编写代码,打开数据文件,一行一行把数据读到数组变量中。
4楼2017-08-02 21:33:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bobo-scut 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +5 1234567wang 2026-05-17 6/300 2026-05-17 21:18 by maxjing
[基金申请] 今年审到国自然15份,谈谈感受 +9 国自然国社科中 2026-05-17 9/450 2026-05-17 21:14 by 雪城001
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-17 19:39 by Equinoxhua
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 9/450 2026-05-17 08:28 by soxiaomucong
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:16 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 5/250 2026-05-17 01:10 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[考博] 现在不知道怎么办,感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯(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,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
信息提示
请填处理意见