Potantial bug with Code_Saturme v8 (atmospheric dispersion in case of imbrication)

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
Julien.Commanay
Posts: 18
Joined: Mon Nov 27, 2023 6:48 pm

Potantial bug with Code_Saturme v8 (atmospheric dispersion in case of imbrication)

Post by Julien.Commanay »

Dear Support team,

I think there is a bug when using atmospheric module of Code_Saturne v8 in case of imbrication with several meteo files.

It seems that Ux interpolated by Cressman on lateral faces is always 0, while values are correct for other fields (Uy, TKE, Eps, ...). I made sure that meteo profiles used for Cresman interpolation are read properly and then I put in place a workaround (maybe not the smartest way to fix the bug) in atimbr.fr to solve this issue. Feel free to ask in case you need additionnal information about this topic (my case is too big to be uploaded there and standard log is not really helpful (I had to add additionnal and very verbose printings to understand what happen).

Kind regards

JCo
Martin FERRAND
Posts: 48
Joined: Wed Mar 14, 2012 10:06 am

Re: Potantial bug with Code_Saturme v8 (atmospheric dispersion in case of imbrication)

Post by Martin FERRAND »

Hello,
Thank you for the report. We recently simplified the use of Cressman in boundary conditions. Could you share with us the way you fix atimbr.f90 ?
Do you have the problem only with multiple méteo profiles or already with one ?

Best regards
Julien.Commanay
Posts: 18
Joined: Mon Nov 27, 2023 6:48 pm

Re: Potantial bug with Code_Saturme v8 (atmospheric dispersion in case of imbrication)

Post by Julien.Commanay »

Hello Martin,

I'll send you the modified routine and set of meteo files by email.

I tested with one and several meteo file and the problem only occur with several meteo (linked to Cresman imbrication, you'll see where I adapted the code).

KR

JCo
Post Reply