It's not quite clear what the three dimensions of your Array N are supposed to mean, so it's a little tricky to figure out what you want plotted here, and how.
There are quite a few problems in your posted Document. As written, the loop code would not run, generating an error about bad index into the Array.
In the plotting command you use `M` but that must surely be `N`.
Inside the inner loop there is `n0` which is not defined, so I suppose that `indexN0` was intended there. I changed this.
Somehow several instances of indexN0 have turned into text (?!) with upright font instead of the italics (usual for general names). I just changed this, though it might not have been necessary.
If I chose wrongly while getting it to run, apologies.
My best guess is that you are trying to make a plot for each "r value" of the second dimension of N, from 1 to nStepsR. And, judging from your intended label on the horizontal axis I figure you want the points of each such curve to be according to the value of the index of the first dimension of N (time, I guess).
Let's clear up something, though. There is no special `surfdata` command that handles multple surfaces. In your previous post the solution I gave wrapped `surfdata` inside a call to `seq`. The `seq` command generates a sequence of calls to its first argument (in that case, `surfdata`). So, you see, `surfdata` was in fact called multiple times.
One way to handle you current Question, if I have understood properly, is to again produce a sequence of individual 2D plots. There are two apects you want, I suspect. Different colors for each individual plot, and proper tickmarks for the horizontal axis. You used Maple 17 before, so I've done so too. I believe that you are trying to plot 1-dimension slices of N, such as N[..,r,1] where r has a fixed value. In other words, the horizontal axis represents the first dimension of Array N (time, I guess). It doesn't matter if I guess the wrong dimension of N. What matters more, I think, is that you have only a single fully explicit dimension of data, which is the value (vertical axis) at the Array entry. The horizontal aspect is just the position of the entries, the first index into N. That's why your label for the horizontal axis mentions scale. Because from one entry to the next the increment (1 to tMax) implicitly represents some unit of time. So you either have to note that unit in the label, or you could create a custon set of tickmarks for the horizontal axis.
Anyway, the command which does 2D plots of an indexed collection of values is plots:-listplot. Unfortunately there is no 2D analogue to `listplot` with a relationship similar to that of the 3D `surfdata` to `matrixplot`. (I'm still rambling about the scale or ticks long the first axis.)
Below, I use a nice feature of newer Maple versions (inclusing Maple 17) to set a palette of colors for the separate 2D plots. This avoids having to generate a list of nStepsR colors.
Another way to handle this whole business is the use the `plot` command. Suppose that we knew that each increment of the first index of N (1 to tMax) represented 3 seconds. In that case you could form a reusable 1D column Vector of time values, and conjoin it with each 1D column Vector of height values (eg, N[..,r,1] for fixed r, changed from 1D Array to column Vector). This way, there is no need to create custom tickmarks.
I hope that at least gives you some ideas for handling this.