There is a parameter in the subroutine posdat:ICLMAP(I)

It is explained in this way:

ICLMAP(I) (PROSTAR)cell number corresponding to STAR cell index I.

So what is the exact meaning of STAR cell index I?

For example,if I get one cell number in PROSTAR which is 10,when I "load STAR transient data " by using cell number 10 from the pstt file by the option load transient post data,I found the results I got are different with the results gotten by using the subroutine posdat. What is the matter?

Is there some difference between the Prostar cell number(10) and the cell number used in posdat(ICLMAP(10)?)?

Thanks a lot for your help!

Lee

> Is there some difference between the Prostar cell number(10) and the cell number used in posdat(ICLMAP(10)?

Yes, there is. ICLMAP accepts STAR cell number and returns PROSTAR cell number. You are using ICLMAP in the wrong way. To get STAR cell number corresponding to PROSTAR cell 10 you have to build the reversed index by looking through all the STAR cell and getting their numbers with ICLMAP.

Thank you,anton!

But how to look through the STAR cell?I mean whether I can do it by a subroutine or other ways.

Whether were the PROSTAR cell number changed into STAR cell number at the moment the model was runing?

Thank you for your help!

Lee

The PROSTAR number is the one you get with ICLMAP(I), where I is the internal STAR cell number (= internal solver cell number). When you write your geometry file, the cell are renumbered to get a small executable as possible. For post-processing, you need to know the PROSTAR cell number attributed to the internal cell after renumbering. This is why you need ICLMAP.

If you want to get values for a particular PROSTAR cell number, say #10, you have to loop over all cells, like:

DO I=1,NCELL

IF (ICLMAP(I).EQ.10) THEN

....

....

END IF END DO