The solvation force of a simple fluid confined between four identical walls is studied with different settings of surface fields. Our numerical results, obtained by densitymatrix renormalization techniques for a strip-like systems of widths up to 500 lattice constants, provide excellent the bulk free-energy extrapolation. The influence of surface fields setting on the pseudo-coexistence line and on the solvation force is presented.