A theoretical analysis of gas-liquid counter-current flow in laminar boundary layers with flat phase boundary based on similarity variables method has been done. The obtained numerical results for the velocity distribution in both phases are compared with analogous results from asymptotic theory and experimental data. The dissipation energy in boundary layer is determined and the results corresponding to counter-current and co-current flows are compared. The comparison shows significant differences in dissipation energy values in gaseous phase.