Based on the now well-known and experimentally observed critical gradient length (LTe) in tokamaks, we present a continuum one-dimensional model for explaining self organized heat transport in tokamaks. The key parameters of this model include a hysteresis parameter, which ensures that the switching of the heat transport coefficient χ upwards and downwards takes place at two different values of R/LTe. Extensive numerical simulations of this model reproduce many features of present day tokamaks, such as sub-marginal temperature profiles, intermittent transport events, 1/f scaling of the frequency spectra, propagating fronts, etc. This model utilizes a minimal set of phenomenological parameters, which may be determined from experiments and/or simulations. Analytical and physical understanding of the observed features has also been attempted.