[Hedge] initialisation, angles, interpolation

[Islam El Salawi]

Hi everybody,

I'd like to set a specific initial  velocity angle at the inlet, and in the 
rest of the field the angle should be 0. Right now I'm using the 
UniformMachflow class, though it sets the same angle for the whole field. I 
tried something like:


T106A = UniformMachFlow(mach=0.59,char_length=0.08592314273,reynolds=50000)

def bound_direction(x,el):
        if x[0]<=min(x):
            return (1,1/sqrt(1-(cos(180-beta_2))**2))
        else:
            return T106A.direction

fields=discr.interpolate_boundary_function(bound_direction,"inflow")

But then an error is raised:

ValueError: non-slice indexing not supported: 0

I followed the example on the scicomp inifileframework site. 
Maybe someone has an idea and can help me.

Regards,
Islam    
-- 
GMX DSL Doppel-Flat ab 19,99 Euro/mtl.! Jetzt mit 
gratis Handy-Flat! http://portal.gmx.net/de/go/dsl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.tiker.net/pipermail/hedge/attachments/20110405/a225b7de/attachment.html>