50 |
|
|
51 |
|
implicit none |
52 |
|
PRIVATE |
53 |
+ |
#define __FORTRAN90 |
54 |
+ |
#include "UseTheForce/DarkSide/fInteractionMap.h" |
55 |
|
|
56 |
|
logical, save :: force_globals_initialized = .false. |
57 |
|
|
65 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: A_Row |
66 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: A_Col |
67 |
|
|
68 |
< |
real( kind = dp ), allocatable, dimension(:), public :: pot_Row |
69 |
< |
real( kind = dp ), allocatable, dimension(:), public :: pot_Col |
70 |
< |
real( kind = dp ), allocatable, dimension(:), public :: pot_Temp |
68 |
> |
real( kind = dp ), allocatable, dimension(:,:), public :: pot_Row |
69 |
> |
real( kind = dp ), allocatable, dimension(:,:), public :: pot_Col |
70 |
> |
real( kind = dp ), allocatable, dimension(:,:), public :: pot_Temp |
71 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: f_Row |
72 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: f_Col |
73 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: f_Temp |
162 |
|
return |
163 |
|
endif |
164 |
|
|
165 |
< |
allocate(pot_row(nAtomsInRow),stat=alloc_stat) |
165 |
> |
allocate(pot_row(LR_POT_TYPES,nAtomsInRow),stat=alloc_stat) |
166 |
|
if (alloc_stat /= 0 ) then |
167 |
|
thisStat = -1 |
168 |
|
return |
169 |
|
endif |
170 |
|
|
171 |
< |
allocate(pot_Col(nAtomsInCol),stat=alloc_stat) |
171 |
> |
allocate(pot_Col(LR_POT_TYPES,nAtomsInCol),stat=alloc_stat) |
172 |
|
if (alloc_stat /= 0 ) then |
173 |
|
thisStat = -1 |
174 |
|
return |
175 |
|
endif |
176 |
|
|
177 |
< |
allocate(pot_Temp(nlocal),stat=alloc_stat) |
177 |
> |
allocate(pot_Temp(LR_POT_TYPES,nlocal),stat=alloc_stat) |
178 |
|
if (alloc_stat /= 0 ) then |
179 |
|
thisStat = -1 |
180 |
|
return |