pressure inlet / subroutines

[koch]

Hello,

I tried hard to find solution for inlet imposed pressure boundary condition. I understand that the only solution is to use subroutines. There is a topic http://code-saturne.org/forum/viewtopic.php?f=9&t=627 where the solution is shown.

When I paste this to my cs_user_boundary_conditions.f90 file, I got an error:

raise RunCaseError('Compile or link error.')
cs_case_domain.RunCaseError: Compile or link error.

I verified log file and there were two types of errors:

- Error: Invalid form of array reference at (1)

- Error: Symbol ‘iel’ at (1) has no IMPLICIT type

I attach log file.

First I found that there was a mistake in subroutine - rcdocl -> rcodcl but it was no difference in output.

I thought that maybe it is a problem with fortran, but I tried to run with empty boundary conditions subroutine, solver started but I got another error:

Error running the calculation.
Check Code_Saturne log (listing) and error* files for details.


In listing I got:

First face with boundary condition definition error
(out of 19)
has boundary condition type 32764, center (1.63787e-316, 6.9396e-310, 6.9396e-310)
cs_boundary_conditions.c:640: Fatal error.
Some boundary condition definitions are incomplete or incorrect.



I got crazy about that. I don't know how to start with subroutines and with this inlet pressure. Can anybody help me?

I use Code_Saturne 4.0.5

Regards
Maciek

[Yvan Fournier]

Hello,

You example is from a very old version, and the syntax has changed. You mostly need to remove the (iphas) references, which have been removed. If iel is not defined, you need to define a local integer variable named iel. Did you check the user boundary condition examples from the code documentation (http://code-saturne.org/doxygen/src/cs_ ... mples.html) ?

Regards,

Yvan