[esp-r] CFD problem

[Andrew Cowie]

Hi all,

I've managed to trace my CFD problem a good way back, but I've reached an impass and I'd be gratefull if anyone had any suggestions.

The problem as described in my previous post occurs because the variable RESORW in CALCW is becoming NaN.  I've traced the NaN's back to BNDARTT in the file bndaryt.F; they first occur in the variable Uqstor(2,3,7).  The root of this is because Tq is calculated as 0 in the loop (8722 i think) where Uq is calculated.  The problem is the array Tf, but for some bizarre reason gdb doesn't recognize it as an array, it thinks it is a scalar of value 15.5.

So this is my impass; i can't debug this in the code because it wont compile if i test Tf as a scalar, and i can't debug it in gdb becuase it won't recognize it as an array and hence i can't inspect the specific value Tf(2,3,7).

I have tried compiling with extra debug options, and nothing is being flagged until the NaN's start cropping up, which obviously isn't helpful.

Any ideas/comments/input-of-any-kind would be gratefully recieved, i have to present my work in a week and it's not gonna look good if i have to make "it doesnt work" last 20 mins.

Thanks in anticipation,
Andy Cowie