[Nek5000-users] Probability Density Function

[nek5000-users at lists.mcs.anl.gov]

Hi Paul,

thanks. It might be that the reason for the negative weights is that zm1(i,j,k+1,e) exceeds the bounds

>> do e=1,nelv
>> do k=1,nz1
>> do j=1,ny1
>> do i=1,nx1
>> wght=(zm1(i,j,k+1,e)-zm1(i,j,k,e)) * area(i,j,f,e)

With a loop k=1,nz1-1 things seem to work fine. I will compare this with bm1(i,j,k,e).

Joerg