44
45
46
47
48
49
50
54 use element_mod , only :nixs,nixc,nixtg
55
56
57
58#include "implicit_f.inc"
59
60
61
62#include "param_c.inc"
63
64
65
66 INTEGER NRT, NINT, NSN, NOINT, NTY
67 INTEGER IRECT(4,*), IXS(NIXS,*), IXC(NIXC,*),
68 . NSV(*), IXTG(NIXTG,*),
69 . KNOD2ELS(*), KNOD2ELC(*), KNOD2ELTG(*), NOD2ELS(*), NOD2ELC(*),
70 . NOD2ELTG(*), IXS10(*), IXS16(*), IXS20(*),
71 . IGEO(*), IWORKSH(*)
73 . x(3,*),geo(npropg,*), pm(npropm,*),thk(*),pm_stack(*)
74 INTEGER ID
75 CHARACTER(LEN=NCHARTITLE) :: TITR
76 TYPE (SURF_) :: IGRSURF
77
78
79
80 INTEGER I, INRT, NELS, NELC, NELTG, J, IG, IL
82
83
84
85 DO i=1,nrt
86 inrt=i
87 CALL inelts(x ,irect,ixs ,nint,nels ,
88 . inrt ,
area ,noint,0 ,igrsurf%ELTYP,
89 . igrsurf%ELEM)
90 CALL ineltc(nelc ,neltg ,inrt ,igrsurf%ELTYP, igrsurf%ELEM)
91 IF(nels+nelc+neltg/=0) cycle
92 CALL insol3(x,irect,ixs,nint,nels,inrt,
93 .
area,noint,knod2els ,nod2els ,0 ,ixs10,
94 . ixs16,ixs20)
95 CALL incoq3(irect,ixc ,ixtg ,nint ,nelc ,
96 . neltg,inrt,geo ,pm ,knod2elc ,
97 . knod2eltg ,nod2elc ,nod2eltg,thk,nty,igeo,
98 . pm_stack , iworksh )
99 IF(nels+nelc+neltg==0)THEN
100 IF(nint>0) THEN
102 . msgtype=msgwarning,
103 . anmode=aninfo_blind_2,
105 . c1=titr,
106 . i2=i)
107 ENDIF
108 IF(nint<0) THEN
110 . msgtype=msgwarning,
111 . anmode=aninfo_blind_2,
113 . c1=titr,
114 . i2=i)
115 ENDIF
116 ENDIF
117 enddo
118
119 DO i=1,nrt
120 DO j=1,4
121 ig=irect(j,i)
123 irect(j,i)=il
124 ENDDO
125 ENDDO
126
127 RETURN
subroutine area(d1, x, x2, y, y2, eint, stif0)
subroutine incoq3(irect, ixc, ixtg, nint, nel, neltg, is, geo, pm, knod2elc, knod2eltg, nod2elc, nod2eltg, thk, nty, igeo, pm_stack, iworksh)
subroutine inelts(x, irect, ixs, nint, nel, i, area, noint, ir, surf_eltyp, surf_elem)
subroutine ineltc(nelc, neltg, is, surf_eltyp, surf_elem)
subroutine insol3(x, irect, ixs, nint, nel, i, area, noint, knod2els, nod2els, ir, ixs10, ixs16, ixs20)
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)