In this paper steady-state synchronous machine model which incorporates magnetic saturation and crossmagnetization is presented. Instead of using synchronous reactances and induced voltage E0 this model uses current dependent flux linkage functions in direct and quadrature axes to model armature reaction. If those functions are modeled as a polynomials of three independent variables which are armature current in direct and quadrature axis and field current, identification of the proposed model's parameters from load operation data becomes a linear problem. Paper also provides a method for calculation arbitrary steady-state load operation point based on the knowledge of the current dependent flux linkage functions. Proposed model, identification and calculation methods are validated and tested on 35 MVA, 10.5 kV, 50 Hz synchronous generator.