A crystallization model appropriate for application in continuum modeling of complex processes is presented. As an extension to the previously developed Schneider equations [W. Schneider, A. Köppel, and J. Berger, "Non-isothermal crystallization of polymers," Int. Polym. Proc. 2, 151 (1988)], the model presented here allows one to account for the growth of crystals of various shapes and to distinguish between one-, two-, and three-dimensional growth, e.g., between rod-like, plate-like, and sphere-like growth. It is explained how a priori knowledge of the shape and growth processes is to be built into the model in a compact form and how experimental data can be used in conjunction with the dynamic model to determine its growth parameters. The model is capable of treating transient processing conditions and permits their straightforward implementation. By using thermodynamic methods, the intimate relation between the crystal shape and the driving forces for phase change is highlighted. All these capabilities and the versatility of the method are made possible by the consistent use of four structural variables to describe the crystal shape and number density, irrespective of the growth dimensionality.