30
31
32
33
34
35
36
37
38
39
41 use element_mod , only : nixs
42
43
44
45#include "implicit_f.inc"
46
47
48
49#include "com04_c.inc"
50#include "units_c.inc"
51#include "scr03_c.inc"
52#include "param_c.inc"
53
54
55
56 INTEGER,INTENT(IN) :: IXS(NIXS,NUMELS),IGEO(NPROPGI,NUMGEO)
57 my_real,
INTENT(IN) :: pm(npropm,nummat), geo(npropg,numgeo)
58 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
59
60
61
62 INTEGER :: IVPR(6), M, IAL, I, IAD1, LGTH
63 INTEGER :: JTHE_FROM_MAT, JALE_FROM_MAT, JALE_FROM_PROP
64
65
66
67 IF(ipri >= 3)THEN
68 WRITE(iout,1200)
69 DO m=1,numels
70 jthe_from_mat = nint(pm(71,iabs(ixs(1,m))))
71 jale_from_mat = nint(pm(72,iabs(ixs(1,m))))
72 jale_from_prop = igeo(62,iabs(ixs(10,m)))
73 ial = jthe_from_mat + jale_from_mat + jale_from_prop
74 IF(ial == 0)cycle
75 iad1 = ale_connectivity%ee_connect%iad_connect(m)
76 lgth = ale_connectivity%ee_connect%iad_connect(m+1)-ale_connectivity%ee_connect%iad_connect(m)
77 DO i=1,lgth
78 ivpr(i)=ale_connectivity%ee_connect%connected(iad1 + i - 1)
79 IF(ivpr(i) /= 0)ivpr(i)=ixs(11,ivpr(i))
80 ENDDO
81 WRITE(iout,'(7I10)')ixs(11,m),(ivpr(i),i=1,6)
82 ENDDO
83 ENDIF
84
85 1200 FORMAT(//
86 & 6x,'A.L.E. ELEMENT CONNECTIVITIES ',/
87 & 6x,'------------------------------ ',/
88 & 6x,'ELEM ELE1 ELE2 ELE3 ELE4 ELE5',
89 & 6x,'ELE6')
90
91 RETURN