Hello,
First I apologize if this is a somewhat trivial question but still getting used to ROMS. I am using it for a river application, mainly looking at velocities in a river meander. I have a couple cases for which I have made orthogonal grids (with gridgen) and have obtained reasonable results. My problem is that I wanted to look at the effects of inertia in the system by comparing runs with without the momentum advection terms. I assumed this could simply be done with #undef UV_ADV but I cannot get a converging solution of the water surface or velocity when I do this.
My case is simple and I only want to include the very basics (no tracers, just trying to essentially solve the nonlinear shallow water equations):
#define CURVGRID
#undef UV_ADV
#define UV_PSOURCE
#define UV_QDRAG
#define UV_VIS2
#define NORTHERN_WALL
#define SOUTHERN_WALL
#define WEST_FSCLAMPED
#define WEST_M2RADIATION
#define EAST_FSRADIATION
#define EAST_M2RADIATION
#define ANA_INITIAL
#define ANA_SMFLUX
#define ANA_PSOURCE
#define ANA_FSOBC
I have tried some other combinations, time steps, etc. to no avail. I had expected it to converge easier in this linear mode than nonlinear but this isn't the case. I've also tried this on a constant depth river and got the same results.
Any guidance is very much appreciated.