Question:

When I tried to make a user subroutine "Postdat", I need to know the global coordinates of cell centers. However, I checked all the data avialable in Postdat file and did not find these data. I don't know if there is no these data in this User Subroutine or I have to use some other method to obtaine these coordinates. If some one has ideas about this, please let me know. Thanks H. Shen

Dear Shen,

Posdat is an extremely powerful routine which can be used for almost any purpose. STAR is designed such a way that you can access almost any data. The best way to achieve it is by using the common block such that

C*****************************

SUBROUTINE POSDAT(KEY,...)

C*********************************

INCLUDE '../parm.inc'

COMMON/DCO21/ CX(3,-NBMAX:NCTMAX)

COMMON /C1C001/ ICLMAP(4)

where CX is the global coordinates of cell center. For example, I want to visit all coordinate:

do ip=1,ncell (ncell - total no. of fluid cells)

x=cx(1,ip)

y=cx(2,ip)

z=cx(3,ip)

enddo

For any star cell no., you can also get its PROSTAR cell no.

IPprostar=ICLMAP(IPstar)

Obviously, with access to all STAR's named commons you can do just about anything you like. However, I can find no documentation that describes the named common blocks: DCO21 and C1C001, so I can't afford to risk using them.

Is there a manual reference for this information?

Dear Mr. Luo: I have got the point. Thanks Shen

Hi, Steve,

All common blocks can be found in $STARDIR/solve/starin.f. It's not quite possible to describe all variables. But if you need some specific information, please let me know.