24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hanye0517

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1700.7
散金: 98
红花: 17
帖子: 144
在线: 116.8小时
虫号: 1253805
注册: 2011-04-02
性别: GG
专业: 核物理

[求助] fortran 如何在子程序中实现动态数组的定义和声明？已有4人参与

我编了一个动态数组的小例子，代码如下

CODE:

program allocatable

implicit none

real:: inteval=0.01,TI=0.03,Bc=0.24,Bdot=4.

real,dimension(:),allocatable :: time

real,dimension(:),allocatable :: B

real,dimension(:,:),allocatable:: TB

integer:: n=1,m=2

integer::i

print *,'The size of the shuzu?'

read *,n

allocate(time(n))

allocate(B(n))

allocate(TB(m,n))

do i=1,n

   time(i)=TI+inteval*i

   B(i)=Bdot*(time(i)-TI)+Bc

enddo

   TB(1,:)=time

   TB(2,:)=B

   print 100,TB

   100 format(f8.5,5x,f8.5)

end program allocatable

但是我觉得声明数组和分配内存的语句在主程序里太碍眼了，很烦躁。
想把他放到subroutine或者module或者function中，让结构更清晰。
但是子程序中的话，能编译，但是显示段错误。
module中干脆不能编译。
function说是只能返回一个数值，不知到怎么能返回一个数组，还是一个动态数组。

其实我想实现的功能就是避开复杂的fortran声明，像matlab那样用一个函数什么的声明一个动态数组。比如上边程序中，如果要是能用time=func（n），生成一个n个元素的一维数组，TB=func（m，n）生成一个2位数组就最好了。不知到这样的程序function要怎么整？求指导！

回复此楼

» 猜你喜欢

球磨粉体时遇到了大的问题，请指教！已经有11人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复

道之所在，虽千万人吾往矣。

1楼 2014-04-21 09:47:12

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

hanye0517

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1700.7
散金: 98
红花: 17
帖子: 144
在线: 116.8小时
虫号: 1253805
注册: 2011-04-02
性别: GG
专业: 核物理

引用回帖:

10楼: Originally posted by jerkwin at 2014-04-21 23:49:01
改改这个代码试试
http://stackoverflow.com/questions/3828094/function-returning-an-array-in-fortran

这个做不到，虽然是函数返回一个数组，但是赋值时仍然要声明一个数组才能让函数赋值。相当与没有用。