Ocean Modeling Discussion

Looks like one of your input has NaN values somewhere. ROMS doesn't allow this even on masked cells. See here: viewtopic.php?f=14&t=4454&view=next

Thanks-- I'll look in to it.

Hmm, it should be u (i-axis) west/east and v (j-axis) north/south. The angle on the grid is close to 0, about 1.5 (the domain is in the mid-latitudes).

The grid is set up just on a mercator curvilinear grid. (It's setup with some matlab code that originally was from ROMS.)

The lon's and lat's ...

hclarke wrote: with no figure. How do I get a figure? Sorry about the basic questions!

Just do pcolor(x_rho, y_rho, h). Matlab has very good documentation.

I've come across something that has stumped me with regards to how ROMS calculates u_eastward and v_northward from the u/v velocities. I took out the data from one of my history files, and I saw this (see below).

My question is, it seems as if U is interpolated in to v_northward, and V interpolated ...

Yes, this is correct, if your parent data array is arranged so that (1,1) is in the southwest corner. Here is a look at how I clip my arrays in matlab:

2-D
zeta_east = S.zeta(end,:);
zeta_west = S.zeta(1, : );
zeta_north = S.zeta(:,end);
zeta_south = S.zeta(:, 1);

3-D
u_east = S.u(end,:,:);
u ...

Yes, indeed, ubar is defined at u-points, vbar at v-points, just as the 3d momentum are at u and v points. This is verified by the script c_initial.m, which defines the dimensions for each of those.