| 1 |
mmeineke |
270 |
subroutine forceFactory(forceName,status,forceSetup) |
| 2 |
|
|
! use TraPPE_ExFF |
| 3 |
|
|
use lj_ff, ONLY: new_lj_atype, init_ljFF,do_lj_ff |
| 4 |
|
|
character(len = *) :: forceName |
| 5 |
|
|
integer, intent(out) :: status |
| 6 |
|
|
external forceSetup |
| 7 |
|
|
|
| 8 |
|
|
status = 0 |
| 9 |
|
|
|
| 10 |
|
|
select case (forceName) |
| 11 |
|
|
|
| 12 |
|
|
! case ("TraPPE_Ex") |
| 13 |
|
|
! call forceSetup(new_TraPPE_ExFF,doTraPPE_ExFF) |
| 14 |
|
|
case("LJ") |
| 15 |
|
|
call forceSetup(new_lj_atype,init_ljFF,do_lj_ff) |
| 16 |
|
|
case default |
| 17 |
|
|
status = -1 |
| 18 |
|
|
return |
| 19 |
|
|
end select |
| 20 |
|
|
|
| 21 |
|
|
|
| 22 |
|
|
|
| 23 |
|
|
|
| 24 |
|
|
|
| 25 |
|
|
|
| 26 |
|
|
end subroutine forceFactory |
