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 "param_c.inc"
51#include "scr03_c.inc"
52
53
54
55 INTEGER, INTENT(IN) :: IXTG(NIXTG, NUMELTG),IGEO(NPROPGI,NUMGEO)
56 my_real,
INTENT(IN) :: pm(npropm, nummat)
57 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
58
59
60
61 INTEGER :: IVPR(4), 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,numeltg
69 jthe_from_mat = nint(pm(71,iabs(ixtg(1,m))))
70 jale_from_mat = nint(pm(72,iabs(ixtg(1,m))))
71 jale_from_prop = igeo(62,iabs(ixtg(5,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,3
77 ivpr(i)=ale_connectivity%ee_connect%connected(iad1 + i - 1)
78 IF(ivpr(i) /= 0)ivpr(i)=ixtg(nixtg,ivpr(i))
79 ENDDO
80 WRITE(iout,'(5I10)')ixtg(nixtg,m),(ivpr(i),i=1,4)
81 ENDDO
82 ENDIF
83
84 RETURN
85
86 1200 FORMAT(//
87 & 6x,'A.L.E. ELEMENT CONNECTIVITIES ',/
88 & 6x,'------------------------------ ',/
89 & 6x,'ELEM ELE1 ELE2 ELE3 ELE4 ')