Hello,
I modified user_properties.f for density, viscosity and thermal conductivity, and I got the following message “Compressible: IJK_P_g remaining undefined.” then the run stopped.
Can you advice me to solve it?
Attach project files
Attach zip file generated from main menu> submit bug report.
1dheatconduction_2022-06-13T225634.793504.zip (222.4 KB)
When you write udfs or modify the code, you should always build the solver with debug flags. This will often point to the issue:
Look at the build output:
source_u_g.f:307:0:
S=(-C_M*(((1-F)**2)/(((F)**3)+E)))* U_g(IJK)
Warning: \u2018f\u2019 may be used uninitialized in this function [-Wmaybe-uninitialized]
The code stops due to out of bound error:
At line 225 of file usr_properties.f
Fortran runtime error: Index '214' of dimension 1 of array 'f' above upper bound of 0
K_g(IJK)=(1-F(IJK))*k_s+(F(IJK)*k_l)
F is defined as a scalar but you use it as an array.