i think i found a bug in the pointinelement routine in the interpolation.F90
THere the distance to quadriterals (case 4) is calculated by
Code: Select all
CASE(4)
sumdist = MAX( ug - 1.0, MAX( -ug -1.0, 0.0 ) )
sumdist = sumdist + MAX( vg - 1.0, MAX( -vg - 1.0, 0.0 ) )
Code: Select all
CASE(4)
sumdist = MAX( ug - 1.0, MAX( -ug , 0.0 ) )
sumdist = sumdist + MAX( vg - 1.0, MAX( -vg , 0.0 ) )
then the term
max(-(-0.5)-1,0)=max(0.5-1,0)=max(-0.5,0)=0
does not give a positive distance as i think it is supposed to.
Also the other therm
max(-0.5-1,0)=max(-1.5,0)=0
does not deliver that. But i think that one should cover the case where u is bigger than 1, and so it is fine.
I think the analogous is true below for the hexagonal case (case 8(?))
best regard
Franz