In other words, a methodology for numerical PDE modeling is presented that is flexible, open ended and readily implemented on modest computers. If the reader is interested in an alternate model, it might possibly be implemented by: (1) modifying and/or extending the current model (for example, by adding terms to the PDEs or adding additional PDEs), or (2) using the reported routines as a prototype for the model of interest.