| 1 |
< |
!! |
| 1 |
> |
!! |
| 2 |
|
!! Copyright (c) 2005 The University of Notre Dame. All Rights Reserved. |
| 3 |
|
!! |
| 4 |
|
!! The University of Notre Dame grants you ("Licensee") a |
| 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 |
| 267 |
|
|
| 268 |
|
vCut = epsilon*((alpha/rCut)**n) |
| 269 |
|
|
| 270 |
< |
call newSpline(MixingMap(i,j)%V, rvals, vvals, & |
| 271 |
< |
0.0d0, 0.0d0, .true.) |
| 272 |
< |
call newSpline(MixingMap(i,j)%phi, rvals, phivals, & |
| 273 |
< |
0.0d0, 0.0d0, .true.) |
| 270 |
> |
call newSpline(MixingMap(i,j)%V, rvals, vvals, .true.) |
| 271 |
> |
call newSpline(MixingMap(i,j)%phi, rvals, phivals, .true.) |
| 272 |
|
|
| 273 |
|
MixingMap(i,j)%epsilon = epsilon |
| 274 |
|
MixingMap(i,j)%m = m |
