How is it possible to get the Star or Prostar cell number of the fluid cell, which is attached to the boundary number "IBP" in the subroutine bcdefw.f ?

Hi Joern,

Maybe this snippet of code from Adapco-Online could help:

do 220 ib=1,nbp23

ip=lx(ib)/7

indx=mod(lx(ib),7)

ibnd=lq(indx,ip)

ibndp=iblmap(-ibnd) 220 continue

Given a certain IB (IBP if passed thru IBLMAP; here it's ibndp), you can get IP with the help of the COMMON LX(IB).

I had this type of coding in mind when I started but it is much easier:


IF(IREG.EQ.2) THEN
WRITE(*,*) ' BCDEFW called for boundary ' , IBP ,
WRITE(*,*) ' BCDEFW called for Prostar CELL ' , IP ,
WRITE(*,*) ' BCDEFW called for Star CELL ' , IPSTAR ,
ENDIF


IP and IPSTAR are already available in the subroutine.

This doesn't work for baffles. IBP is always zero if you are on "side 2" of the baffle The only way round I could find was to write a Prostar subroutine to output a table (using graph registers) of fluid cell and boundary numbers. This table can then be read into BCDEFW and used to find the real IBP number from the IP number by lookup, when the IBP in the subroutine returns zero.

If anyone else has seen this problem or has other solutions I would be interested to know.