code_saturne User's Forum

Hello every body,

I work to measure Y+ on all cells of a profile, like a cylinder, in Code Saturne 1.3.3.

I have wrote the subroutine usproj which I join but i have the error message :
usproj.F gfortran  -D_CS_HAVE_XML  -D_CS_LANG_FR                   -I. -O1 -c usproj.F
In file usproj.F:558
        YPLUS(IEL) = RA(IYPPAR+IEL-1)                                                                                                          1
Error: Unexpected STATEMENT FUNCTION statement at (1)
 In file usproj.F:586
     &              CDGFBO(3, IFAC), YPLUS(IEL)
                                                            1
Error: Function 'yplus' at (1) has no IMPLICIT type
 
Can you help me please? In fact, I have different tests with Code Saturne and verify the function with the users' guide but I have'nt the solution.

Thanks you everybody

Sebastien

You can have a maximum of 72 characters per line in FORTRAN77 so the YPLUS(NFABOR) variable is not recognised. Put it on a new line such as: If you use something like vi to view the file it will highlight when you go over the max number of characters on a line.

Also you define YPLUS to be the size of NFABOR (number boundary faces) but later refer to it with IEL (cell-id), it's quite probable this will lead to errors. I think you should refer to it as YPLUS(IFAC).

Hope that helps,

James

Thank you for your answer.
 
I resolved the problem of implicit and function. But I have a new problem. When I exxecute my case on Code Saturne on 1 proc on my desktop computer, i have a message : ERREUR WHILE RUNNING THE KERNEL (SOLVER). I join my subroutine and the error file and the listing file.
Can you help me please?  
Thanks you everybody
 
Sebastien

Hello Sebastien,
 
FINDPT function isn't properly filled line 565 in your usproj.f90.
 
Best regards
 
Mickaël Hassanaly

Thank you very much
 
Best regards
 
Sebastien Thomas

ello every body,

I work to measure Y+ on all cells of a profile NACA and on a cylinder, in Code Saturne 1.3.3.

I wrote the subroutine USPROJ.F which is joined.

This subroutine run but, I have problem when I interpret my results because I notice the value zero for some points but in paraview they have not this value.


Can you help me please? I attach my subroutine usprof and an exemple of results file, Y+

Hello,
Do you run your case on multiple processors?
Mickaël