Neumann boundary condition

["Si hang" <sh9988@xxxxxxxxxxxxx>]

Dear David and All,

I have encountered a mixed boundary condition (see below) along the 
outer boundary of domain in a 2D time harmonic simulation.

Ht = Jm * exp(-jPhi * x)

Where Jm is the line density of Current, Phi is the phase angle 
and it is a function of x-coordinate. 

This condition can't be set directly to the boundary with the Mixed 
boundary condition in FEMM. The only way I have thought is to split 
the boundary into many small segments. Then manually calculate the 
values on each segment's midpoint, set each segment's condition with 
its midpoint's value. 

Apparently, it's not a good way and cause a lot of work. Is there a 
better way or we can add the feature that support Mixed boundary 
condition with a function of coordinates(like the Prescribed A 
condition)?

Best regards

Si hang