code_saturne User's Forum

Hi Code Saturne develop team:

I am using a subroutine to impose the thermal source term for each cell in a channel flow. In the subroutine, I need to use the cell density and velocity (z direction component) to calculate the heat source term. When I try to compile the subroutine, there is an error comes from the cs_field API, the screen shot of the error is shown, and my heat source term subroutine is attached. Did anyone know what the problem is, thanks!

Regards!
Jundi

Hello,

Yes, cvar_vel is a 2-dimensional array (dim*cells), so you shoumd not declare it with the other scalar arrays on line 49, but as: on a nearby line (note the double instead of single colons).

Regards,

Yvan

Yvan:

OK, thanks! I have another question, if I need to access to the z direction velocity, should I use this:

call field_get_val_v(ivarfl(iw), cvar_vel)
zvelocity=cvar_vel(1,iel)

Regards!
Jundi

Hello,

No, that is the x direction; for the z direction, use:
zvelocity=cvar_vel(3,iel)

Best regards,

Yvan

OK, thanks a lot!