  We present an extension of the previously proposed mean-field renormalization
method to model Hamiltonians which are characterized by more than just one type
of interaction. The method rests on scaling assumptions about the magnetization
of different sublattices of the given lattice and it generates as many flow
equations as coupling constants without arbitrary truncations on the
renormalized Hamiltonian. We obtain good results for the test case of Ising
systems with an additional second-neighbor coupling in two and three
dimensions. An application of the method is also done to a morphological model
of interacting surfaces introduced recenlty by Likos, Mecke and Wagner [J.
Chem. Phys. {\bf{102}}, 9350 (1995)].
  PACS: 64.60.Ak, 64.60.Fr, 05.70.Jk
