Turbulent_heat_flux[Y]

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
C0st4s
Posts: 48
Joined: Fri May 11, 2018 12:21 pm

Turbulent_heat_flux[Y]

Post by C0st4s »

Hi,

I have a question regarding printing the turbulent_heat_flux[Y] of the EBRSM model. Is there a way that I can call this variable and print it in a file using cs_user_extra_operations_1D_profile.f90? I noticed that it is already printed in the monitoring points and I want it printed along a line profile.

Thanks,
Costas
Yvan Fournier
Posts: 4251
Joined: Mon Feb 20, 2012 3:25 pm

Re: Turbulent_heat_flux[Y]

Post by Yvan Fournier »

Hello,

Yes, you just need to adapt the function to use the correct field. You can find the exact field name in the "setup.log" file, but unless you have modified its "output label", it should be the same as for probes (without the [Y] part, which is for the component).

Regards,

Yvan
C0st4s
Posts: 48
Joined: Fri May 11, 2018 12:21 pm

Re: Turbulent_heat_flux[Y]

Post by C0st4s »

Right, thanks so it will be something like call field_get_val_s_by_name() and what is the parallel assignment?

Regards,
Costas
Yvan Fournier
Posts: 4251
Joined: Mon Feb 20, 2012 3:25 pm

Re: Turbulent_heat_flux[Y]

Post by Yvan Fournier »

Hello,

Depending on the field's dimension (easy to check in setup.log), you will need field_get_val_s_by_name() or field_get_val_v_by_name().

I am not sure about what you call the parallel assignment. To avoid parallel issues, using the newer profile system is better. Check the cs_user_postprocess-profiles.c example (new system) instead of cs_user_extra_operations_1D_profile.f90 (old system).

Best regards,

Yvan
C0st4s
Posts: 48
Joined: Fri May 11, 2018 12:21 pm

Re: Turbulent_heat_flux[Y]

Post by C0st4s »

Dear Yvan,

Many thanks, I have sorted it out. The parallel operation is call parall_bcast_r(irangv,uu) for r11 array which I assigned to double precision uu.

Regards,
Costas
Post Reply