In this paper we have studied the problem of an unsteady hydromagnetic forced convective heat transfer slip flow over a porous rotating disk taking into account the temperature dependent density, viscosity and thermal conductivity. The governing non-linear partial differential equations of the flow are transformed into a set of non-linear ordinary differential equations using similarity transformations. The resulting nondimensional equations have been solved numerically by applying Nachtsheim-Swigert shooting iteration technique along with sixth-order Runge-Kutta iteration scheme. Comparison with previously published work for the steady case of the problem is performed and the results are found to be in very good agreement. The results of the numerical solution are presented graphically in the form of velocity (i.e. radial, tangential as well as inward axial), temperature profiles and variable Prandtl number for various values of the model parameters. The corresponding skin friction coefficients (i.e. radial and tangential) and the rate of heat transfer coefficient (i.e. Nusselt number) are also calculated and tabulated. The obtained numerical results show that when modeling a thermal boundary layer with temperature dependent fluid properties, consideration of Prandtl number as constant within the boundary layer produces unrealistic results. Therefore, it must be treated as variable throughout the boundary layer. Results also show that the slip factor significantly controls the flow and heat transfer characteristics.