Y+ problem in Code Saturne 1.3.3

[Sébastien Thomas]

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

[James McNaughton]

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:

Code: Select all

DOUBLE PRECISION YPLUS(NFABOR)

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

[Sébastien Thomas]

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

[Mickael Hassanaly]

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

[Sébastien Thomas]

Thank you very much
 
Best regards
 
Sebastien Thomas

[Sébastien 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+

[Mickael Hassanaly]

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