24小时热门版块排行榜    

查看: 1434  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[基金申请] fileCode有新解读? +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见