posdat.f Problems
 

I uses STAR-CD v3.26

When OPEN function is included, this kind of problem is occuring.

If posdat.f is excluded, it is working very well without any problems

******************* posdat.f *************************

~~~~~~~ abbrev.~~~~~~

IF(LEVEL.EQ.2) THEN

IF(INTFLG(1).EQ.0) THEN

OPEN (84,FILE='MonP.dat',FORM='FORMATTED',STATUS='UNKNO WN') ==> Error position!!!

DO I=1,INUMBER

DO J=1,NCTMXU

IF(IPROSTAR(I) .EQ. ICLMAP(J)) THEN

~~~~~~~ abbrev.~~~~~~

************************************************** *****

******************* Running ********************************

1 5.40E-01 3.13E-01 1.33E+00-1.01E-04 1.40E-05 1.27E-02 1.47E+08 -1.06E-03-1.51E-03 6.03E+00 2.02E+00 6.84E-01 1.92E+01 1.94E+03

-1.97E-07-5.33E-06-3.42E-06 3.58E-07-8.40E-06-7.34E-05 0.00E+00 4.19E-03 6.44E-02 6.35E-05 2.06E-03 9.86E-02 7.18E-01 0.00E+00

-1.60E-05 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.78E-01 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00

[ 1 ] ***Received signal sigsegv - exiting

PNP: Shutdown [2006-02-16-04:13:15] Execution aborted by request (SIGABRT) after 9 seconds (TOTAL ELAPSED TIME).

1060% [kris@node01 transient]#

************************************************** ****

Re: posdat.f Problems
 

It's not easy to see the problem because you only put part of your code, try this IF(LEVEL.EQ.2) THEN

IF(INTFLG(1).EQ.0) THEN OPEN(84,FILE='MonP.dat',FORM='FORMATTED',STATUS='U NKNOWN') INTFLG(1)=1 ENDIF .....

Re: posdat.f Problems
 

in this code

If OPEN command is NOT included, No Error occures.

If OPEN command is included, Error occures and the result is as follow

******************* Running ********************************

1 5.40E-01 3.13E-01 1.33E+00-1.01E-04 1.40E-05 1.27E-02 1.47E+08 -1.06E-03-1.51E-03 6.03E+00 2.02E+00 6.84E-01 1.92E+01 1.94E+03

-1.97E-07-5.33E-06-3.42E-06 3.58E-07-8.40E-06-7.34E-05 0.00E+00 4.19E-03 6.44E-02 6.35E-05 2.06E-03 9.86E-02 7.18E-01 0.00E+00

-1.60E-05 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.78E-01 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00

[ 1 ] ***Received signal sigsegv - exiting

PNP: Shutdown [2006-02-16-04:13:15] Execution aborted by request (SIGABRT) after 9 seconds (TOTAL ELAPSED TIME).

1060% [kris@node01 transient]#

************************************************** ****

POSDAT.F Test Source file

!************************************************* ************************

SUBROUTINE POSDAT(KEY,VOL,U,TE,ED,T,P,VIST,DEN,CP,VISM,CON,

* F,ICLMAP,ICTID,RESOR,VF,FORCB,IRN,PREFM,LEVEL) !C Post-process data !C************************************************ ************************* !C--------------------------------------------------------------------------* !C STAR VERSION 3.20.000 * !C--------------------------------------------------------------------------*

INCLUDE 'comdb.inc'

COMMON/USR001/INTFLG(100)

DIMENSION KEY(-NBMAXU:NCTMXU),VOL(NCTMXU),U(3,-NBMAXU:NCMAXU),

* TE(-NBMAXU:NCMAXU),ED(-NBMAXU:NCMAXU),T(-NBMAXU:NCTMXU,1+NSCU),

* P(-NBMAXU:NCMAXU),VIST(-NBMAXU:NCMAXU),DEN(-NBMAXU:NCTMXU),

* CP(-NBMAXU:NCTMXU),VISM(-NBMXVU:NCMXVU),CON(-NBMXCU:NCMXCU),

* F(3,-NBMAXU:NCMAXU),ICLMAP(NCTMXU),ICTID(NCTMXU),

* RESOR(89,-100:100),VF(NCDMXU),FORCB(3,NWLMX),IRN(NWLMX)

DOUBLE PRECISION P

DOUBLE PRECISION T

DOUBLE PRECISION TE

DOUBLE PRECISION ED

DOUBLE PRECISION U

DIMENSION PREFM(4)

INCLUDE 'usrdat.inc'

DIMENSION IPROSTAR(100),ISTAR(100),IBOOL(100)

DATA IPROSTAR/100*1/

DATA IBOOL/100*0/

DATA ISTAR/100*1/

IF(LEVEL.EQ.2) THEN

IF(INTFLG(1).EQ.0) THEN

OPEN(84,FILE='MonP.dat',FORM='FORMATTED', STATUS='UNKNOWN')

CLOSE(84)

700 FORMAT(A8,1P,I20,I20,I20,I20,I20,I20,I20,I20,I20,I 20,I20

* ,I20,I20,I20,I20,I20,I20,I20,I20,I20,I20,I20) 800 FORMAT(F10.8,1X,E20.13,1X,E20.13,1X,E20.13,1X,E20. 13,

* 1X,E20.13,1X,E20.13,1X,E20.13,1X,E20.13,1X,E20.13, 1X,E20.13,

* 1X,E20.13,1X,E20.13,1X,E20.13,1X,E20.13,1X,E20.13, 1X,E20.13,

* 1X,E20.13,1X,E20.13,1X,E20.13,1X,E20.13,1X,E20.13, 1X,E20.13)

ENDIF

ENDIF

RETURN

end