[esp-r] Re: CFD problem

[Achim Geissler]

Hi Andy

another thought: Do you have FORCHECK available? Looking at the code shows that in some places the variables are explicitly typed, in others not (same Commons). This may be a problem, but not really that likely.

If you have crawled through the code so far, already, looking for possible places where a "divide by zero" could happen would be my next thought. I don't recall what the original difference between your working model and the not working one was. With your new insights, does this difference maybe point to a specific place in the code?

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