linux fortran 大数组,fortran动态数组大小分配问题（allocatable）？ - 程序语言 - 小木虫 - 学术 科研 互动社区...

浔荆

我也常遇到这样的问题，一般我是通过以下方案解决的。

定义一个很大的数组，如sta_info(0:1000)。

之所以定义成sta_info(0:1000)而非sta_info(1:1000)，并不是因为我习惯了C语言里的数组默认下标从0开始。实际上，这里的sta_info(0)和其他元素不一样，它是用来记录整个数组的有效大小的。

例如，我们在第一个循环给元素0:55赋值，其他的元素用不到，可以这么实现：

allocate(sta_info(0:1000))

sta_info=-10000                   !初始化为一个不常见的数，方便调试

sta_info(0)=55                     !确定你需要的数组大小

do i=1,sta_info(0)

sta_info(i)=.....

enddo

简单来说，就是先分配足够的大小，但另外定义一个变量来描述它的实际大小。这样一来避免了数组大小变化的时候，需要频繁的allocate和deallocate。另一方面又可以避免数组太大而带来的速率降低。(排序这样的数组操作，只需要对1：sta_info(0)这些元素进行即可)

pippi6

subroutine DynamicArray

implicit none

real*8, allocatable, dimension( :: sta_info

integer :: i,n=20,j

do i=1,n

allocate(sta_info(i))

do j=1,i

sta_info(j)=j

end do

print*,' i=',i,' sizeof(sta_info)=',sizeof(sta_info)

deallocate(sta_info)

end do

end subroutine DynamicArray

qingkonggzc

引用回帖:

浔荆 at 2015-05-27 08:19:38

我也常遇到这样的问题，一般我是通过以下方案解决的。

定义一个很大的数组，如sta_info(0:1000)。

之所以定义成sta_info(0:1000)而非sta_info(1:1000)，并不是因为我习惯了C语言里的数组默认下标从0开始。实际上， ...

你好，谢谢。

首先这样和直接定义一个较大的数组，然后依次把数据/字符串写进去效果应该是一样的吧，就是有多少数据写多少空间的，剩下的还是空的，然后对其进行排序？这个直接把总共有多少数据计算出来，然后对这1-n个数据进行排序后面的没有影响吗？还是把这1-n个数据提取出来放到一个新的数组里面？

qingkonggzc

引用回帖:

pippi6 at 2015-05-27 13:14:14

subroutine DynamicArray

implicit none

real*8, allocatable, dimension( :: sta_info

integer :: i,n=20,j

do i=1,n

allocate(sta_info(i))

do j=1,i

sta_info(j)=j ...

你好，谢谢你，我的这个是只能在这一个循环结束之后才能读出来第二个循环该有多少个数据的，也不能在开始的时候就给他一个大小值啊？你这还是构造一个中间数组用来保存每一个循环的数据，再把这些数据汇总到一个数组里面去，类似于MATLAB里面的矩阵合并C=[a b]？Fortran里面可以这样合并吗？

浔荆

引用回帖:

qingkonggzc at 2015-05-27 19:23:58

你好，谢谢。

首先这样和直接定义一个较大的数组，然后依次把数据/字符串写进去效果应该是一样的吧，就是有多少数据写多少空间的，剩下的还是空的，然后对其进行排序？这个直接把总共有多少数据计算出来，然后对这 ...

对，是一样的。不一样的只是你可以针对非空元素进行操作而已。

假设有n个有效数据，剩余的全是空值。那么，你排序只需要对这n个数据排序就行了，不需要提取出来也不需要中间数组。直接操作部分数组就行。

你可以把排序子程序写成这样：

subroutine sort(array,n)

integer*4 array(0:1000)

integer*4 n              !只有1~n元素是有效的

......                        !对1~n号元素排序

end subroutine sort

然后调用时，直接call sort(sta_info,sta_info(0))就行了。

这个原理和字符串的原理是类似的，定义一个长度1000的字符串，其实际长度为1001，其中有一位是用于记录字符串有效长度的，有效长度之外全部是空格。

qingkonggzc

引用回帖:

浔荆 at 2015-05-28 08:13:16

对，是一样的。不一样的只是你可以针对非空元素进行操作而已。

假设有n个有效数据，剩余的全是空值。那么，你排序只需要对这n个数据排序就行了，不需要提取出来也不需要中间数组。直接操作部分数组就行。

你 ...

你好，我还是有点不明白只对1~n号元素进行排序的意思，虽然我想要的是这个结果，是对n号元素排序之后新放到一个数组里面还是在原来的数组里面，最后面的还是空？

allocate(sta_info(0:1000))

sta_info=-10000                   !初始化为一个不常见的数，方便调试

sta_info(0)=55                     !确定你需要的数组大小

do i=1,sta_info(0)

sta_info(i)=.....

enddo

sta_info=-10000是相当于给这数组一个unit嘛？不太明白这个初始化的意思。。

sta_info(0)=55                     !确定你需要的数组大小

do i=1,sta_info(0)

sta_info(i)=.....

enddo

这个当把第二个循环的大小确定为100的时候，怎么样能使它接着在sta_info后面写45个数据，这样应该会使它覆盖掉前一循环的数据吧？还是把大小也存放到一个数组里面，每次循环do从sta_num(i-1)~sta_num(i)

，

qingkonggzc

引用回帖:

浔荆 at 2015-05-27 08:19:38

我也常遇到这样的问题，一般我是通过以下方案解决的。

定义一个很大的数组，如sta_info(0:1000)。

之所以定义成sta_info(0:1000)而非sta_info(1:1000)，并不是因为我习惯了C语言里的数组默认下标从0开始。实际上， ...

你好，这是我程序里面的一部分，但是这样在第二个循环的时候就会触发断点，这样的话，存储大小会溢出？如果我直接把大小设置为1000，就不会触发断点

integer::sta_array_size=0

integer::day_sta_num=0

sta_array_size = sta_array_size + day_sol_sum%sta

allocate(sta_info(sta_array_size))

do j=1,day_sol_sum%sta,1

read(101,'(a)',end=200) sum_line

sta_info(day_sta_num+j)%sta_name=sum_line(2:5)

sta_info(day_sta_num+j)%phase_RMS=sum_line(10:13)

sta_info(day_sta_num+j)%weekday=doy-start_doy+1

end do

day_sta_num = day_sol_sum%sta + day_sta_num