14 |
|
#ifdef IS_MPI |
15 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: q_Row |
16 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: q_Col |
17 |
< |
real( kind = dp ), allocatable, dimension(:,:), public :: qcom_Row |
18 |
< |
real( kind = dp ), allocatable, dimension(:,:), public :: qcom_Col |
17 |
> |
real( kind = dp ), allocatable, dimension(:,:), public :: q_group_Row |
18 |
> |
real( kind = dp ), allocatable, dimension(:,:), public :: q_group_Col |
19 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: u_l_Row |
20 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: u_l_Col |
21 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: A_Row |
51 |
|
subroutine InitializeForceGlobals(nlocal, thisStat) |
52 |
|
integer, intent(out) :: thisStat |
53 |
|
integer :: nrow, ncol |
54 |
+ |
integer :: nrow_group, ncol_group |
55 |
|
integer :: nlocal |
56 |
|
integer :: ndim = 3 |
57 |
|
integer :: alloc_stat |
61 |
|
#ifdef IS_MPI |
62 |
|
nrow = getNrow(plan_row) |
63 |
|
ncol = getNcol(plan_col) |
64 |
+ |
nrow_group = getNrowGroup(plan_row) |
65 |
+ |
ncol_group = getNcolGroup(plan_col) |
66 |
+ |
|
67 |
|
#endif |
68 |
|
|
69 |
|
call FreeForceGlobals() |
82 |
|
return |
83 |
|
endif |
84 |
|
|
85 |
< |
allocate(qcom_Row(ndim,nrow),stat=alloc_stat) |
85 |
> |
allocate(q_group_Row(ndim,nrow_group),stat=alloc_stat) |
86 |
|
if (alloc_stat /= 0 ) then |
87 |
|
thisStat = -1 |
88 |
|
return |
89 |
|
endif |
90 |
|
|
91 |
< |
allocate(qcom_Col(ndim,ncol),stat=alloc_stat) |
91 |
> |
allocate(q_group_Col(ndim,ncol_group),stat=alloc_stat) |
92 |
|
if (alloc_stat /= 0 ) then |
93 |
|
thisStat = -1 |
94 |
|
return |
254 |
|
if (allocated(A_Row)) deallocate(A_Row) |
255 |
|
if (allocated(u_l_Col)) deallocate(u_l_Col) |
256 |
|
if (allocated(u_l_Row)) deallocate(u_l_Row) |
257 |
< |
if (allocated(qcom_Col)) deallocate(qcom_Col) |
258 |
< |
if (allocated(qcom_Row)) deallocate(qcom_Row) |
257 |
> |
if (allocated(q_group_Col)) deallocate(q_group_Col) |
258 |
> |
if (allocated(q_group_Row)) deallocate(q_group_Row) |
259 |
|
if (allocated(q_Col)) deallocate(q_Col) |
260 |
|
if (allocated(q_Row)) deallocate(q_Row) |
261 |
|
#else |