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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 请教一个read的用法，谢谢
91/1返回列表
查看: 2447  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

li_clifff

银虫 (正式写手)

[求助] 请教一个read的用法，谢谢

请教一个问题，我要搞清楚程序中这两个语句的结果：

INTEGER:: k, nr,i, nodof=3,nn=96
INTEGER,ALLOCATABLE:: nf(:, : )
ALLOCATE nf(nodof,nn)
READ(10,*)nr,(k,nf(:,k),i=1,nr)

其中文件10为：
81
1 1 0 0   2 1 0 0   3 1 1 0   4 1 0 0   5 1 1 0   6 1 0 0
7 1 1 0   8 1 0 0   9 1 1 0  10 1 0 0  11 1 1 0  12 0 0 0
13 1 0 1  14 1 0 1  15 1 0 1  16 1 0 1  17 0 0 0  18 0 1 0
19 1 0 1  21 1 0 1  23 1 0 1  25 1 0 1  27 1 0 1  28 0 1 0
29 0 0 0  30 1 0 1  31 1 0 1  32 1 0 1  33 1 0 1  34 0 0 0
35 0 1 0  36 1 0 1  38 1 0 1  40 1 0 1  42 1 0 1  44 1 0 1
45 0 1 0  46 0 0 0  47 1 0 1  48 1 0 1  49 1 0 1  50 1 0 1
51 0 0 0  52 0 1 0  53 1 0 1  54 1 1 1  55 1 0 1  57 1 0 1
59 1 0 1  61 1 0 1  62 0 1 0  63 0 0 0  64 1 0 1  65 1 0 1
66 1 0 1  67 1 0 1  68 0 0 0  69 0 1 0  70 1 0 1  72 1 0 1
74 1 0 1  76 1 0 1  78 1 0 1  79 0 1 0  80 0 0 0  81 1 0 1
82 1 0 1  83 1 0 1  84 1 0 1  85 0 0 0  86 0 1 0  87 0 0 0
88 0 1 0  89 0 0 0  90 0 1 0  91 0 0 0  92 0 1 0  93 0 0 0
94 0 1 0  95 0 0 0  96 0 1 0
在debug模式下，nr值为81，k为96（？），nf是3*96的数组，i等于82(?)

能帮我解释下为什么会这样吗？非常感谢
其中nf数组前面几个值为：




[ Last edited by li_clifff on 2011-10-15 at 00:57 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-14 22:50:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)
★ ★
余泽成(金币+2): 谢谢参与应助，欢迎常来程序语言版！ 2011-10-15 19:13:10
ben_ladeng: 专家考核存档 2011-10-17 11:51:35
读入第一个数据， 为nr (number of rows?)
然后 4个一组地读入 到k, nf(1,k), nf(2,k), nf(3,k); 一共读 nr次
赞一下(2人) 回复此楼
更痛苦的是梦没醒路已经不能走了
2楼2011-10-15 03:08:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
余泽成(金币+2): 谢谢参与应助！ 2011-10-15 19:13:21
li_clifff(金币+10): 恩，明白了，谢谢 2011-10-15 23:22:40
ben_ladeng: 专家考核存档 2011-10-17 11:51:44
ALLOCATE nf(nodof,nn) 这句话是错误的吧？
ALLOCATE(nf(nodof,nn))
如果在这之后加一句：
nf = 100000
然后在程序结束之前加上：
CODE:
do i=1,96
   write(*,*) nf(1:3,i)
end do

你就知道你的程序是在做什么了……
赞一下(2人) 回复此楼
3楼2011-10-15 14:30:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
引用回帖:
3楼: Originally posted by snoopyzhao at 2011-10-15 14:30:40:
ALLOCATE nf(nodof,nn) 这句话是错误的吧？
ALLOCATE(nf(nodof,nn))
如果在这之后加一句：
nf = 100000
然后在程序结束之前加上：
CODE:
do i=1,96
   write(*,*) nf(1:3,i)
end do

你就知道你 ...

恩，再请教一个小问题，如果在主程序里定义了变量数组nf是3*9的，即nf（3,9），而且主程序里的 call subroutine（nf，...）,对应的subroutine里面对nf定义为nf（2,9），会有什么问题吗？运行时候到底是nf（3,9），还是nf（2,9），刚想到有这个问题，就顺便问一下，谢谢
赞一下 回复此楼
4楼2011-10-15 23:26:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
引用回帖:
4楼: Originally posted by li_clifff at 2011-10-15 23:26:37:
恩，再请教一个小问题，如果在主程序里定义了变量数组nf是3*9的，即nf（3,9），而且主程序里的 call subroutine（nf，...）,对应的subroutine里面对nf定义为nf（2,9），会有什么问题吗？运行时候到底是nf（3,9） ...

你确定你这种做法可以编译通过吗？呵呵……

如果可以的话，给个示例看看……
赞一下 回复此楼
5楼2011-10-16 13:01:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
引用回帖:
5楼: Originally posted by snoopyzhao at 2011-10-16 13:01:30:
你确定你这种做法可以编译通过吗？呵呵……

如果可以的话，给个示例看看……

program  ex0804
implicit none
  integer :: a(2)
  DATA  a /1,2/
  integer :: b(2)
  data b  /2,3/
  call add(a,b)  ! 把变量a及b交给子程序add来处理
  stop
end


subroutine add(first, second)
implicit none
  integer :: first, second ! first,second的内容会从CALL时得到
  write(*,*) first+second
  return
end
赞一下 回复此楼
6楼2011-10-16 22:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
jjdg(金币+1): 感谢参与 2011-10-17 02:32:35
li_clifff(金币+10): 谢谢你的评论，如果数组参数在主程序和多个子程序都出现，是不是也是这样一步步传递，读入子程序，运行子程序出来新的数组，再读入下一个子程序。。。 2011-10-17 09:44:20
引用回帖:
6楼: Originally posted by li_clifff at 2011-10-16 22:51:30:
program  ex0804
implicit none
  integer :: a(2)
  DATA  a /1,2/
  integer :: b(2)
  data b  /2,3/
  call add(a,b)  ! 把变量a及b交给子程序add来处理
  stop
end


subroutine add(first, se ...

由于 Fortran 是 call by reference，所以，关键还是要看子程中相关变理的声明了……

在你的这个例子中，很显然，只取其中的一个数值。

对于你前面说的，其实，传递给子程的只是数组对应的第一个数值的地址，后面到底如何取值，还要看子程中的定义……这种情况下，一定要注意 Fortran 中是列优先的存贮方式……
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

li_clifff
7楼2011-10-17 00:11:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
送鲜花一朵
引用回帖:
7楼: Originally posted by snoopyzhao at 2011-10-17 00:11:07:
由于 Fortran 是 call by reference，所以，关键还是要看子程中相关变理的声明了……

在你的这个例子中，很显然，只取其中的一个数值。

对于你前面说的，其实，传递给子程的只是数组对应的第一个数值的地址 ...

program  ex0804
implicit none
integer :: a(3)
DATA  a /1,2,3/
integer :: b(3)
data b  /4,5,6/
call add(a,b)  ! 把变量a及b交给子程序add来处理
call add1(a,b)

  stop
end

subroutine add(first, second)
implicit none
  integer :: first(2), second(2) ! first,second的内容会从CALL时得到
  write(*,*) first+second
  return
end

subroutine add1(first, second)
implicit none
  integer :: first, second ! first,second的内容会从CALL时得到
  write(*,*) first+second
  return
end

为什么运行的结果是
     5      7
     5
赞一下 回复此楼
8楼2011-10-17 09:59:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
没有问题，我自己弄错了

program  ex0804
implicit none
!integer:: a=1
integer :: a(3)
DATA  a /1,2,3/
!integer:: b=2
integer :: b(3)
data b  /4,5,6/
call add(a,b)  ! 把变量a及b交给子程序add来处理
call add1(a,b)

stop
end


subroutine add(first, second)
implicit none
  integer :: first(2), second(2) ! first,second的内容会从CALL时得到
! write(*,*)  first+second
  first = first + second
  second =2*second
  return
end


subroutine add1(first, second)
implicit none
  integer :: first(2), second(2) ! first,second的内容会从CALL时得到
  write(*,*) first+second
  return
end

运行结果是

       13               17
赞一下 回复此楼
9楼2011-10-17 10:36:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 li_clifff 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 13/650 2026-08-17 12:33 by gltch
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
关闭
请填处理意见
关闭