A numerical model of mesoscale frontal instabilities and plankton dynamics — I. Model formulation and initial experiments