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