[esp-r] Re: CFD problem

[Achim Geissler]

Hi Andy

it sounds like you have narrowed down the place where problems crop up quite a bit. Have you tried adding some write(*,*) lines e.g. for the Tf(x,y,z)? In addition to other variables that may be involved? This is tedious, but often helps get an idea what may be going wrong (and there should be no "scalar vs array" problem).

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