C READ fragment information!!
C both routines should be in dfa.F file!!
- if (.not. (wdfa_dist.eq.0.0 .and. wdfa_tor.eq.0.0
+ if (.not. (wdfa_dist.eq.0.0 .and. wdfa_tor.eq.0.0 .and.
& wdfa_nei.eq.0.0 .and. wdfa_beta.eq.0.0)) then
call init_dfa_vars
print*, 'init_dfa_vars finished!'