May I use global coordinate of cell centers in subroutine Posdat
 

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

Re: May I use global coordinate of cell centers in subroutine Posdat
 

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)

Re: May I use global coordinate of cell centers in subroutine Posdat
 

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?

Thanks
 

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

Re: May I use global coordinate of cell centers in subroutine Posdat
 

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.