degug of MC method
check for units problems
remove Uij method (too long)
add D+ matrix to extra