Fortran中传递数组，将三维数组传递到二维数组

伟大无敌sl

Fortran中传递数组，将三维数组传递到二维数组

[Fortran] 纯文本查看 复制代码

program ex0819
implicit none
integer,parameter::dim1=2
integer,parameter::dim2=2
integer,parameter::dim3=2
integer::a(dim1,dim2,dim3)
a(:,:,1)=1
a(:,:,2)=2
call getarray1(a(:,:,1),dim1,dim2)
call getarray2(a(1,1,2),dim1)
stop
end program ex0819

subroutine getarray1(a,dim1,dim2)
implicit none
integer::dim1,dim2
integer::a(dim1,dim2)
write(*,*)a
return
end subroutine getarray1

subroutine getarray2(a,dim1)
implicit none
integer::dim1,dim2
integer::a(dim1,*)！最后一维可以不赋值数组大小
integer::i
write(*,*)(a(:,i),i=1,2)
return
end subroutine getarray2

第一个子程序
为什么出来的结果会是1，1，1，1，第二个子程序为什么出来的结果会是2，2，2，2？？？

（彭国伦Fortran95程序设计，p182页例子)

青衣巷

a  :  2x2x2 三维数组。
a(:,:,1)   :   是 2x2 的二维数组，是 a 的前半部分。内容全部等于1
a(:,:,2)   :   也是 2x2 的二维数组，是 a 的后半部分。内容全部等于2

call getarray1(a(:,:,1),dim1,dim2)  这种是推荐使用的方式。没什么好说的。
call getarray2(a(1,1,2),dim1)  这种是不建议的使用方式，容易出错。不易于理解。

a 在内存里是这样存储的
a(1,1,1),a(2,1,1),a(1,2,1),a(2,2,1),a(1,1,2),a(2,1,2),a(1,2,2),a(2,2,2)
1       ,1       ,1       ,1       ,2       ,2       ,2       ,2

call getarray2(a(1,1,2),dim1) 把a(1,1,2) 传递给getarray2函数。
getarray2看到的部分是蓝色的部分：
a(1,1,1),a(2,1,1),a(1,2,1),a(2,2,1),a(1,1,2),a(2,1,2),a(1,2,2),a(2,2,2)
1       ,1       ,1       ,1       ,2       ,2       ,2       ,2      
因此对于getarray2来说，其值是4个2。

这种方式非常不好，不推荐。

参考：http://passarray.w.fcode.cn/

现代Fortran代码，传递数组全部用假定形状+module就好了。