32
33
34
35
36
37
38 USE format_mod , ONLY : fmw_10i
39 use element_mod , only : nixs
40
41
42
43#include "implicit_f.inc"
44
45
46
47#include "units_c.inc"
48#include "com04_c.inc"
49#include "param_c.inc"
50#include "scr03_c.inc"
51#include "inter22.inc"
52
53
54
55 INTEGER , INTENT(IN) :: NBRIC,S_MSV
56 INTEGER , INTENT(IN) :: BRICS(NBRIC), IXS(NIXS,NUMELS), IGEO(NPROPGI,NUMGEO)
57 INTEGER , INTENT(INOUT) :: MSN, MSV(S_MSV)
58 my_real ,
INTENT(IN) :: pm(npropm,nummat)
59
60
61
62 INTEGER :: I, IMAT, IPID, IALE_, IF_ALE, IF_EUL, JALE_FROM_MAT, JALE_FROM_PROP
63
64
65
66 msn = nbric
67 i22_aleul = 0
68 if_ale = 0
69 if_eul = 0
70 DO i=1,nbric
71 msv(i) = brics(i)
72 END DO
73
74 IF(int22/=0)THEN
75 DO i=1,nbric
76 imat = ixs(1,msv(i))
77 ipid = ixs(10,msv(i))
78 jale_from_mat = nint(pm(72,imat))
79 jale_from_prop = igeo(62,ipid)
80 iale_ =
max(jale_from_mat, jale_from_prop)
81 IF(iale_==1)if_ale=1
82 IF(iale_==2)if_eul=1
83 ENDDO
84 IF(if_ale==1)THEN
85 i22_aleul = 1
86 ELSEIF(if_eul==1)THEN
87 i22_aleul = 2
88 ELSE
89 i22_aleul = 0
90 ENDIF
91 ENDIF
92
93 IF(ipri >= 1) THEN
94 WRITE(iout,'(/,A,/)') ' BRICKS USED FOR FLUID DOMAIN DEFINITION'
95 IF (nbric>0)THEN
96 WRITE(iout,fmt=fmw_10i)(ixs(11,msv(i)),i=1,msn)
97 ELSE
98 WRITE(iout,fmt='(A)') ' empty !'
99 ENDIF
100 ENDIF
101
102 RETURN