In this investigation through of Monte Carlo simulations we characterize thermomagnetically two Ising type ferromagnetic models of spins S = ±5/2, ±3/2, ±1/2 and σ = ±7/2, ±5/2, ±3/2, ±1/2 located alternately in a square lattice, formed by two sublattices A and B each with 120 sites containing the spins S and σ, respectively. The ferromagnetic configurations studied are defined by and H1 and H2 Hamiltonians with next neighbors and next nearest neighbors exchange interactions, and single ion anisotropy and external magnetic field interactions. The thermomagnetic analysis of the system defined by H1 Hamiltonian shows discontinuities in the magnetization, as well as a strong influence of magnetic field on the critical and discontinuous transition temperatures. For the case of the ferromagnet defined by H2 Hamiltonian, after analyzing the effects of the next nearest neighbors interactions and magnetic field, we find an almost linear dependence on the critical temperature. The study of hysteresis behavior in both ferromagnets shows the phenomenon of exchange bias.