| 9 |
|
!! @author Matthew Meineke |
| 10 |
|
!! @author Christopher Fennel |
| 11 |
|
!! @author J. Daniel Gezelter |
| 12 |
< |
!! @version $Id: calc_sticky_pair.F90,v 1.11 2003-07-17 20:32:24 gezelter Exp $, $Date: 2003-07-17 20:32:24 $, $Name: not supported by cvs2svn $, $Revision: 1.11 $ |
| 12 |
> |
!! @version $Id: calc_sticky_pair.F90,v 1.12 2003-07-17 20:38:11 gezelter Exp $, $Date: 2003-07-17 20:38:11 $, $Name: not supported by cvs2svn $, $Revision: 1.12 $ |
| 13 |
|
|
| 14 |
|
module sticky_pair |
| 15 |
|
|
| 32 |
|
real( kind = dp ), save :: SSD_ru = 0.0_dp |
| 33 |
|
real( kind = dp ), save :: SSD_rlp = 0.0_dp |
| 34 |
|
real( kind = dp ), save :: SSD_rup = 0.0_dp |
| 35 |
+ |
real( kind = dp ), save :: SSD_rbig = 0.0_dp |
| 36 |
|
|
| 37 |
|
public :: check_sticky_FF |
| 38 |
|
public :: set_sticky_params |
| 63 |
|
SSD_ru = sticky_ru |
| 64 |
|
SSD_rlp = sticky_rlp |
| 65 |
|
SSD_rup = sticky_rup |
| 66 |
+ |
|
| 67 |
+ |
if (SSD_ru .gt. SSD_rup) then |
| 68 |
+ |
SSD_rbig = SSD_ru |
| 69 |
+ |
else |
| 70 |
+ |
SSD_rbig = SSD_rup |
| 71 |
+ |
endif |
| 72 |
|
|
| 73 |
|
sticky_initialized = .true. |
| 74 |
|
return |
| 118 |
|
return |
| 119 |
|
endif |
| 120 |
|
|
| 121 |
< |
if ( rij .LE. SSD_rup ) then |
| 121 |
> |
if ( rij .LE. SSD_rbig ) then |
| 122 |
|
|
| 123 |
|
r3 = r2*rij |
| 124 |
|
r5 = r3*r2 |
