41
42
43
47 USE format_mod , ONLY : fmw_10i
48
49
50
51
52
53
54#include "implicit_f.inc"
55
56
57
58#include "com04_c.inc"
59#include "param_c.inc"
60#include "scr03_c.inc"
61#include "units_c.inc"
62
63
64
65 INTEGER NRT, NINT, NSN, NOINT, NTY, ID, ITIED
66 INTEGER IRECT(4,*), IXS(NIXS,*), IXC(NIXC,*), MWA(*),
67 . NSV(*), IXTG(NIXTG,*), ISEG(*), VARCONVINT,
68 . ITAB(*), FCOUNT(*), KNOD2ELS(*), NOD2ELS(*)
70 . x(3,*), pm(npropm,*)
71 CHARACTER(LEN=NCHARTITLE) :: TITR
72 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
73
74
75
76 INTEGER I, INRT, NELS, J, NN, NF, M, , IL, VARCONV,IAD
79
80
81
82
83 IF(itied/=3)THEN
84 CALL icinv3(ixs,ixc,ixtg,mwa,mwa(1+(2+numnod)))
85 IF(ipri>=1 .AND. nint<0)WRITE(iout,1000)
86 IF(ipri>=1 .AND. nint>0)WRITE(iout,1001)
87 DO i=1,nrt
88 inrt=i
89 CALL i12sol3(x,irect,ixs,nint,nels,inrt,
90 .
area,noint,mwa ,mwa(1+(2+numnod)),nf,itab,
91 . knod2els,nod2els,nty,
id,titr)
92 IF(nels/=0) THEN
93 m=ixs(1,nels)
94 varconv=int(pm(10,m))
95 IF(varconvint==-1)varconvint=varconv
96 IF(varconv==varconvint)THEN
97 segindx = segindx+1
98 iseg(i) = segindx
99 iad = ale_connectivity%ee_connect%iad_connect(nels)
100 ale_connectivity%ee_connect%connected(iad + nf - 1) = -segindx
101 IF(ipri>=1)
102 . WRITE(iout,fmt=fmw_10i)i,ixs(11,nels),nf,segindx,nels
103 ELSE
104 IF(nint>0) WRITE (iout,1101) i, noint
105 IF(nint<0) WRITE (iout,1201) i, noint
106 ierr=ierr+1
107 ENDIF
108 ELSE
109 IF(nint>0) WRITE (iout,1100) i, noint
110 IF(nint<0) WRITE (iout,1200) i, noint
111 ierr=ierr+1
112 ENDIF
113 END DO
115 . msgtype=msgerror,
116 . anmode=aninfo_blind,
118 . c1=titr,
119 . prmod=msg_print)
121 . msgtype=msgerror,
122 . anmode=aninfo_blind,
124 . c1=titr,
125 . prmod=msg_print)
126 END IF
127
128
129
130
131 DO i=1,nrt
132 DO j=1,4
133 ig=irect(j,i)
135 irect(j,i)=il
136 ENDDO
137 ENDDO
138
139
140
141 IF(nint>0)THEN
142 DO i=1,nsn
143 fcount(i)=0
144 ENDDO
145 DO i=1,nrt
146 DO j=1,4
147 nn=irect(j,i)
148 fcount(nn)=fcount(nn)+1
149 ENDDO
150 ENDDO
151 ENDIF
152
153 RETURN
154
155 1000 FORMAT(//' MAIN ELT FACE ')
156 1001 FORMAT(//' SECONDARY ELT FACE ')
157 1100 FORMAT(2x,'** ERROR ** CANNOT FIND FLUID ELT OF MAIN SEGMENT',i10,' OF INTERFACE',i10)
158 1200 FORMAT(2x,'** ERROR ** CANNOT FIND FLUID ELT OF SECONDARY SEGMENT',i10,' OF INTERFACE',i10)
159 1101 FORMAT(2x,'** ERROR ** VARIABLES NOT CONSISTENT IN SECONDARY ELEMENT',i10,' FOR INTERFACE',i10)
160 1201 FORMAT(2x,'** ERROR ** VARIABLES NOT CONSISTENT IN MAIN ELEMENT',i10,' FOR INTERFACE',i10)
subroutine area(d1, x, x2, y, y2, eint, stif0)
subroutine icinv3(ixs, ixc, ixtg, iadd, invc)
subroutine i12sol3(x, irect, ixs, nint, nel, i, area, noint, iadd, invc, nf, itab, knod2els, nod2els, nty, id, titr)
subroutine local_index(il, ig, nodes, n)
integer, parameter nchartitle
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)