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要怎么整？求指导！

回复此楼

» 猜你喜欢

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

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

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

supersunyan

金虫 (正式写手)

应助: 10 (幼儿园)
金币: 1585.9
散金: 1604
红花: 9
帖子: 695
在线: 296小时
虫号: 1145382
注册: 2010-11-13
性别: GG
专业: 流体力学

【答案】应助回帖

感谢参与，应助指数 +1

不知你做什么用，我们计算时会用到多个网格块，每个块都是不定大小的，用时循环定义就是了。如：
module data
  type blocktype
   integer : :n
   real,allocatable:: x(

end blocktype

type(blocktype),alloacatable::blk(

end module data

program main
use data
implicit none
integer :: i

allocate(blk(100))
do i=1,100
blk(i)%n=i
allocate(blk(i)%x(blk(i)%n))
enddo

end program

[ 发自小木虫客户端 ]