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