1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
! { dg-do compile }
! PR fortran/36167
! This used to cause an ICE because of a missing array spec after interface
! mapping.
! Contributed by Frank Muldoon <fmuldoo@me.lsu.edu>
module communication_tools
contains
!*******************************************************************************
function overlap_1(u,lbound_u,ubound_u)
!*******************************************************************************
integer, intent(in), dimension(:) :: lbound_u,ubound_u
real, intent(in), dimension(lbound_u(1):ubound_u(1),lbound_u(2):ubound_u(2),&
lbound_u(3):ubound_u(3)) :: u
real, dimension(&
lbound(u,1):ubound(u,1),&
lbound(u,2):ubound(u,2),&
lbound(u,3):ubound(u,3)) :: overlap_1
return
end function overlap_1
end module communication_tools
!*******************************************************************************
subroutine write_out_particles
!*******************************************************************************
use communication_tools
real, dimension(1:5, 2:4, 3:10) :: vorticityMag
real, allocatable, dimension(:,:,:) :: temp3d
allocate(temp3d( &
lbound(overlap_1(vorticityMag,lbound(vorticityMag),ubound(vorticityMag)),1):&
ubound(overlap_1(vorticityMag,lbound(vorticityMag),ubound(vorticityMag)),1),&
lbound(overlap_1(vorticityMag,lbound(vorticityMag),ubound(vorticityMag)),2):&
ubound(overlap_1(vorticityMag,lbound(vorticityMag),ubound(vorticityMag)),2),&
lbound(overlap_1(vorticityMag,lbound(vorticityMag),ubound(vorticityMag)),3):&
ubound(overlap_1(vorticityMag,lbound(vorticityMag),ubound(vorticityMag)),3)))
return
end subroutine write_out_particles
|
