OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
w_elbuf_str.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| w_elbuf_str ../starter/source/restart/ddsplit/w_elbuf_str.F
25!||--- called by ------------------------------------------------------
26!|| ddsplit ../starter/source/restart/ddsplit/ddsplit.F
27!||--- calls -----------------------------------------------------
28!||--- uses -----------------------------------------------------
29!||====================================================================
30 SUBROUTINE w_elbuf_str(PROC,IPARG,ELBUF_TAB,LEN_AM,FLAG_XFEM)
31C-----------------------------------------------
32C M o d u l e s
33C-----------------------------------------------
34 USE elbufdef_mod
35C-----------------------------------------------
36C I m p l i c i t T y p e s
37C-----------------------------------------------
38#include "implicit_f.inc"
39C-----------------------------------------------
40C C o m m o n B l o c k s
41C-----------------------------------------------
42#include "com01_c.inc"
43#include "param_c.inc"
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 INTEGER PROC,LEN_AM, IPARG(NPARG,*),FLAG_XFEM
48 TYPE(elbuf_struct_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
49C-----------------------------------------------
50C L o c a l V a r i a b l e s
51C-----------------------------------------------
52 INTEGER I,J,K,IL,IR,IS,IT,NG,NEL,NPT,P,L_REL,BUFLEN,NUVAR,
53 . ITY,IGTYP,ERR,IFAIL,IDFAIL,DEBUG_PRINT,IAD,
54 . NFAIL,IXFEM,NPG,NPTT,NPTTOT,INLOC,NONL,
55 . IDAMP_FREQ_RANGE
56 INTEGER ! global variables
57 . G_GAMA,G_SIG,G_OFF,G_NOFF,G_EINT,G_EINS,G_TEMP,
58 . G_RHO,G_PLA,G_VOL,G_EPSD,G_QVIS,G_DELTAX,G_TB,G_RK,G_RE,G_WPLA,
59 . G_SMSTR,G_HOURG,G_BFRAC,G_THK,G_FOR,G_MOM,G_TAG22,G_STRA,G_SIGI,
60 . G_DMG,G_FORPG,G_MOMPG,G_STRPG,G_UELR,G_UELR1,G_EPSQ,G_IERR,
61 . G_DAMDL,G_FORTH,G_EINTTH,G_FILL,G_SEQ,G_STRW,G_STRWPG,G_THK_I,
62 . G_JAC_I,G_AREA,G_SKEW,G_LENGTH,G_TOTDEPL,G_TOTROT,G_FOREP,G_MOMEP,
63 . G_DEP_IN_TENS,G_DEP_IN_COMP,G_ROT_IN_TENS,G_ROT_IN_COMP,
64 . G_POSX,G_POSY,G_POSZ,G_POSXX,G_POSYY,G_POSZZ,G_YIELD,G_LENGTH_ERR,
65 . G_DV,G_DFS,G_SKEW_ERR,G_E6,G_RUPTCRIT,G_MASS,G_V_REPCVT,
66 . G_VR_REPCVT,G_NUVAR,G_NUVARN,G_INIFRIC,G_DT,G_ISMS,G_STRHG,G_ETOTSH,
67 . G_BPRELD,G_ABURN,G_MU,G_PLANL,G_EPSDNL,G_MAXFRAC,G_MAXEPS,G_DT_PITER,
68 . G_TEMPG,G_COR_NF,G_COR_FR,G_COR_XR,G_DEFINI,G_FORINI,G_IDT_TSH,
69 . G_FORPGPINCH,G_MOMPGPINCH,G_EPGPINCHXZ,G_EPGPINCHYZ,G_EPGPINCHZZ,
70 . G_SKEW_ID,G_GAMA_R,G_SLIPRING_ID,G_SLIPRING_STRAND,G_RETRACTOR_ID,G_RINGSLIP,
71 . G_ADD_NODE,G_UPDATE,G_INTVAR,G_BETAORTH,G_AMU,G_SLIPRING_FRAM_ID,G_FRAM_FACTOR,
72 . G_TM_YIELD,G_TM_SEQ,G_TM_EINT,G_TM_DMG,G_TM_SIG,G_TM_STRA,G_TSAIWU,G_DMGSCL,
73 . G_SH_IOFFSET,G_EINT_DISTOR,G_FOR_G,G_FORPG_G,G_YIELD_IN_COMP,G_XXOLD_IN_COMP
74
75 INTEGER ! layer variables
76 . LY_DMG,LY_GAMA,LY_DIRA,LY_DIRB,LY_CRKDIR,LY_PLAPT,LY_SIGPT,
77 . LY_HOURG,LY_UELR,LY_UELR1,Ly_OFFPG,LY_OFF,LEN_PLAPT,LEN_SIGPT
78 INTEGER ! interlayer variables
79 . LY_EINT,LY_COUNT
80 INTEGER ! local variables
81 . L_OFF,L_EINT,L_EINS,L_RHO,L_DP_DRHO,L_VOL,L_PLA,L_TB,L_TEMP,
82 . L_RK,L_RE,L_VK,L_ROB,L_SIG,L_SIGV,L_SIGL,L_SIGA,L_SIGD,L_SIGB,
83 . L_SIGC,L_SIGF,L_STRA,L_EPSA,L_EPSD,L_EPSQ,L_EPSF,L_CRAK,L_WPLA,
84 . L_ANG,L_SF,L_QVIS,L_DELTAX,L_GAMA,
85 . L_DAM,L_DSUM,L_DGLO,L_EPE,L_EPC,L_XST,L_SSP,L_Z,L_FRAC,L_VISC,
86 . L_THK,L_FOR,L_MOM,L_EPS,L_SMSTR,L_BFRAC,L_DMG,L_FORTH,L_EINTTH,
87 . L_SEQ,L_SIGPLY,L_JAC_I,L_FAC_YLD,L_ABURN,L_MU,L_PLANL,L_EPSDNL,
88 . L_DMGSCL,L_PIJ,L_VOL0DP,L_TSAIWU,IDRAPE
89 INTEGER ! failure variables
90 . LF_DAM,LF_DAMMX,LF_DAMINI,LF_TDEL,LF_INDX,LF_OFF
91 INTEGER sum
92 my_real,
93 . ALLOCATABLE, DIMENSION (:) :: rbuf_l
94
95 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: DP_RBUF
96 INTEGER SDP_RBUF,L_DPREL
97
98 TYPE(elbuf_struct_), POINTER :: ELBUF
99 TYPE(g_bufel_) , POINTER :: GBUF
100 TYPE(l_bufel_) , POINTER :: LBUF
101 TYPE(buf_intloc_) , POINTER :: ILBUF
102 TYPE(buf_lay_) , POINTER :: BUFLY
103 TYPE(buf_intlay_) , POINTER :: INTLAY
104 TYPE(fail_loc_) , POINTER :: FLOC
105 TYPE(buf_nloc_) , POINTER :: BUFNL
106 TYPE(buf_nlocts_) , POINTER :: BUFNLTS
107 TYPE(l_bufel_dir_) , POINTER :: LBUF_DIR
108
109C=======================================================================
110 debug_print = 0
111 npg = -huge(npg)
112c
113 DO ng=1,ngroup
114 p = iparg(32,ng)
115 inloc = iparg(78,ng) ! Flag for non-local regularization
116 idamp_freq_range = iparg(93,ng) ! Flag for damping frequency range
117c----
118 IF (p == proc) THEN
119c----
120 ity = iparg(5,ng)
121
122 IF (ity == 1 .OR. ity == 2 .OR. ity == 51 .OR.
123 . ity == 3 .OR. ity == 7 .OR. ity == 101.OR.
124 . ity == 4 .OR. ity == 5 .OR. ity == 6 .OR. ity == 100) THEN
125 elbuf => elbuf_tab(ng)
126c
127 ixfem = iparg(54,ng) ! flag XFEM
128 IF (ixfem == 0 .AND. flag_xfem == 1) cycle ! no "XFEM_TAB" buffer for non xfem PARTS
129c
130 igtyp = elbuf%IGTYP
131 ixfem = elbuf%IXFEM
132 nel = elbuf%NEL
133 idrape = elbuf%IDRAPE
134cc NPT = IPARG(6,NG)
135 npt = elbuf%NLAY*elbuf%NPTR*elbuf%NPTS*elbuf%NPTT
136C
137c---------------------------------------------------
138 IF (ity == 1 .OR. ity == 2 .OR. ity == 51 .OR.
139 . ity == 101) THEN
140 npg = elbuf%NPTR*elbuf%NPTS*elbuf%NPTT
141 npttot = npt
142 ELSEIF (ity == 3 .OR. ity == 7) THEN
143 npg = elbuf%NPTR*elbuf%NPTS
144 npttot = 0
145 DO il=1,elbuf%NLAY
146 npttot = npttot + npg*elbuf%BUFLY(il)%NPTT
147 ENDDO
148 IF (npt == 0) npttot = npt
149 ELSE
150 npg = 0
151 npttot = 0
152 ENDIF
153c---------------------------------------------------
154c allocate and fill the transfer table for restart
155c---------------------------------------------------
156 buflen = iparg(66,ng)
157 ALLOCATE (rbuf_l(buflen), stat=err)
158c
159c---------------------------------------------------
160C IF (DEBUG_PRINT == 1) THEN
161C print*,'**** GROUP,ITY,IGTYP,BUFLEN=',ng,ity,igtyp,BUFLEN
162C ENDIF
163c---------------------------------------------------
164C
165C SMSTR must be written in separate Buffer (Forced in double precision)
166C
167 gbuf => elbuf_tab(ng)%GBUF
168C-- add VOL0DP
169 sdp_rbuf = 0
170 DO il = 1,elbuf%NLAY
171 bufly => elbuf_tab(ng)%BUFLY(il)
172 IF (igtyp == 51 .OR. igtyp == 52) THEN
173 nptt = elbuf_tab(ng)%BUFLY(il)%NPTT
174 ELSE
175 nptt = elbuf_tab(ng)%NPTT
176 ENDIF
177 DO ir = 1,elbuf%NPTR
178 DO is = 1,elbuf%NPTS
179 DO it = 1,nptt
180 sdp_rbuf = sdp_rbuf+bufly%L_VOL0DP*nel
181 ENDDO
182 ENDDO
183 ENDDO
184 ENDDO
185 sdp_rbuf = sdp_rbuf + gbuf%G_SMSTR *nel
186 IF(sdp_rbuf > 0)THEN
187 ALLOCATE(dp_rbuf(sdp_rbuf))
188 ENDIF
189 l_dprel = 0
190C
191
192 rbuf_l(1:buflen) = zero
193 l_rel = 0
194c--------- write top structure parameters: NPAR_TAB = 9
195 rbuf_l(l_rel+1) = elbuf%IGTYP
196 l_rel = l_rel+1
197 rbuf_l(l_rel+1) = elbuf%NEL
198 l_rel = l_rel+1
199 rbuf_l(l_rel+1) = elbuf%NLAY
200 l_rel = l_rel+1
201 rbuf_l(l_rel+1) = elbuf%NINTLAY
202 l_rel = l_rel+1
203 rbuf_l(l_rel+1) = elbuf%NPTR
204 l_rel = l_rel+1
205 rbuf_l(l_rel+1) = elbuf%NPTS
206 l_rel = l_rel+1
207 rbuf_l(l_rel+1) = elbuf%NPTT
208 l_rel = l_rel+1
209 rbuf_l(l_rel+1) = elbuf%IXFEM
210 l_rel = l_rel+1
211 rbuf_l(l_rel+1) = elbuf%NXEL
212 l_rel = l_rel+1
213 rbuf_l(l_rel+1) = elbuf%IDRAPE
214 l_rel = l_rel+1
215c---------------------------------------------------
216c Damping parameters for frequency range
217c---------------------------------------------------
218 IF (idamp_freq_range > 0) THEN
219 rbuf_l(l_rel+1:l_rel+3) = elbuf%DAMP_RANGE%ALPHA(1:3)
220 rbuf_l(l_rel+4:l_rel+6) = elbuf%DAMP_RANGE%TAU(1:3)
221 l_rel = l_rel + 6
222 ENDIF
223c---------------------------------------------------
224c LAYERS : write buffer at layer level
225c---------------------------------------------------
226 sum = l_rel
227
228 IF (debug_print == 1) THEN
229 print*,' NPAR_TAB=',sum
230 ENDIF
231c---------------------------------------------------
232 DO il = 1,elbuf%NLAY
233 bufly => elbuf_tab(ng)%BUFLY(il)
234c
235c--------- write layer parameters: NPAR_LAY = 13
236c
237 rbuf_l(l_rel+1) = bufly%ILAW
238 l_rel = l_rel+1
239 rbuf_l(l_rel+1) = bufly%IMAT
240 l_rel = l_rel+1
241 rbuf_l(l_rel+1) = bufly%IEOS
242 l_rel = l_rel+1
243 rbuf_l(l_rel+1) = bufly%IVISC
244 l_rel = l_rel+1
245 rbuf_l(l_rel+1) = bufly%IPORO
246 l_rel = l_rel+1
247 rbuf_l(l_rel+1) = bufly%NFAIL
248 l_rel = l_rel+1
249 rbuf_l(l_rel+1) = bufly%NVAR_MAT
250 l_rel = l_rel+1
251 rbuf_l(l_rel+1) = bufly%NVAR_EOS
252 l_rel = l_rel+1
253 rbuf_l(l_rel+1) = bufly%NVARTMP
254 l_rel = l_rel+1
255 rbuf_l(l_rel+1) = bufly%NVARTMP_EOS
256 l_rel = l_rel+1
257 rbuf_l(l_rel+1) = bufly%NVAR_VISC
258 l_rel = l_rel+1
259 rbuf_l(l_rel+1) = bufly%NVAR_LOC
260 l_rel = l_rel+1
261 rbuf_l(l_rel+1) = bufly%NVAR_LAY
262 l_rel = l_rel+1
263 rbuf_l(l_rel+1) = bufly%NPTT
264 l_rel = l_rel+1
265c
266c---------------------------------------------------
267 IF (debug_print == 1) THEN
268 print*,' ** LAYER=',il
269 print*,' NPAR_LAY=', l_rel-sum
270 iad = l_rel
271 ENDIF
272c---------------------------------------------------
273c
274c------- write layer variable lengths (BUFLY%NVAR_LAY = 12)
275c
276 ly_dmg = bufly%LY_DMG
277 ly_gama = bufly%LY_GAMA
278 ly_dira = bufly%LY_DIRA
279 ly_dirb = bufly%LY_DIRB
280 ly_crkdir= bufly%LY_CRKDIR
281 ly_plapt = bufly%LY_PLAPT
282 ly_sigpt = bufly%LY_SIGPT
283 ly_hourg = bufly%LY_HOURG
284 ly_uelr = bufly%LY_UELR
285 ly_uelr1 = bufly%LY_UELR1
286 ly_offpg = bufly%LY_OFFPG
287 ly_off = bufly%LY_OFF
288c
289 rbuf_l(l_rel+1) = ly_dmg
290 l_rel = l_rel+1
291 rbuf_l(l_rel+1) = ly_gama
292 l_rel = l_rel+1
293 rbuf_l(l_rel+1) = ly_dira
294 l_rel = l_rel+1
295 rbuf_l(l_rel+1) = ly_dirb
296 l_rel = l_rel+1
297 rbuf_l(l_rel+1) = ly_crkdir
298 l_rel = l_rel+1
299 rbuf_l(l_rel+1) = ly_plapt
300 l_rel = l_rel+1
301 rbuf_l(l_rel+1) = ly_sigpt
302 l_rel = l_rel+1
303 rbuf_l(l_rel+1) = ly_hourg
304 l_rel = l_rel+1
305 rbuf_l(l_rel+1) = ly_uelr
306 l_rel = l_rel+1
307 rbuf_l(l_rel+1) = ly_uelr1
308 l_rel = l_rel+1
309 rbuf_l(l_rel+1) = ly_offpg
310 l_rel = l_rel+1
311 rbuf_l(l_rel+1) = ly_off
312 l_rel = l_rel+1
313c
314c---------------------------------------------------
315 IF (debug_print == 1) THEN
316 print*,' NVAR_LAY=', l_rel-iad
317 iad = l_rel
318 ENDIF
319c---------------------------------------------------
320c------- write local variable lengths
321c
322 l_off = bufly%L_OFF
323 l_gama = bufly%L_GAMA
324 l_stra = bufly%L_STRA
325 l_frac = bufly%L_FRAC
326 l_bfrac = bufly%L_BFRAC
327 l_eint = bufly%L_EINT
328 l_eins = bufly%L_EINS
329 l_rho = bufly%L_RHO
330 l_dp_drho = bufly%L_DP_DRHO
331 l_qvis = bufly%L_QVIS
332 l_deltax = bufly%L_DELTAX
333 l_vol = bufly%L_VOL
334 l_epsa = bufly%L_EPSA
335 l_epsd = bufly%L_EPSD
336 l_epsq = bufly%L_EPSQ
337 l_epsf = bufly%L_EPSF
338 l_pla = bufly%L_PLA
339 l_wpla = bufly%L_WPLA
340 l_temp = bufly%L_TEMP
341 l_tb = bufly%L_TB
342 l_rk = bufly%L_RK
343 l_re = bufly%L_RE
344 l_vk = bufly%L_VK
345 l_sf = bufly%L_SF
346 l_rob = bufly%L_ROB
347 l_dam = bufly%L_DAM
348 l_dsum = bufly%L_DSUM
349 l_dglo = bufly%L_DGLO
350 l_crak = bufly%L_CRAK
351 l_ang = bufly%L_ANG
352 l_epe = bufly%L_EPE
353 l_epc = bufly%L_EPC
354 l_xst = bufly%L_XST
355 l_ssp = bufly%L_SSP
356 l_z = bufly%L_Z
357 l_visc = bufly%L_VISC
358 l_sigl = bufly%L_SIGL
359 l_sigv = bufly%L_SIGV
360 l_siga = bufly%L_SIGA
361 l_sigb = bufly%L_SIGB
362 l_sigc = bufly%L_SIGC
363 l_sigd = bufly%L_SIGD
364 l_sigf = bufly%L_SIGF
365 l_sig = bufly%L_SIG
366 l_sigply = bufly%L_SIGPLY
367 l_for = bufly%L_FOR
368 l_mom = bufly%L_MOM
369 l_thk = bufly%L_THK
370 l_smstr = bufly%L_SMSTR
371 l_dmg = bufly%L_DMG
372 l_forth = bufly%L_FORTH
373 l_eintth = bufly%L_EINTTH
374 l_seq = bufly%L_SEQ
375 l_jac_i = bufly%L_JAC_I
376 l_fac_yld = bufly%L_FAC_YLD
377 l_aburn = bufly%L_ABURN
378 l_mu = bufly%L_MU
379 l_planl = bufly%L_PLANL
380 l_epsdnl = bufly%L_EPSDNL
381 l_dmgscl = bufly%L_DMGSCL
382 l_tsaiwu = bufly%L_TSAIWU
383 l_pij = bufly%L_PIJ
384 l_vol0dp = bufly%L_VOL0DP
385
386
387c---
388 rbuf_l(l_rel+1) = l_off
389 l_rel = l_rel + 1
390 rbuf_l(l_rel+1) = l_gama
391 l_rel = l_rel + 1
392 rbuf_l(l_rel+1) = l_stra
393 l_rel = l_rel + 1
394 rbuf_l(l_rel+1) = l_frac
395 l_rel = l_rel + 1
396 rbuf_l(l_rel+1) = l_bfrac
397 l_rel = l_rel + 1
398 rbuf_l(l_rel+1) = l_eint
399 l_rel = l_rel + 1
400 rbuf_l(l_rel+1) = l_eins
401 l_rel = l_rel + 1
402 rbuf_l(l_rel+1) = l_rho
403 l_rel = l_rel + 1
404 rbuf_l(l_rel+1) = l_dp_drho
405 l_rel = l_rel + 1
406 rbuf_l(l_rel+1) = l_qvis
407 l_rel = l_rel + 1
408 rbuf_l(l_rel+1) = l_deltax
409 l_rel = l_rel + 1
410 rbuf_l(l_rel+1) = l_vol
411 l_rel = l_rel + 1
412 rbuf_l(l_rel+1) = l_epsa
413 l_rel = l_rel + 1
414 rbuf_l(l_rel+1) = l_epsd
415 l_rel = l_rel + 1
416 rbuf_l(l_rel+1) = l_epsq
417 l_rel = l_rel + 1
418 rbuf_l(l_rel+1) = l_epsf
419 l_rel = l_rel + 1
420 rbuf_l(l_rel+1) = l_pla
421 l_rel = l_rel + 1
422 rbuf_l(l_rel+1) = l_wpla
423 l_rel = l_rel + 1
424 rbuf_l(l_rel+1) = l_temp
425 l_rel = l_rel + 1
426 rbuf_l(l_rel+1) = l_tb
427 l_rel = l_rel + 1
428 rbuf_l(l_rel+1) = l_rk
429 l_rel = l_rel + 1
430 rbuf_l(l_rel+1) = l_re
431 l_rel = l_rel + 1
432 rbuf_l(l_rel+1) = l_vk
433 l_rel = l_rel + 1
434 rbuf_l(l_rel+1) = l_sf
435 l_rel = l_rel + 1
436 rbuf_l(l_rel+1) = l_rob
437 l_rel = l_rel + 1
438 rbuf_l(l_rel+1) = l_dam
439 l_rel = l_rel + 1
440 rbuf_l(l_rel+1) = l_dsum
441 l_rel = l_rel + 1
442 rbuf_l(l_rel+1) = l_dglo
443 l_rel = l_rel + 1
444 rbuf_l(l_rel+1) = l_crak
445 l_rel = l_rel + 1
446 rbuf_l(l_rel+1) = l_ang
447 l_rel = l_rel + 1
448 rbuf_l(l_rel+1) = l_epe
449 l_rel = l_rel + 1
450 rbuf_l(l_rel+1) = l_epc
451 l_rel = l_rel + 1
452 rbuf_l(l_rel+1) = l_xst
453 l_rel = l_rel + 1
454 rbuf_l(l_rel+1) = l_ssp
455 l_rel = l_rel + 1
456 rbuf_l(l_rel+1) = l_z
457 l_rel = l_rel + 1
458 rbuf_l(l_rel+1) = l_visc
459 l_rel = l_rel + 1
460 rbuf_l(l_rel+1) = l_sigl
461 l_rel = l_rel + 1
462 rbuf_l(l_rel+1) = l_sigv
463 l_rel = l_rel + 1
464 rbuf_l(l_rel+1) = l_siga
465 l_rel = l_rel + 1
466 rbuf_l(l_rel+1) = l_sigb
467 l_rel = l_rel + 1
468 rbuf_l(l_rel+1) = l_sigc
469 l_rel = l_rel + 1
470 rbuf_l(l_rel+1) = l_sigd
471 l_rel = l_rel + 1
472 rbuf_l(l_rel+1) = l_sigf
473 l_rel = l_rel + 1
474 rbuf_l(l_rel+1) = l_sig
475 l_rel = l_rel + 1
476 rbuf_l(l_rel+1) = l_sigply
477 l_rel = l_rel + 1
478 rbuf_l(l_rel+1) = l_for
479 l_rel = l_rel + 1
480 rbuf_l(l_rel+1) = l_mom
481 l_rel = l_rel + 1
482 rbuf_l(l_rel+1) = l_thk
483 l_rel = l_rel + 1
484 rbuf_l(l_rel+1) = l_smstr
485 l_rel = l_rel + 1
486 rbuf_l(l_rel+1) = l_dmg
487 l_rel = l_rel + 1
488 rbuf_l(l_rel+1) = l_forth
489 l_rel = l_rel + 1
490 rbuf_l(l_rel+1) = l_eintth
491 l_rel = l_rel + 1
492 rbuf_l(l_rel+1) = l_seq
493 l_rel = l_rel + 1
494 rbuf_l(l_rel+1) = l_jac_i
495 l_rel = l_rel + 1
496 rbuf_l(l_rel+1) = l_fac_yld
497 l_rel = l_rel + 1
498 rbuf_l(l_rel+1) = l_aburn
499 l_rel = l_rel + 1
500 rbuf_l(l_rel+1) = l_mu
501 l_rel = l_rel + 1
502 rbuf_l(l_rel+1) = l_planl
503 l_rel = l_rel + 1
504 rbuf_l(l_rel+1) = l_epsdnl
505 l_rel = l_rel + 1
506 rbuf_l(l_rel+1) = l_dmgscl
507 l_rel = l_rel + 1
508 rbuf_l(l_rel+1) = l_tsaiwu
509 l_rel = l_rel + 1
510 rbuf_l(l_rel+1) = l_pij
511 l_rel = l_rel + 1
512 rbuf_l(l_rel+1) = l_vol0dp
513 l_rel = l_rel + 1
514c---------------------------------------------------
515 IF (debug_print == 1) THEN
516 print*,' NVAR_LOC=', l_rel-iad
517 iad = l_rel
518 ENDIF
519c---------------------------------------------------
520c---------
521c ! write layer variables for each element
522c
523 rbuf_l(l_rel+1:l_rel+nel*ly_dmg) = bufly%DMG(1:nel*ly_dmg)
524 l_rel = l_rel+nel*ly_dmg
525 rbuf_l(l_rel+1:l_rel+nel*ly_gama) = bufly%GAMA(1:nel*ly_gama)
526 l_rel = l_rel+nel*ly_gama
527 IF(idrape == 0 .OR.(idrape > 0 .AND. igtyp == 17)) THEN
528 rbuf_l(l_rel+1:l_rel+nel*ly_dira) = bufly%DIRA(1:nel*ly_dira)
529 l_rel = l_rel+nel*ly_dira
530 rbuf_l(l_rel+1:l_rel+nel*ly_dirb) = bufly%DIRB(1:nel*ly_dirb)
531 l_rel = l_rel+nel*ly_dirb
532 ENDIF
533 IF(ly_crkdir > 0) rbuf_l(l_rel+1:l_rel+nel*ly_crkdir) = bufly%CRKDIR(1:nel*ly_crkdir)
534 l_rel = l_rel+nel*ly_crkdir
535
536 IF (npg > 1) THEN
537 IF (elbuf%NLAY > 1) THEN
538 len_plapt = nel*ly_plapt
539 len_sigpt = nel*ly_sigpt
540 ELSE
541 IF (igtyp /= 51 .AND. igtyp /= 52) THEN
542 len_plapt = nel*ly_plapt*npt
543 len_sigpt = nel*ly_sigpt*npt
544 ELSE
545 len_plapt = nel*ly_plapt*bufly%NPTT
546 len_sigpt = nel*ly_sigpt*bufly%NPTT
547 ENDIF
548 ENDIF
549 rbuf_l(l_rel+1:l_rel+len_plapt) = bufly%PLAPT(1:len_plapt)
550 l_rel = l_rel+len_plapt
551 rbuf_l(l_rel+1:l_rel+len_sigpt) = bufly%SIGPT(1:len_sigpt)
552 l_rel = l_rel+len_sigpt
553 ELSEIF (npg == 1) THEN ! nothing to write
554C
555C must point after => LBUF%PLA and LBUF%SIG
556C
557 ENDIF
558c
559 rbuf_l(l_rel+1:l_rel+nel*ly_hourg) = bufly%HOURG(1:nel*ly_hourg)
560 l_rel = l_rel+nel*ly_hourg
561 rbuf_l(l_rel+1:l_rel+nel*ly_uelr) = bufly%UELR(1:nel*ly_uelr)
562 l_rel = l_rel+nel*ly_uelr
563 rbuf_l(l_rel+1:l_rel+nel*ly_uelr1) = bufly%UELR1(1:nel*ly_uelr1)
564 l_rel = l_rel+nel*ly_uelr1
565 rbuf_l(l_rel+1:l_rel+nel*ly_offpg) = bufly%OFFPG(1:nel*ly_offpg)
566 l_rel = l_rel+nel*ly_offpg
567 rbuf_l(l_rel+1:l_rel+nel*ly_off ) = bufly%OFF(1:nel*ly_off)
568 l_rel = l_rel+nel*ly_off
569c
570c---------------------------------------------------
571 IF (debug_print == 1) THEN
572 print*,' LAYER_VARS=', l_rel-iad
573 print*,' total layer vars =', l_rel - sum
574 iad = l_rel
575 ENDIF
576c---------------------------------------------------
577 sum = l_rel
578c
579c---------
580c
581c ! write local layer variables for each element
582c
583
584 IF (igtyp == 51 .OR. igtyp == 52) THEN
585 nptt = elbuf_tab(ng)%BUFLY(il)%NPTT
586 ELSE
587 nptt = elbuf_tab(ng)%NPTT
588 ENDIF
589C
590 DO i = 1,elbuf%NPTR
591 DO j = 1,elbuf%NPTS
592cc DO K = 1,ELBUF%NPTT
593 DO k = 1,nptt
594 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(i,j,k)
595c
596 rbuf_l(l_rel+1) = lbuf%MLAW
597 l_rel = l_rel + 1
598 rbuf_l(l_rel+1) = lbuf%LawID
599 l_rel = l_rel + 1
600c
601 rbuf_l(l_rel+1:l_rel+nel*l_off) = lbuf%OFF(1:nel*l_off)
602 l_rel = l_rel+nel*l_off
603 rbuf_l(l_rel+1:l_rel+nel*l_gama) =lbuf%GAMA(1:nel*l_gama)
604 l_rel = l_rel+nel *l_gama
605 rbuf_l(l_rel+1:l_rel+nel*l_stra) = lbuf%STRA(1:nel*l_stra)
606 l_rel = l_rel+nel*l_stra
607 rbuf_l(l_rel+1:l_rel+nel*l_frac) = lbuf%FRAC(1:nel*l_frac)
608 l_rel = l_rel+nel*l_frac
609 rbuf_l(l_rel+1:l_rel+nel*l_bfrac)= lbuf%BFRAC(1:nel*l_bfrac)
610 l_rel = l_rel+nel*l_bfrac
611 rbuf_l(l_rel+1:l_rel+nel*l_eint)= lbuf%EINT(1:nel*l_eint)
612 l_rel = l_rel+nel*l_eint
613 rbuf_l(l_rel+1:l_rel+nel*l_eins)= lbuf%EINS(1:nel*l_eins)
614 l_rel = l_rel+nel*l_eins
615 rbuf_l(l_rel+1:l_rel+nel*l_rho) = lbuf%RHO(1:nel*l_rho)
616 l_rel = l_rel+nel*l_rho
617 rbuf_l(l_rel+1:l_rel+nel*l_dp_drho)=
618 . lbuf%DP_DRHO(1:nel*l_dp_drho)
619 l_rel = l_rel+nel*l_dp_drho
620 rbuf_l(l_rel+1:l_rel+nel*l_qvis)= lbuf%QVIS(1:nel*l_qvis)
621 l_rel = l_rel+nel*l_qvis
622 rbuf_l(l_rel+1:l_rel+nel*l_deltax)=
623 . lbuf%DELTAX(1:nel*l_deltax)
624 l_rel = l_rel+nel*l_deltax
625 rbuf_l(l_rel+1:l_rel+nel*l_vol) = lbuf%VOL(1:nel*l_vol)
626 l_rel = l_rel+nel*l_vol
627 rbuf_l(l_rel+1:l_rel+nel*l_epsa) = lbuf%EPSA(1:nel*l_epsa)
628 l_rel = l_rel+nel*l_epsa
629 rbuf_l(l_rel+1:l_rel+nel*l_epsd) = lbuf%EPSD(1:nel*l_epsd)
630 l_rel = l_rel+nel*l_epsd
631 rbuf_l(l_rel+1:l_rel+nel*l_epsq) = lbuf%EPSQ(1:nel*l_epsq)
632 l_rel = l_rel+nel*l_epsq
633 rbuf_l(l_rel+1:l_rel+nel*l_epsf) = lbuf%EPSF(1:nel*l_epsf)
634 l_rel = l_rel+nel*l_epsf
635 rbuf_l(l_rel+1:l_rel+nel*l_pla) = lbuf%PLA(1:nel*l_pla)
636 l_rel = l_rel+nel*l_pla
637 rbuf_l(l_rel+1:l_rel+nel*l_wpla) = lbuf%WPLA(1:nel*l_wpla)
638 l_rel = l_rel+nel*l_wpla
639 rbuf_l(l_rel+1:l_rel+nel*l_temp) = lbuf%TEMP(1:nel*l_temp)
640 l_rel = l_rel+nel*l_temp
641 rbuf_l(l_rel+1:l_rel+nel*l_tb) = lbuf%TB(1:nel*l_tb)
642 l_rel = l_rel+nel*l_tb
643 rbuf_l(l_rel+1:l_rel+nel*l_rk) = lbuf%RK(1:nel*l_rk)
644 l_rel = l_rel+nel*l_rk
645 rbuf_l(l_rel+1:l_rel+nel*l_re) = lbuf%RE(1:nel*l_re)
646 l_rel = l_rel+nel*l_re
647 rbuf_l(l_rel+1:l_rel+nel*l_vk) = lbuf%VK(1:nel*l_vk)
648 l_rel = l_rel+nel*l_vk
649 rbuf_l(l_rel+1:l_rel+nel*l_sf) = lbuf%SF(1:nel*l_sf)
650 l_rel = l_rel+nel*l_sf
651 rbuf_l(l_rel+1:l_rel+nel*l_rob) = lbuf%ROB(1:nel*l_rob)
652 l_rel = l_rel+nel*l_rob
653 rbuf_l(l_rel+1:l_rel+nel*l_dam) = lbuf%DAM(1:nel*l_dam)
654 l_rel = l_rel+nel*l_dam
655 rbuf_l(l_rel+1:l_rel+nel*l_dsum) = lbuf%DSUM(1:nel*l_dsum)
656 l_rel = l_rel+nel*l_dsum
657 rbuf_l(l_rel+1:l_rel+nel*l_dglo) = lbuf%DGLO(1:nel*l_dglo)
658 l_rel = l_rel+nel*l_dglo
659 rbuf_l(l_rel+1:l_rel+nel*l_crak) = lbuf%CRAK(1:nel*l_crak)
660 l_rel = l_rel+nel*l_crak
661 rbuf_l(l_rel+1:l_rel+nel*l_ang) = lbuf%ANG(1:nel*l_ang)
662 l_rel = l_rel+nel*l_ang
663 rbuf_l(l_rel+1:l_rel+nel*l_epe) = lbuf%EPE(1:nel*l_epe)
664 l_rel = l_rel+nel*l_epe
665 rbuf_l(l_rel+1:l_rel+nel*l_epc) = lbuf%EPC(1:nel*l_epc)
666 l_rel = l_rel+nel*l_epc
667 rbuf_l(l_rel+1:l_rel+nel*l_xst) = lbuf%XST(1:nel*l_xst)
668 l_rel = l_rel+nel*l_xst
669 rbuf_l(l_rel+1:l_rel+nel*l_ssp) = lbuf%SSP(1:nel*l_ssp)
670 l_rel = l_rel+nel*l_ssp
671 rbuf_l(l_rel+1:l_rel+nel*l_z) = lbuf%Z(1:nel*l_z)
672 l_rel = l_rel+nel*l_z
673 rbuf_l(l_rel+1:l_rel+nel*l_visc) = lbuf%VISC(1:nel*l_visc)
674 l_rel = l_rel+nel*l_visc
675 rbuf_l(l_rel+1:l_rel+nel*l_sigl) = lbuf%SIGL(1:nel*l_sigl)
676 l_rel = l_rel+nel*l_sigl
677 rbuf_l(l_rel+1:l_rel+nel*l_sigv) = lbuf%SIGV(1:nel*l_sigv)
678 l_rel = l_rel+nel*l_sigv
679 rbuf_l(l_rel+1:l_rel+nel*l_siga) = lbuf%SIGA(1:nel*l_siga)
680 l_rel = l_rel+nel*l_siga
681 rbuf_l(l_rel+1:l_rel+nel*l_sigb) = lbuf%SIGB(1:nel*l_sigb)
682 l_rel = l_rel+nel*l_sigb
683 rbuf_l(l_rel+1:l_rel+nel*l_sigc) = lbuf%SIGC(1:nel*l_sigc)
684 l_rel = l_rel+nel*l_sigc
685 rbuf_l(l_rel+1:l_rel+nel*l_sigd) = lbuf%SIGD(1:nel*l_sigd)
686 l_rel = l_rel+nel*l_sigd
687 rbuf_l(l_rel+1:l_rel+nel*l_sigf) = lbuf%SIGF(1:nel*l_sigf)
688 l_rel = l_rel+nel*l_sigf
689 rbuf_l(l_rel+1:l_rel+nel*l_sig) = lbuf%SIG(1:nel*l_sig)
690 l_rel = l_rel+nel*l_sig
691 rbuf_l(l_rel+1:l_rel+nel*l_sigply) = lbuf%SIGPLY(1:nel*l_sigply)
692 l_rel = l_rel+nel*l_sigply
693 rbuf_l(l_rel+1:l_rel+nel*l_for) = lbuf%FOR(1:nel*l_for)
694 l_rel = l_rel+nel*l_for
695 rbuf_l(l_rel+1:l_rel+nel*l_mom) = lbuf%MOM(1:nel*l_mom)
696 l_rel = l_rel+nel*l_mom
697 rbuf_l(l_rel+1:l_rel+nel*l_thk) = lbuf%THK(1:nel*l_thk)
698 l_rel = l_rel+nel*l_thk
699 rbuf_l(l_rel+1:l_rel+nel*l_smstr) = lbuf%SMSTR(1:nel*l_smstr)
700 l_rel = l_rel+nel*l_smstr
701 rbuf_l(l_rel+1:l_rel+nel*l_dmg) = lbuf%DMG(1:nel*l_dmg)
702 l_rel = l_rel+nel*l_dmg
703 rbuf_l(l_rel+1:l_rel+nel*l_forth) = lbuf%FORTH(1:nel*l_forth)
704 l_rel = l_rel+nel*l_forth
705 rbuf_l(l_rel+1:l_rel+nel*l_eintth) = lbuf%EINTTH(1:nel*l_eintth)
706 l_rel = l_rel+nel*l_eintth
707 rbuf_l(l_rel+1:l_rel+nel*l_seq) = lbuf%SEQ(1:nel*l_seq)
708 l_rel = l_rel+nel*l_seq
709 rbuf_l(l_rel+1:l_rel+nel*l_jac_i) = lbuf%JAC_I(1:nel*l_jac_i)
710 l_rel = l_rel+nel*l_jac_i
711 rbuf_l(l_rel+1:l_rel+nel*l_fac_yld) = lbuf%FAC_YLD(1:nel*l_fac_yld)
712 l_rel = l_rel+nel*l_fac_yld
713 rbuf_l(l_rel+1:l_rel+nel*l_aburn)= lbuf%ABURN(1:nel*l_aburn)
714 l_rel = l_rel+nel*l_aburn
715 rbuf_l(l_rel+1:l_rel+nel*l_mu)= lbuf%MU(1:nel*l_mu)
716 l_rel = l_rel+nel*l_mu
717 rbuf_l(l_rel+1:l_rel+nel*l_planl)= lbuf%PLANL(1:nel*l_planl)
718 l_rel = l_rel+nel*l_planl
719 rbuf_l(l_rel+1:l_rel+nel*l_epsdnl)= lbuf%EPSDNL(1:nel*l_epsdnl)
720 l_rel = l_rel+nel*l_epsdnl
721 rbuf_l(l_rel+1:l_rel+nel*l_dmgscl)= lbuf%DMGSCL(1:nel*l_dmgscl)
722 l_rel = l_rel+nel*l_dmgscl
723 rbuf_l(l_rel+1:l_rel+nel*l_tsaiwu)= lbuf%TSAIWU(1:nel*l_tsaiwu)
724 l_rel = l_rel+nel*l_tsaiwu
725 rbuf_l(l_rel+1:l_rel+nel*l_pij) = lbuf%PIJ(1:nel*l_pij)
726 l_rel = l_rel+nel*l_pij
727C--------write first VOL0DP
728 IF(nel*l_vol0dp >0)THEN
729 dp_rbuf(l_dprel+1:l_dprel+nel*l_vol0dp)= lbuf%VOL0DP(1:nel*l_vol0dp)
730 l_dprel = l_dprel + nel*l_vol0dp
731 ENDIF
732 ENDDO
733 ENDDO
734 ENDDO
735 IF(idrape > 0 .AND. (igtyp == 51 .OR. igtyp == 52)) THEN
736 DO k = 1,nptt
737 lbuf_dir => elbuf_tab(ng)%BUFLY(il)%LBUF_DIR(k)
738 rbuf_l(l_rel+1:l_rel+nel*ly_dira) = lbuf_dir%DIRA(1:nel*ly_dira)
739 l_rel = l_rel + nel*ly_dira
740 rbuf_l(l_rel+1:l_rel+nel*ly_dirb) = lbuf_dir%DIRB(1:nel*ly_dirb)
741 l_rel = l_rel + nel*ly_dirb
742 ENDDO
743 ENDIF
744 ENDDO ! ELBUF%NLAY
745c---------------------------------------------------
746 IF (debug_print == 1) THEN
747 print*,' LOCAL VARS (LBUF TABLES)=', l_rel-iad
748 iad = l_rel
749 ENDIF
750c---------------------------------------------------
751c
752c----- Layer Material buffer - UVAR
753c
754 DO il = 1,elbuf%NLAY
755 bufly => elbuf_tab(ng)%BUFLY(il)
756 nuvar = bufly%NVAR_MAT
757 IF (igtyp == 51 .OR. igtyp == 52) THEN
758 nptt = elbuf_tab(ng)%BUFLY(il)%NPTT
759 ELSE
760 nptt = elbuf_tab(ng)%NPTT
761 ENDIF
762 DO ir = 1,elbuf%NPTR
763 DO is = 1,elbuf%NPTS
764cc DO IT = 1,ELBUF%NPTT
765 DO it = 1,nptt
766 rbuf_l(l_rel+1:l_rel+nel*nuvar)=bufly%MAT(ir,is,it)%VAR(1:nel*nuvar)
767 l_rel = l_rel+nel*nuvar
768 ENDDO
769 ENDDO
770 ENDDO
771 ENDDO
772c---------------------------------------------------
773 IF (debug_print == 1) THEN
774 print*,' MAT UVAR=', l_rel-iad
775 iad = l_rel
776 ENDIF
777c---------------------------------------------------
778c
779c----- Layer Fail buffer - UVAR
780c
781 DO il = 1,elbuf%NLAY
782 bufly => elbuf_tab(ng)%BUFLY(il)
783 nfail = bufly%NFAIL
784 IF (igtyp == 51 .OR. igtyp == 52) THEN
785 nptt = elbuf_tab(ng)%BUFLY(il)%NPTT
786 ELSE
787 nptt = elbuf_tab(ng)%NPTT
788 ENDIF
789 DO ir = 1,elbuf%NPTR
790 DO is = 1,elbuf%NPTS
791 DO it = 1,nptt
792 DO k = 1,nfail
793 floc=>elbuf_tab(ng)%BUFLY(il)%FAIL(ir,is,it)%FLOC(k)
794 ifail = floc%ILAWF
795 idfail = floc%IDFAIL
796 nuvar = floc%NVAR
797 lf_dam= floc%LF_DAM
798 lf_dammx= floc%LF_DAMMX
799 lf_damini= floc%LF_DAMINI
800 lf_tdel= floc%LF_TDEL
801 lf_indx= floc%LF_INDX
802 lf_off = floc%LF_OFF
803c
804 rbuf_l(l_rel+1)=ifail
805 l_rel = l_rel+1
806 rbuf_l(l_rel+1)=idfail
807 l_rel = l_rel+1
808 rbuf_l(l_rel+1)=nuvar
809 l_rel = l_rel+1
810 rbuf_l(l_rel+1)=lf_dam
811 l_rel = l_rel+1
812 rbuf_l(l_rel+1)=lf_dammx
813 l_rel = l_rel+1
814 rbuf_l(l_rel+1)=lf_damini
815 l_rel = l_rel+1
816 rbuf_l(l_rel+1)=lf_tdel
817 l_rel = l_rel+1
818 rbuf_l(l_rel+1)=lf_indx
819 l_rel = l_rel+1
820 rbuf_l(l_rel+1)=lf_off
821 l_rel = l_rel+1
822c
823 rbuf_l(l_rel+1:l_rel+nel*nuvar) = floc%VAR(1:nel*nuvar)
824 l_rel = l_rel+nel*nuvar
825 rbuf_l(l_rel+1:l_rel+nel*lf_dam) = floc%DAM(1:nel*lf_dam)
826 l_rel = l_rel+nel*lf_dam
827 rbuf_l(l_rel+1:l_rel+nel*lf_dammx) = floc%DAMMX(1:nel*lf_dammx)
828 l_rel = l_rel+nel*lf_dammx
829 rbuf_l(l_rel+1:l_rel+nel*lf_damini) = floc%DAMINI(1:nel*lf_damini)
830 l_rel = l_rel+nel*lf_damini
831 rbuf_l(l_rel+1:l_rel+nel*lf_tdel) = floc%TDEL(1:nel*lf_tdel)
832 l_rel = l_rel + nel*lf_tdel
833 rbuf_l(l_rel+1:l_rel+nel*lf_indx) = floc%INDX(1:nel*lf_indx)
834 l_rel = l_rel + nel*lf_indx
835 rbuf_l(l_rel+1:l_rel+nel*lf_off) = floc%OFF(1:nel*lf_off)
836 l_rel = l_rel + nel*lf_off
837 ENDDO
838 ENDDO
839 ENDDO
840 ENDDO
841 ENDDO
842c---------------------------------------------------
843 IF (debug_print == 1) THEN
844 print*,' FAIL UVAR=', l_rel-iad
845 iad = l_rel
846 ENDIF
847c---------------------------------------------------
848c
849c----- Layer Visc buffer - UVAR
850c
851 DO il = 1,elbuf%NLAY
852 bufly => elbuf_tab(ng)%BUFLY(il)
853 nuvar = bufly%NVAR_VISC
854 IF(nuvar>0) THEN
855 IF (igtyp == 51 .OR. igtyp == 52) THEN
856 nptt = elbuf_tab(ng)%BUFLY(il)%NPTT
857 ELSE
858 nptt = elbuf_tab(ng)%NPTT
859 ENDIF
860 DO ir = 1,elbuf%NPTR
861 DO is = 1,elbuf%NPTS
862 DO it = 1,nptt
863 rbuf_l(l_rel+1:l_rel+nel*nuvar)=bufly%VISC(ir,is,it)%VAR(1:nel*nuvar)
864 l_rel = l_rel+nel*nuvar
865 ENDDO
866 ENDDO
867 ENDDO
868 ENDIF
869 ENDDO
870c---------------------------------------------------
871 IF (debug_print == 1) THEN
872 print*,' VISC UVAR=', l_rel-iad
873 iad = l_rel
874 ENDIF
875c---------------------------------------------------
876c
877 ! Non-local regularization in the shell thickness
878 IF ((ity==3).OR.(ity==7)) THEN
879 IF ((inloc>0).AND.(elbuf%NPTT>1)) THEN
880 ! Number of non-local nodes in the thickness
881 IF (elbuf%NPTT>2) THEN
882 nonl = elbuf%NPTT+1
883 ELSE
884 nonl = elbuf%NPTT
885 ENDIF
886 ! Copying the buffer
887 DO ir = 1,elbuf%NPTR
888 DO is = 1,elbuf%NPTS
889 ! Pointing the non-local buffer
890 bufnl => elbuf_tab(ng)%NLOC(ir,is)
891 DO i = 1,nel
892 ! Nodal masses
893 rbuf_l(l_rel+1:l_rel+nonl) = bufnl%MASSTH(i,1:nonl)
894 l_rel = l_rel + nonl
895 ! Non-local in-thickness cumulated variables
896 rbuf_l(l_rel+1:l_rel+nonl) = bufnl%UNLTH(i,1:nonl)
897 l_rel = l_rel + nonl
898 ! Non-local in-thickness velocities
899 rbuf_l(l_rel+1:l_rel+nonl) = bufnl%VNLTH(i,1:nonl)
900 l_rel = l_rel + nonl
901 ! Non-local in-thickness forces
902 rbuf_l(l_rel+1:l_rel+nonl) = bufnl%FNLTH(i,1:nonl)
903 l_rel = l_rel + nonl
904 ENDDO
905 ENDDO
906 ENDDO
907 ENDIF
908 ELSEIF ((ity == 1).AND.(elbuf_tab(ng)%NLAY>1)) THEN
909 IF (inloc>0) THEN
910 ! Number of non-local nodes in the thickshell thickness
911 nonl = elbuf%NLAY+1
912 ! Copying the buffer
913 DO ir = 1,elbuf%NPTR
914 DO is = 1,elbuf%NPTS
915 ! Pointing the non-local buffer
916 bufnlts => elbuf_tab(ng)%NLOCTS(ir,is)
917 DO i = 1,nel
918 ! Nodal masses
919 rbuf_l(l_rel+1:l_rel+nonl) = bufnlts%MASSTH(i,1:nonl)
920 l_rel = l_rel + nonl
921 ! Non-local in-thickness cumulated variables
922 rbuf_l(l_rel+1:l_rel+nonl) = bufnlts%UNLTH(i,1:nonl)
923 l_rel = l_rel + nonl
924 ! Non-local in-thickness velocities
925 rbuf_l(l_rel+1:l_rel+nonl) = bufnlts%VNLTH(i,1:nonl)
926 l_rel = l_rel + nonl
927 ! Non-local in-thickness forces
928 rbuf_l(l_rel+1:l_rel+nonl) = bufnlts%FNLTH(i,1:nonl)
929 l_rel = l_rel + nonl
930 ENDDO
931 ENDDO
932 ENDDO
933 ENDIF
934 ENDIF
935c---------------------------------------------------
936 IF (debug_print == 1) THEN
937 print*,' NLOC UVAR=', l_rel-iad
938 iad = l_rel
939 ENDIF
940c---------------------------------------------------
941
942c---------------------------------------------------
943c
944c----- EoS
945c
946 DO il = 1,elbuf%NLAY
947 bufly => elbuf_tab(ng)%BUFLY(il)
948 nuvar = bufly%NVAR_EOS
949 IF(ASSOCIATED(bufly%EOS))THEN
950 IF (igtyp == 51 .OR. igtyp == 52) THEN
951 nptt = elbuf_tab(ng)%BUFLY(il)%NPTT
952 ELSE
953 nptt = elbuf_tab(ng)%NPTT
954 ENDIF
955 DO ir = 1,elbuf%NPTR
956 DO is = 1,elbuf%NPTS
957 DO it = 1,nptt
958 rbuf_l(l_rel+1:l_rel+nel*nuvar)=bufly%EOS(ir,is,it)%VAR(1:nel*nuvar)
959 l_rel = l_rel+nel*nuvar
960 ENDDO
961 ENDDO
962 ENDDO
963 endif!IF(ASSOCIATED(BUFLY%EOS))
964 ENDDO
965c---------------------------------------------------
966 IF (debug_print == 1) THEN
967 print*,' EOS UVAR=', l_rel-iad
968 iad = l_rel
969 ENDIF
970c---------------------------------------------------
971 sum = l_rel
972c
973c---------------------------------------------------
974c INTER LAYERS (Ply Xfem only)
975c---------------------------------------------------
976 l_eps = 3
977 l_sig = 3
978 ly_eint = 1
979 ly_count= 1
980c
981 DO il = 1,elbuf%NINTLAY
982 intlay => elbuf_tab(ng)%INTLAY(il)
983c
984c--------- write inter layer parameters: = 5
985c
986 rbuf_l(l_rel+1) = intlay%ILAW
987 l_rel = l_rel+1
988 rbuf_l(l_rel+1) = intlay%IMAT
989 l_rel = l_rel+1
990 rbuf_l(l_rel+1) = intlay%NFAIL
991 l_rel = l_rel+1
992 rbuf_l(l_rel+1) = intlay%NVAR_MAT
993 l_rel = l_rel+1
994 rbuf_l(l_rel+1) = intlay%NVARTMP
995 l_rel = l_rel+1
996c
997c ! write inter layer variables for each layer
998 rbuf_l(l_rel+1:l_rel+nel*ly_eint) = intlay%EINT(1:nel*ly_eint)
999 l_rel = l_rel+nel*ly_eint
1000 rbuf_l(l_rel+1:l_rel+nel*ly_count) = intlay%COUNT(1:nel*ly_count)
1001 l_rel = l_rel+nel*ly_count
1002
1003c ! write local inter layer variables for each element
1004 DO i = 1,elbuf%NPTR
1005 DO j = 1,elbuf%NPTS
1006 ilbuf => elbuf_tab(ng)%INTLAY(il)%ILBUF(i,j)
1007c
1008 rbuf_l(l_rel+1:l_rel+nel*l_eps) =ilbuf%EPS(1:nel*l_eps)
1009 l_rel = l_rel+nel*l_eps
1010 rbuf_l(l_rel+1:l_rel+nel*l_sig) =ilbuf%SIG(1:nel*l_sig)
1011 l_rel = l_rel+nel*l_sig
1012 ENDDO
1013 ENDDO
1014 ENDDO
1015c
1016c----- Inter Layer Material buffer - UVAR
1017c
1018 DO il = 1,elbuf%NINTLAY
1019 intlay => elbuf_tab(ng)%INTLAY(il)
1020 nuvar = intlay%NVAR_MAT
1021 DO ir = 1,elbuf%NPTR
1022 DO is = 1,elbuf%NPTS
1023 rbuf_l(l_rel+1:l_rel+nel*nuvar)=intlay%MAT(ir,is)%VAR(1:nel*nuvar)
1024 l_rel = l_rel+nel*nuvar
1025 ENDDO
1026 ENDDO
1027 ENDDO
1028c
1029c----- Inter Layer Fail buffer - UVAR
1030c
1031 DO il = 1,elbuf%NINTLAY
1032 intlay => elbuf_tab(ng)%INTLAY(il)
1033 nfail = intlay%NFAIL
1034 DO ir = 1,elbuf%NPTR
1035 DO is = 1,elbuf%NPTS
1036 DO k = 1,nfail
1037 floc=>elbuf_tab(ng)%INTLAY(il)%FAIL(ir,is)%FLOC(k)
1038 ifail = floc%ILAWF
1039 idfail = floc%IDFAIL
1040 nuvar = floc%NVAR
1041 rbuf_l(l_rel+1)=ifail
1042 l_rel = l_rel+1
1043 rbuf_l(l_rel+1)=idfail
1044 l_rel = l_rel+1
1045 rbuf_l(l_rel+1)=nuvar
1046 l_rel = l_rel+1
1047 rbuf_l(l_rel+1:l_rel+nel*nuvar)=floc%VAR(1:nel*nuvar)
1048 l_rel = l_rel+nel*nuvar
1049 ENDDO
1050 ENDDO
1051 ENDDO
1052 ENDDO
1053c---------------------------------------------------
1054 IF (debug_print == 1) THEN
1055 print*,' total interlayer=', L_REL-SUM
1056 IAD = L_REL
1057 ENDIF
1058c---------------------------------------------------
1059c
1060c-----------------------------------------------------------------------
1061c GLOBAL BUFFER
1062c-----------------------------------------------------------------------
1063 GBUF => ELBUF_TAB(NG)%GBUF
1064c
1065 G_NOFF = GBUF%G_NOFF
1066 G_IERR = GBUF%G_IERR
1067 G_OFF = GBUF%G_OFF
1068 G_GAMA = GBUF%G_GAMA
1069 G_SMSTR = GBUF%G_SMSTR
1070 G_HOURG = GBUF%G_HOURG
1071 G_BFRAC = GBUF%G_BFRAC
1072 G_EINT = GBUF%G_EINT
1073 G_EINS = GBUF%G_EINS
1074 G_RHO = GBUF%G_RHO
1075 G_QVIS = GBUF%G_QVIS
1076 G_DELTAX= GBUF%G_DELTAX
1077 G_VOL = GBUF%G_VOL
1078 G_EPSD = GBUF%G_EPSD
1079 G_EPSQ = GBUF%G_EPSQ
1080 G_PLA = GBUF%G_PLA
1081 G_WPLA = GBUF%G_WPLA
1082 G_TEMP = GBUF%G_TEMP
1083 G_TB = GBUF%G_TB
1084 G_RK = GBUF%G_RK
1085 G_RE = GBUF%G_RE
1086 G_SIG = GBUF%G_SIG
1087 G_FOR = GBUF%G_FOR
1088 G_MOM = GBUF%G_MOM
1089 G_THK = GBUF%G_THK
1090 G_TAG22 = GBUF%G_TAG22
1091 G_STRA = GBUF%G_STRA
1092 G_SIGI = GBUF%G_SIGI
1093 G_DMG = GBUF%G_DMG
1094 G_FORPG = GBUF%G_FORPG
1095 G_MOMPG = GBUF%G_MOMPG
1096 G_MAXEPS= GBUF%G_MAXEPS
1097 G_AMU = GBUF%G_AMU
1098 G_MAXFRAC = GBUF%G_MAXFRAC
1099 G_GAMA_R = GBUF%G_GAMA_R
1100C
1101 G_FORPGPINCH = GBUF%G_FORPGPINCH
1102 G_MOMPGPINCH = GBUF%G_MOMPGPINCH
1103 G_EPGPINCHXZ = GBUF%G_EPGPINCHXZ
1104 G_EPGPINCHYZ = GBUF%G_EPGPINCHYZ
1105 G_EPGPINCHZZ = GBUF%G_EPGPINCHZZ
1106C
1107 G_STRPG = GBUF%G_STRPG
1108 G_UELR = GBUF%G_UELR
1109 G_UELR1 = GBUF%G_UELR1
1110 G_DAMDL = GBUF%G_DAMDL
1111 G_FORTH = GBUF%G_FORTH
1112 G_EINTTH = GBUF%G_EINTTH
1113 G_FILL = GBUF%G_FILL
1114 G_SEQ = GBUF%G_SEQ
1115 G_STRW = GBUF%G_STRW
1116 G_STRWPG= GBUF%G_STRWPG
1117 G_THK_I = GBUF%G_THK_I
1118 G_JAC_I = GBUF%G_JAC_I
1119 G_DT = GBUF%G_DT
1120 G_ISMS = GBUF%G_ISMS
1121 G_BPRELD= GBUF%G_BPRELD
1122 G_ABURN = GBUF%G_ABURN
1123 G_MU = GBUF%G_MU
1124 G_PLANL = GBUF%G_PLANL
1125 G_EPSDNL = GBUF%G_EPSDNL
1126 G_DMGSCL = GBUF%G_DMGSCL
1127 G_TSAIWU = GBUF%G_TSAIWU
1128 G_COR_NF = GBUF%G_COR_NF
1129 G_COR_FR = GBUF%G_COR_FR
1130 G_COR_XR = GBUF%G_COR_XR
1131 G_SH_IOFFSET = GBUF%G_SH_IOFFSET
1132 G_EINT_DISTOR = GBUF%G_EINT_DISTOR
1133 G_FOR_G = GBUF%G_FOR_G
1134 G_FORPG_G = GBUF%G_FORPG_G
1135
1136c---
1137c 1-D elements (TRUSS, BEAM, SPRING)
1138c---
1139 G_AREA = GBUF%G_AREA
1140 G_SKEW = GBUF%G_SKEW
1141 G_LENGTH = GBUF%G_LENGTH
1142 G_TOTDEPL = GBUF%G_TOTDEPL
1143 G_TOTROT = GBUF%G_TOTROT
1144 G_FOREP = GBUF%G_FOREP
1145 G_MOMEP = GBUF%G_MOMEP
1146 G_DEP_IN_TENS = GBUF%G_DEP_IN_TENS
1147 G_DEP_IN_COMP = GBUF%G_DEP_IN_COMP
1148 G_ROT_IN_TENS = GBUF%G_ROT_IN_TENS
1149 G_ROT_IN_COMP = GBUF%G_ROT_IN_COMP
1150 G_POSX = GBUF%G_POSX
1151 G_POSY = GBUF%G_POSY
1152 G_POSZ = GBUF%G_POSZ
1153 G_POSXX = GBUF%G_POSXX
1154 G_POSYY = GBUF%G_POSYY
1155 G_POSZZ = GBUF%G_POSZZ
1156 G_YIELD = GBUF%G_YIELD
1157 G_LENGTH_ERR = GBUF%G_LENGTH_ERR
1158 G_DV = GBUF%G_DV
1159 G_DFS = GBUF%G_DFS
1160 G_SKEW_ERR = GBUF%G_SKEW_ERR
1161 G_E6 = GBUF%G_E6
1162 G_RUPTCRIT = GBUF%G_RUPTCRIT
1163 G_MASS = GBUF%G_MASS
1164 G_V_REPCVT = GBUF%G_V_REPCVT
1165 G_VR_REPCVT = GBUF%G_VR_REPCVT
1166 G_NUVAR = GBUF%G_NUVAR
1167 G_NUVARN = GBUF%G_NUVARN
1168 G_DEFINI = GBUF%G_DEFINI
1169 G_FORINI = GBUF%G_FORINI
1170 G_INIFRIC = GBUF%G_INIFRIC
1171 G_STRHG = GBUF%G_STRHG
1172 G_ETOTSH = GBUF%G_ETOTSH
1173 G_DT_PITER = GBUF%G_DT_PITER
1174 G_TEMPG = GBUF%G_TEMPG
1175 G_SKEW_ID = GBUF%G_SKEW_ID
1176 G_IDT_TSH = GBUF%G_IDT_TSH
1177 G_YIELD_IN_COMP = GBUF%G_YIELD_IN_COMP
1178 G_XXOLD_IN_COMP = GBUF%G_XXOLD_IN_COMP
1179c---
1180c seatbelts elements
1181c---
1182 G_SLIPRING_ID = GBUF%G_SLIPRING_ID
1183 G_SLIPRING_FRAM_ID = GBUF%G_SLIPRING_FRAM_ID
1184 G_SLIPRING_STRAND = GBUF%G_SLIPRING_STRAND
1185 G_RETRACTOR_ID= GBUF%G_RETRACTOR_ID
1186 G_RINGSLIP = GBUF%G_RINGSLIP
1187 G_ADD_NODE = GBUF%G_ADD_NODE
1188 G_UPDATE = GBUF%G_UPDATE
1189 G_FRAM_FACTOR = GBUF%G_FRAM_FACTOR
1190 G_INTVAR = GBUF%G_INTVAR
1191C
1192 G_BETAORTH = GBUF%G_BETAORTH
1193C ---H3D/TMAX in BUFLEN
1194 G_TM_YIELD = GBUF%G_TM_YIELD
1195 G_TM_SEQ = GBUF%G_TM_SEQ
1196 G_TM_EINT = GBUF%G_TM_EINT
1197 G_TM_DMG = GBUF%G_TM_DMG
1198 G_TM_SIG = GBUF%G_TM_SIG
1199 G_TM_STRA = GBUF%G_TM_STRA
1200c
1201c--------------------
1202 sum = l_rel
1203c--------------------
1204c
1205 RBUF_L(L_REL+1) = GBUF%NVAR_GLOB
1206 L_REL = L_REL + 1
1207c
1208c ! write global variable lengths (GBUF%NVAR_GLOB)
1209c
1210 RBUF_L(L_REL+1) = G_NOFF
1211 L_REL = L_REL + 1
1212 RBUF_L(L_REL+1) = G_IERR
1213 L_REL = L_REL + 1
1214 RBUF_L(L_REL+1) = G_OFF
1215 L_REL = L_REL + 1
1216 RBUF_L(L_REL+1) = G_GAMA
1217 L_REL = L_REL + 1
1218 RBUF_L(L_REL+1) = G_SMSTR
1219 L_REL = L_REL + 1
1220 RBUF_L(L_REL+1) = G_HOURG
1221 L_REL = L_REL + 1
1222 RBUF_L(L_REL+1) = G_BFRAC
1223 L_REL = L_REL + 1
1224 RBUF_L(L_REL+1) = G_EINT
1225 L_REL = L_REL + 1
1226 RBUF_L(L_REL+1) = G_EINS
1227 L_REL = L_REL + 1
1228 RBUF_L(L_REL+1) = G_RHO
1229 L_REL = L_REL + 1
1230 RBUF_L(L_REL+1) = G_QVIS
1231 L_REL = L_REL + 1
1232 RBUF_L(L_REL+1) = G_DELTAX
1233 L_REL = L_REL + 1
1234 RBUF_L(L_REL+1) = G_VOL
1235 L_REL = L_REL + 1
1236 RBUF_L(L_REL+1) = G_EPSD
1237 L_REL = L_REL + 1
1238 RBUF_L(L_REL+1) = G_EPSQ
1239 L_REL = L_REL + 1
1240 RBUF_L(L_REL+1) = G_PLA
1241 L_REL = L_REL + 1
1242 RBUF_L(L_REL+1) = G_WPLA
1243 L_REL = L_REL + 1
1244 RBUF_L(L_REL+1) = G_TEMP
1245 L_REL = L_REL + 1
1246 RBUF_L(L_REL+1) = G_TB
1247 L_REL = L_REL + 1
1248 RBUF_L(L_REL+1) = G_RK
1249 L_REL = L_REL + 1
1250 RBUF_L(L_REL+1) = G_RE
1251 L_REL = L_REL + 1
1252 RBUF_L(L_REL+1) = G_SIG
1253 L_REL = L_REL + 1
1254 RBUF_L(L_REL+1) = G_FOR
1255 L_REL = L_REL + 1
1256 RBUF_L(L_REL+1) = G_MOM
1257 L_REL = L_REL + 1
1258 RBUF_L(L_REL+1) = G_THK
1259 L_REL = L_REL + 1
1260 RBUF_L(L_REL+1) = G_TAG22
1261 L_REL = L_REL + 1
1262 RBUF_L(L_REL+1) = G_STRA
1263 L_REL = L_REL + 1
1264 RBUF_L(L_REL+1) = G_SIGI
1265 L_REL = L_REL + 1
1266 RBUF_L(L_REL+1) = G_DMG
1267 L_REL = L_REL + 1
1268 RBUF_L(L_REL+1) = G_FORPG
1269 L_REL = L_REL + 1
1270 RBUF_L(L_REL+1) = G_MOMPG
1271 L_REL = L_REL + 1
1272C
1273 RBUF_L(L_REL+1) = G_FORPGPINCH
1274 L_REL = L_REL + 1
1275 RBUF_L(L_REL+1) = G_MOMPGPINCH
1276 L_REL = L_REL + 1
1277 RBUF_L(L_REL+1) = G_EPGPINCHXZ
1278 L_REL = L_REL + 1
1279 RBUF_L(L_REL+1) = G_EPGPINCHYZ
1280 L_REL = L_REL + 1
1281 RBUF_L(L_REL+1) = G_EPGPINCHZZ
1282 L_REL = L_REL + 1
1283C
1284 RBUF_L(L_REL+1) = G_STRPG
1285 L_REL = L_REL + 1
1286 RBUF_L(L_REL+1) = G_UELR
1287 L_REL = L_REL + 1
1288 RBUF_L(L_REL+1) = G_UELR1
1289 L_REL = L_REL + 1
1290 RBUF_L(L_REL+1) = G_DAMDL
1291 L_REL = L_REL + 1
1292 RBUF_L(L_REL+1) = G_FORTH
1293 L_REL = L_REL + 1
1294 RBUF_L(L_REL+1) = G_EINTTH
1295 L_REL = L_REL + 1
1296 RBUF_L(L_REL+1) = G_FILL
1297 L_REL = L_REL + 1
1298 RBUF_L(L_REL+1) = G_SEQ
1299 L_REL = L_REL + 1
1300 RBUF_L(L_REL+1) = G_STRW
1301 L_REL = L_REL + 1
1302 RBUF_L(L_REL+1) = G_STRWPG
1303 L_REL = L_REL + 1
1304 RBUF_L(L_REL+1) = G_THK_I
1305 L_REL = L_REL + 1
1306 RBUF_L(L_REL+1) = G_JAC_I
1307 L_REL = L_REL + 1
1308 RBUF_L(L_REL+1) = G_DT
1309 L_REL = L_REL + 1
1310 RBUF_L(L_REL+1) = G_ISMS
1311 L_REL = L_REL + 1
1312 RBUF_L(L_REL+1) = G_BPRELD
1313 L_REL = L_REL + 1
1314 RBUF_L(L_REL+1) = G_ABURN
1315 L_REL = L_REL + 1
1316 RBUF_L(L_REL+1) = G_MU
1317 L_REL = L_REL + 1
1318 RBUF_L(L_REL+1) = G_PLANL
1319 L_REL = L_REL + 1
1320 RBUF_L(L_REL+1) = G_EPSDNL
1321 L_REL = L_REL + 1
1322 RBUF_L(L_REL+1) = G_DMGSCL
1323 L_REL = L_REL + 1
1324 RBUF_L(L_REL+1) = G_TSAIWU
1325 L_REL = L_REL + 1
1326 RBUF_L(L_REL+1) = G_TEMPG
1327 L_REL = L_REL + 1
1328 RBUF_L(L_REL+1) = G_COR_NF
1329 L_REL = L_REL + 1
1330 RBUF_L(L_REL+1) = G_COR_FR
1331 L_REL = L_REL + 1
1332 RBUF_L(L_REL+1) = G_COR_XR
1333 L_REL = L_REL + 1
1334 RBUF_L(L_REL+1) = G_MAXFRAC
1335 L_REL = L_REL + 1
1336 RBUF_L(L_REL+1) = G_MAXEPS
1337 L_REL = L_REL + 1
1338C
1339 RBUF_L(L_REL+1) = G_BETAORTH
1340 L_REL = L_REL + 1
1341 RBUF_L(L_REL+1) = G_AMU
1342 L_REL = L_REL + 1
1343 RBUF_L(L_REL+1) = G_SH_IOFFSET
1344 L_REL = L_REL + 1
1345 RBUF_L(L_REL+1) = G_EINT_DISTOR
1346 L_REL = L_REL + 1
1347 RBUF_L(L_REL+1) = G_FOR_G
1348 L_REL = L_REL + 1
1349 RBUF_L(L_REL+1) = G_FORPG_G
1350 L_REL = L_REL + 1
1351c---
1352c 1-D elements (TRUSS, BEAM, SPRING)
1353c---
1354 RBUF_L(L_REL+1) = G_AREA
1355 L_REL = L_REL + 1
1356 RBUF_L(L_REL+1) = G_SKEW
1357 L_REL = L_REL + 1
1358 RBUF_L(L_REL+1) = G_LENGTH
1359 L_REL = L_REL + 1
1360 RBUF_L(L_REL+1) = G_TOTDEPL
1361 L_REL = L_REL + 1
1362 RBUF_L(L_REL+1) = G_TOTROT
1363 L_REL = L_REL + 1
1364 RBUF_L(L_REL+1) = G_FOREP
1365 L_REL = L_REL + 1
1366 RBUF_L(L_REL+1) = G_MOMEP
1367 L_REL = L_REL + 1
1368 RBUF_L(L_REL+1) = G_DEP_IN_TENS
1369 L_REL = L_REL + 1
1370 RBUF_L(L_REL+1) = G_DEP_IN_COMP
1371 L_REL = L_REL + 1
1372 RBUF_L(L_REL+1) = G_ROT_IN_TENS
1373 L_REL = L_REL + 1
1374 RBUF_L(L_REL+1) = G_ROT_IN_COMP
1375 L_REL = L_REL + 1
1376 RBUF_L(L_REL+1) = G_POSX
1377 L_REL = L_REL + 1
1378 RBUF_L(L_REL+1) = G_POSY
1379 L_REL = L_REL + 1
1380 RBUF_L(L_REL+1) = G_POSZ
1381 L_REL = L_REL + 1
1382 RBUF_L(L_REL+1) = G_POSXX
1383 L_REL = L_REL + 1
1384 RBUF_L(L_REL+1) = G_POSYY
1385 L_REL = L_REL + 1
1386 RBUF_L(L_REL+1) = G_POSZZ
1387 L_REL = L_REL + 1
1388 RBUF_L(L_REL+1) = G_YIELD
1389 L_REL = L_REL + 1
1390 RBUF_L(L_REL+1) = G_LENGTH_ERR
1391 L_REL = L_REL + 1
1392 RBUF_L(L_REL+1) = G_DV
1393 L_REL = L_REL + 1
1394 RBUF_L(L_REL+1) = G_DFS
1395 L_REL = L_REL + 1
1396 RBUF_L(L_REL+1) = G_SKEW_ERR
1397 L_REL = L_REL + 1
1398 RBUF_L(L_REL+1) = G_E6
1399 L_REL = L_REL + 1
1400 RBUF_L(L_REL+1) = G_RUPTCRIT
1401 L_REL = L_REL + 1
1402 RBUF_L(L_REL+1) = G_MASS
1403 L_REL = L_REL + 1
1404 RBUF_L(L_REL+1) = G_V_REPCVT
1405 L_REL = L_REL + 1
1406 RBUF_L(L_REL+1) = G_VR_REPCVT
1407 L_REL = L_REL + 1
1408 RBUF_L(L_REL+1) = G_NUVAR
1409 L_REL = L_REL + 1
1410 RBUF_L(L_REL+1) = G_NUVARN
1411 L_REL = L_REL + 1
1412 RBUF_L(L_REL+1) = G_DEFINI
1413 L_REL = L_REL + 1
1414 RBUF_L(L_REL+1) = G_FORINI
1415 L_REL = L_REL + 1
1416 RBUF_L(L_REL+1) = G_INIFRIC
1417 L_REL = L_REL + 1
1418 RBUF_L(L_REL+1) = G_STRHG
1419 L_REL = L_REL + 1
1420 RBUF_L(L_REL+1) = G_ETOTSH
1421 L_REL = L_REL + 1
1422 RBUF_L(L_REL+1) = G_DT_PITER
1423 L_REL = L_REL + 1
1424 RBUF_L(L_REL+1) = G_SKEW_ID
1425 L_REL = L_REL + 1
1426 RBUF_L(L_REL+1) = G_GAMA_R
1427 L_REL = L_REL + 1
1428 RBUF_L(L_REL+1) = G_YIELD_IN_COMP
1429 L_REL = L_REL + 1
1430 RBUF_L(L_REL+1) = G_XXOLD_IN_COMP
1431 L_REL = L_REL + 1
1432c---
1433c seatbelts elements
1434c---
1435 RBUF_L(L_REL+1) = G_SLIPRING_ID
1436 L_REL = L_REL + 1
1437 RBUF_L(L_REL+1) = G_SLIPRING_FRAM_ID
1438 L_REL = L_REL + 1
1439 RBUF_L(L_REL+1) = G_SLIPRING_STRAND
1440 L_REL = L_REL + 1
1441 RBUF_L(L_REL+1) = G_RETRACTOR_ID
1442 L_REL = L_REL + 1
1443 RBUF_L(L_REL+1) = G_RINGSLIP
1444 L_REL = L_REL + 1
1445 RBUF_L(L_REL+1) = G_ADD_NODE
1446 L_REL = L_REL + 1
1447 RBUF_L(L_REL+1) = G_UPDATE
1448 L_REL = L_REL + 1
1449 RBUF_L(L_REL+1) = G_FRAM_FACTOR
1450 L_REL = L_REL + 1
1451 RBUF_L(L_REL+1) = G_INTVAR
1452 L_REL = L_REL + 1
1453 RBUF_L(L_REL+1) = G_IDT_TSH
1454 L_REL = L_REL + 1
1455c---
1456c H3D/MAX
1457c---
1458 RBUF_L(L_REL+1) = G_TM_YIELD
1459 L_REL = L_REL + 1
1460 RBUF_L(L_REL+1) = G_TM_SEQ
1461 L_REL = L_REL + 1
1462 RBUF_L(L_REL+1) = G_TM_EINT
1463 L_REL = L_REL + 1
1464 RBUF_L(L_REL+1) = G_TM_DMG
1465 L_REL = L_REL + 1
1466 RBUF_L(L_REL+1) = G_TM_SIG
1467 L_REL = L_REL + 1
1468 RBUF_L(L_REL+1) = G_TM_STRA
1469 L_REL = L_REL + 1
1470c---------------------------------
1471 IF (DEBUG_PRINT == 1) THEN
1472 print*,' nvar_glob=',GBUF%NVAR_GLOB,L_REL-SUM
1473 IAD = L_REL
1474 ENDIF
1475c---------------------------------
1476c
1477 sum = l_rel
1478c
1479c ! write global variables for each element
1480c
1481 RBUF_L(L_REL+1:L_REL+NEL*G_NOFF) = GBUF%NOFF(1:NEL*G_NOFF)
1482 L_REL = L_REL + NEL*G_NOFF
1483 RBUF_L(L_REL+1:L_REL+NEL*G_IERR) = GBUF%IERR(1:NEL*G_IERR)
1484 L_REL = L_REL + NEL*G_IERR
1485 RBUF_L(L_REL+1:L_REL+NEL*G_GAMA) = GBUF%GAMA(1:NEL*G_GAMA)
1486 L_REL = L_REL + NEL*G_GAMA
1487 IF(NEL*G_SMSTR >0)THEN
1488 DP_RBUF(L_DPREL+1:L_DPREL+NEL*G_SMSTR)= GBUF%SMSTR(1:NEL*G_SMSTR)
1489 L_DPREL = L_DPREL + NEL*G_SMSTR
1490 ENDIF
1491 RBUF_L(L_REL+1:L_REL+NEL*G_HOURG)= GBUF%HOURG(1:NEL*G_HOURG)
1492 L_REL = L_REL + NEL*G_HOURG
1493 RBUF_L(L_REL+1:L_REL+NEL*G_TAG22)= GBUF%TAG22(1:NEL*G_TAG22)
1494 L_REL = L_REL+NEL*G_TAG22
1495 RBUF_L(L_REL+1:L_REL+NEL*G_STRA) = GBUF%STRA(1:NEL*G_STRA)
1496 L_REL = L_REL+NEL*G_STRA
1497 RBUF_L(L_REL+1:L_REL+NEL*G_SIGI) = GBUF%SIGI(1:NEL*G_SIGI)
1498 L_REL = L_REL+NEL*G_SIGI
1499 RBUF_L(L_REL+1:L_REL+NEL*G_DMG) = GBUF%DMG(1:NEL*G_DMG)
1500 L_REL = L_REL+NEL*G_DMG
1501 RBUF_L(L_REL+1:L_REL+NEL*G_UELR) = GBUF%UELR(1:NEL*G_UELR)
1502 L_REL = L_REL+NEL*G_UELR
1503 RBUF_L(L_REL+1:L_REL+NEL*G_UELR1) = GBUF%UELR1(1:NEL*G_UELR1)
1504 L_REL = L_REL+NEL*G_UELR1
1505 RBUF_L(L_REL+1:L_REL+NEL*G_DAMDL) = GBUF%DAMDL(1:NEL*G_DAMDL)
1506 L_REL = L_REL+NEL*G_DAMDL
1507 RBUF_L(L_REL+1:L_REL+NEL*G_FOR) = GBUF%FOR(1:NEL*G_FOR)
1508 L_REL = L_REL+NEL*G_FOR
1509 RBUF_L(L_REL+1:L_REL+NEL*G_MOM) = GBUF%MOM(1:NEL*G_MOM)
1510 L_REL = L_REL+NEL*G_MOM
1511 RBUF_L(L_REL+1:L_REL+NEL*G_THK) = GBUF%THK(1:NEL*G_THK)
1512 L_REL = L_REL+NEL*G_THK
1513 RBUF_L(L_REL+1:L_REL+NEL*G_STRW)=GBUF%STRW(1:NEL*G_STRW)
1514 L_REL = L_REL+NEL*G_STRW
1515 RBUF_L(L_REL+1:L_REL+NEL*G_THK_I) = GBUF%THK_I(1:NEL*G_THK_I)
1516 L_REL = L_REL+NEL*G_THK_I
1517 RBUF_L(L_REL+1:L_REL+NEL*G_JAC_I) = GBUF%JAC_I(1:NEL*G_JAC_I)
1518 L_REL = L_REL+NEL*G_JAC_I
1519 RBUF_L(L_REL+1:L_REL+NEL*G_DT) = GBUF%DT(1:NEL*G_DT)
1520 L_REL = L_REL+NEL*G_DT
1521 RBUF_L(L_REL+1:L_REL+NEL*G_ISMS) = GBUF%ISMS(1:NEL*G_ISMS)
1522 L_REL = L_REL+NEL*G_ISMS
1523 RBUF_L(L_REL+1:L_REL+NEL*G_BPRELD) = GBUF%BPRELD(1:NEL*G_BPRELD)
1524 L_REL = L_REL+NEL*G_BPRELD
1525 RBUF_L(L_REL+1:L_REL+NEL*G_COR_NF) = GBUF%COR_NF(1:NEL*G_COR_NF)
1526 L_REL = L_REL+NEL*G_COR_NF
1527 RBUF_L(L_REL+1:L_REL+NEL*G_COR_FR) = GBUF%COR_FR(1:NEL*G_COR_FR)
1528 L_REL = L_REL+NEL*G_COR_FR
1529 RBUF_L(L_REL+1:L_REL+NEL*G_COR_XR) = GBUF%COR_XR(1:NEL*G_COR_XR)
1530 L_REL = L_REL+NEL*G_COR_XR
1531 RBUF_L(L_REL+1:L_REL+NEL*G_MAXFRAC)= GBUF%MAXFRAC(1:NEL*G_MAXFRAC)
1532 L_REL = L_REL + NEL*G_MAXFRAC
1533 RBUF_L(L_REL+1:L_REL+NEL*G_MAXEPS)= GBUF%MAXEPS(1:NEL*G_MAXEPS)
1534 L_REL = L_REL + NEL*G_MAXEPS
1535C
1536 RBUF_L(L_REL+1:L_REL+NEL*G_BETAORTH) = GBUF%BETAORTH(1:NEL*G_BETAORTH)
1537 L_REL = L_REL+NEL*G_BETAORTH
1538 RBUF_L(L_REL+1:L_REL+NEL*G_AMU) = GBUF%AMU(1:NEL*G_AMU)
1539 L_REL = L_REL + NEL*G_AMU
1540 RBUF_L(L_REL+1:L_REL+NEL*G_SH_IOFFSET) = GBUF%SH_IOFFSET(1:NEL*G_SH_IOFFSET)
1541 L_REL = L_REL + NEL*G_SH_IOFFSET
1542 RBUF_L(L_REL+1:L_REL+NEL*G_EINT_DISTOR) = GBUF%EINT_DISTOR(1:NEL*G_EINT_DISTOR)
1543 L_REL = L_REL + NEL*G_EINT_DISTOR
1544 RBUF_L(L_REL+1:L_REL+NEL*G_FOR_G) = GBUF%FOR_G(1:NEL*G_FOR_G)
1545 L_REL = L_REL + NEL*G_FOR_G
1546 RBUF_L(L_REL+1:L_REL+NEL*G_WPLA) = GBUF%WPLA(1:NEL*G_WPLA)
1547 L_REL = L_REL+NEL*G_WPLA
1548c---------------------------------------------------
1549 IF (DEBUG_PRINT == 1) THEN
1550 print*,' COMMON glob=', L_REL-SUM
1551 IAD = L_REL
1552 ENDIF
1553c---------------------------------------------------
1554c
1555 IF (NPG == 1) THEN
1556 IF (G_STRPG>G_STRA) THEN
1557 RBUF_L(L_REL+1:L_REL+NEL*G_STRPG) = GBUF%STRPG(1:NEL*G_STRPG)
1558 L_REL = L_REL+NEL*G_STRPG
1559 END IF
1560 ELSEIF (NPG > 1) THEN
1561 RBUF_L(L_REL+1:L_REL+NEL*G_FORPG) = GBUF%FORPG(1:NEL*G_FORPG)
1562 L_REL = L_REL+NEL*G_FORPG
1563 RBUF_L(L_REL+1:L_REL+NEL*G_MOMPG) = GBUF%MOMPG(1:NEL*G_MOMPG)
1564 L_REL = L_REL+NEL*G_MOMPG
1565C
1566 RBUF_L(L_REL+1:L_REL+NEL*G_FORPGPINCH) = GBUF%FORPGPINCH(1:NEL*G_FORPGPINCH)
1567 L_REL = L_REL+NEL*G_FORPGPINCH
1568 RBUF_L(L_REL+1:L_REL+NEL*G_MOMPGPINCH) = GBUF%MOMPGPINCH(1:NEL*G_MOMPGPINCH)
1569 L_REL = L_REL+NEL*G_MOMPGPINCH
1570 RBUF_L(L_REL+1:L_REL+NEL*G_EPGPINCHXZ) = GBUF%EPGPINCHXZ(1:NEL*G_EPGPINCHXZ)
1571 L_REL = L_REL+NEL*G_EPGPINCHXZ
1572 RBUF_L(L_REL+1:L_REL+NEL*G_EPGPINCHYZ) = GBUF%EPGPINCHYZ(1:NEL*G_EPGPINCHYZ)
1573 L_REL = L_REL+NEL*G_EPGPINCHYZ
1574 RBUF_L(L_REL+1:L_REL+NEL*G_EPGPINCHZZ) = GBUF%EPGPINCHZZ(1:NEL*G_EPGPINCHZZ)
1575 L_REL = L_REL+NEL*G_EPGPINCHZZ
1576C
1577 RBUF_L(L_REL+1:L_REL+NEL*G_STRPG) = GBUF%STRPG(1:NEL*G_STRPG)
1578 L_REL = L_REL+NEL*G_STRPG
1579 RBUF_L(L_REL+1:L_REL+NEL*G_STRWPG)= GBUF%STRWPG(1:NEL*G_STRWPG)
1580 L_REL = L_REL+NEL*G_STRWPG
1581 RBUF_L(L_REL+1:L_REL+NEL*G_FORPG_G)= GBUF%FORPG_G(1:NEL*G_FORPG_G)
1582 L_REL = L_REL+NEL*G_FORPG_G
1583 ENDIF
1584c---------------------------------------------------
1585 IF (DEBUG_PRINT == 1) THEN
1586 print*,' glob npg > 1 =', L_REL-IAD
1587 IAD = L_REL
1588 ENDIF
1589c---------------------------------------------------
1590c-------
1591 sum=L_REL
1592 IF (NPTTOT /= 1) THEN
1593 RBUF_L(L_REL+1:L_REL+NEL*G_BFRAC)= GBUF%BFRAC(1:NEL*G_BFRAC)
1594 L_REL = L_REL + NEL*G_BFRAC
1595 RBUF_L(L_REL+1:L_REL+NEL*G_OFF) = GBUF%OFF(1:NEL*G_OFF)
1596 L_REL = L_REL+NEL*G_OFF
1597 RBUF_L(L_REL+1:L_REL+NEL*G_EINT) = GBUF%EINT(1:NEL*G_EINT)
1598 L_REL = L_REL+NEL*G_EINT
1599 RBUF_L(L_REL+1:L_REL+NEL*G_EINS) = GBUF%EINS(1:NEL*G_EINS)
1600 L_REL = L_REL+NEL*G_EINS
1601 RBUF_L(L_REL+1:L_REL+NEL*G_RHO) = GBUF%RHO(1:NEL*G_RHO)
1602 L_REL = L_REL+NEL*G_RHO
1603 RBUF_L(L_REL+1:L_REL+NEL*G_QVIS) = GBUF%QVIS(1:NEL*G_QVIS)
1604 L_REL = L_REL+NEL*G_QVIS
1605 RBUF_L(L_REL+1:L_REL+NEL*G_DELTAX)=GBUF%DELTAX(1:NEL*G_DELTAX)
1606 L_REL = L_REL+NEL*G_DELTAX
1607 RBUF_L(L_REL+1:L_REL+NEL*G_VOL) = GBUF%VOL(1:NEL*G_VOL)
1608 L_REL = L_REL+NEL*G_VOL
1609 RBUF_L(L_REL+1:L_REL+NEL*G_EPSD) = GBUF%EPSD(1:NEL*G_EPSD)
1610 L_REL = L_REL+NEL*G_EPSD
1611 RBUF_L(L_REL+1:L_REL+NEL*G_EPSQ) = GBUF%EPSQ(1:NEL*G_EPSQ)
1612 L_REL = L_REL+NEL*G_EPSQ
1613 RBUF_L(L_REL+1:L_REL+NEL*G_PLA) = GBUF%PLA(1:NEL*G_PLA)
1614 L_REL = L_REL+NEL*G_PLA
1615 RBUF_L(L_REL+1:L_REL+NEL*G_TEMP) = GBUF%TEMP(1:NEL*G_TEMP)
1616 L_REL = L_REL+NEL*G_TEMP
1617 RBUF_L(L_REL+1:L_REL+NEL*G_TB) = GBUF%TB(1:NEL*G_TB)
1618 L_REL = L_REL+NEL*G_TB
1619 RBUF_L(L_REL+1:L_REL+NEL*G_RK) = GBUF%RK(1:NEL*G_RK)
1620 L_REL = L_REL+NEL*G_RK
1621 RBUF_L(L_REL+1:L_REL+NEL*G_RE) = GBUF%RE(1:NEL*G_RE)
1622 L_REL = L_REL+NEL*G_RE
1623 RBUF_L(L_REL+1:L_REL+NEL*G_SIG) = GBUF%SIG(1:NEL*G_SIG)
1624 L_REL = L_REL+NEL*G_SIG
1625 RBUF_L(L_REL+1:L_REL+NEL*G_FORTH)= GBUF%FORTH(1:NEL*G_FORTH)
1626 L_REL = L_REL+NEL*G_FORTH
1627 RBUF_L(L_REL+1:L_REL+NEL*G_EINTTH)= GBUF%EINTTH(1:NEL*G_EINTTH)
1628 L_REL = L_REL+NEL*G_EINTTH
1629 RBUF_L(L_REL+1:L_REL+NEL*G_SEQ) = GBUF%SEQ(1:NEL*G_SEQ)
1630 L_REL = L_REL+NEL*G_SEQ
1631 RBUF_L(L_REL+1:L_REL+NEL*G_ABURN)= GBUF%ABURN(1:NEL*G_ABURN)
1632 L_REL = L_REL + NEL*G_ABURN
1633 RBUF_L(L_REL+1:L_REL+NEL*G_MU)= GBUF%MU(1:NEL*G_MU)
1634 L_REL = L_REL + NEL*G_MU
1635 ENDIF
1636C---------
1637 RBUF_L(L_REL+1:L_REL+NEL*G_FILL) = GBUF%FILL(1:NEL*G_FILL)
1638 L_REL = L_REL+NEL*G_FILL
1639c---
1640c 1-D elements (TRUSS, BEAM, SPRING)
1641c---
1642.and. IF(IGTYP == 3 ITY == 5) THEN
1643 RBUF_L(L_REL+1:L_REL+NEL*G_DMGSCL) = GBUF%DMGSCL(1:NEL*G_DMGSCL)
1644 L_REL = L_REL+NEL*G_DMGSCL
1645 ENDIF
1646 RBUF_L(L_REL+1:L_REL+NEL*G_AREA) = GBUF%AREA(1:NEL*G_AREA)
1647 L_REL = L_REL +NEL*G_AREA
1648 RBUF_L(L_REL+1:L_REL+NEL*G_SKEW) = GBUF%SKEW(1:NEL*G_SKEW)
1649 L_REL = L_REL +NEL*G_SKEW
1650 RBUF_L(L_REL+1:L_REL+NEL*G_LENGTH) = GBUF%LENGTH(1:NEL*G_LENGTH)
1651 L_REL = L_REL +NEL*G_LENGTH
1652 RBUF_L(L_REL+1:L_REL+NEL*G_TOTDEPL) = GBUF%TOTDEPL(1:NEL*G_TOTDEPL)
1653 L_REL = L_REL +NEL*G_TOTDEPL
1654 RBUF_L(L_REL+1:L_REL+NEL*G_TOTROT) = GBUF%TOTROT(1:NEL*G_TOTROT)
1655 L_REL = L_REL +NEL*G_TOTROT
1656 RBUF_L(L_REL+1:L_REL+NEL*G_FOREP) = GBUF%FOREP(1:NEL*G_FOREP)
1657 L_REL = L_REL +NEL*G_FOREP
1658 RBUF_L(L_REL+1:L_REL+NEL*G_MOMEP) = GBUF%MOMEP(1:NEL*G_MOMEP)
1659 L_REL = L_REL +NEL*G_MOMEP
1660 RBUF_L(L_REL+1:L_REL+NEL*G_DEP_IN_TENS) = GBUF%DEP_IN_TENS(1:NEL*G_DEP_IN_TENS)
1661 L_REL = L_REL + NEL*G_DEP_IN_TENS
1662 RBUF_L(L_REL+1:L_REL+NEL*G_DEP_IN_COMP) = GBUF%DEP_IN_COMP(1:NEL*G_DEP_IN_COMP)
1663 L_REL = L_REL +NEL*G_DEP_IN_COMP
1664 RBUF_L(L_REL+1:L_REL+NEL*G_ROT_IN_TENS) = GBUF%ROT_IN_TENS(1:NEL*G_ROT_IN_TENS)
1665 L_REL = L_REL +NEL*G_ROT_IN_TENS
1666 RBUF_L(L_REL+1:L_REL+NEL*G_ROT_IN_COMP) = GBUF%ROT_IN_COMP(1:NEL*G_ROT_IN_COMP)
1667 L_REL = L_REL +NEL*G_ROT_IN_COMP
1668 RBUF_L(L_REL+1:L_REL+NEL*G_POSX) = GBUF%POSX(1:NEL*G_POSX)
1669 L_REL = L_REL +NEL*G_POSX
1670 RBUF_L(L_REL+1:L_REL+NEL*G_POSY) = GBUF%POSY(1:NEL*G_POSY)
1671 L_REL = L_REL +NEL*G_POSY
1672 RBUF_L(L_REL+1:L_REL+NEL*G_POSZ) = GBUF%POSZ(1:NEL*G_POSZ)
1673 L_REL = L_REL +NEL*G_POSZ
1674 RBUF_L(L_REL+1:L_REL+NEL*G_POSXX) = GBUF%POSXX(1:NEL*G_POSXX)
1675 L_REL = L_REL +NEL*G_POSXX
1676 RBUF_L(L_REL+1:L_REL+NEL*G_POSYY) = GBUF%POSYY(1:NEL*G_POSYY)
1677 L_REL = L_REL +NEL*G_POSYY
1678 RBUF_L(L_REL+1:L_REL+NEL*G_POSZZ) = GBUF%POSZZ(1:NEL*G_POSZZ)
1679 L_REL = L_REL +NEL*G_POSZZ
1680 RBUF_L(L_REL+1:L_REL+NEL*G_YIELD) = GBUF%YIELD(1:NEL*G_YIELD)
1681 L_REL = L_REL +NEL*G_YIELD
1682 RBUF_L(L_REL+1:L_REL+NEL*G_LENGTH_ERR) = GBUF%LENGTH_ERR(1:NEL*G_LENGTH_ERR)
1683 L_REL = L_REL +NEL*G_LENGTH_ERR
1684 RBUF_L(L_REL+1:L_REL+NEL*G_DV) = GBUF%DV(1:NEL*G_DV)
1685 L_REL = L_REL +NEL*G_DV
1686 RBUF_L(L_REL+1:L_REL+NEL*G_DFS) = GBUF%DFS(1:NEL*G_DFS)
1687 L_REL = L_REL +NEL*G_DFS
1688 RBUF_L(L_REL+1:L_REL+NEL*G_SKEW_ERR) = GBUF%SKEW_ERR(1:NEL*G_SKEW_ERR)
1689 L_REL = L_REL +NEL*G_SKEW_ERR
1690 RBUF_L(L_REL+1:L_REL+NEL*G_E6) = GBUF%E6(1:NEL*G_E6)
1691 L_REL = L_REL +NEL*G_E6
1692 RBUF_L(L_REL+1:L_REL+NEL*G_RUPTCRIT) = GBUF%RUPTCRIT(1:NEL*G_RUPTCRIT)
1693 L_REL = L_REL +NEL*G_RUPTCRIT
1694 RBUF_L(L_REL+1:L_REL+NEL*G_MASS) = GBUF%MASS(1:NEL*G_MASS)
1695 L_REL = L_REL +NEL*G_MASS
1696 RBUF_L(L_REL+1:L_REL+NEL*G_V_REPCVT) = GBUF%V_REPCVT(1:NEL*G_V_REPCVT)
1697 L_REL = L_REL +NEL*G_V_REPCVT
1698 RBUF_L(L_REL+1:L_REL+NEL*G_VR_REPCVT) = GBUF%VR_REPCVT(1:NEL*G_VR_REPCVT)
1699 L_REL = L_REL +NEL*G_VR_REPCVT
1700 RBUF_L(L_REL+1:L_REL+NEL*G_NUVAR) = GBUF%VAR(1:NEL*G_NUVAR)
1701 L_REL = L_REL +NEL*G_NUVAR
1702 RBUF_L(L_REL+1:L_REL+NEL*G_NUVARN) = GBUF%VARN(1:NEL*G_NUVARN)
1703 L_REL = L_REL +NEL*G_NUVARN
1704 RBUF_L(L_REL+1:L_REL+NEL*G_DEFINI) = GBUF%DEFINI(1:NEL*G_DEFINI)
1705 L_REL = L_REL +NEL*G_DEFINI
1706 RBUF_L(L_REL+1:L_REL+NEL*G_FORINI) = GBUF%FORINI(1:NEL*G_FORINI)
1707 L_REL = L_REL +NEL*G_FORINI
1708 RBUF_L(L_REL+1:L_REL+NEL*G_INIFRIC) = GBUF%INIFRIC(1:NEL*G_INIFRIC)
1709 L_REL = L_REL +NEL*G_INIFRIC
1710 RBUF_L(L_REL+1:L_REL+NEL*G_STRHG) = GBUF%STRHG(1:NEL*G_STRHG)
1711 L_REL = L_REL +NEL*G_STRHG
1712 RBUF_L(L_REL+1:L_REL+NEL*G_DT_PITER) = GBUF%DT_PITER(1:NEL*G_DT_PITER)
1713 L_REL = L_REL +NEL*G_DT_PITER
1714 RBUF_L(L_REL+1:L_REL+NEL*G_SKEW_ID) = GBUF%SKEW_ID(1:NEL*G_SKEW_ID)
1715 L_REL = L_REL +NEL*G_SKEW_ID
1716 RBUF_L(L_REL+1:L_REL+NEL*G_YIELD_IN_COMP) = GBUF%YIELD_IN_COMP(1:NEL*G_YIELD_IN_COMP)
1717 L_REL = L_REL +NEL*G_YIELD_IN_COMP
1718 RBUF_L(L_REL+1:L_REL+NEL*G_XXOLD_IN_COMP) = GBUF%XXOLD_IN_COMP(1:NEL*G_XXOLD_IN_COMP)
1719 L_REL = L_REL +NEL*G_XXOLD_IN_COMP
1720c
1721c failure model for beams type 3
1722.and. IF (IGTYP == 3 G_NOFF == 1) THEN
1723 IFAIL = ELBUF_TAB(NG)%GBUF%FAIL(1)%ILAWF
1724 IDFAIL = ELBUF_TAB(NG)%GBUF%FAIL(1)%IDFAIL
1725 NUVAR = ELBUF_TAB(NG)%GBUF%FAIL(1)%NVAR
1726 LF_DAM = ELBUF_TAB(NG)%GBUF%FAIL(1)%LF_DAM
1727 LF_DAMMX = ELBUF_TAB(NG)%GBUF%FAIL(1)%LF_DAMMX
1728 LF_DAMINI= ELBUF_TAB(NG)%GBUF%FAIL(1)%LF_DAMINI
1729 LF_TDEL = ELBUF_TAB(NG)%GBUF%FAIL(1)%LF_TDEL
1730 LF_INDX = ELBUF_TAB(NG)%GBUF%FAIL(1)%LF_INDX
1731 LF_OFF = ELBUF_TAB(NG)%GBUF%FAIL(1)%LF_OFF
1732c
1733 RBUF_L(L_REL+1)=IFAIL
1734 L_REL = L_REL+1
1735 RBUF_L(L_REL+1)=IDFAIL
1736 L_REL = L_REL+1
1737 RBUF_L(L_REL+1)=NUVAR
1738 L_REL = L_REL+1
1739 RBUF_L(L_REL+1)=LF_DAM
1740 L_REL = L_REL+1
1741 RBUF_L(L_REL+1)=LF_DAMMX
1742 L_REL = L_REL+1
1743 RBUF_L(L_REL+1)=LF_DAMINI
1744 L_REL = L_REL+1
1745 RBUF_L(L_REL+1)=LF_TDEL
1746 L_REL = L_REL+1
1747 RBUF_L(L_REL+1)=LF_INDX
1748 L_REL = L_REL+1
1749 RBUF_L(L_REL+1)=LF_OFF
1750 L_REL = L_REL+1
1751c
1752 RBUF_L(L_REL+1:L_REL+NEL*NUVAR) = ELBUF_TAB(NG)%GBUF%FAIL(1)%VAR(1:NEL*NUVAR)
1753 L_REL = L_REL+NEL*NUVAR
1754 RBUF_L(L_REL+1:L_REL+NEL*LF_DAM) = ELBUF_TAB(NG)%GBUF%FAIL(1)%DAM(1:NEL*LF_DAM)
1755 L_REL = L_REL+NEL*LF_DAM
1756 RBUF_L(L_REL+1:L_REL+NEL*LF_DAMMX) = ELBUF_TAB(NG)%GBUF%FAIL(1)%DAMMX(1:NEL*LF_DAMMX)
1757 L_REL = L_REL+NEL*LF_DAMMX
1758 RBUF_L(L_REL+1:L_REL+NEL*LF_DAMINI) = ELBUF_TAB(NG)%GBUF%FAIL(1)%DAMINI(1:NEL*LF_DAMINI)
1759 L_REL = L_REL+NEL*LF_DAMINI
1760 RBUF_L(L_REL+1:L_REL+NEL*LF_TDEL) = ELBUF_TAB(NG)%GBUF%FAIL(1)%TDEL(1:NEL*LF_TDEL)
1761 L_REL = L_REL + NEL*LF_TDEL
1762 RBUF_L(L_REL+1:L_REL+NEL*LF_INDX) = ELBUF_TAB(NG)%GBUF%FAIL(1)%INDX(1:NEL*LF_INDX)
1763 L_REL = L_REL + NEL*LF_INDX
1764 RBUF_L(L_REL+1:L_REL+NEL*LF_OFF) = ELBUF_TAB(NG)%GBUF%FAIL(1)%OFF(1:NEL*LF_OFF)
1765 L_REL = L_REL + NEL*LF_OFF
1766 END IF
1767c---
1768c seatbelts elements
1769c---
1770 RBUF_L(L_REL+1:L_REL+NEL*G_SLIPRING_ID) = GBUF%SLIPRING_ID(1:NEL*G_SLIPRING_ID)
1771 L_REL = L_REL +NEL*G_SLIPRING_ID
1772 RBUF_L(L_REL+1:L_REL+NEL*G_SLIPRING_FRAM_ID) = GBUF%SLIPRING_FRAM_ID(1:NEL*G_SLIPRING_FRAM_ID)
1773 L_REL = L_REL +NEL*G_SLIPRING_FRAM_ID
1774 RBUF_L(L_REL+1:L_REL+NEL*G_SLIPRING_STRAND) = GBUF%SLIPRING_STRAND(1:NEL*G_SLIPRING_STRAND)
1775 L_REL = L_REL +NEL*G_SLIPRING_STRAND
1776 RBUF_L(L_REL+1:L_REL+NEL*G_RETRACTOR_ID)= GBUF%RETRACTOR_ID(1:NEL*G_RETRACTOR_ID)
1777 L_REL = L_REL +NEL*G_RETRACTOR_ID
1778 RBUF_L(L_REL+1:L_REL+NEL*G_RINGSLIP)= GBUF%RINGSLIP(1:NEL*G_RINGSLIP)
1779 L_REL = L_REL +NEL*G_RINGSLIP
1780 RBUF_L(L_REL+1:L_REL+NEL*G_ADD_NODE)= GBUF%ADD_NODE(1:NEL*G_ADD_NODE)
1781 L_REL = L_REL +NEL*G_ADD_NODE
1782 RBUF_L(L_REL+1:L_REL+NEL*G_UPDATE) = GBUF%UPDATE(1:NEL*G_UPDATE)
1783 L_REL = L_REL +NEL*G_UPDATE
1784 RBUF_L(L_REL+1:L_REL+NEL*G_FRAM_FACTOR) = GBUF%FRAM_FACTOR(1:NEL*G_FRAM_FACTOR)
1785 L_REL = L_REL +NEL*G_FRAM_FACTOR
1786 RBUF_L(L_REL+1:L_REL+NEL*G_INTVAR) = GBUF%INTVAR(1:NEL*G_INTVAR)
1787 L_REL = L_REL +NEL*G_INTVAR
1788C----
1789 RBUF_L(L_REL+1:L_REL+NEL*G_IDT_TSH) = GBUF%IDT_TSH(1:NEL*G_IDT_TSH)
1790 L_REL = L_REL +NEL*G_IDT_TSH
1791 IF (NPG > 1) THEN
1792 RBUF_L(L_REL+1:L_REL + NEL*GBUF%G_TEMPG)= GBUF%TEMPG(1:NEL*GBUF%G_TEMPG)
1793 L_REL = L_REL + NEL*GBUF%G_TEMPG
1794 ENDIF
1795c---
1796c TMAX(H3D)
1797c---
1798 RBUF_L(L_REL+1:L_REL+NEL*G_TM_YIELD) = GBUF%TM_YIELD(1:NEL*G_TM_YIELD)
1799 L_REL = L_REL +NEL*G_TM_YIELD
1800 RBUF_L(L_REL+1:L_REL+NEL*G_TM_SEQ) = GBUF%TM_SEQ(1:NEL*G_TM_SEQ)
1801 L_REL = L_REL +NEL*G_TM_SEQ
1802 RBUF_L(L_REL+1:L_REL+NEL*G_TM_EINT) = GBUF%TM_EINT(1:NEL*G_TM_EINT)
1803 L_REL = L_REL +NEL*G_TM_EINT
1804 RBUF_L(L_REL+1:L_REL+NEL*G_TM_DMG) = GBUF%TM_DMG(1:NEL*G_TM_DMG)
1805 L_REL = L_REL +NEL*G_TM_DMG
1806 RBUF_L(L_REL+1:L_REL+NEL*G_TM_SIG) = GBUF%TM_SIG1(1:NEL*G_TM_SIG)
1807 L_REL = L_REL +NEL*G_TM_SIG
1808 RBUF_L(L_REL+1:L_REL+NEL*G_TM_SIG) = GBUF%TM_SIG3(1:NEL*G_TM_SIG)
1809 L_REL = L_REL +NEL*G_TM_SIG
1810 RBUF_L(L_REL+1:L_REL+NEL*G_TM_STRA) = GBUF%TM_STRA1(1:NEL*G_TM_STRA)
1811 L_REL = L_REL +NEL*G_TM_STRA
1812 RBUF_L(L_REL+1:L_REL+NEL*G_TM_STRA) = GBUF%TM_STRA3(1:NEL*G_TM_STRA)
1813 L_REL = L_REL +NEL*G_TM_STRA
1814c---------------------------------------------------
1815 IF (DEBUG_PRINT == 1) THEN
1816 print*,' total glob vars =', L_REL-SUM
1817 IAD = L_REL
1818 ENDIF
1819c---------------------------------------------------
1820C---------
1821 CALL WRITE_DB (RBUF_L,L_REL)
1822 IF( SDP_RBUF >0)THEN
1823 CALL WRITE_DPDB(DP_RBUF,SDP_RBUF)
1824 DEALLOCATE(DP_RBUF)
1825 ENDIF
1826C---------
1827 DEALLOCATE (RBUF_L)
1828c---------
1829 LEN_AM = LEN_AM + L_REL
1830 ENDIF ! ity
1831 ENDIF ! p = proc
1832 ENDDO ! NG=1,NGROUP
1833C-----------
1834 RETURN
1835 END
my_real
#define my_real
Definition cppsort.cpp:32
w_elbuf_str
subroutine w_elbuf_str(proc, iparg, elbuf_tab, len_am, flag_xfem)
Definition w_elbuf_str.F:31