wvelocity

[jcwarner]

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)

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(i,j,N).

[arango]

Yes, the units for the bottom and surface boundary conditions do not match. The boundary conditions can be computed in different ways. This corrected by a linear extrapolation.