ready to do serious programming for Nh = 2 case, (phase3 correctly reading in Moperator2)