put the for loop over ih OUTSIDE MoperatorGeneralFill, created dummy Matrix for getting ownership ranges to set rows efficiency, also no longer need the "if i < Nxyzcr  lines that were preventing the spurious nonzeros in Moperator when constructing multimode J