[esp-r] Re: CFD problem

[Achim Geissler]

Hi Andy

another thing, Tf is actually a scalar in the common /climi/ - again, could the difference between the working and non-working model lead to additional modules in the "name space"? I'm afraid, though, that my Fortran knowledge (which common is defined where at run time) is limited and this could absolutely be a wild goose chase (FORCHECK should report on such possible name clashes, though) - possibly the compiler, also, however.

Best
Achim



On Sep 20, 2012, at 6:49 PM, Andrew Cowie wrote:

> 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
> 
> _______________________________________________
> esp-r mailing list
> esp-r at lists.strath.ac.uk
> http://lists.strath.ac.uk/mailman/listinfo/esp-r
> 

achim.geissler at intergga.ch