Fortranでのポインタ配列操作 - 下限エラー

Question

3D配列ポインタを4D配列ポインタに格納したいが、何か問題がある。例えば：functionはp(:,:,:)などのポインタを返す関数へのポインタであり、これは正常に動作

real,pointer :: p(:,:,:,:) => null() 
integer :: ndims,d 

ndims = 3 
do d=1,ndims 
    p(d,:,:,:) => function 
end do 

。私が得ているエラーは、この行のLower bound has to be presentです：p(d,:,:,:) => functionしかし、私はこの種の操作を（ポインタではなく）実行していて、それは問題なく動作しました。私はここで何が欠けていますか？

Answer 1

私は、これは実際にはFortranでやっていること

real,pointer :: p(:,:,:,:) => null() 

の一般的な誤解から来ていると思います。

ポインタの4D配列ではなく、4D配列へのポインタ 'p'です。あなたは何ができるか

は、ちょっとしたトリックを再生し、内部のポインタで構造体を作成することです：

type pointer4D 
    real, pointer :: p 
end type pointer4D 
type(pointer4D), dimension(:) :: arr 

あなたのコードは、このようになる：「関数」を仮定し

do d=1,ndims 
    arr(d)%p => function 
end do 

は、ある種のです3Dオブジェクトを使用すると、次元Nの要素I、j、kにアクセスできます。

arr(N)%p(I,j,k)