We investigate numerically the global pressure-drop of fluids with temperature dependent viscosity, flowing through a porous medium channel bounded by two parallel isoflux surfaces. By reviewing the development of the Hazen-Dupuit-Darcy (HDD) equation we bring to light the inappropriateness of the model in estimating the global pressure-drop of fluids with temperature dependent viscosity. Albeit this observation, we tested the accuracy of the HDD model in comparison with numerical results by using three alternatives, namely (1) fluid viscosity determined at the average bulk temperature, (2) fluid viscosity determined at the log-mean bulk temperature and (3) fluid viscosity replaced by a channel-length averaged fluid viscosity. The HDD model is inadequate because the temperature dependent fluid viscosity surprisingly affects both, viscous and form, global drag terms. We propose and validate a new global model, which accounts for the effects of temperature dependent viscosity in both drag terms of the original HDD model. Based on our new model, two regimes are discovered as the surface heat flux increases. In the first regime both drag terms are affected, while in the second regime only the form drag term is affected, prior to the model reaching an inviscid limit. Predictive empirical relations correcting the viscous and form drag terms, complementing the new model, are obtained as functions of the surface heat flux.