40
41
42
43 USE elbufdef_mod
49 use element_mod , only : nixr
50
51
52
53#include "implicit_f.inc"
54
55
56
57#include "com01_c.inc"
58#include "com04_c.inc"
59#include "com10_c.inc"
60#include "param_c.inc"
61#include "tabsiz_c.inc"
62
63
64
65 INTEGER IPARG(,*),IGEO(NPROPGI,*),IXR(NIXR,*)
66 INTEGER, DIMENSION(NUMNOD), INTENT(in) :: WEIGHT
67
68 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP) :: ELBUF_TAB
69 TYPE (TH_), INTENT(IN),TARGET :: TH
70
71
72
73 LOGICAL :: BOOL
74 INTEGER :: ID,LOCAL_ID
75 INTEGER, POINTER :: ITHGR,ITHB
76 INTEGER, DIMENSION(10) :: ALREADY_DONE
77 INTEGER :: SITHBUF
78
79
80
81
82
83
84 already_done(1:10) = 0
85
87 IF(nthgrp01(
id)/=0)
THEN
89 ithgr => th%ITHGRPA(1)
90 ithb=> th%ITHBUFA(1)
91 sithbuf = th%SITHBUFA
92 ENDIF
94 ithgr => th%ITHGRPB(1)
95 ithb=> th%ITHBUFB(1)
96 sithbuf = th%SITHBUFB
97 ENDIF
99 ithgr => th%ITHGRPC(1)
100 ithb =>th%ITHBUFC(1)
101 sithbuf = th%SITHBUFC
102 ENDIF
104 ithgr => th%ITHGRPD(1)
105 ithb => th%ITHBUFD(1)
106 sithbuf = th%SITHBUFD
107 ENDIF
109 ithgr => th%ITHGRPE(1)
110 ithb => th%ITHBUFE(1)
111 sithbuf = th%SITHBUFE
112 ENDIF
114 ithgr => th%ITHGRPF(1)
115 ithb => th%ITHBUFF(1)
116 sithbuf = th%SITHBUFF
117 ENDIF
119 ithgr => th%ITHGRPG(1)
120 ithb => th%ITHBUFG(1)
121 sithbuf = th%SITHBUFG
122 ENDIF
124 ithgr => th%ITHGRPH(1)
125 ithb => th%ITHBUFH(1)
126 sithbuf = th%SITHBUFH
127 ENDIF
129 ithgr => th%ITHGRPI(1)
130 ithb => th%ITHBUFI(1)
131 sithbuf = th%SITHBUFI
132 ENDIF
134 CALL init_th(iparg,ithb,elbuf_tab,igeo,ixr,
135 . ithgr,nthgrp1(
id),
id,weight,sithbuf)
136 ENDIF
137 ENDDO
138
141 CALL init_th(iparg,th%ITHBUF,elbuf_tab,igeo,ixr,
142 . th%ITHGRP,nthgrp,
id,weight,th%SITHBUF)
143
145 bool=.false.
146 IF(abfile(
id)/=0)
THEN
147 IF(
id==2.OR.
id==4.OR.
id==5.OR.
id==6.OR.
id==7.OR.
id==8.OR.
id==9)
THEN
148 local_id = 1
149 IF(already_done(local_id)==0) THEN
150 ithgr => th%ITHGRPA(1)
151 ithb=> th%ITHBUFA(1)
152 bool=.true.
153 ENDIF
154
156 local_id = 2
157 IF(already_done(local_id)==0) THEN
158 ithgr => th%ITHGRPB(1)
159 ithb=> th%ITHBUFB(1)
160 bool=.true.
161 ENDIF
162 ENDIF
163 ENDIF
164 IF(bool.EQV..true.) THEN
165 already_done(local_id)=1
166 CALL init_th(iparg,ithb,elbuf_tab,igeo,ixr,
167 . ithgr,nthgrp1(
id),
id,weight,th%SITHBUF)
168 ENDIF
169 ENDDO
170
171
172 RETURN
subroutine init_th(iparg, ithbuf, elbuf_tab, igeo, ixr, ithgrp, nthgrp2, id, weight, sithbuf)