code_saturne User's Forum

Thank you very much Yvan for your answer. It helped me updating the typecl.f90 file to version 5.3.4. The turbulence model is in version 5.1.6 and 5.3.4, but this particular modified non-user-defined subroutine, which is used for very specific purpose was not moved to version 5.3.4, which I just did...

Hi all, I tried to use SRC files of version 5.1.6 with the version 5.3.4 of Code Saturne but I am experiencing a SIGSEGV error. I isolated the issue. It is coming from my file typecl.f90 (see attached document) that is used to impose a pressure difference between 2 channels that are linked with a pi...

Hi Yvan,

I am using the 5.3.4 version. Is the syntax adapted to this version? If not, could you give me the appropriate one? There is no example in the MEG dialog.

Thanks again
Best regards,
Martin

Hi all, I would like to save outputs with the writer to make animation on Paraview. I need to save the data every 20 timesteps but only after the 10000th timestep. I write in the setup file: <frequency period="formula"> if ( niter &gt; 10000 ) { if (mod((niter - 10000), 20) == 0) { iac...

Hi Yvan and thanks for your answer. Practically, I should create a file named cs_restart.c, then write in it: /* Creation of the restart file */ cs_restart_create(const char * my_quantity1_file, const char * save_additional_quantities_directory, cs_restart_mode_t write) /* Declaration of the used su...

Hi all, I have created a new variable ("r") that is used in a home made model. I would like to store it at the end of a computation to use it to initialise a restart. For the already existing quantities there is no issue. For instance for the streamwise velocity, in the cs_user_initialisat...

Ok thanks!

Best regards

Indeed I forgot to close the loop. Thanks for your advices! I just check Code_Saturne results on Paraview and the the point data ID of the cells vertices. I noticed that two vertices of the same face are not sequential. So the point data ID numerotation in Paraview is different from the cell vertice...

Hi Yvan and thanks a lot for your answer. Here is the code for internal faces, may you just have a look to it and tell me if you see major issues or clumsy loops? ! FOR INTERNAL FACES allocate(Lf(nfac)) ! Vector containing the maximum length of each face allocate(DDmax(ncel)) ! Vector containing the...

Hi all, I develop a new turbulence model using Code_Saturne 5.3.4 (the moving to the lastest version will be performed in the near future). I use the fortran routines and I need to compute the longest edge of each cell of my mesh to use it in my turbulence model (the diagonales should not be taken i...