It seems that DDM method can greatly improve the throughput. Just to include an option as `DDM=Y` in the control file. Will implement once I have time.