obtaining nodal (local) vector solution

[suren]

Hi.

I've coded Stokes equation (2-D, Newtonian fluid).. its working well and yielding results (as can be seen in the ElmerPost) same as the ones from the standard NS solver (FlowSolver)... no problem until here...

Now I am trying to include the material nonlinearity, which demands the local vector (velocity) solution... I've used CALL GetVectorLocalSolution (Velocity)

It gives correct velocity in second axis (vertical z-axis, max_velocity_w ~120 m/a), but gives some funny numbers (e.g. 6874589.03 m/a) in the first axis (x-axis or direction of flow... max_velocity_u should be around 270 m/a )... Couldnt understand what went wrong...

The solver section in .SIF looks like

Solver 1
...
Variable = [Velocity:2 Pressure:1]
...
FYI: CALL GetScalarLocalSolution (Pressure) is also giving the right values...

hope to hearing soon,
many thanks
~suren.

[suren]

I've got it, now... but with the sort of coding:
DO i=1, dim
Velocity(i,1:n) = FlowSolution((dim+1)*(FlowPerm(NodeIndexes(1:n))-1) + i)
END DO

... getting the nodal solution values directly, rather than using default subroutine GetVectorLocalSolution(Velocity)...

~suren