OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
lec_inistate_yfile.F File Reference
#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr16_c.inc"
#include "scr17_c.inc"
#include "scry_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine lec_inistate_yfile (nsigsh, nsigs, nsigsph, nsigrs, nsigi, sigsh, sigi, sigsph, sigrs, sigsp, isigsh, iuser, id_sigsh, id_solid_sigi, id_quad_sigi)

Function/Subroutine Documentation

◆ lec_inistate_yfile()

subroutine lec_inistate_yfile ( integer nsigsh,
integer nsigs,
integer nsigsph,
integer nsigrs,
integer nsigi,
sigsh,
sigi,
sigsph,
sigrs,
sigsp,
integer isigsh,
integer iuser,
integer, dimension(*) id_sigsh,
integer, dimension(*) id_solid_sigi,
integer, dimension(*) id_quad_sigi )

Definition at line 31 of file lec_inistate_yfile.F.

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
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21