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)

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).