'Error: "The leftmost part-ref in a data-ref can not be a function reference."

I am trying to write a FORTRAN subroutine for ABAQUS that would modify the seepage coefficient and thus the flow depending on whether there is contact on the surface. In order to do so I need 2 subroutines URDFIL to retrieve the node data and FLOW to modify the seepage coefficient.

However, when I compile my subroutine I get the following errors:

flow_us.for(81): error #6837: The leftmost part-ref in a data-ref can not be a function reference. [K_ELE_DETAILS]
       IF(K_ELE_DETAILS(E_INDX)%IS_CONT(N_INDX).EQ.0)THEN
----------^
flow_us.for(81): error #6158: The structure-name is invalid or is missing. [K_ELE_DETAILS]
       IF(K_ELE_DETAILS(E_INDX)%IS_CONT(N_INDX).EQ.0)THEN

Obviously it is repeated for the 3 lines (if's) that contain such structure (81,8, and 89).

Please find the code below and hopefully someone will be able to help

****************************************************************************************
***SUBROUTINE FOR ADAPTIVE FLUID FLOW
****************************************************************************************
****************************************************************************************
**
**
**
*USER SUBROUTINE
    SUBROUTINE URDFIL(LSTOP,LOVRWRT, KSTEP, KINC, DTIME, TIME)
    INCLUDE 'ABA_PARAM.INC'
C
    DIMENSION ARRAY(513),JRRAY(NPRECD,513),TIME(2)
    EQUIVALENCE (ARRAY(1),JRRAY(1,1))
C DECLARATIONS
    TYPE ELE_DATA
        SEQUENCE
        DOUBLE PRECISION :: NODE_COORD(9)
        DOUBLE PRECISION :: OPP_NODE_COORD(9)
        DOUBLE PRECISION :: IPT_COORD(9)
        DOUBLE PRECISION :: POR(3)
        DOUBLE PRECISION :: OPP_POR(3)
        INTEGER :: ELE_NUM
        INTEGER :: OPP_ELE_NUM
        INTEGER :: NODE_NUM(3)
        INTEGER :: OPP_NODE_NUM(3)
        INTEGER :: IPT_NUM(3)
        INTEGER :: OPP_IS_CONT(3)           
    END TYPE ELE_DATA

    TYPE(ELE_DATA)::K_ELE_DETAILS(500)
    COMMON K_ELE_DETAILS
    PARAMETER (THRESHOLD_CSTRESS=1.0E-6)
*******************************************************
    INTEGER :: NO_OF_NODES
    INTEGER :: NO_OF_ELEMENTS
    INTEGER :: NO_OF_DIM
    COMMON NO_OF_DIM, NO_OF_NODES, NO_OF_ELEMENTS
********************************************************
C INITIALIZE   
    LSTOP=0
    LOVRWRT=1
    LOP=0
    NO_OF_NODES=10000
    NO_OF_ELEMENTS=10000
    NO_OF_DIM=2

        DO K1=1,999999
            CALL DBFILE(0, ARRAY,JRCD)
            IF (JCRD.NE.0) GO TO 110
            KEY=JRRAY(1,2)
*******************************************************
C THE KEYS USED IN THE FOLLOWING LINES REFER
C TO INFORMATION ON THE SURFACE, NODES, CONTACT ETC
            IF (KEY.EQ.1501) THEN
            ELSE IF (KEY.EQ.1502)THEN
            ELSE IF(KEY.EQ.1911) THEN
            ELSE IF(KEY.EQ.108.AND.SURFACE_N_SET.EQ.'N_TOP') THEN
            ELSE IF(KEY.EQ.107.AND.SURFACE_N_SET.EQ.'N_TOP') THEN
            ELSE IF(KEY.EQ.1503)THEN
            ELSE IF(KEY.EQ.1504.AND.K_NODE_SET.EQ.'N_TOP')THEN
            ELSE IF(KEY.EQ.1511.AND.K_NODE_SET.EQ.'N_TOP')THEN
C IS THE NODE IN CONTACT?
120  CONTINUE
        END IF
        END DO 
110  CONTINUE
            RETURN
            END

**********************************************************
**********************************************************
*USER SUBROUTINE
    SUBROUTINE FLOW(H, SINK, KSTEP, KINC, TIME, NOEL, NPT, COORDS,
    1 JLTYP,SNAME)
    INCLUDE 'ABA_PARAM.INC'
C
    DIMENSION TIME(2), COORDS(3)
    CHARACTER*80 SNAME

    MIN_DIST=10E-20
    DO K25=1,NO_OF_ELEMENTS
        DO K26=1,NO_OF_NODES
C FINDS THE CLOSEST NODE TO THE INTEGRATION POINT NPT
        END DO
    END DO
C NOT IN CONTACT
    IF(K_ELE_DETAILS(E_INDX)%IS_CONT(N_INDX).EQ.0)THEN
        IF(K_ELE_DETAILS(E_INDX)%POR(N_INDX).GE.0) THEN
            SINK=0
            H=0.001
        ELSE
            SINK=0
            H=1
        END IF
    ELSE IF(K_ELE_DETAILS(E_INDX)%IS_CONT(N_INDX).EQ.1)THEN
C IF THERE IS CONTACT
        SINK=0
        H=0
    END IF
    RETURN
    END
fortransubroutineabaqus


Solution 1:[1]

I think you meant to write

IF(K_ELE_DETAILS(E_INDX)%OPP_IS_CONT(N_INDX).EQ.0) THEN

instead of

IF(K_ELE_DETAILS(E_INDX)%IS_CONT(N_INDX).EQ.0)THEN

Also, if you intent to extend this code with type-bound procedures, you cannot reference a component of a function result of derived data type without an appropriate interface. Such issues are rarely encountered if the type declaration is confined to a module. The module will automatically generate the correct interface for type-bound procedures.

Sources

This article follows the attribution requirements of Stack Overflow and is licensed under CC BY-SA 3.0.

Source: Stack Overflow

Solution Source
Solution 1 jlokimlin

Related Questions