| 1 |
|
module reaction_field |
| 2 |
< |
use simulation |
| 2 |
> |
use force_globals |
| 3 |
|
use definitions |
| 4 |
|
use atype_module |
| 5 |
|
use vector_class |
| 17 |
|
|
| 18 |
|
contains |
| 19 |
|
|
| 20 |
< |
subroutine initialize_rf() |
| 21 |
< |
rrf = getRrf() |
| 22 |
< |
rt = getRt() |
| 23 |
< |
dielect = getDielect() |
| 20 |
> |
subroutine initialize_rf(this_rrf, this_rt, this_dielect) |
| 21 |
> |
real(kind=dp), intent(in) :: this_rrf, this_rt, this_dielect |
| 22 |
> |
|
| 23 |
> |
rrf = this_rrf |
| 24 |
> |
rt = this_rt |
| 25 |
> |
dielect = this_dielect |
| 26 |
|
|
| 27 |
|
rrfsq = rrf * rrf |
| 28 |
|
pre = 14.38362d0*2.0d0*(dielect-1.0d0)/((2.0d0*dielect+1.0d0)*rrfsq*rrf) |
