#1 .../src/futils.F90

[Anonymous]

I think the calculation of the physical vertical velocity is wrong

around line 149:
this is the old version:
:
indx=i+(j-jmin)*(imax-imin+1)
:
do k=0,kmin(i,j)-2
indx = indx+l
end do
:
do k=kmax,kmin(i,j),-1
ws(indx)=0.5*(ws(indx)+ws(indx-l))
indx = indx-l
end do
I'm not sure that indx-1 is always the point just below indx!

Presently in the newest version the calculation of the destaggered vertical velocity is commented out.

Karsten