44 |
|
|
45 |
|
|
46 |
|
module suttonchen |
47 |
+ |
use definitions |
48 |
|
use simulation |
49 |
|
use force_globals |
50 |
|
use status |
60 |
|
#define __FORTRAN90 |
61 |
|
#include "UseTheForce/DarkSide/fInteractionMap.h" |
62 |
|
|
62 |
– |
INTEGER, PARAMETER :: DP = selected_real_kind(15) |
63 |
|
!! number of points for the spline approximations |
64 |
|
INTEGER, PARAMETER :: np = 3000 |
65 |
|
|
254 |
|
n = 0.5_dp*(n1+n2) |
255 |
|
|
256 |
|
dr = (rCut) / dble(np-1) |
257 |
< |
rvals(1) = 0.0d0 |
258 |
< |
vvals(1) = 0.0d0 |
259 |
< |
phivals(1) = 0.0d0 |
257 |
> |
rvals(1) = 0.0_dp |
258 |
> |
vvals(1) = 0.0_dp |
259 |
> |
phivals(1) = 0.0_dp |
260 |
|
|
261 |
|
do k = 2, np |
262 |
|
r = dble(k-1)*dr |