hm_read_initemp.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    hm_read_initemp            ../starter/source/initial_conditions/thermic/hm_read_initemp.F
!||--- called by ------------------------------------------------------
!||    lectur                     ../starter/source/starter/lectur.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                     ../starter/source/output/message/message.F
!||    hm_get_float_array_index   ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
!||    hm_get_floatv              ../starter/source/devtools/hm_reader/hm_get_floatv.f
!||    hm_get_int_array_index     ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
!||    hm_get_intv                ../starter/source/devtools/hm_reader/hm_get_intv.F
!||    hm_option_read_key         ../starter/source/devtools/hm_reader/hm_option_read_key.f
!||    hm_option_start            ../starter/source/devtools/hm_reader/hm_option_start.F
!||    udouble                    ../starter/source/system/sysfus.F
!||    usr2sys                    ../starter/source/system/sysfus.F
!||    usrtos                     ../starter/source/system/sysfus.f
!||--- uses       -----------------------------------------------------
!||    hm_option_read_mod         ../starter/share/modules1/hm_option_read_mod.F
!||    message_mod                ../starter/share/message_module/message_mod.F
!||    submodel_mod               ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE hm_read_initemp(TEMP    ,NINTEMP ,ITHERM_FE,ITAB     ,ITABM1  ,
     .                           IGRNOD  ,INITIDS ,UNITAB   ,LSUBMODEL)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE message_mod
      USE groupdef_mod
      USE submodel_mod
      USE hm_option_read_mod
      USE names_and_titles_mod , ONLY : nchartitle, ncharkey
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      "com01_c.inc"
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
      INTEGER ITAB(*), ITABM1(*),INITIDS(*)
      INTEGER ,INTENT(IN) :: NINTEMP
      INTEGER ,INTENT(IN) :: ITHERM_FE
      my_real :: temp(*)
      TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)
C-----------------------------------------------
      TYPE (GROUP_)  , DIMENSION(NGRNOD)  :: IGRNOD
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,N,NNOD,NOSYS,ITYPE,ID,ISK,IGR,IGRS,NBTEMP,BID,
     .        flag_fmt,flag_fmt_tmp,ifix_tmp,uid,typ,nod,nodsys,nb_line
      my_real
     .  temp0
      CHARACTER MESS*40
      CHARACTER(LEN=NCHARTITLE) :: TITR
      CHARACTER(LEN=NCHARKEY) :: KEY
      LOGICAL IS_AVAILABLE
      my_real, DIMENSION(:), ALLOCATABLE :: list_temp0
      INTEGER, DIMENSION(:), ALLOCATABLE :: LIST_NOD
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER USR2SYS, USRTOS
      DATA mess/'INITIAL TEMPERATURE DEFINITION           '/
C=======================================================================
      bid =0
      nbtemp = 0
      is_available = .false.   
      nb_line = 0                      
C--------------------------------------------------
C START BROWSING MODEL INITEMP
C--------------------------------------------------
      CALL hm_option_start('/INITEMP')
C--------------------------------------------------
C BROWSING /INITEMP OPTIONS 1->NRADIA 
C--------------------------------------------------
      DO i=1,nintemp
        titr = ''
        CALL hm_option_read_key(lsubmodel,
     .                      unit_id = uid,
     .                      option_id = id,
     .                      option_titr = titr)
 
!        IF (ITHERM_FE == 0) THEN
!          CALL ANCMSG(MSGID=858,
!     .                MSGTYPE=MSGERROR,
!     .                ANMODE=ANINFO)
!        ENDIF
C--------------------------------------------------
C EXTRACT DATAS (INTEGER VALUES)
C--------------------------------------------------
        CALL hm_get_intv('distribution',typ,is_available,lsubmodel)
        CALL hm_get_intv('entityid',igr,is_available,lsubmodel)
        IF (typ == 1) THEN
          CALL hm_get_intv('grnd_ID',igr,is_available,lsubmodel)
          CALL hm_get_intv('distribution_table_count',nb_line,is_available,lsubmodel)
          IF(.NOT.ALLOCATED(list_nod)) ALLOCATE(list_nod(nb_line))
          list_nod(1:nb_line) = 0
          IF(.NOT.ALLOCATED(list_temp0)) ALLOCATE(list_temp0(nb_line))
          list_temp0(1:nb_line) = zero
          DO j=1,nb_line
            CALL hm_get_int_array_index('location_unit_node',list_nod(j),j,is_available,lsubmodel)
          ENDDO
        ENDIF
C--------------------------------------------------
C EXTRACT DATAS (REAL VALUES)
C--------------------------------------------------
        IF (typ == 0) THEN
          CALL hm_get_floatv('magnitude',temp0,is_available,lsubmodel,unitab)
        ELSEIF (typ == 1) THEN
          CALL hm_get_float_array_index('T0',temp0,j,is_available,lsubmodel,unitab)
          DO j=1,nb_line
            CALL hm_get_float_array_index('magnitude',list_temp0(j),j,is_available,lsubmodel,unitab)
          ENDDO
        ENDIF
C--------------------------------------------------  
        nbtemp = nbtemp+1
        initids(nbtemp)=id              
        igrs=0     
        IF (typ == 0) THEN                                           
          IF (igr == 0)THEN 
            CALL ancmsg(msgid=668,
     .                  msgtype=msgerror,
     .                  anmode=aninfo,
     .                  c1='/INITEM',
     .                  c2='/INITEM',
     .                  c3=titr,
     .                  i1=id)
          ENDIF                                                       
          DO j=1,ngrnod                                               
            IF (igr == igrnod(j)%ID) igrs=j
          ENDDO                                                       
          IF(igrs /= 0)THEN                                           
            DO j=1,igrnod(igrs)%NENTITY
              nosys=igrnod(igrs)%ENTITY(j)
              temp(nosys)= temp0
            ENDDO                                                     
            nnod=igrnod(igrs)%NENTITY
          ELSE
            CALL ancmsg(msgid=53,
     .                  msgtype=msgerror,
     .                  anmode=aninfo,
     .                  c1='IN /INITEM OPTION',
     .                  i1=igr)
          ENDIF   
        ELSE                                                   
          DO j=1,ngrnod                                               
            IF (igr == igrnod(j)%ID) igrs=j
          ENDDO                                                       
          IF(igrs /= 0)THEN                                           
            DO j=1,igrnod(igrs)%NENTITY
              nosys=igrnod(igrs)%ENTITY(j)
              temp(nosys)= temp0
            ENDDO                                                     
            nnod=igrnod(igrs)%NENTITY
          ENDIF      
          DO j=1,nb_line
            nodsys=usr2sys(list_nod(j),itabm1,mess,id)
            IF (list_nod(j) == 0) THEN
              CALL ancmsg(msgid=78,
     .                    msgtype=msgerror,
     .                    anmode=aninfo,
     .                    c1='/INITEM',
     .                    i1=id,  
     .                    i2=nod)
            ENDIF
            IF (nodsys /= 0) temp(nodsys)= list_temp0(j)
          ENDDO 
        ENDIF 
        IF(ALLOCATED(list_temp0)) DEALLOCATE(list_temp0) 
        IF(ALLOCATED(list_nod)) DEALLOCATE(list_nod) 
       ENDDO  
C---                                                           
      CALL udouble(initids,1,nbtemp,mess,0,bid)
C--------------------------------------------------
C     PRINT
C--------------------------------------------------
      j=0
      RETURN
C
2000  FORMAT(/, ' INITIAL TEMPERATURE ',/' -------------------',//
     + 6x,'NODE',17x,'TEMP ' )
      RETURN
      END