call flush(iout)
call molread(*10)
call flush(iout)
- if (constr_dist.gt.0) call read_dist_constr
+ if ((constr_dist.gt.0).and.(genconstr.eq.0)) call read_dist_constr
#ifdef MPI
write (iout,*) "Calling proc_groups"
call proc_groups
call promienie
call read_compar
call read_ref_structure(*10)
+ if ((constr_dist.gt.0).and.(genconstr.gt.0)) call gen_dist_constr
call proc_cont
call fragment_list
endif
+
C if (constr_dist.gt.0) call read_dist_constr
write (iout,*) "Begin read_database"
call flush(iout)