OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
lec_inistate_yfile.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!|| lec_inistate_yfile ../starter/source/initial_conditions/inista/lec_inistate_yfile.F
25!||--- called by ------------------------------------------------------
26!|| lec_inistate ../starter/source/elements/initia/lec_inistate.f
27!||--- uses -----------------------------------------------------
28!|| format_mod ../starter/share/modules1/format_mod.F90
29!|| reader_old_mod ../starter/share/modules1/reader_old_mod.F90
30!||====================================================================
31 SUBROUTINE lec_inistate_yfile(
32 1 NSIGSH ,NSIGS ,NSIGSPH ,NSIGRS ,NSIGI ,
33 2 SIGSH ,SIGI ,SIGSPH ,SIGRS ,SIGSP ,
34 3 ISIGSH ,IUSER ,
35 4 ID_SIGSH , ID_SOLID_SIGI, ID_QUAD_SIGI )
36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE format_mod
40 USE reader_old_mod , ONLY : line
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45C-----------------------------------------------
46C C o m m o n B l o c k s
47C-----------------------------------------------
48#include "com01_c.inc"
49#include "com04_c.inc"
50#include "scr16_c.inc"
51#include "scr17_c.inc"
52#include "scry_c.inc"
53#include "units_c.inc"
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 INTEGER NSIGI, NSIGSH, NSIGS, NSIGSPH, NSIGRS, ISIGSH, IUSER
58 INTEGER ID_SIGSH(*), ID_SOLID_SIGI(*), ID_QUAD_SIGI(*)
59 my_real SIGI(NSIGS,*),SIGSH(MAX(1,NSIGSH),*),SIGSP(NSIGI,*),SIGSPH(NSIGSPH,*),SIGRS(NSIGRS,*)
60C-----------------------------------------------
61C L o c a l V a r i a b l e s
62C-----------------------------------------------
63 INTEGER K, N, I,J, IGTYP,NUMS,IHBE,NE, ISH3N,IIS,NIP,IPG,NPG,PT,J1,JJ1,NU,IP,JJ,NUVAR,NVARSH,NUMR,NUMSOLID,KK,JJHBE
64 INTEGER KN
65 INTEGER IS_SOLID,IS_QUAD
66 my_real THKI, EM, EB, H1, H2, H3, R0 , EIN, SMOOTHL
67 my_real s(6), sigloc(6)
68C=======================================================================
69C
70C -- LECTURE OF INITIAL STATE DATA - EXTRACTED FROM INITIA.F
71C
72C=======================================================================
73
74C------------------
75C FICHIER Y000
76C------------------
77 is_solid = 0
78 is_quad = 0
79 110 READ(iin5,fmt='(A)',END=199,ERR=199)line
80 IF(line(1:7)/='/SOLIDE'.AND.n2d == 0)GOTO 110
81 IF(line(1:5)/='/QUAD'.AND.n2d/=0) GOTO 110
82 IF(line(1:7) == '/SOLIDE')is_solid=1
83 IF(line(1:5) == '/QUAD')is_quad=1
84C
85 READ(iin5,fmt='(A)',END=199,ERR=199)line
86 115 READ(iin5,fmt='(A)',END=120,ERR=199)line
87 IF(line(1:1) == '#')GOTO 115
88 IF(ioutp_fmt == 2) THEN
89 IF(line(1:8) == ' ')GOTO 115
90 ELSE
91 IF(line(1:10) == ' ')GOTO 115
92 ENDIF
93 IF(line(1:1) == '/')GOTO 117
94 IF(ioutp_fmt == 2)THEN
95 READ(line,'(2I8)')i,n
96 ELSE
97 READ(line,'(2I10)')i,n
98 ENDIF
99 IF(is_solid == 1) id_solid_sigi(i) = n
100 IF(is_quad == 1) id_quad_sigi(i) = n
101 GOTO 115
102 117 IF(line(1:5)/='/TETRA'.AND.numels10 == 0) GOTO 120
103 READ(iin5,fmt='(A)',END=199,ERR=199)line
104 118 READ(iin5,fmt='(A)',END=120,ERR=199)line
105 IF(line(1:1) == '#')GOTO 118
106 IF(ioutp_fmt == 2) THEN
107 IF(line(1:8) == ' ')GOTO 118
108 ELSE
109 IF(line(1:10) == ' ')GOTO 118
110 ENDIF
111 IF(line(1:1) == '/')GOTO 120
112 IF (ioutp_fmt == 2) THEN
113 READ(line,'(2I8)')i,n
114 ELSE
115 READ(line,'(2I10)')i,n
116 ENDIF
117 id_solid_sigi(i) = n
118 GOTO 118
119C------------------
120C FICHIER Ynnn
121C EPS PLASTIC
122C------------------
123 120 READ(iin4,fmt='(A)',END=199,ERR=199)line
124 IF(line(1:8) == '/ENDDATA')THEN
125 rewind(iin4)
126 GOTO 130
127 ENDIF
128 122 IF(line(1:27)/='/SOLID /SCALAR /EPSP')GOTO 120
129 READ(iin4,fmt='(A)',END=199,ERR=199)line
130 i = 0
131 125 READ(iin4,fmt='(A)',END=130,ERR=199)line
132 IF(line(1:1) == '#')GOTO 125
133 IF(line(1:1) == '/')GOTO 132
134 IF (ioutp_fmt == 2) THEN
135 READ(line,'(6F12.0)')(s(j),j=1,6)
136 ELSE
137 READ(line,'(6F20.0)')(s(j),j=1,6)
138 ENDIF
139 DO j=1,6
140 i=i+1
141 IF(i>numsol+numquad) GOTO 130
142 sigi(10,i) = s(j)
143 ENDDO
144 GOTO 125
145C------------------
146C ENERGY
147C------------------
148 130 READ(iin4,fmt='(A)',END=199,ERR=199)line
149 IF(line(1:8) == '/ENDDATA')THEN
150 rewind(iin4)
151 GOTO 140
152 ENDIF
153 132 IF(line(1:27)/='/SOLID /SCALAR /ENER')GOTO 130
154 READ(iin4,fmt='(A)',END=199,ERR=199)line
155C
156 i = 0
157 135 READ(iin4,fmt='(A)',END=140,ERR=199)line
158 IF(line(1:1) == '#')GOTO 135
159 IF(line(1:1) == '/')GOTO 142
160 IF (ioutp_fmt == 2) THEN
161 READ(line,'(6F12.0)')(s(j),j=1,6)
162 ELSE
163 READ(line,'(6F20.0)')(s(j),j=1,6)
164 ENDIF
165 DO j=1,6
166 i=i+1
167 IF(i>numsol+numquad) GOTO 140
168 sigi(9,i) = s(j)
169 ENDDO
170 GOTO 135
171C------------------
172C RHO
173C------------------
174 140 READ(iin4,fmt='(A)',END=199,ERR=199)line
175 IF(line(1:8) == '/ENDDATA')THEN
176 rewind(iin4)
177 GOTO 144
178 ENDIF
179 142 IF(line(1:27)/='/SOLID /SCALAR /DENS')GOTO 140
180 READ(iin4,fmt='(A)',END=199,ERR=199)line
181C
182 i = 0
183 145 READ(iin4,fmt='(A)',END=144,ERR=199)line
184 IF(line(1:1) == '#')GOTO 145
185 IF(line(1:1) == '/')GOTO 146
186 IF (ioutp_fmt == 2) THEN
187 READ(line,'(6F12.0)')(s(j),j=1,6)
188 ELSE
189 READ(line,'(6F20.0)')(s(j),j=1,6)
190 ENDIF
191 DO j=1,6
192 i=i+1
193 IF(i>numsol+numquad) GOTO 144
194 sigi(8,i) = s(j)
195 ENDDO
196 GOTO 145
197C
198 144 READ(iin4,fmt='(A)',END=199,ERR=199)line
199 IF(line(1:8) == '/ENDDATA')THEN
200 rewind(iin4)
201 GOTO 150
202 ENDIF
203 146 IF(line(1:28)/='/SOLID /SCALAR /USERS')GOTO 144
204 READ(iin4,fmt='(A)',END=199,ERR=199)line
205 iuser = 1
206C
207 i = 0
208 147 READ(iin4,fmt='(A)',END=150,ERR=199)line
209 IF(line(1:1) == '#')GOTO 147
210 IF(line(1:1) == '/')GOTO 152
211 IF(ioutp_fmt == 2)THEN
212 READ(line,'(4I8)')nums,nip,nuvar,jjhbe
213 ELSE
214 READ(line,'(4I10)')nums,nip,nuvar,jjhbe
215 ENDIF
216C
217 i=i+1
218 sigsp(nvsolid1 + nvsolid2 +1 , i) = nums
219 sigsp(nvsolid1 + nvsolid2 +2 , i) = nip
220 sigsp(nvsolid1 + nvsolid2 +3 , i) = nuvar
221 sigsp(nvsolid1 + nvsolid2 +4 , i) = jjhbe
222 nu = (nuvar - mod(nuvar,6))/6
223 kn = 6
224 IF (nuvar == 0) THEN
225 GOTO 147
226 ENDIF
227 IF(nuvar < 6 )THEN
228 nu = 1
229 kn = nuvar
230 ENDIF
231 IF(i>numsol+numquad) GOTO 150
232 IF(nums == 10.OR.nums == 16.OR.nums == 20.OR.
233 . (nums == 8.AND.jjhbe == 14).OR.(nums == 8.AND.jjhbe == 17).OR.
234 . ((nums == 6.OR.nums == 8).AND.jjhbe == 15).OR.
235 . (nums == 8.AND.nip == 1.AND.
236 . jjhbe/=14.AND.jjhbe/=15).OR.nums == 4)THEN
237 DO k=1,nip
238 iis = nvsolid1 + nvsolid2 + 4 + nuvar*(k - 1)
239 DO j1 = 1,nu
240 jj1 = (j1 - 1)*6
241 149 READ(iin4,fmt='(A)',END=150,ERR=199)line
242 IF(line(1:1) == '#')GOTO 149
243 IF(line(1:1) == '/')GOTO 152
244 IF(ioutp_fmt == 2)THEN
245 READ(line,fmt='(6F12.5)')(sigsp(iis + jj1 + j,i),j=1,kn)
246 ELSE
247 READ(line,fmt='(6F20.0)')(sigsp(iis+jj1 + j,i),j=1,kn)
248 ENDIF
249 ENDDO
250 IF(nuvar > nu*6)THEN
251 151 READ(iin4,fmt='(A)',END=150,ERR=199)line
252 IF(line(1:1) == '#')GOTO 151
253 IF(line(1:1) == '/')GOTO 152
254 IF(ioutp_fmt == 2)THEN
255 READ(line,fmt='(6F12.5)') (sigsp(iis+ nu*6 + j,i),j=1,nuvar - 6*nu)
256 ELSE
257 READ(line,fmt='(6F20.0)') (sigsp(iis +6*nu + j,i),j=1,nuvar - 6*nu)
258 ENDIF
259 ENDIF
260 ENDDO ! K=1,NIP
261 ENDIF
262 GO TO 147
263C------------------
264C STRESS
265C------------------
266 150 READ(iin4,fmt='(A)',END=199,ERR=199)line
267 IF(line(1:8) == '/ENDDATA')THEN
268 rewind(iin4)
269 GOTO 170
270 ENDIF
271 152 IF(line(1:29)/='/SOLID /TENSOR /STRESS')GOTO 150
272 READ(iin4,fmt='(A)',END=199,ERR=199)line
273C
274 i = 0
275 155 READ(iin4,fmt='(A)',END=170,ERR=199)line
276 IF(line(1:1) == '#')GOTO 155
277 IF(line(1:1) == '/')GOTO 172
278 IF (ioutp_fmt == 2) THEN
279 READ(line,'(6F12.0)')(s(j),j=1,6)
280 ELSE
281 READ(line,'(6F20.0)')(s(j),j=1,6)
282 ENDIF
283 i=i+1
284 IF(i>numsol+numquad) GOTO 170
285 DO j=1,6
286 sigi(j,i) = s(j)
287 ENDDO
288 GOTO 155
289C------------------------------------------------------------------------
290C STRESS FULL + Elastoplastic + Energy + RHO in each point integration
291C------------------------------------------------------------------------
292 170 READ(iin4,fmt='(A)',END=199,ERR=199)line
293 IF(line(1:8) == '/ENDDATA')THEN
294 rewind(iin4)
295 GOTO 160
296 ENDIF
297 172 IF(line(1:30)/='/SOLID /TENSOR /STR_FUL')GOTO 170
298 READ(iin4,fmt='(A)',END=199,ERR=199)line
299C
300 i = 0
301 j = 0
302 174 READ(iin4,fmt='(A)',END=160,ERR=199)line
303 IF(line(1:1) == '#')GOTO 174
304 IF(line(1:1) == '/')GOTO 160
305 j=j+1
306 i=i+1
307 IF(i>numsol+numquad) GOTO 170
308 IF (ioutp_fmt == 2) THEN
309 READ(line,'(3I8)') nip,nums,jjhbe
310 ELSE
311 READ(line,'(3I10)') nip,nums, jjhbe
312 ENDIF
313 IF((nums == 8.OR.nums == 4).AND.jjhbe == 0)THEN
314 175 READ(iin4,fmt='(A)',END=160,ERR=199)line
315 IF(line(1:1) == '#')GOTO 175
316 IF(line(1:1) == '/')GOTO 160
317 IF (ioutp_fmt == 2) THEN
318 READ(line,fmt='(2F12.0)') ein,r0
319 ELSE
320 READ(line,fmt='(2F20.0)') ein,r0
321 ENDIF
322 sigsp(2,i) = nip
323 IF(nip == 8)THEN
324 sigsp(1,i) = 1
325 sigi(8,i) = r0
326 sigi(9,i) = ein
327 ELSEIF(nip == 1) THEN
328 sigi(8,i) = r0
329 sigi(9,i) = ein
330 ENDIF
331 sigsp(3,i) = ein
332 sigsp(4,i) = r0
333C
334 IF (nip == 1) THEN
335 186 READ(iin4,fmt='(A)',END=160,ERR=199)line
336 IF(line(1:1) == '#')GOTO 186
337 IF(line(1:1) == '/')GOTO 160
338 IF (ioutp_fmt == 2) THEN
339 READ(line,fmt='(6F12.0)') (sigi(j,i), j=1,6)
340 ELSE
341 READ(line,fmt='(6F20.0)') (sigi(j,i), j=1,6)
342 ENDIF
343 187 READ(iin4,fmt='(A)',END=160,ERR=199)line
344 IF(line(1:1) == '#')GOTO 187
345 IF(line(1:1) == '/')GOTO 160
346 IF (ioutp_fmt == 2) THEN
347 READ(line,fmt='(1F12.0)') sigi(10,i)
348 ELSE
349 READ(line,fmt='(1F20.0)') sigi(10,i)
350 ENDIF
351 ELSE
352 DO j=1,6
353 sigi(j,i) = zero
354 ENDDO
355 sigi(10,i) = zero
356C
357 DO k=1,nip
358 iis= 4 + (k-1)*7
359 176 READ(iin4,fmt='(A)',END=160,ERR=199)line
360 IF(line(1:1) == '#')GOTO 176
361 IF(line(1:1) == '/')GOTO 160
362 IF (ioutp_fmt == 2) THEN
363 READ(line,fmt='(6F12.0)')
364 . sigsp(iis+1,i),sigsp(iis+2,i),sigsp(iis+3,i),
365 . sigsp(iis+4,i),sigsp(iis+5,i),sigsp(iis+6,i)
366 ELSE
367 READ(line,fmt='(6F20.0)')
368 . sigsp(iis+1,i),sigsp(iis+2,i),sigsp(iis+3,i),
369 . sigsp(iis+4,i),sigsp(iis+5,i),sigsp(iis+6,i)
370 ENDIF
371 177 READ(iin4,fmt='(A)',END=160,ERR=199)line
372 IF(line(1:1) == '#')GOTO 177
373 IF(line(1:1) == '/')GOTO 160
374 IF (ioutp_fmt == 2) THEN
375 READ(line,fmt='(1F12.0)') sigsp(iis+7,i)
376 ELSE
377 READ(line,fmt='(1F20.0)') sigsp(iis+7,i)
378 ENDIF
379C
380 DO j=1,6
381 sigi(j,i) = sigi(j,i) + 0.125*sigsp(iis+j,i)
382 ENDDO
383 sigi(10,i)= sigi(10,i) + 0.125*sigsp(iis+7,i)
384 ENDDO
385 ENDIF
386c
387 ELSEIF(nums == 10.OR.nums == 16.OR.nums == 20.OR.
388 . (nums == 8.AND.jjhbe == 14).OR.(nums == 8.AND.jjhbe == 17).OR.
389 . ((nums==6.OR.nums==8).AND.(jjhbe==15.or.jjhbe==12)))THEN
390c
391 DO k=1,nip
392 iis=1 + (k-1)*9
393 178 READ(iin4,fmt='(A)',END=160,ERR=199)line
394 IF(line(1:1) == '#')GOTO 178
395 IF(line(1:1) == '/')GOTO 160
396 READ(line,fmt='(6F12.0)')
397 IF (ioutp_fmt == 2) THEN
398 READ(line,fmt='(6F12.0)')
399 . sigsp(iis+1,i),sigsp(iis+2,i),sigsp(iis+3,i),
400 . sigsp(iis+4,i),sigsp(iis+5,i),sigsp(iis+6,i)
401 ELSE
402 READ(line,fmt='(6F20.0)')
403 . sigsp(iis+1,i),sigsp(iis+2,i),sigsp(iis+3,i),
404 . sigsp(iis+4,i),sigsp(iis+5,i),sigsp(iis+6,i)
405 ENDIF
406
407 179 READ(iin4,fmt='(A)',END=160,ERR=199)line
408 IF(line(1:1) == '#')GOTO 179
409 IF(line(1:1) == '/')GOTO 160
410 IF (ioutp_fmt == 2) THEN
411 READ(line,fmt='(3F12.0)') sigsp(iis+7,i),sigsp(iis+8,i), sigsp(iis+9,i)
412 ELSE
413 READ(line,fmt='(3F20.0)') sigsp(iis+7,i),sigsp(iis+8,i), sigsp(iis+9,i)
414 ENDIF
415 ENDDO
416 ENDIF
417 GOTO 174
418 160 CONTINUE
419 rewind(iin4)
420C--------------------------------------
421C STRAIN FULL in each point integration
422C--------------------------------------
423 180 READ(iin4,fmt='(A)',END=199,ERR=199)line
424 IF(line(1:8) == '/ENDDATA')THEN
425 rewind(iin4)
426 GOTO 199
427 ENDIF
428 182 IF(line(1:33)/='/SOLID /TENSOR /STRAIN_FUL')GOTO 180
429 READ(iin4,fmt='(A)',END=199,ERR=199)line
430C
431 i = 0
432 j = 0
433 184 READ(iin4,fmt='(A)',END=199,ERR=199)line
434 IF(line(1:1) == '#')GOTO 184
435 IF(line(1:1) == '/')GOTO 199
436 j=j+1
437 IF(i>numsol+numquad) GOTO 199
438 IF (ioutp_fmt == 2) THEN
439 READ(line,'(3I8)')nip,nums,numsolid
440 ELSE
441 READ(line,'(3I10)')nip,nums,numsolid
442 ENDIF
443C
444 DO kk = 1,numsolid
445 i = i + 1
446 DO k=1,nip
447 iis= nvsolid1 + (k-1)*6
448 185 READ(iin4,fmt='(A)',END=199,ERR=199)line
449 IF(line(1:1) == '#')GOTO 185
450 IF(line(1:1) == '/')GOTO 199
451 IF (ioutp_fmt == 2) THEN
452 READ(line,fmt='(6F12.0)')
453 . sigsp(iis+1,i),sigsp(iis+2,i),sigsp(iis+3,i),
454 . sigsp(iis+4,i),sigsp(iis+5,i),sigsp(iis+6,i)
455 ELSE
456 READ(line,fmt='(6F20.0)')
457 . sigsp(iis+1,i),sigsp(iis+2,i),sigsp(iis+3,i),
458 . sigsp(iis+4,i),sigsp(iis+5,i),sigsp(iis+6,i)
459 ENDIF
460 ENDDO
461 ENDDO
462 GOTO 184
463 199 CONTINUE
464C
465C--------------------------------------------------------
466C CONTRAINTES INITIALES + energies densites eps-plast
467C SHELL 4N/3N
468C--------------------------------------------------------
469C
470C FICHIER Y000
471C
472 rewind(iin4)
473 rewind(iin5)
474C
475 310 READ(iin5,fmt='(A)',END=399,ERR=399)line
476 IF(line(1:6)/='/SHELL')GOTO 310
477 READ(iin5,fmt='(A)',END=399,ERR=399)line
478C
479 312 READ(iin5,fmt='(A)',END=314,ERR=399)line
480 IF(line(1:1) == '#')GOTO 312
481 IF(ioutp_fmt == 2) THEN
482 IF(line(1:8) == ' ')GOTO 312
483 ELSE
484 IF(line(1:10) == ' ')GOTO 312
485 ENDIF
486 IF(line(1:1) == '/')GOTO 314
487 IF (ioutp_fmt == 2) THEN
488 READ(line,'(2I8)')i,n
489 ELSE
490 READ(line,'(2I10)')i,n
491 ENDIF
492 IF(i>numshel) GOTO 314
493 sigsh(1,i) = n
494 id_sigsh(i) = n
495 GOTO 312
496C
497 314 rewind(iin5)
498 315 READ(iin5,fmt='(A)',END=399,ERR=399)line
499 IF(line(1:8)/='/SHELL3N')GOTO 315
500 READ(iin5,fmt='(A)',END=399,ERR=399)line
501C
502 316 READ(iin5,fmt='(A)',END=320,ERR=399)line
503 IF(line(1:1) == '#')GOTO 316
504 IF(ioutp_fmt == 2) THEN
505 IF(line(1:8) == ' ')GOTO 316
506 ELSE
507 IF(line(1:10) == ' ')GOTO 316
508 ENDIF
509 IF(line(1:1) == '/')GOTO 320
510 IF (ioutp_fmt == 2) THEN
511 READ(line,'(2I8)')i,n
512 ELSE
513 READ(line,'(2I10)')i,n
514 ENDIF
515 sigsh(1,i+numshel) = n
516 id_sigsh(i+numshel) = n
517 IF(i>numsh3n) GOTO 320
518 GOTO 316
519C----------------------------------------------
520C FICHIER Ynnn
521C NIP + THK + ENER + STRESS + EPS PLASTIC
522C----------------------------------------------
523 320 READ(iin4,fmt='(A)',END=399,ERR=399)line
524 IF(line(1:8) == '/ENDDATA')THEN
525 rewind(iin4)
526 GOTO 350
527 ENDIF
528 322 IF(line(1:33)/='/SHELL /TENSOR /STRESS_FUL')GOTO 320
529 READ(iin4,fmt='(A)',END=399,ERR=399)line
530 isigsh =1
531C
532 i = 0
533 j = 0
534 325 READ(iin4,fmt='(A)',END=350,ERR=399)line
535 IF(line(1:1) == '#')GOTO 325
536 IF(line(1:1) == '/')GOTO 352
537 j=j+1
538 i=i+1
539 IF(j>numshel+numsh3n) GOTO 350
540 IF (ioutp_fmt == 2) THEN
541 READ(line,'(2I8)')nip,npg
542 READ(iin4,fmt='(6F12.0)',END=399,ERR=399)THKI,EM,EB,H1,H2,h3
543 ELSE
544 READ(line,fmt=fmt_2i) nip,npg
545 READ(iin4,fmt='(6F20.0)',END=399,ERR=399)THKI,EM,EB,H1,H2,h3
546 ENDIF
547 sigsh(2,i) = nip
548 sigsh(3,i) = thki
549 sigsh(4,i) = em
550 sigsh(5,i) = eb
551 sigsh(nvshell - 1,i) = one
552 IF (npg == 0.OR.npg == 1)THEN
553 sigsh(14,i) = h1
554 sigsh(15,i) = h2
555 sigsh(16,i) = h3
556 IF(nip == 0)THEN
557 IF (ioutp_fmt == 2) THEN
558 READ(iin4,fmt='(6F12.0)',END=399,ERR=399)
559 . sigsh(22,i),sigsh(23,i),sigsh(24,i),
560 . sigsh(25,i),sigsh(26,i),sigsh(27,i)
561 READ(iin4,fmt='(3F12.0)',END=399,ERR=399)
562 . sigsh(28,i),sigsh(29,i),sigsh(30,i)
563 ELSE
564 READ(iin4,fmt='(6F20.0)',END=399,ERR=399)
565 . sigsh(22,i),sigsh(23,i),sigsh(24,i),
566 . sigsh(25,i),sigsh(26,i),sigsh(27,i)
567 READ(iin4,fmt='(3F20.0)',END=399,ERR=399)
568 . sigsh(28,i),sigsh(29,i),sigsh(30,i)
569 ENDIF
570 ELSEIF(nip/=0)THEN
571 IF (ioutp_fmt == 2) THEN
572 DO k=1,nip
573 pt = 22 + 6*(k - 1)
574 READ(iin4,fmt='(6F12.0)',END=399,ERR=399)
575 . sigsh(pt,i ),sigsh(pt + 1,i),sigsh(pt + 2,i),
576 . sigsh(pt + 3,i),sigsh(pt + 4,i),sigsh(pt + 5,i)
577 ENDDO
578 ELSE
579 DO k=1,nip
580 pt = 22 + 6*(k - 1)
581 READ(iin4,fmt='(6F20.0)',END=399,ERR=399)
582 . sigsh(pt,i),sigsh(pt + 1,i),sigsh(pt + 2,i),
583 . sigsh(pt + 3,i),sigsh(pt + 4,i),sigsh(pt + 5,i)
584 ENDDO
585 ENDIF
586 ENDIF
587 ELSEIF(npg>1) THEN
588 sigsh(nvshell,i) = npg
589 IF(nip == 0)THEN
590 IF (ioutp_fmt == 2) THEN
591 DO n=1,npg
592 pt= 22 + (n-1)*9
593 READ(iin4,fmt='(6F12.0)',END=399,ERR=399)
594 . sigsh(pt,i),sigsh(pt+1,i),sigsh(pt+2,i),
595 . sigsh(pt+3,i),sigsh(pt+4,i),sigsh(pt+5,i)
596 READ(iin4,fmt='(3F12.0)',END=399,ERR=399)
597 . sigsh(pt+6,i),sigsh(pt+7,i),sigsh(pt+8,i)
598 ENDDO
599 ELSE
600 DO n=1,npg
601 pt= 22 + (n-1)*9
602 READ(iin4,fmt='(6F20.0)',END=399,ERR=399)
603 . sigsh(pt,i),sigsh(pt+1,i),sigsh(pt+2,i),
604 . sigsh(pt+3,i),sigsh(pt+4,i),sigsh(pt+5,i)
605 READ(iin4,fmt='(3F20.0)',END=399,ERR=399)
606 . sigsh(pt+6,i),sigsh(pt+7,i),sigsh(pt+8,i)
607 ENDDO
608 ENDIF
609 ELSE
610 pt = 22
611 IF (ioutp_fmt == 2) THEN
612 DO k=1,nip
613 DO n=1,npg
614 READ(iin4,fmt='(6F12.0)',END=399,ERR=399)
615 . sigsh(pt,i),sigsh(pt+1,i),sigsh(pt+2,i),
616 . sigsh(pt+3,i),sigsh(pt+4,i),sigsh(pt+5,i)
617 pt = pt + 6
618 ENDDO
619 ENDDO
620 ELSE
621 DO k=1,nip
622 DO n=1,npg
623 READ(iin4,fmt='(6F20.0)',END=399,ERR=399)
624 . sigsh(pt,i),sigsh(pt+1,i),sigsh(pt+2,i),
625 . sigsh(pt+3,i),sigsh(pt+4,i),sigsh(pt+5,i)
626 pt = pt + 6
627 ENDDO
628 ENDDO
629 ENDIF
630 ENDIF
631 ELSE
632c CALL ANCERR(58,ANINFO_BLIND_2)
633 ENDIF
634 GOTO 325
635C------------------
636C------------------
637 350 READ(iin4,fmt='(A)',END=360,ERR=399)line
638 IF(line(1:8) == '/ENDDATA')THEN
639 rewind(iin4)
640 GOTO 360
641 ENDIF
642 352 IF(line(1:33)/='/SHELL /TENSOR /STRAIN_FUL')GOTO 350
643 READ(iin4,fmt='(A)',END=399,ERR=399)line
644 isigsh =1
645C
646 i = 0
647 j = 0
648 355 READ(iin4,fmt='(A)',END=360,ERR=399)line
649 IF(line(1:1) == '#')GOTO 355
650 IF(line(1:1) == '/')GOTO 360
651 j=j+1
652 i=i+1
653 IF(j> numshel+numsh3n) GOTO 360
654 sigsh(nvshell - 1,i) = one
655 IF (ioutp_fmt == 2) THEN
656 READ(line,'(6F12.0)')sigsh(6,i),sigsh(7,i),sigsh(8,i),
657 . sigsh(9,i),sigsh(10,i),sigsh(27,i)
658 READ(iin4,fmt='(A)',END=399,ERR=399)line
659 READ(line,'(3F12.0)')sigsh(11,i),sigsh(12,i),sigsh(13,i)
660 ELSE
661 READ(line,'(6F20.0)')sigsh(6,i),sigsh(7,i),sigsh(8,i),
662 . sigsh(9,i),sigsh(10,i),sigsh(27,i)
663 READ(iin4,fmt='(A)',END=399,ERR=399)line
664 READ(line,'(3F20.0)')sigsh(11,i),sigsh(12,i),sigsh(13,i)
665 ENDIF
666 GOTO 355
667C
668 360 CONTINUE
669C--------------------------------
670C /SHELL /SCALAR /EPSP
671C--------------------------------
672 rewind(iin4)
673 420 READ(iin4,fmt='(A)',END=399,ERR=399)line
674 IF(line(1:8) == '/ENDDATA')THEN
675 rewind(iin4)
676 GO TO 520
677 ENDIF
678 422 IF(line(1:27)/='/SHELL /SCALAR /EPSP')GOTO 420
679 READ(iin4,fmt='(A)',END=399,ERR=399)line
680 isigsh =1
681C
682 i = 0
683 j = 0
684 425 READ(iin4,fmt='(A)',END=520,ERR=399)line
685 IF(line(1:1) == '#')GOTO 425
686 IF(line(1:1) == '/')GOTO 522
687 j=j+1
688 i=i+1
689 IF(j>numshel+numsh3n) GOTO 520
690 sigsh(nvshell - 1,i) = one
691 IF (ioutp_fmt == 2) THEN
692 READ(line,fmt='(6F12.0)')sigloc
693 ELSE
694 READ(line,fmt='(6F20.0)')sigloc
695 ENDIF
696 sigsh(27,i)=sigloc(1)
697 j=j+1
698 i=i+1
699 IF(j>numshel+numsh3n) GOTO 520
700 sigsh(27,i)=sigloc(2)
701 j=j+1
702 i=i+1
703 IF(j>numshel+numsh3n) GOTO 520
704 sigsh(27,i)=sigloc(3)
705 j=j+1
706 i=i+1
707 IF(j>numshel+numsh3n) GOTO 520
708 sigsh(27,i)=sigloc(4)
709 j=j+1
710 i=i+1
711 IF(j>numshel+numsh3n) GOTO 520
712 sigsh(27,i)=sigloc(5)
713 j=j+1
714 i=i+1
715 IF(j>numshel+numsh3n) GOTO 520
716 sigsh(27,i)=sigloc(6)
717 GOTO 425
718 520 READ(iin4,fmt='(A)',END=399,ERR=399)line
719 IF(line(1:8) == '/ENDDATA')THEN
720 rewind(iin4)
721 GOTO 450
722 ENDIF
723 522 IF(line(1:28)/='/shell /scalar /users')GOTO 520
724 READ(IIN4,FMT='(a)',END=399,ERR=399)LINE
725C
726 IUSER = 1
727C
728 I = 0
729 J = 0
730 525 READ(IIN4,FMT='(a)',END=450,ERR=399)LINE
731 IF(LINE(1:1) == '#')GOTO 525
732 IF(line(1:1) == '/')GOTO 450
733 j=j+1
734 i=i+1
735 IF(j>numshel+numsh3n) GOTO 450
736 IF(ioutp_fmt == 2)THEN
737 READ(line,fmt='(4I8)') ihbe,nip,npg,nuvar
738 ELSE
739 READ(line,fmt='(4I10)')ihbe,nip,npg,nuvar
740 ENDIF
741 sigsh(nvshell + 1 ,i) = ihbe
742 sigsh(nvshell + 2 ,i) = nip
743 sigsh(nvshell + 3 ,i) = npg
744 sigsh(nvshell + 4 ,i) = nuvar
745 nvarsh = nvshell + 4
746 nu = (nuvar - mod(nuvar,6))/6
747 kn = 6
748 IF( nuvar < 6)THEN
749 nu = 1
750 kn = nuvar
751 ENDIF
752C-----------------------
753C COQUE STANDARD
754C IHBE different de 11
755C-----------------------
756 IF(npg == 0)THEN
757 IF(nip == 0)THEN
758 DO j1 = 1,nu
759 jj1 = (j1 - 1)*6
760 526 READ(iin4,fmt='(A)',END=450,ERR=399)line
761 IF(line(1:1) == '#')GOTO 526
762 IF(line(1:1) == '/')GOTO 450
763 IF(ioutp_fmt == 2)THEN
764 READ(line,'(6F12.5)')(sigsh(nvarsh+k+jj1,i),k=1,kn)
765 ELSE
766 READ(line,'(6F20.0)')(sigsh(nvarsh+k+jj1,i),k=1,kn)
767 ENDIF
768 ENDDO
769 IF(nuvar > 6*nu)THEN
770 536 READ(iin4,fmt='(A)',END=450,ERR=399)line
771 IF(line(1:1) == '#')GOTO 536
772 IF(line(1:1) == '/')GOTO 450
773 IF(ioutp_fmt == 2)THEN
774 READ(line,'(6F12.5)')(sigsh(nvarsh+k+6*nu,i),k=1,nuvar-6*nu)
775 ELSE
776 READ(line,'(6F20.0)')(sigsh(nvarsh+k+6*nu,i),k=1,nuvar-6*nu)
777 ENDIF
778 ENDIF
779 ELSE
780 pt = 0
781 DO ip = 1,nip
782 DO j1 = 1, nu
783 jj1 = (j1 - 1)*6
784 527 READ(iin4,fmt='(A)',END=450,ERR=399)line
785 IF(line(1:1) == '#')GOTO 527
786 IF(line(1:1) == '/')GOTO 450
787 IF(ioutp_fmt == 2)THEN
788 READ(line,'(6F12.5)')(sigsh(nvarsh+pt+jj1+k,i),k=1,kn)
789 ELSE
790 READ(line,'(6F20.0)')(sigsh(nvarsh+pt+jj1+k,i),k=1,kn)
791 ENDIF
792 ENDDO
793 IF(nuvar > nu*6)THEN
794 537 READ(iin4,fmt='(A)',END=450,ERR=399)line
795 IF(line(1:1) == '#')GOTO 537
796 IF(line(1:1) == '/')GOTO 450
797 IF(ioutp_fmt == 2)THEN
798 READ(line,'(6F12.5)')
799 . (sigsh(nvarsh+pt+nu*6+k,i),k=1,nuvar-6*nu)
800 ELSE
801 READ(line,'(6F20.0)')
802 . (sigsh(nvarsh+pt+nu*6+k,i),k=1,nuvar-6*nu)
803 ENDIF
804 ENDIF
805 pt = pt + nuvar
806 ENDDO
807 ENDIF
808 ELSE
809 IF(nip == 0)THEN
810 pt = 0
811 DO n=1,npg
812 DO j1 = 1,nu
813 jj1 = (j1 - 1)*6
814 528 READ(iin4,fmt='(A)',END=450,ERR=399)line
815 IF(line(1:1) == '#')GOTO 528
816 IF(line(1:1) == '/')GOTO 450
817 IF(ioutp_fmt == 2)THEN
818 READ(line,'(6F12.5)')(sigsh(nvarsh+pt+jj1+k,i),k=1,kn)
819 ELSE
820 READ(line,'(6F20.0)')(sigsh(nvarsh+pt+jj1+k,i),k=1,kn)
821 ENDIF
822 ENDDO
823 IF(nuvar > nu*6)THEN
824 538 READ(iin4,fmt='(A)',END=450,ERR=399)line
825 IF(line(1:1) == '#')GOTO 538
826 IF(line(1:1) == '/')GOTO 450
827 IF(ioutp_fmt == 2)THEN
828 READ(line,'(6F12.5)')
829 . (sigsh(nvarsh+pt+nu*6+k,i),k=1,nuvar-6*nu)
830 ELSE
831 READ(line,'(6F20.0)')
832 . (sigsh(nvarsh+pt+nu*6+k,i),k=1,nuvar-6*nu)
833 ENDIF
834 ENDIF
835 pt = pt + nuvar
836 ENDDO
837 ELSE
838 pt = 0
839 DO n = 1,npg
840 DO ip = 1,nip
841 DO j1 = 1, nu
842 jj1 = (j1 -1)*6
843 529 READ(iin4,fmt='(A)',END=450,ERR=399)line
844 IF(line(1:1) == '#')GOTO 529
845 IF(line(1:1) == '/')GOTO 450
846 IF(ioutp_fmt == 2)THEN
847 READ(line,'(6F12.5)')(sigsh(nvarsh+k+pt+jj1,i),k=1,kn)
848 ELSE
849 READ(line,'(6F20.0)')(sigsh(nvarsh+k+pt+jj1,i),k=1,kn)
850 ENDIF
851 ENDDO
852 IF(nuvar > 6*nu)THEN
853 539 READ(iin4,fmt='(A)',END=450,ERR=399)line
854 IF(line(1:1) == '#')GOTO 539
855 IF(line(1:1) == '/')GOTO 450
856 IF(ioutp_fmt == 2)THEN
857 READ(line,'(6F12.5)')
858 . (sigsh(nvarsh+k+pt+6*nu,i),k=1,nuvar - 6*nu)
859 ELSE
860 READ(line,'(6F20.0)')
861 . ( sigsh( nvarsh+k+pt+6*nu,i ),k=1,nuvar - 6*nu )
862 ENDIF
863 ENDIF
864 pt = pt + nuvar
865 ENDDO
866 ENDDO
867 ENDIF
868 ENDIF
869 GO TO 525
870 450 CONTINUE
871 399 CONTINUE
872C
873C--------------------------------------------------------
874C CONTRAINTES INITIALES + ENERGIES DENSITES EPS-PLAST
875C SPH particles
876C--------------------------------------------------------
877C
878C FICHIER Y000
879C
880 rewind(iin4)
881 rewind(iin5)
882 610 READ(iin5,fmt='(A)',END=699,ERR=699)line
883 IF(line(1:7)/='/SPHCEL')GOTO 610
884C
885 READ(iin5,fmt='(A)',END=699,ERR=699)line
886C
887 615 READ(iin5,fmt='(A)',END=620,ERR=699)line
888 IF(line(1:1) == '#')GOTO 615
889 IF(line(1:8) == ' ')GOTO 615
890 IF(line(1:1) == '/')GOTO 620
891 IF(ioutp_fmt == 2)THEN
892 READ(line,'(2I8)')i,n
893 ELSE
894 READ(line,'(2I10)')i,n
895 END IF
896 sigsph(7,i) = n
897 GOTO 615
898C
899C FICHIER Ynnn
900C EPS PLASTIC
901C
902 620 READ(iin4,fmt='(A)',END=699,ERR=699)line
903 IF(line(1:8) == '/ENDDATA')THEN
904 rewind(iin4)
905 GOTO 630
906 ENDIF
907 622 IF(line(1:27)/='/SPHCEL /SCALAR /EPSP')GOTO 620
908 READ(iin4,fmt='(A)',END=699,ERR=699)line
909C
910 i = 0
911 625 READ(iin4,fmt='(A)',END=630,ERR=699)line
912 IF(line(1:1) == '#')GOTO 625
913 IF(line(1:1) == '/')GOTO 632
914 IF(ioutp_fmt == 2)THEN
915 READ(line,'(6F12.0)')(s(j),j=1,6)
916 ELSE
917 READ(line,'(6F20.0)')(s(j),j=1,6)
918 END IF
919 DO j=1,6
920 i=i+1
921 IF(i>numsphy) GOTO 630
922 sigsph(10,i) = s(j)
923 ENDDO
924 GOTO 625
925C------------------
926C ENERGY
927C------------------
928 630 READ(iin4,fmt='(A)',END=699,ERR=699)line
929 IF(line(1:8) == '/ENDDATA')THEN
930 rewind(iin4)
931 GOTO 640
932 ENDIF
933 632 IF(line(1:27)/='/SPHCEL /SCALAR /ENER')GOTO 630
934 READ(iin4,fmt='(A)',END=699,ERR=699)line
935C
936 i = 0
937 635 READ(iin4,fmt='(A)',END=640,ERR=699)line
938 IF(line(1:1) == '#')GOTO 635
939 IF(line(1:1) == '/')GOTO 642
940 IF(ioutp_fmt == 2)THEN
941 READ(line,'(6F12.0)')(s(j),j=1,6)
942 ELSE
943 READ(line,'(6F20.0)')(s(j),j=1,6)
944 END IF
945 DO j=1,6
946 i=i+1
947 IF(i>numsphy) GOTO 640
948 sigsph(9,i) = s(j)
949 ENDDO
950 GOTO 635
951C------------------
952C RHO
953C------------------
954 640 READ(iin4,fmt='(A)',END=699,ERR=699)line
955 IF(line(1:8) == '/ENDDATA')THEN
956 rewind(iin4)
957 GOTO 650
958 ENDIF
959 642 IF(line(1:27)/='/SPHCEL /SCALAR /DENS')GOTO 640
960 READ(iin4,fmt='(A)',END=699,ERR=699)line
961C
962 i = 0
963 645 READ(iin4,fmt='(A)',END=650,ERR=699)line
964 IF(line(1:1) == '#')GOTO 645
965 IF(line(1:1) == '/')GOTO 652
966 IF(ioutp_fmt == 2)THEN
967 READ(line,'(6F12.0)')(s(j),j=1,6)
968 ELSE
969 READ(line,'(6F20.0)')(s(j),j=1,6)
970 END IF
971 DO j=1,6
972 i=i+1
973 IF(i>numsphy) GOTO 650
974 sigsph(8,i) = s(j)
975 ENDDO
976 GOTO 645
977C------------------
978C STRESS
979C------------------
980 650 READ(iin4,fmt='(A)',END=699,ERR=699)line
981 IF(line(1:8) == '/ENDDATA')THEN
982 rewind(iin4)
983 GOTO 660
984 ENDIF
985 652 IF(line(1:27)/='/SPHCEL /TENSOR /STRESS')GOTO 650
986 READ(iin4,fmt='(A)',END=699,ERR=699)line
987C
988 i = 0
989 655 READ(iin4,fmt='(A)',END=660,ERR=699)line
990 IF(line(1:1) == '#')GOTO 655
991 IF(line(1:1) == '/')GOTO 662
992 IF(ioutp_fmt == 2)THEN
993 READ(line,'(6F12.0)')(s(j),j=1,6)
994 ELSE
995 READ(line,'(6F20.0)')(s(j),j=1,6)
996 END IF
997 i=i+1
998 IF(i>numsphy) GOTO 660
999 DO j=1,6
1000 sigsph(j,i) = s(j)
1001 ENDDO
1002 GOTO 655
1003C-----------------------
1004C Smoothing length
1005C-----------------------
1006 660 READ(iin4,fmt='(A)',END=699,ERR=699)line
1007 IF(line(1:8) == '/ENDDATA')THEN
1008 rewind(iin4)
1009 GOTO 670
1010 ENDIF
1011 662 IF(line(1:27)/='/SPHCEL /SCALAR /SLEN')GOTO 660
1012 READ(iin4,fmt='(A)',END=699,ERR=699)line
1013C
1014 i = 0
1015 665 READ(iin4,fmt='(A)',END=670,ERR=699)line
1016 IF(line(1:1) == '#')GOTO 665
1017 IF(line(1:1) == '/')GOTO 672
1018 IF(ioutp_fmt == 2)THEN
1019 READ(line,'(6F12.0)')(s(j),j=1,6)
1020 ELSE
1021 READ(line,'(6F20.0)')(s(j),j=1,6)
1022 END IF
1023 DO j=1,6
1024 i=i+1
1025 IF(i>numsphy) GOTO 670
1026 sigsph(11,i) = s(j)
1027 ENDDO
1028 GOTO 665
1029C-------------------------------------------------------------------
1030C STRESS FULL + Elastoplastic + Energy + RHO + Smoothing length
1031C-------------------------------------------------------------------
1032 670 READ(iin4,fmt='(A)',END=699,ERR=699)line
1033 IF(line(1:8) == '/ENDDATA')THEN
1034 rewind(iin4)
1035 GOTO 680
1036 ENDIF
1037 672 IF(line(1:28)/='/SPHCEL /TENSOR /STR_FUL')GOTO 670
1038 READ(iin4,fmt='(A)',END=699,ERR=699)line
1039C
1040 i = 0
1041 675 READ(iin4,fmt='(A)',END=680,ERR=699)line
1042 IF(line(1:1) == '#')GOTO 675
1043 IF(line(1:1) == '/')GOTO 680
1044 i=i+1
1045 IF(i>numsphy) GOTO 670
1046 IF(ioutp_fmt == 2)THEN
1047 READ(line,fmt='(3F12.0)') ein,r0,smoothl
1048 ELSE
1049 READ(line,fmt='(3F20.0)') ein,r0,smoothl
1050 END IF
1051 sigsph(8,i) = r0
1052 sigsph(9,i) = ein
1053 sigsph(11,i) = smoothl
1054
1055 686 READ(iin4,fmt='(A)',END=680,ERR=699)line
1056 IF(line(1:1) == '#')GOTO 686
1057 IF(line(1:1) == '/')GOTO 680
1058 IF(ioutp_fmt == 2)THEN
1059 READ(line,fmt='(6F12.0)')(sigsph(j,i), j=1,6)
1060 ELSE
1061 READ(line,fmt='(6F20.0)')(sigsph(j,i), j=1,6)
1062 END IF
1063 687 READ(iin4,fmt='(A)',END=680,ERR=699)line
1064 IF(line(1:1) == '#')GOTO 687
1065 IF(line(1:1) == '/')GOTO 680
1066 IF(ioutp_fmt == 2)THEN
1067 READ(line,fmt='(1F12.0)') sigsph(10,i)
1068 ELSE
1069 READ(line,fmt='(1F20.0)') sigsph(10,i)
1070 END IF
1071 GOTO 675
1072 680 CONTINUE
1073 699 CONTINUE
1074C
1075C--------------------------------------------------------------
1076C EFFORT et MOMENTS INITIALES + energies densites eps-plast
1077C Ressort
1078C--------------------------------------------------------------
1079C
1080C /INISPRI - *STY file -
1081C
1082C
1083C FICHIER Y000
1084C
1085 rewind(iin4)
1086 rewind(iin5)
1087 710 READ(iin5,fmt='(A)',END=899,ERR=899)line
1088 IF (line(1:7) /= '/SPRING') GOTO 710
1089C
1090 READ(iin5,fmt='(A)',END=899,ERR=899)line
1091C
1092 715 READ(iin5,fmt='(A)',END=720,ERR=899)line
1093 IF (line(1:1) == '#') GOTO 715
1094 IF(ioutp_fmt == 2) THEN
1095 IF(line(1:8) == ' ')GOTO 715
1096 ELSE
1097 IF(line(1:10) == ' ')GOTO 715
1098 ENDIF
1099 IF (line(1:1) == '/') GOTO 720
1100 READ(line,fmt_2i) i,n ! SYSSPRI, USRSPRI
1101 sigrs(1,i) = n
1102 GOTO 715
1103C
1104C FICHIER Ynnn
1105C FORCE + MOMENT + EPSP + ENERGY
1106C
1107 720 READ(iin4,fmt='(A)',END=899,ERR=899)line
1108 IF (line(1:8) == '/ENDDATA') THEN
1109 rewind(iin4)
1110 GOTO 740
1111 ENDIF
1112 722 IF (line(1:16) /= '/SPRING /FULL') GOTO 720
1113 READ(iin4,fmt='(A)',END=740,ERR=899)line
1114C
1115 i = 0
1116 725 READ(iin4,fmt='(A)',END=740,ERR=899)line
1117 IF (line(1:1) == '#') GOTO 725
1118 IF (line(1:1) == '/') GOTO 740
1119 IF (ioutp_fmt == 2) THEN
1120 READ(line,'(3I8)')igtyp,numr,nuvar
1121 ELSE
1122 READ(line,'(3I10)')igtyp,numr,nuvar
1123 ENDIF ! IF (IOUTP_FMT == 2)
1124C------
1125 IF (igtyp == 4 .OR. igtyp == 12) THEN
1126C------
1127 DO k=1,numr
1128 i = i + 1
1129 IF (i > numelr) GOTO 899
1130 726 READ(iin4,fmt='(A)',END=740,ERR=899)line
1131 IF (line(1:1) == '#') GOTO 726
1132 IF (line(1:1) == '/') GOTO 740
1133 READ(line,fmt=fmt_5f) (sigrs(1+j,i),j=1,5)
1134 READ(iin4,fmt='(A)',END=899,ERR=899)line
1135 READ(line,fmt=fmt_3f) (sigrs(6+j,i),j=1,3)
1136 ENDDO ! DO K=1,NUMR
1137C------
1138 ELSEIF (igtyp == 8 .OR. igtyp == 13 .OR.
1139 . igtyp == 23 .OR. igtyp == 25) THEN
1140C------
1141 DO k=1,numr
1142 i = i + 1
1143 IF (i > numelr) GOTO 899
1144 728 READ(iin4,fmt='(A)',END=740,ERR=899)line
1145 IF (line(1:1) == '#') GOTO 728
1146 IF (line(1:1) == '/') GOTO 740
1147 READ(line,fmt=fmt_5f) (sigrs(1+j,i),j=1,5)
1148 READ(iin4,fmt='(A)',END=899,ERR=899)line
1149 READ(line,fmt=fmt_5f) (sigrs(6+j,i),j=1,5)
1150 READ(iin4,fmt='(A)',END=899,ERR=899)line
1151 READ(line,fmt=fmt_5f) (sigrs(11+j,i),j=1,5)
1152 READ(iin4,fmt='(A)',END=899,ERR=899)line
1153 READ(line,fmt=fmt_5f) (sigrs(16+j,i),j=1,5)
1154 READ(iin4,fmt='(A)',END=899,ERR=899)line
1155 READ(line,fmt=fmt_5f) (sigrs(21+j,i),j=1,5)
1156 READ(iin4,fmt='(A)',END=899,ERR=899)line
1157 READ(line,fmt=fmt_5f) (sigrs(26+j,i),j=1,5)
1158 READ(iin4,fmt='(A)',END=899,ERR=899)line
1159 READ(line,fmt=fmt_5f) (sigrs(31+j,i),j=1,5)
1160 READ(iin4,fmt='(A)',END=899,ERR=899)line
1161 READ(line,fmt=fmt_5f) (sigrs(36+j,i),j=1,5)
1162 ENDDO ! DO K=1,NUMR
1163C------
1164 ELSEIF (igtyp == 26) THEN
1165C------
1166 DO k=1,numr
1167 i = i + 1
1168 IF (i > numelr) GOTO 899
1169 729 READ(iin4,fmt='(A)',END=740,ERR=899)line
1170 IF (line(1:1) == '#') GOTO 729
1171 IF (line(1:1) == '/') GOTO 740
1172 READ(line,fmt=fmt_3f) (sigrs(1+j,i),j=1,3)
1173 READ(iin4,fmt='(A)',END=899,ERR=899)line
1174 READ(line,fmt=fmt_3f) (sigrs(6+j,i),j=1,3)
1175 ENDDO ! DO K=1,NUMR
1176C------
1177 ELSEIF (igtyp == 29 .OR. igtyp == 30 .OR. igtyp == 31 .OR.
1178 . igtyp == 32 .OR. igtyp == 33 .OR. igtyp == 35 .OR.
1179 . igtyp == 36 .OR. igtyp == 44 .OR. igtyp == 45 .OR.
1180 . igtyp == 46) THEN
1181C
1182 DO k= 1,numr
1183 i = i + 1
1184 IF (i > numelr) GOTO 899
1185 731 READ(iin4,fmt='(A)',END=740,ERR=899)line
1186 IF (line(1:1) == '#') GOTO 731
1187 IF (line(1:1) == '/') GOTO 740
1188 READ(line,fmt=fmt_2f) (sigrs(1+j,i),j=1,2)
1189 READ(iin4,fmt='(A)',END=899,ERR=899)line
1190 READ(line,fmt=fmt_2f) (sigrs(3+j,i),j=1,2)
1191 READ(iin4,fmt='(a)',END=899,ERR=899)LINE
1192 READ(LINE,FMT=FMT_2F) (SIGRS(5+J,I),J=1,2)
1193 READ(IIN4,FMT='(a)',END=899,ERR=899)LINE
1194 READ(LINE,FMT=FMT_2F) (SIGRS(7+J,I),J=1,2)
1195 READ(IIN4,FMT='(a)',END=899,ERR=899)LINE
1196 READ(LINE,FMT=FMT_2F) (SIGRS(9+J,I),J=1,2)
1197 READ(IIN4,FMT='(a)',END=899,ERR=899)LINE
1198 READ(LINE,FMT=FMT_2F) (SIGRS(11+J,I),J=1,2)
1199 READ(IIN4,FMT='(a)',END=899,ERR=899)LINE
1200 READ(LINE,FMT=FMT_F) SIGRS(14,I)
1201C
1202 PT = 14
1203C
1204C -- UVAR --
1205C
1206 NU = (NUVAR - MOD(NUVAR,5))/5
1207 KN = 5
1208 IF (NUVAR < 5) THEN
1209 NU = 1
1210 KN = NUVAR
1211 ENDIF
1212C
1213 DO J1= 1,NU
1214 JJ1 = (J1 -1)*5
1215 READ(IIN4,FMT='(a)',END=899,ERR=899)LINE
1216 READ(LINE,FMT=FMT_5F) (SIGRS(PT+J+JJ1,I),J=1,KN)
1217 ENDDO
1218 IF (NUVAR > 5*NU) THEN
1219 READ(IIN4,FMT='(a)',END=899,ERR=899)LINE
1220 READ(LINE,FMT=FMT_5F) (SIGSH(PT+NU*5+J,I),J=1,NUVAR-5*NU)
1221 ENDIF ! IF (NUVAR > 5*NU)
1222C
1223 ENDDO ! DO K= 1,NUMR
1224 ENDIF ! IF (IGTYP)
1225 GOTO 725
1226C
1227 740 CONTINUE
1228 899 CONTINUE
1229C
1230 RETURN
1231 END
initia
subroutine initia(iparg, elbuf, ms, in, v, x, ixs, ixq, ixc, ixt, ixp, ixr, detonators, geo, pm, rby, npby, lpby, npc, npts, pld, veul, ale_connectivity, skew, fill, ipart, itab, sensors, skvol, ixtg, thk, nloc_dmg, group_param_tab, glob_therm, igrnod, igrsurf, bufsf, vr, bufmat, xlas, las, dtelem, mss, msq, msc, mst, msp, msr, mstg, ptg, inc, nod2eltg, knod2eltg, inp, inr, intg, index, itri, kxx, ixx, xelemwa, iwa, nod2elq, knod2elq, nod2els, knod2els, kxsp, ixsp, nod2sp, ispcond, icode, iskew, iskn, ispsym, xframe, isptag, spbuf, mssx, nsigi, npbyl, lpbyl, rbyl, msnf, mssf, nsigsh, igeo, ipm, nsigs, nsigsph, vns, vnsx, stc, stt, stp, str, sttg, stur, bns, bnsx, volnod, bvolnod, etnod, nshnod, stifint, fxbdep, fxbvit, fxbacc, fxbipm, fxbrpm, fxbelm, fxbsig, fxbmod, ins, ptshel, ptsh3n, ptsol, ptquad, wma, ptsph, fxbnod, mbufel, mdepl, fxani, numel, nsigrs, sh4tree, sh3tree, mcp, temp, imerge2, iadmerge2, slnrbm, nslnrbm, rmstifn, rmstifr, ms_layer, zi_layer, itag, itagel, mcpc, mcptg, xrefc, xreftg, xrefs, mssa, msrt, irbe2, lrbe2, inivol, kvol, nbsubmat, ixs10, ixs16, ixs20, totaddmas, ipmas, stifn, msz2, itagn, sitage, itage, ixr_kj, elbuf_tab, nom_opt, ptr_nopt_rbe2, ptr_nopt_adm, ptr_nopt_fun, sol2sph, irst, sh3trim, xfem_tab, kxig3d, ixig3d, msig3d, knot, nctrlmax, wige, stack, rnoise, drape, sh4ang, sh3ang, geo_stack, igeo_stack, stifintr, strc, strp, strr, strtg, perturb, itagnd, nativ_sms, iloadp, facload, ptspri, nsigbeam, ptbeam, nsigtruss, pttruss, multi_fvm, sigi, sigsh, sigsp, sigsph, sigrs, sigbeam, sigtruss, strsglob, straglob, orthoglob, isigsh, iyldini, ksigsh3, fail_ini, iusolyld, iuser, iddlevel, inimap1d, inimap2d, func2d, fvm_inivel, tagprt_sms, igrbric, igrquad, igrsh4n, igrsh3n, igrpart, totmas, knotlocpc, knotlocel, vnige, bnige, fxbglm, fxbcpm, fxbcps, fxblm, fxbfls, fxbdls, fxb_matrix, fxb_matrix_add, fxb_last_adress, ptr_nopt_fxb, r_skew, knod2el1d, nod2el1d, ebcs_tab, rby_iniaxis, alea, knod2elc, nod2elc, dr, slrbody, drapeg, ipari, intbuf_tab, interfaces, mat_param, npreload_a, preload_a, fail_fractal, fail_brokmann, defaults, ndamp_freq_range, dampr, ibeam_vector, rbeam_vector, ikine)
Definition initia.F:188
lec_inistate
subroutine lec_inistate(ixs, ixq, ixc, ixt, ixp, ixr, geo, pm, kxsp, ixtg, index, itri, nsigsh, igeo, ipm, nsigs, nsigsph, ksysusr, ptshel, ptsh3n, ptsol, ptquad, ptsph, numel, nsigrs, unitab, isolnodd00, lsubmodel, rtrans, idrape, nsigi, ptspri, nsigbeam, ptbeam, nsigtruss, pttruss, sigi, sigsh, sigsp, sigsph, sigrs, sigbeam, sigtruss, strsglob, straglob, orthoglob, isigsh, iyldini, ksigsh3, fail_ini, iusolyld, iuser, igrbric, map_tables, iparg, stack, iworksh, mat_param, numsph, nisp)
Definition lec_inistate.F:49
lec_inistate_yfile
subroutine lec_inistate_yfile(nsigsh, nsigs, nsigsph, nsigrs, nsigi, sigsh, sigi, sigsph, sigrs, sigsp, isigsh, iuser, id_sigsh, id_solid_sigi, id_quad_sigi)
Definition lec_inistate_yfile.F:36
starter
program starter
Definition starter.F:39