id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
95	wvelocity	jcwarner	arango	"wvelocity.F, near the bottom we have:

          wvel(i,j,0)=cff1*(vert(i,j,1)-
     &                      (vert(i,j,2)-vert(i,j,1))*
     &                      2.0_r8*(z_r(i,j,1)-z_w(i,j,0)))+
     &                cff2*vert(i,j,1)-
     &                cff3*vert(i,j,2)

written out this is:

wvel(i,j,0)=cff1*
(vert(i,j,1)-(vert(i,j,2)-vert(i,j,1))*2.0_r8*(z_r(i,j,1)-z_w(i,j,0)))+
cff2*vert(i,j,1)-cff3*vert(i,j,2)

cff1, cff2, and cff3 are coefficients for the interpolation (just numbers).
vert has dimensions of m/s.  so this is:

( m/s - ( m/s - m/s)* 2 * (m))   +   m/s  +   m/s
This does not look dimensionally correct.
Same issue for wvel (N).

"	bug	new	major	Release ROMS/TOMS 3.1	Nonlinear	3.1			WigginsW@…