In our numerical simulations, a nondimensional form of our previous model (17) was used to reduce significantly the number of parameters and to prevent the occurrence of numerical errors in the computation.