call exec_softreg
else if (modecalc.eq.12) then
call exec_MD
+ call exec_checkgrad
else if (modecalc.eq.14) then
call exec_MREMD
else
subroutine exec_MD
use MPI_data !include 'COMMON.SETUP'
use control_data !include 'COMMON.CONTROL'
- use geometry, only:chainbuild
+ use geometry, only:chainbuild,chainbuild_cart
use MDyn
use io_units !include 'COMMON.IOUNITS'
+ use compare, only:alloc_compare_arrays
! use io_common
implicit none
! include 'DIMENSIONS'
#endif
print *,'Start MD'
call alloc_MD_arrays
+ call alloc_compare_arrays
print *,'After MD alloc'
if (me.eq.king .or. .not. out1file) &
write (iout,*) "Calling chainbuild"
+ if (extconf) then
call chainbuild
+ else
+ call chainbuild_cart
+ endif
call MD
return
end subroutine exec_MD
use REMD_data !include 'COMMON.REMD'
use geometry, only:chainbuild
use MREMDyn
+ use compare, only:alloc_compare_arrays
implicit none
! include 'DIMENSIONS'
integer :: i
call alloc_MD_arrays
call alloc_MREMD_arrays
-
+ call alloc_compare_arrays
! if (me.eq.king .or. .not. out1file) &
! write (iout,*) "Calling chainbuild"
! call chainbuild