Dear XERF!
I read the your answers in this forum to get state variables. I copied the method to receive the coordinates by SDVINI then to give them to HETVAL. I have inserted the "*INITIAL CONDITIONS, TYPE=SOLUTION, USER" into the Model, but I cant run the Subroutine SDVINI. I want get the coordinates of nodes as well node number to transfer to Subrountine HETVAL. I cant get the Coords(3). What is wrong?
My Subroutines:
----------------------
SUBROUTINE SDVINI(STATEV, COORDS, NSTATV, NCRDS, NOEL, NPT,
1 LAYER, KSPT)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME
DIMENSION STATEV(NSTATV), COORDS(NCRDS)
C
STATEV(1)=NOEL
STATEV(2)=NPT
STATEV(3)=COORDS(1)
STATEV(4)=COORDS(2)
STATEV(5)=COORDS(3)
RETURN
END
SUBROUTINE HETVAL(CMNAME, TEMP, TIME, DTIME, STATEV, FLUX,
1 PREDEF, DPRED)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION COORDS(3)
INTEGER NPT, NOEL, OpenStatus
REAL A
CHARACTER*80 CMNAME
C
DIMENSION TEMP(2), STATEV(5), PREDEF(*), TIME(*), FLUX(2),
1 DPRED(*)
COORDS(1) = STATEV(3)
COORDS(2) = STATEV(4)
COORDS(3) = STATEV(5)
IF (COORDS(1)>10) FLUX(1)=1000
RETURN
END
-----------
I think It is good, but I cant take the coord(3) to analysis. Can U help me?