24小时热门版块排行榜

>论坛更新日志 (4579)
>文献求助 (470)
>虫友互识 (376)
>导师招生 (311)
>考博 (185)
>招聘信息布告栏 (164)
>硕博家园 (157)
>休闲灌水 (115)
>论文道贺祈福 (114)
>博后之家 (100)
>考研 (88)
>论文投稿 (78)
>基金申请 (67)
>教师之家 (67)
>公派出国 (62)
>找工作 (55)

返回列表

lijunjiexm

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 22.8
红花: 1
帖子: 182
在线: 87.1小时
虫号: 1198843
注册: 2011-02-03
专业: 常微分方程与动力系统

[求助] 求助简单的一个参数传递问题，谢谢

请教一个简单程序MID,
   PROGRAM MID

implicit real*8(a-h,o-z)
   real*8 ia(5)

open(10,file='inid.dat')
open(21,file='mesh1.dat')
   read(21,*) numnp,numel

do 300 i=1,numnp
      read(21,2100) idum,xreal,yreal,ia(i)
      write(10,*) ia(i)
c       write(10,2100)  idum,xreal,yreal,ia(i)
300 continue

2100  format(i5,2x,e15.5,2x,e12.5,2x,i2)

   end
C---------------------------------------------------------------------------------------
mesh1.dat为
      5       2 numnp , numel
1    0.00000E+00 0.00000E+00 1
2    0.00000E+00 0.35803E-03 1
3    0.00000E+00 0.70174E-03 1
4    0.00000E+00 0.10317E-02 1
5    0.00000E+00 0.13485E-02 1
但是用 write(10,*)  ia(i)语句后，写进去的inid.dat文件为：
  4.940656458412465E-324
  4.940656458412465E-324
  4.940656458412465E-324
  4.940656458412465E-324
  4.940656458412465E-324
而选择write(10,2100)  idum,xreal,yreal,ia(i) ，则inid.dat为
1    0.00000E+00 0.00000E+00 1
2    0.00000E+00 0.35803E-03 1
3    0.00000E+00 0.70174E-03 1
4    0.00000E+00 0.10317E-02 1
5    0.00000E+00 0.13485E-02 1
但在在debug过程ia的数据居然是 4.940656458412465E-324，这个怎么解释呢？怎么使得ia数组等于 1 1 1 1 1？请多多指教，谢谢

回复此楼

» 猜你喜欢

真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
论文投稿，期刊推荐已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2011-09-01 21:42:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

【答案】应助回帖

★ ★
xzhdty(金币+2): 欢迎常来程序语言看看 2011-09-02 09:20:12

把 real*8 ia(5) 改成 integer ia(5) 就可以了（gfortran 验证通过）。如果不修改，在 gfortran 下会报如下错误：

Fortran runtime error: Expected INTEGER for item 5 in formatted transfer, got RE
AL
(i5,2x,e15.5,2x,e12.5,2x,i2)
^

赞一下(1人)

回复此楼

2楼2011-09-01 22:07:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lijunjiexm

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 22.8
红花: 1
帖子: 182
在线: 87.1小时
虫号: 1198843
注册: 2011-02-03
专业: 常微分方程与动力系统

引用回帖:

2楼: Originally posted by snoopyzhao at 2011-09-01 22:07:11:
把 real*8 ia(5) 改成 integer ia(5) 就可以了（gfortran 验证通过）。如果不修改，在 gfortran 下会报如下错误：

Fortran runtime error: Expected INTEGER for item 5 in formatted transfer, got RE
AL
( ...

为什么real不行呢？试了下把real*8改成dimension好像也可以，这个能稍微解释下吗，谢谢

赞一下

回复此楼

3楼2011-09-01 22:19:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖