Hello I have some code that
 

Hello
I have some code that worked in OpenFoam-1.4.1 , but i can't seem to get it working on rhoCentralFoam.

I am running a very simple test
1)Geometry is just a rectangular channel.
2)If i specify all variables at inlet(as in supersonic inflow), and extrapolate all values(as in supersonic outlet), everything is fine. As expected there is no change in any flow parameters
3)If i change the bc's to subsonic inflow and outflow, all hell breaks lose.
I specify p,T on inlet. U is computed at inlet using the keyword pressureInletVelocity,
At outlet i specify p and extrapolate T,U using the zeroGradient keyword.

However even in such a simple case i get a floating point exception from sqrt

Could some one tell me what i am doing wrong?
Or is this a bug?

I have been looking at this for the last 3 days, with no results

Thanks
Srinath