hm__read__radiation_8F_source.html

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_radiation      ../starter/source/loads/thermic/hm_read_radiation.F

!||--- called by ------------------------------------------------------

!||    lectur                 ../starter/source/starter/lectur.F

!||--- calls      -----------------------------------------------------

!||    ancmsg                 ../starter/source/output/message/message.F

!||    hm_get_floatv          ../starter/source/devtools/hm_reader/hm_get_floatv.F

!||    hm_get_floatv_dim      ../starter/source/devtools/hm_reader/hm_get_floatv_dim.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

!||    usr2sys                ../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_radiation(IB, FAC, ITAB, IXS ,IGRSURF, UNITAB ,LSUBMODEL,

     .                             NRADIA,NUMRADIA,NIRADIA,LFACTHER)

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

      use element_mod , only : nixs

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      "param_c.inc"

#include      "units_c.inc"

#include      "com04_c.inc"

C-----------------------------------------------

C   D u m m y   A r g u m e n t s

C-----------------------------------------------

      INTEGER ,INTENT(IN) :: NRADIA

      INTEGER ,INTENT(IN) :: NUMRADIA

      INTEGER ,INTENT(IN) :: NIRADIA

      INTEGER ,INTENT(IN) :: LFACTHER

      TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB

      INTEGER IB(NIRADIA,*), ITAB(*), IXS(NIXS,*)

      my_real fac(lfacther,*)


      TYPE (SURF_)   , DIMENSION(NSURF)   :: IGRSURF

      TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER K, M, I1, I2, I3, I4, IFU, I, ISENS,NPR0,NN,ISU,IS,

     .        id,j,uid,iflagunit,ity

      my_real

     .   fcx, fcy, fac_l, fac_t, fac_m, emi, sigma, tstart, tstop,

     .   emiss(numradia),fcx_dim,fcy_dim,tstop_dim

      CHARACTER MESS*40

      CHARACTER(LEN=NCHARTITLE) :: TITR

      LOGICAL IS_AVAILABLE

C-----------------------------------------------

C   E x t e r n a l   F u n c t i o n s

C-----------------------------------------------

      INTEGER USR2SYS

      DATA mess/'RADIATIVE FLUX DEFINITION               '/

C=======================================================================

      is_available = .false.

      k = 0

C--------------------------------------------------

C START BROWSING MODEL RADIATION

C--------------------------------------------------

      CALL hm_option_start('/RADIATION')

C--------------------------------------------------

C BROWSING /RADIATION OPTIONS 1->NRADIA

C--------------------------------------------------

      DO i=1,nradia

        titr = ''

        CALL hm_option_read_key(lsubmodel,

     .                       unit_id = uid,

     .                       option_id = id,

     .                       option_titr = titr)

        iflagunit = 0

        DO j=1,unitab%NUNITS

          IF (unitab%UNIT_ID(j) == uid) THEN

            iflagunit = 1

            EXIT

          ENDIF

        ENDDO

        IF (uid /= 0.AND.iflagunit == 0) THEN

          CALL ancmsg(msgid=659,anmode=aninfo,msgtype=msgerror,

     .                i2=uid,i1=id,c1='CONVECTION HEAT',

     .                 c2='CONVECTION HEAT',

     .                 c3=titr)

        ENDIF

C--------------------------------------------------

C EXTRACT DATAS (INTEGER VALUES)

C--------------------------------------------------

        CALL hm_get_intv('entityid',isu,is_available,lsubmodel)

        CALL hm_get_intv('curveid',ifu,is_available,lsubmodel)

        CALL hm_get_intv('rad_sensor_id',isens,is_available,lsubmodel)

C--------------------------------------------------

C EXTRACT DATAS (REAL VALUES)

C--------------------------------------------------

        CALL hm_get_floatv('xscale',fcx,is_available,lsubmodel,unitab)

        CALL hm_get_floatv_dim('xscale',fcx_dim,is_available,lsubmodel,unitab)

        CALL hm_get_floatv('magnitude',fcy,is_available,lsubmodel,unitab)

        CALL hm_get_floatv_dim('magnitude',fcy_dim,is_available,lsubmodel,unitab)

        CALL hm_get_floatv('rad_tstart',tstart,is_available,lsubmodel,unitab)

        CALL hm_get_floatv('rad_tstop',tstop,is_available,lsubmodel,unitab)

        CALL hm_get_floatv_dim('rad_tstop',tstop_dim,is_available,lsubmodel,unitab)

        CALL hm_get_floatv('flux',emi,is_available,lsubmodel,unitab)

C--------------------------------------------------

        IF (fcx == zero) fcx = fcx_dim

        IF (fcy == zero) fcy = fcy_dim

        IF(tstop == zero) tstop= ep30 * tstop_dim

C

        fac_m = unitab%FAC_M_WORK

        fac_t = unitab%FAC_T_WORK

        sigma=stefboltz*(fac_t*fac_t*fac_t)/fac_m

C

        is=0

        DO j=1,nsurf

          IF (isu == igrsurf(j)%ID) is=j

        ENDDO

        IF(is > 0)THEN

          nn=igrsurf(is)%NSEG

          DO j=1,nn

            k=k+1

            ib(1,k)=igrsurf(is)%NODES(j,1)

            ib(2,k)=igrsurf(is)%NODES(j,2)

            ib(3,k)=igrsurf(is)%NODES(j,3)

            ity    =igrsurf(is)%ELTYP(j)

            IF(ity==7)THEN

C          true triangles (not segments built from 3 nodes).

             ib(4,k)=0

            ELSE

             ib(4,k)=igrsurf(is)%NODES(j,4)

            ENDIF

            ib(5,k)  = ifu

            ib(6,k)  = isens

            ib(7,k)  = igrsurf(is)%ELTYP(j)

            ib(8,k)  = igrsurf(is)%ELEM(j)

            IF(ity == 1) THEN

              ib(9,k)  = ixs(11,igrsurf(is)%ELEM(j))

            ELSE

              ib(9,k)  = 0

            ENDIF

C

            fac(1,k) = fcy

            fac(2,k) = one/fcx

            fac(3,k) = emi*sigma

            fac(4,k) = tstart

            fac(5,k) = tstop

            fac(6,k) = one

C

C temporary storage for print out.

            emiss(k)=emi

          ENDDO

        ENDIF

      ENDDO

C

      i1=1

      i2=min0(50,numradia)

C

   90  WRITE (iout,2000)

       WRITE (iout,2001)

       DO i=i1,i2

         WRITE (iout,'(5(1X,I10),1X,1F10.3,2(1X,I10),1X,4G20.13)') i,

     .         itab(ib(1,i)),itab(ib(2,i)),itab(ib(3,i)),itab(ib(4,i)),

     .         emiss(i),ib(5,i),ib(6,i),fac(4,i),fac(5,i),one/fac(2,i),

     .         fac(1,i)

       ENDDO

       IF(i2 == numradia)GOTO 200

       i1=i1+50

       i2=min0(i2+50,numradia)

       GOTO 90

 200  RETURN

 300  CALL ancmsg(msgid=157,

     .            msgtype=msgerror,

     .            anmode=aninfo,

     .            i1=k)

C---

 2000 FORMAT(//

     .'     RADIATION   HEAT  '/

     .'     ----------------  ')

 2001 FORMAT(/

     .'    SEGMENT      NODE1      NODE2      NODE3      NODE4  EMISSIVITY',

     .'     CURVE     SENSOR          T-START              T-STOP', 8x,

     .'     SCALE-X             SCALE-Y')

C-----------

      RETURN


      END


my_real
#define my_real
Definition cppsort.cpp:32

hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418

hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:161

hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528

hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:253

hm_read_radiation
subroutine hm_read_radiation(ib, fac, itab, ixs, igrsurf, unitab, lsubmodel, nradia, numradia, niradia, lfacther)
Definition hm_read_radiation.F:42

hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:414

id
initmumps id
Definition multiplerhs_example.m:4

groupdef_mod
Definition groupdef_mod.F:668

hm_option_read_mod
Definition hm_option_read_mod.F:412

message_mod
Definition message_mod.F:1257

names_and_titles_mod
Definition names_and_titles_mod.F:1002

names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008

submodel_mod
Definition submodel_mod.F:717

unitab_mod
Definition unitab_mod.F:517

ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895

submodel_mod::submodel_data
Definition submodel_mod.F:728