user prostprocessing output

[frpageot]

Hello,

I tried to add a user output for postprocessing and then as an example tried to use one of the cases exposed in he usvpst.f90 file. I launched the calculation and an error orred whan it was supposed to write the data :
Signal SIGSEGV (accès à une zone mémoire interdite) intercepté !

Pile d'appels :
1: 0x7fe08459b76d <fvm_convert_array+0x344d> (libfvm.so.0)
2: 0x7fe08460fef3 <fvm_writer_field_helper_step_e+0xb23> (libfvm.so.0)
3: 0x7fe08461ac2d <fvm_to_med_export_field+0x151d> (libfvm.so.0)
4: 0x7fe08460d767 <fvm_writer_export_field+0xb7> (libfvm.so.0)
5: 0x7fe08497563d <cs_post_write_var+0x16d> (libsaturne.so.0)
6: 0x7fe08497871a <pstev1_+0x10a> (libsaturne.so.0)
7: 0x7fe08497b949 <psteva_+0x5b> (libsaturne.so.0)
8: 0x40b532 <usvpst_+0x16e> (cs_solver)
9: 0x7fe0849799db <pstvar_+0x10cb> (libsaturne.so.0)
10: 0x7fe08487dedc <caltri_+0x76ec> (libsaturne.so.0)
11: 0x7fe0848569e3 <cs_run+0x843> (libsaturne.so.0)
12: 0x7fe084856cc5 <main+0x1f5> (libsaturne.so.0)
13: 0x7fe080ff5c4d <__libc_start_main+0xfd> (libc.so.6)
14: 0x40b309 <> (cs_solver)
Fin de la pile

I just used the first example with the enclosed file. Can someone hlep me ? Thanks.

[Yvan Fournier]

Hello,

Unless you only have output on the volume mesh (i.e. the surface mesh is deactivated, either through the GUI or using ichrbo = 0), you should not set ipart = -1. The routine is called for all output meshes, with a different value of ipart.

Also, what turbulence model are you using ? ir11, ir22, and ir33 are defined correctly only when you are using a Reynolds stress model (Rij SSG or LRR) ?

Regards,

Yvan

[Fabien852]

Hello,

I wish to use the postprocessing output like in the GUI.

Could you please give me some step to do it ? I read the user's guide but i don't understand.

I just want a postprocessing at each 10 iterations on a MED file.

I use the version 2.0.2.


Best regards,

Fabien

[Yvan Fournier]

Hello,

If you do not need anything more advanced than what is available using the GUI, why do you need user subroutines (unless you want simultaneous MED and EnSight output) ?

If you only need an additional mesh or sub-mesh or output format with the same variables, you only need to use usdpst (check for pstcat in the examples), not usvpst.

You did not answer my previous question (what turbulence model are you using), so I cannot be sure what caused your crash.

Regards,

Yvan

[Fabien852]

Hello Yvan,

I'm not the one who made this topic, this is why I don't answer .

I use user subroutines because I use the combustion module and I want to see the flame evolution. So I will check usdpst file.


Best regards,


Fabien

[Fabien852]

Hello,

Can you help me to find the line where i can set the "Post-processing every "n" time steps" in the file usdpst.f90 like in the GUI ?

I just want to see the evolution of temperature and velocity in a burner.

Best regards,


Fabien

[Jacques Fontaine]

Hello,
The output frequency for post-pro is fixed with ntchr = n_iter in usini1.f90.

Best regards,
JF

[frpageot]

Hello,

I am sorry I did not answer to Yvan, my apologizes. I found what was the problem, but they were several. Actually if I do no rememeber well how I solved the problem : iit seems that I was saving my user routines in a folder named with the same name CS is using for temporay files, and it seems that there was a conflict. Is it possible ? nevertheless i succed to have the appropriate syntax for the user routine, and the user output is now running OK. Many thanks for this wonderfull tool which are Code_Saturne, and this user forum.

Fred

[Fabien852]

Thank you, it works !