A computational method is developed for predicting the steady-state temperature field in an LMR core. Detailed core-wide coolant temperature profiles are efficiently calculated using the simplified energy equation mixing model[1] and the subchannel analysis method. The $ heta$-method is employed for discretizing the energy equations in the axial direction. The interassembly coupling is achieved by interassembly gap flow. Cladding and fuel temperatures are calculated with the one-dimensional conduction model and temperature integrals of conductivities. The accuracy of the method is tested by performing several benchmark calculations for too LMR problems. The results indicate that the accuracy is comparable to the other methods based on ENERGY model. It is also shown that the implicit scheme for the axial discretization is more efficient than the explicit scheme.