update 8/11

corrected error with coh sub elem ordering when nsub=4, now matrix crack
can propagate as expected for angle from 30 to 90

however, when angle smaller than 30, error occurs, dont know why