hm_read_mat_user29_31.F File Reference

#include "implicit_f.inc"
#include "scr15_c.inc"
#include "scr17_c.inc"
#include "param_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_mat29_31 (ilaw, iuser_key, userl_avail, uparam, maxuparam, nuparam, nuvar, ifunc, maxfunc, nfunc, parmat, lsubmodel, pm, matparam)

Function/Subroutine Documentation

hm_read_mat29_31()

subroutine hm_read_mat29_31	(	integer	ilaw,
		character(len=ncharline)	iuser_key,
		integer	userl_avail,
		dimension(maxuparam)	uparam,
		integer	maxuparam,
		integer	nuparam,
		integer	nuvar,
		integer, dimension(maxfunc)	ifunc,
		integer	maxfunc,
		integer	nfunc,
		dimension(100)	parmat,
		type(submodel_data), dimension(nsubmod), intent(in)	lsubmodel,
		intent(inout)	pm,
		type (matparam_struct_), intent(inout)	matparam )

Definition at line 40 of file hm_read_mat_user29_31.F.

C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   Read user material 29 to 31 - generic routine
C------------------------------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C
C     NAME          DESCRIPTION                         
C
C     USERL_AVAIL   Flag if userlibrary was load
C     MATNUM        Number of Material 29, 30, 31
C     IS_AVAILABLE  Bool / Result of HM_interface
C     LSUBMODEL     SUBMODEL Structure.
C------------------------------------------------------------------
C
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE submodel_mod , ONLY : nsubmod, submodel_data
      USE matparam_def_mod
      USE names_and_titles_mod , ONLY : nchartitle, ncharline
      USE format_mod , ONLY : fmt_2f
      USE reader_old_mod , ONLY : kcur, kline, line,key0
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      "scr15_c.inc"
#include      "scr17_c.inc"
#include      "param_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ILAW
      INTEGER USERL_AVAIL
      INTEGER MAXUPARAM
      INTEGER NUPARAM 
      INTEGER NUVAR
      INTEGER NFUNC,MAXFUNC
      TYPE(SUBMODEL_DATA),INTENT(IN)  :: LSUBMODEL(NSUBMOD)
      INTEGER, DIMENSION(MAXFUNC)     :: IFUNC
      my_real, DIMENSION(100)         :: parmat
      my_real, DIMENSION(MAXUPARAM)   :: uparam
      my_real, DIMENSION(NPROPM) ,INTENT(INOUT)   :: pm
      CHARACTER(LEN=NCHARLINE) :: IUSER_KEY
      TYPE (MATPARAM_STRUCT_) ,INTENT(INOUT) :: MATPARAM
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      CHARACTER(LEN=4096) :: SCR_FILE_NAME
      CHARACTER(LEN=NCHARLINE) :: RLINE
      CHARACTER (LEN=4) :: CLAW
      LOGICAL :: IS_AVAILABLE
      INTEGER NLINES,J
      INTEGER SCR_FILE_NAME_LEN
      my_real rho0,rhor
!
      CHARACTER OPTION*256
      INTEGER SIZE
C-----------------------------------------------
      is_available = .false.
!
      IF (userl_avail == 0)THEN
      ! ERROR to be printed & exit
        option='/MAT/'//iuser_key
        size=len_trim(option)
        CALL ancmsg(msgid=1130,
     .              msgtype=msgerror,c1=option(1:size),anmode=aninfo)
        CALL arret(2)
      ENDIF
 
     
      CALL hm_get_intv ('Number_of_datalines' ,nlines ,is_available, lsubmodel)
 
      IF(nlines > 1)THEN
 
        ! Rho & Rho0_r are read out of user laws
        j=1
        CALL hm_get_string_index('arraydatalines', rline, 1, ncharline, is_available)
        READ(rline,err=999,fmt=fmt_2f)rho0,rhor
 
        pm(1)=rhor
        pm(89)=rho0
 
        ! Create tempo file  
        WRITE(claw,'(I4.4)')ilaw
        scr_file_name='SI'//rootnam(1:rootlen)//'_'//claw//'.scr'
        scr_file_name_len=len_trim(scr_file_name)
        OPEN(unit=30,file=trim(scr_file_name),form='FORMATTED',recl=ncharline)
 
        ! Read & Dump in scratch file
        DO j=2,nlines   
         CALL hm_get_string_index('arraydatalines', rline, j, ncharline, is_available)
         WRITE(30,fmt='(A)')trim(rline)
        ENDDO
        CLOSE(unit=30)
  
        CALL st_userlib_usermat(ilaw,rootnam,rootlen,
     .                          uparam,maxuparam,nuparam,
     .                          nuvar ,ifunc,maxfunc,nfunc   ,parmat)
 
        CALL user_output(iout,ilaw,rootnam,rootlen,1)
      ELSE
         ! Error MESSAGE_Mess
      ENDIF
 
      ! Properties compatibility
      CALL init_mat_keyword(matparam,"SOLID_ALL")
      CALL init_mat_keyword(matparam,"SHELL_ALL")
      CALL init_mat_keyword(matparam,"SPH")
 
      RETURN
 999  CALL ancmsg(msgid=55,anmode=aninfo,msgtype=msgerror,c1=key0(kcur),c2=kline,c3=line)
      CALL arret(2)