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