logical function IsMomAllowedDet(nJ) integer :: nJ(nEl) IsMomAllowedDet = IsMomAllowedDetAnyParent(nJ, HFSym%Sym) end function IsMomAllowedDet