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