law112_upd.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	law112_upd (titr, uparam, npc, pld, numtabl, itable, table, mat_id)

Function/Subroutine Documentation

law112_upd()

subroutine law112_upd	(	character(len=nchartitle)	titr,
			uparam,
		integer, dimension(*)	npc,
			pld,
		integer	numtabl,
		integer, dimension(numtabl)	itable,
		type(ttable), dimension(ntable)	table,
		integer	mat_id )

Definition at line 33 of file law112_upd.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE table_mod
      USE names_and_titles_mod , ONLY : nchartitle
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      CHARACTER(LEN=NCHARTITLE)  :: TITR
      INTEGER MAT_ID,NUMTABL,ITABLE(NUMTABL)
      INTEGER NPC(*)
      my_real uparam(*),pld(*)
      TYPE(TTABLE), DIMENSION(NTABLE) ::  TABLE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,K,ITAB
      my_real slopemin,dx,dy
C=======================================================================
      ! Recovering tabulated flag
      itab = nint(uparam(14))
c
      ! Checking yield stress functions slope                                                         
      IF (itab > 0) THEN    
        DO i = 1,numtabl                                          
          slopemin = infinity
          DO k = 2,SIZE(table(itable(i))%X(1)%VALUES)
            dx = table(itable(i))%X(1)%VALUES(k) - table(itable(i))%X(1)%VALUES(k-1)
            dy = table(itable(i))%Y%VALUES(k) - table(itable(i))%Y%VALUES(k-1)
            slopemin = min(slopemin,dy/dx)
          ENDDO    
          IF (slopemin < zero) THEN
            CALL ancmsg(msgid=2059 ,
     .               msgtype=msgerror,
     .               anmode=aninfo_blind_2,
     .               i1=mat_id,
     .               i2=table(itable(i))%NOTABLE,
     .               c1=titr)
          ENDIF
        ENDDO  
      ENDIF
c-----------
      RETURN