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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » Fortran 77语言的数组定义问题
41/1返回列表
查看: 1481  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

bobo-scut

新虫 (初入文坛)

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

我编写了一个abaqus（有限元）子程序，使用Fortran 77语言调用一个数组（180*180），因为Fortran 77 单个数组最多只能有512行，一行6个数，所以需要平均分为12段赋值给新数组输出，为什么输出的数组混乱了？第10、11、12段数组跑到前面来了，导致输出的数组顺序颠倒。。。求大神指教
回复此楼

» 猜你喜欢
1楼 2017-08-02 19:01:55
已阅   回复此楼   关注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 ...

你的问题我也不知道是怎么回事。
给你一个小建议：把数据保存到单独的数据文件中；然后在程序中编写代码，打开数据文件，一行一行把数据读到数组变量中。
赞一下(1人) 回复此楼
4楼2017-08-02 21:33:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bobo-scut 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-06 04:49 by qL1cU46YIdK0
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +7 ncTVBM6AOJHF 2026-11-04 19/950 2026-11-06 04:38 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 PYZmcHr5kSaq 2026-11-05 7/350 2026-11-06 02:59 by qL1cU46YIdK0
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 4/200 2026-11-06 02:37 by qL1cU46YIdK0
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 02:24 by qL1cU46YIdK0
[基金申请] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 13/650 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 15/750 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:31 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 01:20 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 13/650 2026-11-05 22:55 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-05 22:44 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
[基金申请] 各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
关闭
请填处理意见
关闭