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