zero velocity at the wall

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
Ionut G
Posts: 43
Joined: Fri Apr 20, 2018 2:43 pm

zero velocity at the wall

Post by Ionut G »

Hello,

I am trying to set from "cs_user_boundary_conditions.f90" zero velocity at the wall.
I am using this code:

Code: Select all

call getfbr('PLATE', nlelt, lstelt)
do ilelt = 1, nlelt

  ifac = lstelt(ilelt)
  
  rcodcl(ifac, iu, 1) = 0.d0
  rcodcl(ifac, iv, 1) = 0.d0
  rcodcl(ifac, iw, 1) = 0.d0
  
enddo
I get this error: SIGSEGV signal (forbidden memory area access) intercepted!

What I'm doing wrong?

Thank you,
ionut
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: zero velocity at the wall

Post by Yvan Fournier »

Hello,

Did you allocate the lstelt array?

Regards,

Yvan
Ionut G
Posts: 43
Joined: Fri Apr 20, 2018 2:43 pm

Re: zero velocity at the wall

Post by Ionut G »

Hi Yvan,

yes, i've allocated it.
I attached the file too.

Thank you,
Ionut
Attachments
cs_user_boundary_conditions.f90
(20.88 KiB) Downloaded 191 times
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: zero velocity at the wall

Post by Yvan Fournier »

Hello,

Why did you comment some arguments from the cs_f_user_boundary_conditions subroutine ?

This is probably the cause of the crash.

Regards,

Yvan
Ionut G
Posts: 43
Joined: Fri Apr 20, 2018 2:43 pm

Re: zero velocity at the wall

Post by Ionut G »

Hi,

I thought that if I don't need them I could comment them. I didn't want to use more memory if is not needed for the software.

I uncommented them and now is working.

Thank you,
Ionut
Post Reply