I am trying to access pressure and position of a boundary region using posdat.f I identify mine region with the following code.

do 220 ib=1,NBP23

irnb=IRN(ib)

if(irnb.eq.27) then

ENDIF 220 continue

what I need to access pressure and position of the boundary? Can someone help me?

Star3.26 by the way

stop using CONTINUE from fortran60, use do/enddo and indent your code. this is rubish

I think you can access pressure simply by the variable P(ib).

For accessing the centroid of a boundary in postdat you need to add the following pointer to your coding:

POINTER (P_CX,CX(3,-NBMAX:NCTMAX))

COMMON /DC021/ P_CX

CX(I,IB) returns you the coordinates of the cell center were 'I' can take the Values 1,2,3 for X,Y or Z coordinate.

Regards

andreas

Hi all

It appears to not work properly.

The boundaries in the CX(I,IB) array are not in the IB position.

Can someone clarify this please?

thanks

jmcp

I need to know simple aircraft design rules for homebuilet aircrafts.