A finite difference scheme for solving the problem of natural transport of heat, mass, momentum and species concentration along vertical porous plates is presented. Several drying related problems are numerically solved, by including a gas-injection boundary condition directly into the governing equations. The effect of variable physical properties is investigated by means of direct comparison against experimental data obtained through holographic interferometry. The relative importance of wall diffusive and convective fluxes is examined. Sherwood and Nusselt numbers can be accurately obtained by means of the proposed techniques.