comfiling problem ..
 

Hi all I set the density based on USER-f(P.T) in morecular property,and i made simple fortran code below. but densit.f doesn't comfiled so that the output of density shows constant value wich i gave as a initial value. what is problem in this case.

regards

mspark

fortran code

SUBROUTINE DENSIT(DEN,DENDP)

INCLUDE 'comdb.inc'

COMMON/USR001/INTFLG(100)

INCLUDE 'usrdat.inc'

DIMENSION SCALAR(50)

EQUIVALENCE( UDAT12(001), ICTID )

EQUIVALENCE( UDAT11(001), CP )

EQUIVALENCE( UDAT11(005), PR )

EQUIVALENCE( UDAT11(006), P )

EQUIVALENCE( UDAT11(007), T )

EQUIVALENCE( UDAT11(009), SCALAR(01) )

EQUIVALENCE( UDAT11(059), U )

EQUIVALENCE( UDAT11(060), V )

EQUIVALENCE( UDAT11(061), W )

EQUIVALENCE( UDAT11(067), X )

EQUIVALENCE( UDAT11(068), Y )

EQUIVALENCE( UDAT11(069), Z )

DENSITYL1=96693.*T**(-0.7684)

DENSITYV1=4.*10.**(-22.)*T**9.2727

DENSITYL=3589.7*(P)**-0.081+0.001/DENSITYL1

DENSITYV=7E-05*(P)**0.9772+0.001/DENSITYV1

SV_F=1/DENSITYL

SV_G=1/DENSITYV

SV_FG=abs(SV_G-SV_F)

REALH=CP*T

VT=72.424*(P)**0.1054

HFLUID=7.2051*VT**1.8239

HLATENT=20000.*10000*VT**(-1.2023)

QUALITY=(REALH-HLATENT)/HFLUID

DEN=1/(SV_F+QUALITY*SV_FG)

RETURN

END