Hello,
As you declare that there is no tidal flats in your mesh, it could be that somewhere in your initial conditions the depth is 0 or negative (case 1, however you specify a constant depth, so this could be only explained by a specific programming of CONDIM.F). If you do a specific spacing of planes in subroutine CONDIM it could also be that you need to add the keyword :
MESH TRANSFORMATION : 2 or 3
in your parameter file (case 2). In that case please send your fortran file so that I can check what you do. Another current problem is old versions of subroutine CALCOT remaining in some test cases, where truncation errors may cause an intercrossing of planes (at some time we had to change a formula to avoid this, though the two versions look identical in theory) (case 3, but needs tidal flats).
If you have activated an example given in subroutine condim, a recent mistake that occured to a user is also that the lines :
ZSTAR%R(1)=0.D0
ZSTAR%R(NPLAN)=1.D0
are missing (case 4).
So this is 4 possibilities to check (that can be resolved by seing your fortran file), and we'll see,
With best regards,
Jean-Michel Hervouet