OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
anim_dcod_key_0.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!|| anim_dcod_key_0 ../engine/source/output/anim/reader/anim_dcod_key_0.F
25!||--- called by ------------------------------------------------------
26!|| freanim ../engine/source/output/anim/reader/freanim.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!||--- uses -----------------------------------------------------
30!|| message_mod ../engine/share/message_module/message_mod.F
31!|| names_and_titles_mod ../common_source/modules/names_and_titles_mod.F
32!|| output_mod ../common_source/modules/output/output_mod.F90
33!|| stack_mod ../engine/share/modules/stack_mod.F
34!||====================================================================
35 SUBROUTINE anim_dcod_key_0(OUTPUT,KEY2,KEY3,KEY4,KEY5,KEY6,JERR,IXITKEY)
36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE message_mod
40 USE stack_mod
41 USE names_and_titles_mod, ONLY: ncharkey
42 USE output_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "scr06_c.inc"
51#include "scr14_c.inc"
52#include "scr25_c.inc"
53#include "sms_c.inc"
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 TYPE(output_), intent(inout) :: OUTPUT
58 INTEGER JERR, IXITKEY
59 CHARACTER(LEN=NCHARKEY) :: KEY2
60 CHARACTER(LEN=NCHARKEY) :: KEY3
61 CHARACTER(LEN=NCHARKEY) :: KEY4
62 CHARACTER(LEN=NCHARKEY) :: KEY5
63 CHARACTER(LEN=NCHARKEY) :: KEY6
64C-----------------------------------------------
65C L o c a l V a r i a b l e s
66C-----------------------------------------------
67 INTEGER I, IDX,IDX0,IDX1,IDX2, N1, N2, N3, ITYP,IADG, J,NTITLE,IUS,ILAY,IPT
68 my_real x0,y0,z0,vnx,vny,vnz,v0
69 LOGICAL :: LAW51_PHASE1, LAW51_PHASE2, LAW51_PHASE3, LAW51_PHASE4
70C-----------------------------------------------
71C Output parameter for LAW51 multiphase data
72C By default set to FALSE
73C Activated in engine file with
74C /ANIM/BRIC(ELEM)/LAW51/ALL
75C /ANIM/BRIC(ELEM)/LAW51/1
76C and so on
77 law51_phase1 = .false.
78 law51_phase2 = .false.
79 law51_phase3 = .false.
80 law51_phase4 = .false.
81C-----------------------------------------------
82 IF(key2(1:4)=='COMP')THEN
83 izip=1
84 ELSEIF(key2(1:4)=='GZIP')THEN
85 izip=2
86 ELSEIF(key2(1:4)=='VERS')THEN
87 READ(key3,'(I2)')anim_vers
88 ELSEIF(key2(1:3)=='MAT')THEN
89 anim_mat = 1
90C
91 ELSEIF(key2(1:4)=='VECT')THEN
92C--------------------------
93C VECTEUR
94C--------------------------
95 IF (key3(1:3)=='VEL')THEN
96 anim_v(1) = 1
97 ELSEIF(key3(1:3)=='DIS')THEN
98 anim_v(2) = 1
99 ELSEIF(key3(1:3)=='ACC')THEN
100 anim_v(3) = 1
101 ELSEIF(key3(1:4)=='CONT' .AND. key3(1:5)/='CONT2'.AND.key4(1:3)/='MAX')THEN
102 anim_v(4) = 1
103 animcont = anim_v(4)
104 ELSEIF(key3(1:4)=='FINT')THEN
105 anim_v(5) = 1
106 ELSEIF(key3(1:4)=='FEXT')THEN
107 anim_v(6) = 1
108 ELSEIF(key3(1:4)=='FOPT')THEN
109 anim_v(7) = 1
110 anim_v(8) = 1
111C rotation velocities
112 ELSEIF(key3(1:4)=='VROT')THEN
113 anim_v(9) = 1
114C fluid velocities for alebag and /bem/flow
115 ELSEIF(key3(1:4)=='FVEL'.OR.key3(1:4)=='VFLU')THEN
116C fluid velocities for alebag and /bem/flow
117 anim_v(10) = 1
118 ELSEIF(key3(1:5)=='FRES')THEN
119 IF ((anim_v(5)/=1.OR.anim_v(6)/=1).AND.idtmins==0.AND.idtmins_int==0)THEN
120 CALL ancmsg(msgid=190,msgtype=msgerror,anmode=aninfo)
121 jerr =jerr+1
122 ixitkey=ixitkey+1
123 ENDIF
124 anim_v(11) = 1
125 ELSEIF(key3(1:5)=='PCONT'.AND. key3(1:6)/='PCONT2')THEN
126 anim_v(12) = 1
127 ELSEIF(key3(1:5)=='CONT2')THEN
128 anim_v(13) = 1
129 ELSEIF(key3(1:4)=='DROT')THEN
130 anim_v(14) = 1
131 ELSEIF(key3(1:5)=='DXANC')THEN
132 anim_v(15) = 1
133 ELSEIF(key3(1:4)=='GVEL')THEN
134 anim_v(16) = 1
135 ELSEIF(key3(1:5)=='FREAC')THEN
136 anim_v(17) = 1
137 ELSEIF(key3(1:5)=='MREAC')THEN
138 anim_v(18) = 1
139 ELSEIF (key3(1:5)=='CLUST')THEN
140 IF(key4(1:5)=='FORCE')THEN
141 anim_v(19) = 1
142 ELSEIF (key4(1:3)=='MOM')THEN
143 anim_v(20) = 1
144 ENDIF
145 ELSEIF (key3(1:4)=='ZVEL')THEN !inter22
146 anim_v(21) = 1
147 ELSEIF (key3(1:5)=='ZFVEL')THEN !inter22
148 anim_v(22) = 1
149 ELSEIF (key3(1:4)=='ZMOM')THEN !inter22
150 anim_v(23) = 1
151 ELSEIF (key3(1:3)=='ZFP')THEN !inter22
152 anim_v(24) = 1
153 ELSEIF (key3(1:5)=='ZFINT')THEN !inter22
154 anim_v(25) = 1
155 ELSEIF (key3(1:4)=='CONT'.AND.key4(1:3)=='MAX')THEN
156 anim_v(26) = 1
157 animcont = anim_v(4)
158 anim_v(4) = 1
159 ELSEIF(key3(1:6)=='PCONT2')THEN
160 anim_v(27) = 1
161 ELSE
162 ixitkey=ixitkey+1
163 ENDIF
164C
165 ELSEIF(key2(1:4)=='TENS')THEN
166C--------------------------
167C TENSEUR OLD
168C--------------------------
169 IF (key3(1:6)=='STRESS')THEN
170 IF (key4(1:4)=='MEMB')THEN
171 anim_t(1) = 1
172 anim_ct(1) = 1
173 ELSEIF(key4(1:4)=='BEND')THEN
174 anim_t(2) = 1
175 anim_ct(2) = 1
176 ELSEIF(key4(1:5)=='UPPER')THEN
177 anim_t(3) = 1
178 anim_ct(3) = 1
179 ELSEIF(key4(1:5)=='LOWER')THEN
180 anim_t(4) = 1
181 anim_ct(4) = 1
182 ELSEIF(key4(1:3)=='ALL')THEN
183 istresall = 1
184 ELSE
185 READ(key4(1:5),'(I5)')n1
186 IF(n1>0.AND.n1<=100)THEN
187 anim_t(100+n1) = 1
188 anim_ct(100+n1) = 1
189 ELSE
190 ixitkey=ixitkey+1
191 ENDIF
192 ENDIF
193
194 ELSEIF(key3(1:5)=='STRAI')THEN
195 IF (key4(1:4)=='MEMB')THEN
196 anim_t(5) = 1
197 anim_ct(5) = 1
198 ELSEIF(key4(1:4)=='BEND')THEN
199 anim_t(6) = 1
200 anim_ct(6) = 1
201 ELSEIF(key4(1:5)=='UPPER')THEN
202 anim_t(7) = 1
203 anim_ct(7) = 1
204 ELSEIF(key4(1:5)=='LOWER')THEN
205 anim_t(8) = 1
206 anim_ct(8) = 1
207 ELSEIF(key4(1:3)=='ALL')THEN
208 istraiall = 1
209 ELSE
210 READ(key4(1:5),'(I5)')n1
211 IF(n1>0.AND.n1<=100)THEN
212 anim_t(200+n1) = 1
213 anim_ct(200+n1) = 1
214 ELSE
215 ixitkey=ixitkey+1
216 ENDIF
217 ENDIF
218 ELSEIF(key3(1:5)=='EPSDO')THEN
219 iepsdot = 1
220 IF (key4(1:4)=='MEMB')THEN
221 anim_t(91) = 1
222 anim_ct(91) = 1
223 ELSEIF(key4(1:4)=='BEND')THEN
224 anim_t(92) = 1
225 anim_ct(92) = 1
226 ELSEIF(key4(1:5)=='UPPER')THEN
227 anim_t(93) = 1
228 anim_ct(93) = 1
229 ELSEIF(key4(1:5)=='LOWER')THEN
230 anim_t(94) = 1
231 anim_ct(94) = 1
232 ELSEIF(key4(1:4)=='ALL')THEN
233 iepsdoall = 1
234 ELSE
235 READ(key4(1:5),'(I5)')n1
236 IF(n1>0.AND.n1<=100)THEN
237 anim_t(300+n1) = 1
238 anim_ct(300+n1) = 1
239 ELSE
240 ixitkey=ixitkey+1
241 ENDIF
242 ENDIF
243 ELSE
244 ixitkey=ixitkey+1
245 ENDIF
246C--------------------------
247C SHELL
248C--------------------------
249 ELSEIF(key2(1:4)=='SHEL')THEN
250C
251 fmt_ani = 4
252C-----------------------------
253C ply shell for composite with xfem ply formulation
254C-------------------------------
255 IF(key3(1:3)=='PLY') THEN
256 anim_ply = 1
257 !READ(KEY4,'(I5)')N1
258 ! IF(N1>0.AND.N1<=200)THEN
259 ! ANIM_PLIES(N1) = 1
260 ! ELSE
261 ! IXITKEY=IXITKEY+1
262 ! ENDIF
263C----------------------------------------------------
264C outp for inter ply composites shell
265C---------------------------------------
266 ELSEIF(key3(1:4)=='IPLY') THEN
267 IF(key4(1:3)=='DAM' .OR. key4(1:5)=='MINDA') THEN
268 anim_ce(10240) = 1
269 ELSEIF(key4(1:5)=='SIGZZ') THEN
270 anim_ce(10241) = 1
271 ELSEIF(key4(1:5)=='SIGYZ') THEN
272 anim_ce(10242) = 1
273 ELSEIF(key4(1:5)=='SIGXZ') THEN
274 anim_ce(10243) = 1
275 ELSEIF(key4(1:5)=='EPSZZ') THEN
276 anim_ce(10244) = 1
277 ELSEIF(key4(1:5)=='EPSYZ') THEN
278 anim_ce(10245) = 1
279 ELSEIF(key4(1:5)=='EPSXZ') THEN
280 anim_ce(10246) = 1
281 ELSEIF(key4(1:4)=='EINT') THEN
282 anim_ce(10247)= 1
283 ELSEIF(key4(1:5)=='MAXDA') THEN
284 anim_ce(10669)= 1
285 ENDIF
286 ELSE IF(key3(1:3)=='CRK') THEN
287 ! XFEM for layered shell (real crack path visualisation)
288C ANIM_CRK = 1
289C--------------------------
290C SHELL:TENSOR
291C--------------------------
292 ELSEIF (key3(1:4) == 'TENS') THEN
293c
294 IF (key4(1:5) == 'STRES') THEN
295 n1 = 0
296 n2 = 0
297 IF (key5(1:4) == 'MEMB') THEN
298 anim_ct(1) = 1
299 ELSEIF (key5(1:4) == 'BEND') THEN
300 anim_ct(2) = 1
301 ELSEIF (key5(1:5) == 'UPPER') THEN
302 anim_ct(3) = 1
303 ELSEIF (key5(1:5) == 'LOWER') THEN
304 anim_ct(4) = 1
305 ELSEIF (key5(1:3) == 'ALL') THEN
306 istresall = 1 ! mean layer value
307 IF (key6(1:3) == 'ALL') istresfull = 2 ! all NPTT
308 ELSE
309
310 READ(key5(1:5),'(I5)') n1
311 IF (n1 > 0 .AND. n1 <= 100) THEN
312 ! case of PID51 (read NPTT, or, UPPER or LOWER within NPTT)
313 IF (key6(1:5) == 'UPPER') THEN
314 anim_ct(400+n1) = 1
315 ELSEIF (key6(1:5) == 'LOWER') THEN
316 anim_ct(500+n1) = 1
317 ELSEIF (key6(1:3) == 'ALL') THEN
318 istresfull = 1
319 anim_stress(n1) = 1
320 ELSE
321 READ(key6(1:5),'(I5)') n2 ! NPTT (max=10) through layer N1 (max=100)
322 IF (n2 > 0 .AND. n2 <= 10) THEN
323 anim_ct(600 + 10*n1 + n2) = 1 ! NPTT through N1
324 ELSE
325 anim_ct(100+n1) = 1 ! classic layer out (mean value)
326 ENDIF ! IF (N2 > 0 .AND. N2 <= 10)
327 ENDIF ! IF (KEY6(1:5) == 'UPPER')
328 ELSE
329 ixitkey=ixitkey+1
330 ENDIF
331 ENDIF
332c---
333 ELSE IF (key4(1:6) == 'MSTRES') THEN
334 ! IGTYP = 1,9,10,11,16
335 ilay = 0
336 ipt = 0
337 idx = 3120 + 3*mx_ply_anim
338 IF (key5(1:4) == 'MEMB') THEN
339 anim_ct(idx+1) = 1
340 ELSEIF (key5(1:5) == 'UPPER') THEN
341 anim_ct(idx+2) = 1
342 ELSEIF (key5(1:5) == 'LOWER') THEN
343 anim_ct(idx+3) = 1
344 ELSE
345 ! /TESNSOR/MSTRESS/ILAY
346 READ(key5(1:5),'(I5)') ilay
347 IF (ilay > 0 .AND. ilay <= 100) THEN
348 anim_ct(idx+ilay) = 1
349 anim_stress(ilay) = 1
350 ELSE
351 ixitkey=ixitkey+1
352 ENDIF
353 ENDIF
354C----
355 ELSEIF (key4(1:5) == 'STRAI') THEN
356 n1 = 0
357 n2 = 0
358 idx = 1610 + 3*mx_ply_anim
359 IF (key5(1:4) == 'MEMB') THEN
360 anim_ct(5) = 1
361 ELSEIF (key5(1:4) == 'BEND') THEN
362 anim_ct(6) = 1
363 ELSEIF (key5(1:5) == 'UPPER') THEN
364 anim_ct(7) = 1
365 ELSEIF (key5(1:5) == 'LOWER') THEN
366 anim_ct(8) = 1
367 ELSEIF (key5(1:3) == 'ALL') THEN
368 istraiall = 1 ! mean layer value
369 IF (key6(1:3) == 'ALL') istrainfull = 2 ! all NPTT
370 ELSE
371 READ(key5(1:5),'(I5)') n1
372 IF (n1 > 0 .AND. n1 <= 100) THEN
373 ! case of PID_51, 52 (read NPTT, or, UPPER or LOWER, within NPTT)
374 IF (key6(1:5) == 'UPPER') THEN
375 anim_ct(idx + n1) = 1
376 ELSEIF (key6(1:5) == 'LOWER') THEN
377 anim_ct(idx + 100 + n1) = 1
378 ELSEIF (key6(1:3) == 'ALL') THEN
379 istrainfull = 1
380 anim_strain(n1) = 1
381 ELSE
382 READ(key6(1:5),'(I5)') n2 ! NPTT (max=10) through layer N1 (max=100)
383 IF (n2 > 0 .AND. n2 <= 10) THEN
384 anim_ct(idx + 200 + 10*n1 + n2) = 1 ! NPTT through N1
385 ELSE
386 anim_ct(200+n1) = 1 ! classic layer out (mean value)
387 ENDIF ! IF (N2 > 0 .AND. N2 <= 10)
388 ENDIF ! IF (KEY6(1:5) == 'UPPER')
389 ELSE
390 ixitkey=ixitkey+1
391 ENDIF ! IF (N1 > 0 .AND. N1 <= 100)
392 ENDIF ! IF (KEY5(1:4)
393C---
394 ELSEIF(key4(1:5)=='EPSDO')THEN
395 iepsdot = 1
396 n1 = 0
397 n2 = 0
398 idx = 2820 + 3*mx_ply_anim
399 IF (key5(1:4) == 'MEMB') THEN
400 anim_ct(91) = 1
401 ELSEIF (key5(1:4) == 'BEND') THEN
402 anim_ct(92) = 1
403 ELSEIF (key5(1:5) == 'UPPER') THEN
404 anim_ct(93) = 1
405 ELSEIF (key5(1:5) == 'LOWER') THEN
406 anim_ct(94) = 1
407 ELSEIF(key5(1:5) == 'ALL') THEN
408 iepsdoall = 1 ! mean layer value
409 IF (key6(1:3) == 'ALL') iepsdofull = 2 ! all NPTT
410 ELSE
411 READ(key5(1:5),'(I5)') n1
412 IF (n1 > 0 .AND. n1 <= 100) THEN
413 ! case of PID_51, 52 (read NPTT, or, UPPER or LOWER, within NPTT)
414 IF (key6(1:5) == 'UPPER') THEN
415 anim_ct(idx + n1) = 1
416 ELSEIF (key6(1:5) == 'LOWER') THEN
417 anim_ct(idx + 100 + n1) = 1
418 ELSEIF (key6(1:3) == 'ALL') THEN
419 iepsdofull = 1
420 anim_epsdot(n1) = 1
421 ELSE
422 READ(key6(1:5),'(I5)') n2 ! NPTT (max=10) through layer N1 (max=100)
423 IF (n2 > 0 .AND. n2 <= 10) THEN
424 anim_ct(idx + 200 + 10*n1 + n2) = 1 ! NPTT through N1
425 ELSE
426 anim_ct(300+n1) = 1 ! classic layer out (mean value)
427 ENDIF ! IF (N2 > 0 .AND. N2 <= 10)
428 ENDIF ! IF (KEY6(1:5) == 'UPPER')
429 ELSE
430 ixitkey=ixitkey+1
431 ENDIF ! IF (N1 > 0 .AND. N1 <= 100)
432 ENDIF ! IF (KEY5(1:4)
433!
434 ENDIF
435C--------------------------
436C COQUES:SCALAIRES
437C--------------------------
438 ELSEIF(key3(1:4)=='EPSP')THEN
439 n1 = 0
440 n2 = 0
441 IF(key4(1:5)=='UPPER')THEN
442 anim_ce(2040)=1
443 ELSEIF(key4(1:5)=='LOWER')THEN
444 anim_ce(2041)=1
445 ELSEIF(key4(1:3)=='ALL')THEN
446 iepspall = 1 ! mean layer value
447 IF(key5(1:3)=='ALL') iepspfull = 2 ! all NPTT
448 ELSE
449 READ(key4(1:5),'(I5)')n1 ! layer
450 IF (n1 > 0 .AND. n1 <= 100) THEN
451 ! case of PID51 (read NPTT, or, UPPER or LOWER within NPTT)
452 IF(key5(1:5)=='UPPER')THEN
453 anim_ce(10677+n1) = 1
454 ELSEIF(key5(1:5)=='LOWER')THEN
455 anim_ce(10777+n1) = 1
456 ELSEIF(key5(1:3)=='ALL')THEN
457 iepspfull = 1
458 anim_epsp(n1) = 1
459 ELSE
460 READ(key5(1:5),'(I5)')n2 ! NPTT (max=10) through layer N1 (max=100)
461 IF (n2 > 0 .AND. n2 <= 10) THEN
462 anim_ce(10877 + 10*n1 + n2) = 1 ! NPTT through N1
463 ELSE
464 anim_ce(2041+n1) = 1 ! classic layer out (mean value)
465 ENDIF
466 ENDIF
467 ELSE
468 anim_ce(1) = 1
469 ENDIF
470 ENDIF
471C-----
472 ELSEIF(key3(1:4)=='DENS')THEN
473 anim_ce(2) = 1
474 ELSEIF(key3(1:4)=='ENER')THEN
475 anim_ce(3) = 1
476 ELSEIF(key3(1:4)=='TEMP')THEN
477 anim_ce(4) = 1
478 ELSEIF(key3(1:4)=='THIC')THEN
479 anim_ce(5) = 1
480 ELSEIF(key3(1:1)=='P' .AND. key3(1:3) /= 'PHI' .AND. key3(1:5) /= 'PLA51')THEN
481 anim_ce(6) = 1
482 ELSEIF(key3(1:4)=='VONM')THEN
483 anim_ce(7) = 1
484 ELSEIF(key3(1:1)=='K')THEN
485 anim_ce(8) = 1
486 ELSEIF(key3(1:4)=='TVIS')THEN
487 anim_ce(9) = 1
488 ELSEIF(key3(1:5)=='VORTX')THEN
489 anim_ce(10) = 1
490 ELSEIF(key3(1:5)=='DAM1')THEN
491 anim_ce(11) = 1
492 ELSEIF(key3(1:5)=='DAM2')THEN
493 anim_ce(12) = 1
494 ELSEIF(key3(1:5)=='DAM3')THEN
495 anim_ce(13) = 1
496 ELSEIF(key3(1:5)=='SIGX ')THEN
497 anim_ce(14) = 1
498 ELSEIF(key3(1:5)=='SIGY ')THEN
499 anim_ce(15) = 1
500 ELSEIF(key3(1:5)=='SIGZ ')THEN
501 anim_ce(16) = 1
502 ELSEIF(key3(1:5)=='SIGXY'.OR.key3(1:5)=='SIGYX')THEN
503 anim_ce(17) = 1
504 ELSEIF(key3(1:5)=='SIGYZ'.OR.key3(1:5)=='sigzy')THEN
505 ANIM_CE(18) = 1
506 ELSEIF(KEY3(1:5)=='sigzx.OR.'KEY3(1:5)=='sigxz')THEN
507 ANIM_CE(19) = 1
508 ELSEIF(KEY3(1:5)=='usr1 .OR.'KEY3(1:5)=='user1')THEN
509 ANIM_CE(20) = 1
510 READ(KEY4(1:5),'(i5)')N1
511 IF(N1/=0)THEN
512 ANIM_CE(39+N1)=1
513 ANIM_CE(20) = 0
514 ENDIF
515 ELSEIF(KEY3(1:5)=='usr2 .OR.'KEY3(1:5)=='user2')THEN
516 ANIM_CE(21) = 1
517 READ(KEY4(1:5),'(i5)')N1
518 IF(N1/=0)THEN
519 ANIM_CE(139+N1)=1
520 ANIM_CE(21) = 0
521 ENDIF
522 ELSEIF(KEY3(1:5)=='usr3 .OR.'KEY3(1:5)=='user3')THEN
523 ANIM_CE(22) = 1
524 READ(KEY4(1:5),'(i5)')N1
525 IF(N1/=0)THEN
526 ANIM_CE(239+N1)=1
527 ANIM_CE(22) = 0
528 ENDIF
529 ELSEIF(KEY3(1:5)=='usr4 .OR.'KEY3(1:5)=='user4')THEN
530 ANIM_CE(23) = 1
531 READ(KEY4(1:5),'(i5)')N1
532 IF(N1/=0)THEN
533 ANIM_CE(339+N1)=1
534 ANIM_CE(23) = 0
535 ENDIF
536 ELSEIF(KEY3(1:5)=='usr5 .OR.'KEY3(1:5)=='user5')THEN
537 ANIM_CE(24) = 1
538 READ(KEY4,'(i5)')N1
539 IF(N1/=0)THEN
540 ANIM_CE(439+N1)=1
541 ANIM_CE(24) = 0
542 ENDIF
543 ELSEIF(KEY3(1:4)=='hour')THEN
544 ANIM_CE(25) = 1
545 ELSEIF (KEY3(1:4)=='epsd')THEN
546 anim_ce(26) = 1
547 ELSEIF(key3(1:5)=='USR6 '.OR.key3(1:5)=='USER6')THEN
548 anim_ce(27) = 1
549 READ(key4(1:5),'(I5)')n1
550 IF(n1/=0)THEN
551 anim_ce(539+n1)=1
552 anim_ce(27) = 0
553 ENDIF
554 ELSEIF(key3(1:5)=='USR7 '.OR.key3(1:5)=='USER7')THEN
555 anim_ce(28) = 1
556 READ(key4(1:5),'(I5)')n1
557 IF(n1/=0)THEN
558 anim_ce(639+n1)=1
559 anim_ce(28) = 0
560 ENDIF
561 ELSEIF(key3(1:5)=='USR8 '.OR.key3(1:5)=='USER8')THEN
562 anim_ce(29) = 1
563 READ(key4(1:5),'(I5)')n1
564 IF(n1/=0)THEN
565 anim_ce(739+n1)=1
566 anim_ce(29) = 0
567 ENDIF
568 ELSEIF(key3(1:5)=='USR9 '.OR.key3(1:5)=='USER9')THEN
569 anim_ce(30) = 1
570 READ(key4(1:5),'(I5)')n1
571 IF(n1/=0)THEN
572 anim_ce(839+n1)=1
573 anim_ce(30) = 0
574 ENDIF
575 ELSEIF(key3(1:5)=='USR10')THEN
576 anim_ce(31) = 1
577 READ(key4(1:5),'(I5)')n1
578 IF(n1/=0)THEN
579 anim_ce(939+n1)=1
580 anim_ce(31) = 0
581 ENDIF
582 ELSEIF(key3(1:5)=='USR11')THEN
583 anim_ce(32) = 1
584 READ(key4(1:5),'(I5)')n1
585 IF(n1/=0)THEN
586 anim_ce(1039+n1)=1
587 anim_ce(32) = 0
588 ENDIF
589 ELSEIF(key3(1:5)=='USR12')THEN
590 anim_ce(33) = 1
591 READ(key4(1:5),'(I5)')n1
592 IF(n1/=0)THEN
593 anim_ce(1139+n1)=1
594 anim_ce(33) = 0
595 ENDIF
596 ELSEIF(key3(1:5)=='USR13')THEN
597 anim_ce(34) = 1
598 READ(key4(1:5),'(I5)')n1
599 IF(n1/=0)THEN
600 anim_ce(1239+n1)=1
601 anim_ce(34) = 0
602 ENDIF
603 ELSEIF(key3(1:5)=='USR14')THEN
604 anim_ce(35) = 1
605 READ(key4(1:5),'(I5)')n1
606 IF(n1/=0)THEN
607 anim_ce(1339+n1)=1
608 anim_ce(35) = 0
609 ENDIF
610 ELSEIF(key3(1:5)=='USR15')THEN
611 anim_ce(36) = 1
612 READ(key4(1:5),'(I5)')n1
613 IF(n1/=0)THEN
614 anim_ce(1439+n1)=1
615 anim_ce(36) = 0
616 ENDIF
617 ELSEIF(key3(1:5)=='usr16')THEN
618 ANIM_CE(37) = 1
619 READ(KEY4(1:5),'(i5)')N1
620 IF(N1/=0)THEN
621 ANIM_CE(1539+N1)=1
622 ANIM_CE(37) = 0
623 ENDIF
624 ELSEIF(KEY3(1:5)=='usr17')THEN
625 ANIM_CE(38) = 1
626 READ(KEY4(1:5),'(i5)')N1
627 IF(N1/=0)THEN
628 ANIM_CE(1639+N1)=1
629 ANIM_CE(38) = 0
630 ENDIF
631 ELSEIF(KEY3(1:5)=='usr18')THEN
632 ANIM_CE(39) = 1
633 READ(KEY4(1:5),'(i5)')N1
634 IF(N1/=0)THEN
635 ANIM_CE(1739+N1)=1
636 ANIM_CE(39)=0
637 ENDIF
638 ELSEIF(KEY3(1:5)=='usr19')THEN
639 READ(KEY4(1:5),'(i5)')N1
640 IF(N1/=0)THEN
641 ANIM_CE(1839+N1)=1
642 ENDIF
643 ELSEIF(KEY3(1:5)=='usr20')THEN
644 READ(KEY4(1:5),'(i5)')N1
645 IF(N1/=0)THEN
646 ANIM_CE(1939+N1)=1
647 ENDIF
648 ELSEIF(KEY3(1:4)=='fail')THEN
649 ANIM_CE(2142) = 1
650 READ(KEY4(1:5),'(i5)')IFAILA
651 IF(IFAILA == 0) IFAILA = 1
652 ELSEIF(KEY3(1:5)=='crmas')THEN
653 ANIM_CE(2143) = 1
654 ELSEIF(KEY3(1:5)=='crvel')THEN
655 ANIM_CE(2144) = 1
656 ELSEIF(KEY3(1:4)=='thin')THEN
657 ANIM_CE(2155) = 1
658 ELSEIF(KEY3(1:5)=='error')THEN
659 IF(KEY4(1:5)=='thick')THEN
660 ANIM_CE(2156) = 1
661 ELSE
662 ENDIF
663 ELSEIF(KEY3(1:5)=='usr21')THEN
664 READ(KEY4(1:5),'(i5)')N1
665 IF(N1/=0)THEN
666 ANIM_CE(2239+N1)=1
667 ENDIF
668 ELSEIF(KEY3(1:5)=='usr22')THEN
669 READ(KEY4(1:5),'(i5)')N1
670 IF(N1/=0)THEN
671 ANIM_CE(2339+N1)=1
672 ENDIF
673 ELSEIF(KEY3(1:5)=='usr23')THEN
674 READ(KEY4(1:5),'(i5)')N1
675 IF(N1/=0)THEN
676 ANIM_CE(2439+N1)=1
677 ENDIF
678 ELSEIF(KEY3(1:5)=='usr24')THEN
679 READ(KEY4(1:5),'(i5)')N1
680 IF(N1/=0)THEN
681 ANIM_CE(2539+N1)=1
682 ENDIF
683 ELSEIF(KEY3(1:5)=='usr25')THEN
684 READ(KEY4(1:5),'(i5)')N1
685 IF(N1/=0)THEN
686 ANIM_CE(2639+N1)=1
687 ENDIF
688 ELSEIF(KEY3(1:5)=='usr26')THEN
689 READ(KEY4(1:5),'(i5)')N1
690 IF(N1/=0)THEN
691 ANIM_CE(2739+N1)=1
692 ENDIF
693 ELSEIF(KEY3(1:5)=='usr27')THEN
694 READ(KEY4(1:5),'(i5)')N1
695 IF(N1/=0)THEN
696 ANIM_CE(2839+N1)=1
697 ENDIF
698 ELSEIF(KEY3(1:5)=='usr28')THEN
699 READ(KEY4(1:5),'(i5)')N1
700 IF(N1/=0)THEN
701 ANIM_CE(2939+N1)=1
702 ENDIF
703 ELSEIF(KEY3(1:5)=='usr29')THEN
704 READ(KEY4(1:5),'(i5)')N1
705 IF(N1/=0)THEN
706 ANIM_CE(3039+N1)=1
707 ENDIF
708 ELSEIF(KEY3(1:5)=='usr30')THEN
709 READ(KEY4,'(i5)')N1
710 IF(N1/=0)THEN
711 ANIM_CE(3139+N1)=1
712 ENDIF
713 ELSEIF(KEY3(1:5)=='usr31')THEN
714 READ(KEY4,'(i5)')N1
715 IF(N1/=0)THEN
716 ANIM_CE(3239+N1)=1
717 ENDIF
718 ELSEIF(KEY3(1:5)=='usr32')THEN
719 READ(KEY4(1:5),'(i5)')N1
720 IF(N1/=0)THEN
721 ANIM_CE(3339+N1)=1
722 ENDIF
723 ELSEIF(KEY3(1:5)=='usr33')THEN
724 READ(KEY4(1:5),'(i5)')N1
725 IF(N1/=0)THEN
726 ANIM_CE(3439+N1)=1
727 ENDIF
728 ELSEIF(KEY3(1:5)=='usr34')THEN
729 READ(KEY4(1:5),'(i5)')N1
730 IF(N1/=0)THEN
731 ANIM_CE(3539+N1)=1
732 ENDIF
733 ELSEIF(KEY3(1:5)=='usr35')THEN
734 READ(KEY4(1:5),'(i5)')n1
735 IF(n1/=0)THEN
736 anim_ce(3639+n1)=1
737 ENDIF
738 ELSEIF(key3(1:5)=='USR36')THEN
739 READ(key4(1:5),'(I5)')n1
740 IF(n1/=0)THEN
741 anim_ce(3739+n1)=1
742 ENDIF
743 ELSEIF(key3(1:5)=='USR37')THEN
744 READ(key4(1:5),'(I5)')n1
745 IF(n1/=0)THEN
746 anim_ce(3839+n1)=1
747 ENDIF
748 ELSEIF(key3(1:5)=='USR38')THEN
749 READ(key4(1:5),'(I5)')n1
750 IF(n1/=0)THEN
751 anim_ce(3939+n1)=1
752 ENDIF
753 ELSEIF(key3(1:5)=='USR39')THEN
754 READ(key4(1:5),'(I5)')n1
755 IF(n1/=0)THEN
756 anim_ce(4039+n1)=1
757 ENDIF
758 ELSEIF(key3(1:5)=='USR40')THEN
759 READ(key4(1:5),'(I5)')n1
760 IF(n1/=0)THEN
761 anim_ce(4139+n1)=1
762 ENDIF
763 ELSEIF(key3(1:5)=='USR41')THEN
764 READ(key4(1:5),'(I5)')n1
765 IF(n1/=0)THEN
766 anim_ce(4239+n1)=1
767 ENDIF
768 ELSEIF(key3(1:5)=='USR42')THEN
769 READ(key4(1:5),'(I5)')n1
770 IF(n1/=0)THEN
771 anim_ce(4339+n1)=1
772 ENDIF
773 ELSEIF(key3(1:5)=='USR43')THEN
774 READ(key4(1:5),'(I5)')n1
775 IF(n1/=0)THEN
776 anim_ce(4439+n1)=1
777 ENDIF
778 ELSEIF(key3(1:5)=='USR44')THEN
779 READ(key4(1:5),'(I5)')n1
780 IF(n1/=0)THEN
781 anim_ce(4539+n1)=1
782 ENDIF
783 ELSEIF(key3(1:5)=='USR45')THEN
784 READ(key4(1:5),'(I5)')n1
785 IF(n1/=0)THEN
786 anim_ce(4639+n1)=1
787 ENDIF
788 ELSEIF(key3(1:5)=='USR46')THEN
789 READ(key4(1:5),'(I5)')n1
790 IF(n1/=0)THEN
791 anim_ce(4739+n1)=1
792 ENDIF
793 ELSEIF(key3(1:5)=='USR47')THEN
794 READ(key4(1:5),'(I5)')n1
795 IF(n1/=0)THEN
796 anim_ce(4839+n1)=1
797 ENDIF
798 ELSEIF(key3(1:5)=='USR48')THEN
799 READ(key4(1:5),'(I5)')n1
800 IF(n1/=0)THEN
801 anim_ce(4939+n1)=1
802 ENDIF
803 ELSEIF(key3(1:5)=='USR49')THEN
804 READ(key4(1:5),'(I5)')n1
805 IF(n1/=0)THEN
806 anim_ce(5039+n1)=1
807 ENDIF
808 ELSEIF(key3(1:5)=='USR50')THEN
809 READ(key4(1:5),'(I5)')n1
810 IF(n1/=0)THEN
811 anim_ce(5139+n1)=1
812 ENDIF
813 ELSEIF(key3(1:5)=='USR51')THEN
814 READ(key4(1:5),'(I5)')n1
815 IF(n1/=0)THEN
816 anim_ce(5239+n1)=1
817 ENDIF
818 ELSEIF(key3(1:5)=='USR52')THEN
819 READ(key4(1:5),'(I5)')n1
820 IF(n1/=0)THEN
821 anim_ce(5339+n1)=1
822 ENDIF
823 ELSEIF(key3(1:5)=='USR53')THEN
824 READ(key4(1:5),'(I5)')n1
825 IF(n1/=0)THEN
826 anim_ce(5439+n1)=1
827 ENDIF
828 ELSEIF(key3(1:5)=='USR54')THEN
829 READ(key4(1:5),'(I5)')n1
830 IF(n1/=0)THEN
831 anim_ce(5539+n1)=1
832 ENDIF
833 ELSEIF(key3(1:5)=='USR55')THEN
834 READ(key4(1:5),'(I5)')n1
835 IF(n1/=0)THEN
836 anim_ce(5639+n1)=1
837 ENDIF
838 ELSEIF(key3(1:5)=='USR56')THEN
839 READ(key4(1:5),'(I5)')n1
840 IF(n1/=0)THEN
841 anim_ce(5739+n1)=1
842 ENDIF
843 ELSEIF(key3(1:5)=='USR57')THEN
844 READ(key4(1:5),'(I5)')n1
845 IF(n1/=0)THEN
846 anim_ce(5839+n1)=1
847 ENDIF
848 ELSEIF(key3(1:5)=='USR58')THEN
849 READ(key4(1:5),'(I5)')n1
850 IF(n1/=0)THEN
851 anim_ce(5939+n1)=1
852 ENDIF
853 ELSEIF(key3(1:5)=='USR59')THEN
854 READ(key4(1:5),'(I5)')n1
855 IF(n1/=0)THEN
856 anim_ce(6039+n1)=1
857 ENDIF
858 ELSEIF(key3(1:5)=='USR60')THEN
859 READ(key4(1:5),'(I5)')n1
860 IF(n1/=0)THEN
861 anim_ce(6139+n1)=1
862 ENDIF
863 ELSEIF(key3(1:5)=='USR61')THEN
864 READ(key4(1:5),'(I5)')n1
865 IF(n1/=0)THEN
866 anim_ce(6239+n1)=1
867 ENDIF
868 ELSEIF(key3(1:5)=='USR62')THEN
869 READ(key4(1:5),'(I5)')n1
870 IF(n1/=0)THEN
871 anim_ce(6339+n1)=1
872 ENDIF
873 ELSEIF(key3(1:5)=='USR63')THEN
874 READ(key4(1:5),'(I5)')n1
875 IF(n1/=0)THEN
876 anim_ce(6439+n1)=1
877 ENDIF
878 ELSEIF(key3(1:5)=='USR64')THEN
879 READ(key4(1:5),'(I5)')n1
880 IF(n1/=0)THEN
881 anim_ce(6539+n1)=1
882 ENDIF
883 ELSEIF(key3(1:5)=='USR65')THEN
884 READ(key4(1:5),'(I5)')n1
885 IF(n1/=0)THEN
886 anim_ce(6639+n1)=1
887 ENDIF
888 ELSEIF(key3(1:5)=='USR66')THEN
889 READ(key4(1:5),'(I5)')n1
890 IF(n1/=0)THEN
891 anim_ce(6739+n1)=1
892 ENDIF
893 ELSEIF(key3(1:5)=='USR67')THEN
894 READ(key4(1:5),'(I5)')n1
895 IF(n1/=0)THEN
896 anim_ce(6839+n1)=1
897 ENDIF
898 ELSEIF(key3(1:5)=='USR68')THEN
899 READ(key4(1:5),'(I5)')n1
900 IF(n1/=0)THEN
901 anim_ce(6939+n1)=1
902 ENDIF
903 ELSEIF(key3(1:5)=='USR69')THEN
904 READ(key4(1:5),'(I5)')n1
905 IF(n1/=0)THEN
906 anim_ce(7039+n1)=1
907 ENDIF
908 ELSEIF(key3(1:5)=='USR70')THEN
909 READ(key4(1:5),'(I5)')n1
910 IF(n1/=0)THEN
911 anim_ce(7139+n1)=1
912 ENDIF
913 ELSEIF(key3(1:5)=='USR71')THEN
914 READ(key4(1:5),'(I5)')n1
915 IF(n1/=0)THEN
916 anim_ce(7239+n1)=1
917 ENDIF
918 ELSEIF(key3(1:5)=='USR72')THEN
919 READ(key4(1:5),'(I5)')n1
920 IF(n1/=0)THEN
921 anim_ce(7339+n1)=1
922 ENDIF
923 ELSEIF(key3(1:5)=='USR73')THEN
924 READ(key4(1:5),'(I5)')n1
925 IF(n1/=0)THEN
926 anim_ce(7439+n1)=1
927 ENDIF
928 ELSEIF(key3(1:5)=='USR74')THEN
929 READ(key4(1:5),'(I5)')n1
930 IF(n1/=0)THEN
931 anim_ce(7539+n1)=1
932 ENDIF
933 ELSEIF(key3(1:5)=='USR75')THEN
934 READ(key4(1:5),'(I5)')n1
935 IF(n1/=0)THEN
936 anim_ce(7639+n1)=1
937 ENDIF
938 ELSEIF(key3(1:5)=='USR76')THEN
939 READ(key4(1:5),'(I5)')n1
940 IF(n1/=0)THEN
941 anim_ce(7739+n1)=1
942 ENDIF
943 ELSEIF(key3(1:5)=='USR77')THEN
944 READ(key4,'(I5)')n1
945 IF(n1/=0)THEN
946 anim_ce(7839+n1)=1
947 ENDIF
948 ELSEIF(key3(1:5)=='USR78')THEN
949 READ(key4(1:5),'(I5)')n1
950 IF(n1/=0)THEN
951 anim_ce(7939+n1)=1
952 ENDIF
953 ELSEIF(key3(1:5)=='USR79')THEN
954 READ(key4(1:5),'(I5)')n1
955 IF(n1/=0)THEN
956 anim_ce(8039+n1)=1
957 ENDIF
958 ELSEIF(key3(1:5)=='USR80')THEN
959 READ(key4(1:5),'(I5)')n1
960 IF(n1/=0)THEN
961 anim_ce(8139+n1)=1
962 ENDIF
963 ELSEIF(key3(1:5)=='USR81')THEN
964 READ(key4(1:5),'(I5)')n1
965 IF(n1/=0)THEN
966 anim_ce(8239+n1)=1
967 ENDIF
968 ELSEIF(key3(1:5)=='USR82')THEN
969 READ(key4(1:5),'(I5)')n1
970 IF(n1/=0)THEN
971 anim_ce(8339+n1)=1
972 ENDIF
973 ELSEIF(key3(1:5)=='USR83')THEN
974 READ(key4(1:5),'(I5)')n1
975 IF(n1/=0)THEN
976 anim_ce(8439+n1)=1
977 ENDIF
978 ELSEIF(key3(1:5)=='USR84')THEN
979 READ(key4,'(I5)')n1
980 IF(n1/=0)THEN
981 anim_ce(8539+n1)=1
982 ENDIF
983 ELSEIF(key3(1:5)=='USR85')THEN
984 READ(key4(1:5),'(I5)')n1
985 IF(n1/=0)THEN
986 anim_ce(8639+n1)=1
987 ENDIF
988 ELSEIF(key3(1:5)=='USR86')THEN
989 READ(key4(1:5),'(I5)')n1
990 IF(n1/=0)THEN
991 anim_ce(8739+n1)=1
992 ENDIF
993 ELSEIF(key3(1:5)=='USR87')THEN
994 READ(key4(1:5),'(I5)')n1
995 IF(n1/=0)THEN
996 anim_ce(8839+n1)=1
997 ENDIF
998 ELSEIF(key3(1:5)=='USR88')THEN
999 READ(key4(1:5),'(I5)')n1
1000 IF(n1/=0)THEN
1001 anim_ce(8939+n1)=1
1002 ENDIF
1003 ELSEIF(key3(1:5)=='USR89')THEN
1004 READ(key4(1:5),'(I5)')n1
1005 IF(n1/=0)THEN
1006 anim_ce(9039+n1)=1
1007 ENDIF
1008 ELSEIF(key3(1:5)=='USR90')THEN
1009 READ(key4(1:5),'(I5)')n1
1010 IF(n1/=0)THEN
1011 anim_ce(9139+n1)=1
1012 ENDIF
1013 ELSEIF(key3(1:5)=='USR91')THEN
1014 READ(key4(1:5),'(I5)')n1
1015 IF(n1/=0)THEN
1016 anim_ce(9239+n1)=1
1017 ENDIF
1018 ELSEIF(key3(1:5)=='USR92')THEN
1019 READ(key4(1:5),'(I5)')n1
1020 IF(n1/=0)THEN
1021 anim_ce(9339+n1)=1
1022 ENDIF
1023 ELSEIF(key3(1:5)=='USR93')THEN
1024 READ(key4(1:5),'(I5)')n1
1025 IF(n1/=0)THEN
1026 anim_ce(9439+n1)=1
1027 ENDIF
1028 ELSEIF(key3(1:5)=='USR94')THEN
1029 READ(key4(1:5),'(I5)')n1
1030 IF(n1/=0)THEN
1031 anim_ce(9539+n1)=1
1032 ENDIF
1033 ELSEIF(key3(1:5)=='USR95')THEN
1034 READ(key4(1:5),'(I5)')n1
1035 IF(n1/=0)THEN
1036 anim_ce(9639+n1)=1
1037 ENDIF
1038 ELSEIF(key3(1:5)=='USR96')THEN
1039 READ(key4(1:5),'(I5)')n1
1040 IF(n1/=0)THEN
1041 anim_ce(9739+n1)=1
1042 ENDIF
1043 ELSEIF(key3(1:5)=='USR97')THEN
1044 READ(key4(1:5),'(I5)')n1
1045 IF(n1/=0)THEN
1046 anim_ce(9839+n1)=1
1047 ENDIF
1048 ELSEIF(key3(1:5)=='USR98')THEN
1049 READ(key4(1:5),'(I5)')n1
1050 IF(n1/=0)THEN
1051 anim_ce(9939+n1)=1
1052 ENDIF
1053 ELSEIF(key3(1:5)=='USR99')THEN
1054 READ(key4(1:5),'(I5)')n1
1055 IF(n1/=0)THEN
1056 anim_ce(10039+n1)=1
1057 ENDIF
1058 ELSEIF(key3(1:3)=='PHI')THEN
1059 IF(key4(1:3)=='ALL')THEN
1060 iphiall = 1
1061 ELSEIF(key4(1:4)=='memb')THEN
1062 ANIM_CE(10673) = 1
1063 ELSEIF(KEY4(1:5)=='upper')THEN
1064 ANIM_CE(10674) = 1
1065 ELSEIF(KEY4(1:5)=='lower')THEN
1066 anim_ce(10675) = 1
1067 ELSE
1068 READ(key4(1:5),'(I5)')n1
1069 IF(n1/=0)THEN
1070 anim_ce(10139+n1)=1
1071 ENDIF
1072 ENDIF
1073 ELSEIF(key3(1:4)=='DAMA')THEN
1074 n1 = 0
1075 n2 = 0
1076 idx = 11925 + 4*mx_ply_anim + 6
1077 ! in addition DAMA for PID_51 and PID_52 - shell composite -
1078 IF(key4(1:5)=='UPPER')THEN
1079 anim_ce(10257)=1
1080 ELSEIF(key4(1:5)=='LOWER')THEN
1081 anim_ce(10258)=1
1082 ELSEIF(key4(1:4)=='MEMB')THEN
1083 anim_ce(10259)=1
1084 ELSEIF(key4(1:3)=='ALL')THEN
1085 idamaall = 1 ! mean layer value
1086 IF(key5(1:3)=='ALL') idamafull = 2 ! all NPTT
1087 ELSE
1088 READ(key4(1:5),'(I5)')n1 ! layer
1089 IF (n1 > 0 .AND. n1 <= 100) THEN
1090 ! case of PID_51, 52 (read NPTT, or, UPPER or LOWER, MEMB within NPTT)
1091 IF (key5(1:5) == 'UPPER') THEN
1092 anim_ce(idx + n1) = 1
1093 ELSEIF (key5(1:5) == 'LOWER') THEN
1094 anim_ce(idx + 100 + n1) = 1
1095 ELSEIF (key5(1:4) == 'MEMB') THEN
1096 anim_ce(idx + 200 + n1) = 1
1097 ELSEIF(key5(1:3) == 'ALL')THEN
1098 idamafull = 1
1099 anim_dama(n1) = 1
1100 ELSE
1101 READ(key5(1:5),'(I5)')n2 ! NPTT (max=10) through layer N1 (max=100)
1102 IF (n2 > 0 .AND. n2 <= 10) THEN
1103 anim_ce(idx + 300 + 10*n1 + n2) = 1 ! NPTT through N1
1104 ELSE
1105 anim_ce(10259+n1) = 1 ! classic layer out (mean value)
1106 ENDIF
1107 ENDIF
1108 ELSE
1109 anim_ce(10256) = 1
1110 ENDIF
1111 ENDIF ! IF(KEY4=='UPPER')THEN
1112 ELSEIF(key3(1:4)=='NXTF')THEN
1113 IF(key4(1:5)=='UPPER')THEN
1114 anim_ce(10360)=1
1115 ELSEIF(key4(1:5)=='LOWER')THEN
1116 anim_ce(10361)=1
1117 ELSEIF(key4(1:4)=='MEMB')THEN
1118 anim_ce(10362)=1
1119 ELSEIF(key4(1:3)=='ALL')THEN
1120 inxtfall = 1
1121 ELSE
1122 READ(key4(1:5),'(I5)')n1
1123 IF (n1/=0) THEN
1124 anim_ce(10362+n1) = 1
1125 ENDIF
1126 ENDIF
1127 ELSEIF(key3(1:5)=='SIG1H')THEN
1128 IF(key4(1:5)=='UPPER')THEN
1129 anim_ce(10463)=1
1130 ELSEIF(key4(1:5)=='LOWER')THEN
1131 anim_ce(10464)=1
1132 ELSEIF(key4(1:4)=='MEMB')THEN
1133 anim_ce(10465)=1
1134 ELSEIF(key4(1:3)=='ALL')THEN
1135 sigh1all = 1
1136 ELSE
1137 READ(key4(1:5),'(I5)')n1
1138 IF (n1/=0) THEN
1139 anim_ce(10465+n1) = 1
1140 ENDIF
1141 ENDIF
1142 ELSEIF(key3(1:5)=='SIG2H')THEN
1143 IF(key4(1:5)=='UPPER')THEN
1144 anim_ce(10566)=1
1145 ELSEIF(key4(1:5)=='LOWER')THEN
1146 anim_ce(10567)=1
1147 ELSEIF(key4(1:4)=='MEMB')THEN
1148 anim_ce(10568)=1
1149 ELSEIF(key4(1:3)=='ALL')THEN
1150 sigh2all = 1
1151 ELSE
1152 READ(key4(1:5),'(I5)')n1
1153 IF (n1/=0) THEN
1154 anim_ce(10568+n1) = 1
1155 ENDIF
1156 ENDIF
1157 ELSEIF(key3(1:5)=='TDEL')THEN
1158 anim_ce(10670) = 1
1159 ELSEIF(key3(1:3)=='SSP')THEN
1160 anim_ce(10671) = 1
1161 ELSEIF(key3(1:5)=='SCHLI')THEN !.OR. 'SCHLIEREN'
1162 anim_ce(10672) = 1
1163 ELSEIF(key3(1:5)=='DOMAI')THEN
1164 anim_ce(10676) = 1
1165 !C element equivalent stress criteria
1166 ELSEIF(key3(1:5)=='SIGEQ')THEN
1167 anim_ce(10677) = 1
1168 ELSEIF(key3(1:4)=='BULK' .OR. key3(1:4)=='QVIS')THEN
1169 anim_ce(11888) = 1 !PID15 UNIM_CE (10678: 11887)
1170 ELSEIF(key3(1:4)=='TDET')THEN
1171 anim_ce(11889) = 1
1172 ELSEIF(key3(1:5)=='LAW20')THEN
1173 anim_ce(11890) = 1 !dens
1174 anim_ce(11891) = 1
1175 anim_ce(11892) = 0
1176 anim_ce(11893) = 0
1177 anim_ce(11894) = 1 !ener
1178 anim_ce(11895) = 1
1179 anim_ce(11896) = 0
1180 anim_ce(11897) = 0
1181 anim_ce(11898) = 1 !temp
1182 anim_ce(11899) = 1
1183 anim_ce(11900) = 0
1184 anim_ce(11901) = 0
1185 anim_ce(11902) = 1 !pres
1186 anim_ce(11903) = 1
1187 anim_ce(11904) = 0
1188 anim_ce(11905) = 0
1189 anim_ce(11906) = 1 !plas
1190 anim_ce(11907) = 1
1191 anim_ce(11908) = 0
1192 anim_ce(11909) = 0
1193 anim_ce(11910) = 1 !ssp
1194 anim_ce(11911) = 1
1195 anim_ce(11912) = 0
1196 anim_ce(11913) = 0
1197 anim_ce(11914) = 1 !volume
1198 anim_ce(11915) = 1
1199 anim_ce(11916) = 0
1200 anim_ce(11917) = 0
1201 anim_ce(11918) = 1 !mass
1202 anim_ce(11919) = 1
1203 anim_ce(11920) = 0
1204 anim_ce(11921) = 0
1205 anim_ce(11922) = 1 !Qvis
1206 anim_ce(11923) = 1
1207 anim_ce(11924) = 0
1208 anim_ce(11925) = 0
1209
1210 ELSEIF(key3(1:5)=='IDPLY')THEN
1211 !ANIM_CE(11925) -> ANIM_CE(11925+3*MX_PLY_ANIM)
1212 !ANIM_CT(1610) -> ANIM_CT(1610+2*MX_PLY_ANIM)
1213 ! === /ANIM/SHELL/IDPLY/STRESS ===
1214 IF (key4(1:5)=='STRES') THEN
1215 IF(key5(1:3) == 'ALL')THEN
1216 istresall_ply = 1
1217 ELSE
1218 READ(key5(1:5),'(I5)')n1
1219 READ(key6(1:5),'(I5)')ipt
1220 nbply_anim_stress = nbply_anim_stress + 1
1221 ply_anim_stress(3*(nbply_anim_stress-1)+2) = 2
1222 ply_anim_stress(3*(nbply_anim_stress-1)+1) = n1
1223 ply_anim_stress(3*(nbply_anim_stress-1)+3) = ipt
1224 anim_ct(1610 + nbply_anim_stress) = 1
1225 IF(nbply_anim_stress > mx_ply_anim)
1226 . CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1227 ENDIF
1228c
1229 ELSE IF (key4(1:6)=='MSTRES') THEN
1230 READ(key5(1:5),'(I5)') ilay
1231 READ(key6(1:5),'(I5)') ipt
1232c
1233 idx = 3120 + 3*mx_ply_anim + 103
1234 nbply_anim_stress = nbply_anim_stress + 1
1235 ply_anim_stress(3*(nbply_anim_stress-1)+2) = 3
1236 ply_anim_stress(3*(nbply_anim_stress-1)+1) = ilay
1237 ply_anim_stress(3*(nbply_anim_stress-1)+3) = ipt
1238 anim_ct(idx + nbply_anim_stress) = 1
1239 IF (nbply_anim_stress > mx_ply_anim) THEN
1240 CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1241 ENDIF
1242c
1243 ! /ANIM/SHELL/IDPLY/STRAIN
1244 ELSEIF(key4(1:5)=='STRAI')THEN
1245 IF(key5(1:3) == 'ALL')THEN
1246 istrainall_ply = 1
1247 ELSE
1248 READ(key5(1:5),'(I5)')n1
1249 READ(key6(1:5),'(I5)')ipt
1250 nbply_anim_strain = nbply_anim_strain + 1
1251 ply_anim_strain(3*(nbply_anim_strain-1)+1) = n1
1252 ply_anim_strain(3*(nbply_anim_strain-1)+2) = 3
1253 ply_anim_strain(3*(nbply_anim_strain-1)+3) = ipt
1254 anim_ct( (1610+ mx_ply_anim) + nbply_anim_strain) = 1
1255 IF(nbply_anim_strain > mx_ply_anim)
1256 . CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1257 ENDIF
1258 ! /ANIM/SHELL/IDPLY/EPSDOT
1259 ELSEIF(key4(1:5)=='EPSDO')THEN
1260 iepsdot = 1
1261 IF(key5(1:3) == 'ALL')THEN
1262 iepsdotall_ply = 1
1263 ELSE
1264 READ(key5(1:5),'(I5)')n1
1265 READ(key6(1:5),'(I5)')ipt
1266 nbply_anim_epsdot = nbply_anim_epsdot + 1
1267 ply_anim_epsdot(3*(nbply_anim_epsdot-1)+1) = n1
1268 ply_anim_epsdot(3*(nbply_anim_epsdot-1)+2) = 6
1269 ply_anim_epsdot(3*(nbply_anim_epsdot-1)+3) = ipt
1270 anim_ct( (1610+ 2*mx_ply_anim) + nbply_anim_epsdot) = 1
1271 IF(nbply_anim_epsdot > mx_ply_anim)
1272 . CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1273 ENDIF
1274 ELSEIF(key4(1:3)=='PHI')THEN
1275 IF(key5(1:3) == 'ALL')THEN
1276 ! /ANIM/SHELL/IDPLY/PHI/ALL
1277 iphiall_ply = 1
1278 ELSE
1279 ! /ANIM/SHELL/IDPLY/PHI/Id
1280 READ(key5(1:5),'(I5)')n1
1281 nbply_anim_phi = nbply_anim_phi + 1
1282 IF(nbply_anim_phi > mx_ply_anim)
1283 . CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1284 ply_anim_phi(3*(nbply_anim_phi-1)+1) = n1
1285 ply_anim_phi(3*(nbply_anim_phi-1)+2) = 4
1286 ply_anim_phi(3*(nbply_anim_phi-1)+3) = 0
1287 anim_ce( (11925 + mx_ply_anim) + nbply_anim_phi ) = 1
1288 ENDIF
1289 ELSEIF(key4(1:4)=='EPSP')THEN
1290 IF(key5(1:3) == 'ALL')THEN
1291 ! /ANIM/SHELL/IDPLY/EPSP/ALL
1292 iepspall_ply = 1
1293 ELSE
1294 ! /ANIM/SHELL/IDPLY/EPSP/Id_ply_id_point
1295 READ(key5(1:5),'(I5)')n1
1296 READ(key6(1:5),'(I5)')ipt
1297 nbply_anim_epsp = nbply_anim_epsp + 1
1298 IF(nbply_anim_epsp > mx_ply_anim)
1299 . CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1300 ply_anim_epsp(3*(nbply_anim_epsp-1)+1) = n1
1301 ply_anim_epsp(3*(nbply_anim_epsp-1)+2) = 5
1302 ply_anim_epsp(3*(nbply_anim_epsp-1)+3) = ipt
1303 anim_ce( (11925 + 2 * mx_ply_anim) + nbply_anim_epsp ) = 1
1304 ENDIF
1305 ELSEIF(key4(1:4)=='DAMA')THEN
1306 IF(key5(1:3) == 'ALL')THEN
1307 ! /ANIM/SHELL/IDPLY/DAMA/ALL
1308 idamaall_ply = 1
1309 ELSE
1310 ! /ANIM/SHELL/IDPLY/DAMA/Id_ply_id_point
1311 READ(key5(1:5),'(I5)')n1
1312 READ(key6(1:5),'(I5)')ipt
1313 nbply_anim_dama = nbply_anim_dama + 1
1314 IF(nbply_anim_dama > mx_ply_anim)
1315 . CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1316 ply_anim_dama(3*(nbply_anim_dama-1)+1) = n1
1317 ply_anim_dama(3*(nbply_anim_dama-1)+2) = 7
1318 ply_anim_dama(3*(nbply_anim_dama-1)+3) = ipt
1319 anim_ce( (11925 + 3 * mx_ply_anim) + nbply_anim_dama ) = 1
1320 ENDIF
1321 ELSEIF(key4(1:3)=='ALL')THEN
1322 ! /ANIM/SHELL/IDPLY/ALL
1323 iplyall = 1
1324 ELSE
1325 ! /ANIM/SHELL/IDPLY/Id
1326 READ(key4(1:5),'(I5)')n1
1327 nbply_anim = nbply_anim + 1
1328 IF(nbply_anim > mx_ply_anim)
1329 . CALL ancmsg(msgid=268,msgtype=msgwarning,i1=mx_ply_anim,anmode=aninfo)
1330 ply_anim(3*(nbply_anim-1)+1) = n1
1331 ply_anim(3*(nbply_anim-1)+2) = 1
1332 ply_anim(3*(nbply_anim-1)+3) = 0
1333 anim_ce(11925 + nbply_anim) = 1
1334 ENDIF
1335 ELSEIF(key3(1:4) == 'FLDF')THEN
1336 idx = 11925+4*mx_ply_anim
1337 IF (key4(1:5) == 'UPPER') THEN
1338 anim_ce(idx+1)=1
1339 ELSEIF (key4(1:5) == 'LOWER') THEN
1340 anim_ce(idx+2)=1
1341 ELSEIF (key4(1:4) == 'MEMB') THEN
1342 anim_ce(idx+3)=1
1343 ENDIF
1344 ELSEIF(key3(1:4) == 'FLDZ')THEN
1345 idx = 11925+4*mx_ply_anim + 3
1346 IF (key4(1:5) == 'UPPER') THEN
1347 anim_ce(idx+1)=1
1348 ELSEIF (key4(1:5) == 'LOWER') THEN
1349 anim_ce(idx+2)=1
1350 ELSEIF (key4(1:4) == 'MEMB') THEN
1351 anim_ce(idx+3)=1
1352 ENDIF
1353 ELSEIF(key3(1:5) == 'DT ')THEN
1354 idx = 13242 + 4*mx_ply_anim
1355 anim_ce(idx) = 1
1356 ELSEIF(key3(1:3) == 'AMS')THEN
1357 idx = 13242 + 4*mx_ply_anim + 1
1358 anim_ce(idx) = 1
1359 !NIM/SHELL/WPLA/KEY4/KEY5
1360 ELSEIF(key3(1:4) == 'WPLA')THEN
1361 n1 = 0
1362 n2 = 0
1363 idx = 13245 + 4*mx_ply_anim
1364 IF(key4(1:5)=='UPPER')THEN
1365 anim_ce(idx + 1)=1
1366 ELSEIF(key4(1:5)=='LOWER')THEN
1367 anim_ce(idx + 2)=1
1368 ELSEIF(key4(1:5)=='ALL')THEN
1369 iwplaall = 1 ! mean layer value
1370 IF(key5(1:3)=='ALL') iwplafull = 2 ! all NPTT
1371 ELSE
1372 READ(key4(1:5),'(I5)')n1 ! layer
1373 IF (n1 > 0 .AND. n1 <= 100) THEN
1374 ! case of PID51 (read NPTT, or, UPPER or LOWER within NPTT)
1375 idx = 13347 + 4*mx_ply_anim
1376 IF(key5(1:5)=='UPPER')THEN
1377 anim_ce(idx + n1) = 1
1378 ELSEIF(key5(1:5)=='LOWER')THEN
1379 anim_ce(idx + 100 + n1) = 1
1380 ELSEIF(key5(1:3)=='ALL')THEN
1381 iwplafull = 1
1382 anim_wpla(n1) = 1
1383 ELSE
1384 idx = 13547 + 4*mx_ply_anim
1385 READ(key5(1:5),'(I5)')n2 ! NPTT (max=10) through layer N1 (max=100)
1386 IF (n2 > 0 .AND. n2 <= 10) THEN
1387 anim_ce(idx + 10*(n1-1) + n2) = 1 ! NPTT through N1
1388 ELSE
1389 idx = 13247 + 4*mx_ply_anim !!! N1 <= 100
1390 anim_ce(idx + n1) = 1 ! classic layer out (mean value)
1391 ENDIF
1392 ENDIF
1393 ELSE
1394 idx = 13245 + 4*mx_ply_anim
1395 anim_ce(idx) = 1
1396 ENDIF
1397 ENDIF
1398 ! IDX = 13547 + 4*MX_PLY_ANIM + 10*100 limit
1399 ELSEIF(key3(1:3) == 'OFF')THEN
1400 idx = 13547 + 4*mx_ply_anim + 1000 + 1
1401 anim_ce(idx) = 1
1402 ELSEIF(key3(1:4) == 'MACH')THEN
1403 idx = 13547 + 4*mx_ply_anim + 1000 + 2
1404 anim_ce(idx) = 1
1405 ELSEIF(key3(1:5) == 'COLOR')THEN
1406 idx = 13547 + 4*mx_ply_anim + 1000 + 3
1407 anim_ce(idx) = 1
1408 ELSEIF(key3(1:4) == 'DAMG') THEN
1409 idx = 13547 + 4*mx_ply_anim + 1000 + 4
1410 n1 = 0
1411 IF(key4(1:4)=='MEAN')THEN
1412 anim_ce(idx) = 1
1413 ELSEIF(key4(1:5)=='UPPER')THEN
1414 anim_ce(idx+1)=1
1415 ELSEIF(key4(1:5)=='LOWER')THEN
1416 anim_ce(idx+2)=1
1417 ELSEIF(key4(1:4)=='MEMB')THEN
1418 anim_ce(idx+3)=1
1419 ELSE
1420 READ(key4(1:5),'(I5)')n1 ! NPTT
1421 IF (n1 > 0 .AND. n1 <= 11) THEN
1422 anim_ce(idx+3+n1)=1
1423 ELSE
1424 anim_ce(idx) = 1
1425 ENDIF
1426 ENDIF
1427 ELSEIF(key3(1:7)=='NL_EPSP')THEN
1428 idx = 4*mx_ply_anim + 14567
1429 n1 = 0
1430 IF(key4(1:5)=='UPPER')THEN
1431 anim_ce(idx+1)=1
1432 ELSEIF(key4(1:5)=='LOWER')THEN
1433 anim_ce(idx+2)=1
1434 ELSEIF(key4(1:3)=='ALL')THEN
1435 iepspnlall = 1
1436 ELSE
1437 READ(key4(1:5),'(I5)') n1
1438 IF (n1 > 0 .AND. n1 <= 11) THEN
1439 anim_ce(idx+2+n1) = 1
1440 ELSE
1441 anim_ce(idx) = 1
1442 ENDIF
1443 ENDIF
1444 ELSEIF(key3(1:7)=='NL_EPSD')THEN
1445 idx = 4*mx_ply_anim + 14581
1446 n1 = 0
1447 IF(key4(1:5)=='UPPER')THEN
1448 anim_ce(idx+1)=1
1449 ELSEIF(key4(1:5)=='LOWER')THEN
1450 anim_ce(idx+2)=1
1451 ELSEIF(key4(1:3)=='ALL')THEN
1452 iepsdnlall = 1
1453 ELSE
1454 READ(key4(1:5),'(I5)') n1
1455 IF (n1 > 0 .AND. n1 <= 11) THEN
1456 anim_ce(idx+2+n1) = 1
1457 ELSE
1458 anim_ce(idx) = 1
1459 ENDIF
1460 ENDIF
1461 !ANIM/SHELL/TSAIWU/KEY4/KEY5
1462 ELSEIF (key3(1:6) == 'TSAIWU') THEN
1463 n1 = 0
1464 n2 = 0
1465 idx = 14595 + 4*mx_ply_anim
1466 IF(key4(1:5)=='UPPER')THEN
1467 anim_ce(idx + 1)=1
1468 ELSEIF(key4(1:5)=='LOWER')THEN
1469 anim_ce(idx + 2)=1
1470 ELSEIF(key4(1:5)=='ALL')THEN
1471 itsaiwuall = 1 ! mean layer value
1472 IF(key5(1:3)=='ALL') itsaiwuall = 2 ! all NPTT
1473 ELSE
1474 READ(key4(1:5),'(I5)') n1 ! layer
1475 IF (n1 > 0 .AND. n1 <= 100) THEN
1476 ! case of PID51 (read NPTT, or, UPPER or LOWER within NPTT)
1477 idx = 14697 + 4*mx_ply_anim
1478 IF (key5(1:5)=='UPPER') THEN
1479 anim_ce(idx + n1) = 1
1480 ELSEIF(key5(1:5)=='LOWER') THEN
1481 anim_ce(idx + 100 + n1) = 1
1482 ELSEIF(key5(1:3)=='ALL')THEN
1483 itsaiwufull = n1
1484 ELSE
1485 idx = 14897 + 4*mx_ply_anim
1486 READ(key5(1:5),'(I5)') n2 ! NPTT (max=10) through layer N1 (max=100)
1487 IF (n2 > 0 .AND. n2 <= 10) THEN
1488 anim_ce(idx + 10*(n1-1) + n2) = 1 ! NPTT through N1
1489 ELSE
1490 idx = 14597 + 4*mx_ply_anim ! N1 <= 100
1491 anim_ce(idx + n1) = 1 ! classic layer out (mean value)
1492 ENDIF
1493 ENDIF
1494 ELSE
1495 idx = 14595 + 4*mx_ply_anim
1496 anim_ce(idx) = 1
1497 ENDIF
1498 ENDIF
1499 ELSEIF(key3(1:9) == 'TILLOTSON')THEN
1500 idx = 15898 + 4*mx_ply_anim
1501 anim_ce(idx) = 1
1502 ELSEIF(key3(1:7) == 'VSTRAIN')THEN
1503 idx = 15898 + 4*mx_ply_anim
1504 read(key4, '(I2)') ilay
1505 if(ilay >= 1 .and. ilay <= 10)then
1506 anim_ce(idx+ilay) = 1
1507 else
1508 anim_ce(idx) = 1
1509 endif
1510 !ELSEIF(KEY3(:) == '')THEN
1511 ! IDX = 15022 + 4*MX_PLY_ANIM ?
1512 ! ANIM_CE(IDX) = 1
1513 ELSE
1514 ixitkey=ixitkey+1
1515c
1516 ENDIF ! KEY3 SHELLS
1517c
1518C--------------------------
1519C BRIQUES
1520C--------------------------
1521 ELSEIF(key2(1:4)=='BRIC')THEN
1522C
1523 fmt_ani = 4
1524C--------------------------
1525C BRIQUES:TENSEUR
1526C--------------------------
1527 IF(key3(1:4)=='TENS')THEN !ANIM/BRICK/TENS
1528 IF(key4(1:5)=='STRES')THEN
1529 IF(key5(1:3) == 'ALL')THEN
1530 ibrick_stressall = 1
1531 ELSE
1532 READ(key5(1:5),'(I5)')n1
1533 IF(n1/=0)THEN
1534 IF ( mod(abs(n1)/10,10) == 0) THEN
1535 READ(key6(1:5),'(I5)')n2
1536 anim_st(2010+mod(abs(n1),100)+
1537 . n2*10+(abs(n1)/100)*2010)=1
1538 ELSE
1539 anim_st(10+n1)=1
1540 ENDIF
1541 ELSE
1542 anim_st(1) = 1
1543 ENDIF
1544 ENDIF
1545 ELSEIF(key4(1:5)=='STRAI')THEN
1546 IF(key5(1:3) == 'ALL')THEN
1547 ibrick_strainall = 1
1548 ELSE
1549 READ(key5(1:5),'(I5)')n1
1550 IF(n1/=0)THEN
1551 IF ( mod(abs(n1)/10,10) == 0) THEN
1552 READ(key6(1:5),'(I5)')n2
1553 anim_st(22110+mod(abs(n1),100)+
1554 . n2*10+(abs(n1)/100)*2010)=1
1555 ELSE
1556 anim_st(1010+n1)=1
1557 ENDIF
1558 ELSE
1559 anim_st(2) = 1
1560 ENDIF
1561 ENDIF
1562 ELSEIF(key4(1:5)=='EPSDO')THEN
1563C IEPSDOT = 1
1564C ANIM_ST(3) = 1
1565 ELSEIF(key4(1:4)=='DAMA')THEN
1566 anim_st(4) = 1
1567 ELSEIF(key4(1:4)=='EPSP')THEN
1568 IF(key5(1:3) == 'ALL')THEN
1569 ibrick_epspall = 1
1570 ELSE
1571 READ(key5(1:5),'(I5)')n1
1572 IF(n1/=0)THEN
1573 IF ( mod(abs(n1)/10,10) == 0) THEN
1574 READ(key6(1:5),'(I5)')n2
1575 anim_st(43210+mod(abs(n1),100)+
1576 . n2*10+(abs(n1)/100)*2010)=1
1577 ELSE
1578 anim_st(42210+n1)=1
1579 ENDIF
1580 ELSE
1581 anim_st(5) = 1
1582 ENDIF
1583 ENDIF
1584 ELSE
1585 ixitkey=ixitkey+1
1586 ENDIF
1587 ELSE IF (key3(1:5) == 'LAW51') THEN
1588 IF (key4(1:3) == 'ALL') THEN
1589 law51_phase1 = .true.
1590 law51_phase2 = .true.
1591 law51_phase3 = .true.
1592 law51_phase4 = .true.
1593 ELSE IF (key4(1:1) == '1') THEN
1594 law51_phase1 = .true.
1595 ELSE IF (key4(1:1) == '2') THEN
1596 law51_phase2 = .true.
1597 ELSE IF (key4(1:1) == '3') THEN
1598 law51_phase3 = .true.
1599 ELSE IF (key4(1:1) == '4') THEN
1600 law51_phase4 = .true.
1601 ENDIF
1602C--------------------------
1603C BRIQUES:SCALAIRES
1604C--------------------------
1605 ELSEIF(key3(1:4)=='EPSP')THEN
1606 anim_se(1) = 1
1607 IF (law51_phase1) anim_se(4913) = 1 !EPSP_law51-phase1
1608 IF (law51_phase2) anim_se(4914) = 1 !EPSP_law51-phase2
1609 IF (law51_phase3) anim_se(4915) = 1 !EPSP_law51-phase3
1610 IF (law51_phase4) anim_se(4916) = 1 !EPSP_law51-phase4
1611 ELSEIF(key3(1:5)=='PLA51')THEN !obsolete
1612 anim_se(1) = 1
1613 ELSEIF(key3(1:5)=='DENS3')THEN
1614 !/ANIM/ELEM/DENS37
1615 anim_se(4935) = 1 !law37 -rho_liquid
1616 anim_se(4936) = 1 !law37 - rho_gas
1617 ELSEIF(key3(1:4)=='DENS')THEN
1618 anim_se(2) = 1
1619 IF (law51_phase1) anim_se(4897) = 1 !DENS_law51-phase1
1620 IF (law51_phase2) anim_se(4898) = 1 !DENS_law51-phase2
1621 IF (law51_phase3) anim_se(4899) = 1 !DENS_law51-phase3
1622 IF (law51_phase4) anim_se(4900) = 1 !DENS_law51-phase4
1623 !Anim_se (4935) = 1! Law37 -rho_liquid
1624 !ANIM_SE(4936) = 1 !law37 - rho_gas
1625 ELSEIF(key3(1:4)=='ENER')THEN
1626 anim_se(3) = 1
1627 IF (law51_phase1) anim_se(4901) = 1 !ENER_law51-phase1
1628 IF (law51_phase2) anim_se(4902) = 1 !ENER_law51-phase2
1629 IF (law51_phase3) anim_se(4903) = 1 !ENER_law51-phase3
1630 IF (law51_phase4) anim_se(4904) = 1 !ENER_law51-phase4
1631 ELSEIF(key3(1:4)=='TEMP')THEN
1632 anim_se(4) = 1
1633 IF (law51_phase1) anim_se(4905) = 1 !TEMP_law51-phase1
1634 IF (law51_phase2) anim_se(4906) = 1 !TEMP_law51-phase2
1635 IF (law51_phase3) anim_se(4907) = 1 !TEMP_law51-phase3
1636 IF (law51_phase4) anim_se(4908) = 1 !TEMP_law51-phase4
1637 ELSEIF(key3(1:5)=='TEM51')THEN
1638 anim_se(4) = 1
1639 ELSEIF(key3(1:1)=='P'.AND.key3(1:5)/='PLA51')THEN
1640 anim_se(6) = 1
1641 IF (law51_phase1) anim_se(4909) = 1 !PRES_law51-phase1
1642 IF (law51_phase2) anim_se(4910) = 1 !PRES_law51-phase2
1643 IF (law51_phase3) anim_se(4911) = 1 !PRES_law51-phase3
1644 IF (law51_phase4) anim_se(4912) = 1 !PRES_law51-phase4
1645 ELSEIF(key3(1:4)=='VONM')THEN
1646 anim_se(7) = 1
1647 ELSEIF(key3(1:1)=='K')THEN
1648 anim_se(8) = 1
1649 ELSEIF(key3(1:4)=='TVIS')THEN
1650 anim_se(9) = 1
1651 ELSEIF(key3(1:5)=='VORTX')THEN
1652 anim_se(10) = 1
1653 ELSEIF(key3(1:4)=='DAM1')THEN
1654 anim_se(11) = 1
1655 ELSEIF(key3(1:4)=='DAM2')THEN
1656 anim_se(12) = 1
1657 ELSEIF(key3(1:4)=='DAM3')THEN
1658 anim_se(13) = 1
1659 ELSEIF(key3(1:5)=='SIGX ')THEN
1660 anim_se(14) = 1
1661 ELSEIF(key3(1:5)=='SIGY ')THEN
1662 anim_se(15) = 1
1663 ELSEIF(key3(1:5)=='SIGZ ')THEN
1664 anim_se(16) = 1
1665 ELSEIF(key3(1:5)=='SIGXY'.OR.key3(1:5)=='SIGYX')THEN
1666 anim_se(17) = 1
1667 ELSEIF(key3(1:5)=='SIGYZ'.OR.key3(1:5)=='SIGZY')THEN
1668 anim_se(18) = 1
1669 ELSEIF(key3(1:5)=='SIGZX'.OR.key3(1:5)=='SIGXZ')THEN
1670 anim_se(19) = 1
1671 ELSEIF(key3(1:5)=='USR1 '.OR.key3(1:5)=='USER1')THEN
1672 anim_se(20) = 1
1673 ELSEIF(key3(1:5)=='USR2 '.OR.key3(1:5)=='USER2')THEN
1674 anim_se(21) = 1
1675 ELSEIF(key3(1:5)=='usr3 .OR.'KEY3(1:5)=='user3')THEN
1676 ANIM_SE(22) = 1
1677 ELSEIF(KEY3(1:5)=='usr4 .OR.'KEY3(1:5)=='user4')THEN
1678 ANIM_SE(23) = 1
1679 ELSEIF(KEY3(1:5)=='usr5 .OR.'KEY3(1:5)=='user5')THEN
1680 ANIM_SE(24) = 1
1681 ELSEIF(KEY3(1:4)=='hour')THEN
1682 ANIM_SE(25) = 1
1683 ELSEIF(KEY3(1:4)=='epsd')THEN
1684 ANIM_SE(26) = 1
1685 ELSEIF(KEY3(1:5)=='usr6 .OR.'KEY3(1:5)=='user6')THEN
1686 ANIM_SE(27) = 1
1687 ELSEIF(KEY3(1:5)=='usr7 .OR.'KEY3(1:5)=='user7')THEN
1688 ANIM_SE(28) = 1
1689 ELSEIF(KEY3(1:5)=='usr8 .OR.'KEY3(1:5)=='user8')THEN
1690 ANIM_SE(29) = 1
1691 ELSEIF(KEY3(1:5)=='usr9 .OR.'KEY3(1:5)=='user9')THEN
1692 ANIM_SE(30) = 1
1693 ELSEIF(KEY3(1:5)=='usr10')THEN
1694 ANIM_SE(31) = 1
1695 ELSEIF(KEY3(1:5)=='usr11')THEN
1696 ANIM_SE(32) = 1
1697 ELSEIF(KEY3(1:5)=='usr12')THEN
1698 ANIM_SE(33) = 1
1699 ELSEIF(KEY3(1:5)=='usr13')THEN
1700 ANIM_SE(34) = 1
1701 ELSEIF(KEY3(1:5)=='usr14')THEN
1702 ANIM_SE(35) = 1
1703 ELSEIF(KEY3(1:5)=='usr15')THEN
1704 ANIM_SE(36) = 1
1705 ELSEIF(KEY3(1:5)=='usr16')THEN
1706 ANIM_SE(37) = 1
1707 ELSEIF(KEY3(1:5)=='usr17')THEN
1708 ANIM_SE(38) = 1
1709 ELSEIF(KEY3(1:5)=='usr18')THEN
1710 ANIM_SE(39) = 1
1711 ELSEIF(KEY3(1:5)=='usr19')THEN
1712 ANIM_SE(40) = 1
1713 ELSEIF(KEY3(1:5)=='usr20')THEN
1714 ANIM_SE(41) = 1
1715 ELSEIF(KEY3(1:5)=='usr21')THEN
1716 ANIM_SE(42) = 1
1717 ELSEIF(KEY3(1:5)=='usr22')THEN
1718 ANIM_SE(43) = 1
1719 ELSEIF(KEY3(1:5)=='usr23')THEN
1720 ANIM_SE(44) = 1
1721 ELSEIF(KEY3(1:5)=='usr24')THEN
1722 ANIM_SE(45) = 1
1723 ELSEIF(KEY3(1:5)=='usr25')THEN
1724 ANIM_SE(46) = 1
1725 ELSEIF(KEY3(1:5)=='usr26')THEN
1726 ANIM_SE(47) = 1
1727 ELSEIF(KEY3(1:5)=='usr27')THEN
1728 ANIM_SE(48) = 1
1729 ELSEIF(KEY3(1:5)=='usr28')THEN
1730 ANIM_SE(49) = 1
1731 ELSEIF(KEY3(1:5)=='usr29')THEN
1732 ANIM_SE(50) = 1
1733 ELSEIF(KEY3(1:5)=='usr30')THEN
1734 ANIM_SE(51) = 1
1735 ELSEIF(KEY3(1:5)=='usr31')THEN
1736 ANIM_SE(52) = 1
1737 ELSEIF(KEY3(1:5)=='usr32')THEN
1738 ANIM_SE(52) = 1
1739 ELSEIF(KEY3(1:5)=='usr33')THEN
1740 ANIM_SE(54) = 1
1741 ELSEIF(KEY3(1:5)=='usr34')THEN
1742 ANIM_SE(55) = 1
1743 ELSEIF(KEY3(1:5)=='usr35')THEN
1744 ANIM_SE(56) = 1
1745 ELSEIF(KEY3(1:5)=='usr36')THEN
1746 ANIM_SE(57) = 1
1747 ELSEIF(KEY3(1:5)=='usr37')THEN
1748 ANIM_SE(58) = 1
1749 ELSEIF(KEY3(1:5)=='usr38')THEN
1750 ANIM_SE(59) = 1
1751 ELSEIF(KEY3(1:5)=='usr39')THEN
1752 ANIM_SE(60) = 1
1753 ELSEIF(KEY3(1:5)=='usr40')THEN
1754 ANIM_SE(61) = 1
1755 ELSEIF(KEY3(1:5)=='usr41')THEN
1756 ANIM_SE(62) = 1
1757 ELSEIF(KEY3(1:5)=='usr42')THEN
1758 ANIM_SE(63) = 1
1759 ELSEIF(KEY3(1:5)=='usr43')THEN
1760 ANIM_SE(64) = 1
1761 ELSEIF(KEY3(1:5)=='usr44')THEN
1762 ANIM_SE(65) = 1
1763 ELSEIF(KEY3(1:5)=='usr45')THEN
1764 ANIM_SE(66) = 1
1765 ELSEIF(KEY3(1:5)=='usr46')THEN
1766 ANIM_SE(67) = 1
1767 ELSEIF(KEY3(1:5)=='usr47')THEN
1768 ANIM_SE(68) = 1
1769 ELSEIF(KEY3(1:5)=='usr48')THEN
1770 ANIM_SE(69) = 1
1771 ELSEIF(KEY3(1:5)=='usr49')THEN
1772 ANIM_SE(70) = 1
1773 ELSEIF(KEY3(1:5)=='usr50')THEN
1774 ANIM_SE(71) = 1
1775 ELSEIF(KEY3(1:5)=='usr51')THEN
1776 ANIM_SE(72) = 1
1777 ELSEIF(KEY3(1:5)=='usr52')THEN
1778 ANIM_SE(73) = 1
1779 ELSEIF(KEY3(1:5)=='usr53')THEN
1780 ANIM_SE(74) = 1
1781 ELSEIF(KEY3(1:5)=='usr54')THEN
1782 ANIM_SE(75) = 1
1783 ELSEIF(KEY3(1:5)=='usr55')THEN
1784 ANIM_SE(76) = 1
1785 ELSEIF(KEY3(1:5)=='usr56')THEN
1786 ANIM_SE(77) = 1
1787 ELSEIF(KEY3(1:5)=='usr57')THEN
1788 ANIM_SE(78) = 1
1789 ELSEIF(KEY3(1:5)=='usr58')THEN
1790 ANIM_SE(79) = 1
1791 ELSEIF(KEY3(1:5)=='usr59')THEN
1792 ANIM_SE(80) = 1
1793 ELSEIF(KEY3(1:5)=='usr60')THEN
1794 ANIM_SE(81) = 1
1795 ELSEIF(KEY3(1:4)=='wpla')THEN
1796 READ(KEY4(1:5),'(i5)')N1
1797 IF (N1 /= 0 ) THEN
1798 ANIM_SE(81+N1) = 1
1799 ELSE
1800 ANIM_SE(4964) = 1
1801 ENDIF
1802 ELSEIF(KEY3(1:4)=='flay')THEN
1803 ANIM_SE(282) = 1
1804 ELSEIF(KEY3(1:5)=='vfrac')THEN
1805 !BIMAT(law20)
1806 !BIPHAS(law37) : TRIMAT = -2
1807 !MULTIMAT(law51) : TRIMAT = 3|4
1808 ! law51 + law37 : TRIMAT = 3|4
1809 !consistency with material law checked in radioss2.F since TRIMAT common var not yet read
1810 ANIM_SE(283) = 1 !vfrac phase1 USR4 (law51 & 37 & 20)
1811 ANIM_SE(284) = 1 !vfrac phase2 USR4 (law51 & 37 & 20)
1812 ANIM_SE(285) = 1 !vfrac phase3 USR22 (law51)
1813 ANIM_SE(286) = 1 !vfrac phase4 USR48 (law51)
1814 ELSEIF(KEY3(1:5)=='bfrac.OR.' KEY3(1:4)=='bf51')THEN
1815 ANIM_SE(887) = 1
1816 ELSEIF(KEY3(1:5)=='orthd')THEN
1817 IF(KEY4=='all')THEN
1818 IORTHDALL = 1
1819 ELSE
1820 READ(KEY4(1:5),'(i5)')N1
1821 IF(N1/=0)THEN
1822 ANIM_SE(286+3*(N1-1)+1)=1
1823 ANIM_SE(286+3*(N1-1)+2)=1
1824 ANIM_SE(286+3*(N1-1)+3)=1
1825 ENDIF
1826 ENDIF
1827 ELSEIF(KEY3(1:5)=='vdam1')THEN
1828 READ(KEY4(1:5),'(i5)')N1 !/ijk
1829 IF(N1 > 0)THEN
1830 ANIM_SE(890+N1)=1
1831 ELSE
1832 ANIM_SE(888) = 1
1833 ENDIF
1834 ELSEIF(KEY3(1:5)=='vdam2')THEN
1835 READ(KEY4(1:5),'(i5)')N1 !/ijk
1836 IF(N1 > 0)THEN
1837 ANIM_SE(1890+N1)=1
1838 ELSE
1839 ANIM_SE(889) = 1
1840 ENDIF
1841 ELSEIF(KEY3(1:5)=='vdam3')THEN
1842 READ(KEY4(1:5),'(i5)')N1 !/ijk
1843 IF(N1 > 0)THEN
1844 ANIM_SE(2890+N1)=1
1845 ELSE
1846 ANIM_SE(890) = 1
1847 ENDIF
1848 ELSEIF(KEY3(1:4)=='dama')THEN
1849 READ(KEY4(1:5),'(i5)')N1 !/ijk
1850 IF(N1 > 0)THEN
1851 IF ( MOD(ABS(N1)/10,10) == 0) THEN
1852 READ(KEY5,'(i5)')N2
1853 ANIM_SE(3890+MOD(ABS(N1),100)+ N2*10+(ABS(N1)/100)*2010)=1
1854 ELSE
1855 ANIM_SE(3890+N1)=1
1856 ENDIF
1857 ELSE
1858 ANIM_SE(3890) = 1
1859 ENDIF
1860 ELSEIF(KEY3(1:4)=='tdel')THEN
1861 ANIM_SE(4890)=1
1862 ELSEIF(KEY3(1:3)=='ssp')THEN
1863 ANIM_SE(4891) = 1
1864 IF (LAW51_PHASE1) ANIM_SE(4917) = 1 !SSP_law51-phase1
1865 IF (LAW51_PHASE2) ANIM_SE(4918) = 1 !SSP_law51-phase2
1866 IF (LAW51_PHASE3) ANIM_SE(4919) = 1 !SPP_law51-phase3
1867 IF (LAW51_PHASE4) ANIM_SE(4920) = 1 !SSP_law51-phase4
1868 ELSEIF(KEY3(1:5)=='schli')THEN
1869 ANIM_SE(4892) = 1
1870 ELSEIF(KEY3(1:5)=='domai')THEN
1871 ANIM_SE(4893) = 1
1872 ELSEIF(KEY3(1:4)=='fill')THEN
1873 ANIM_SE(4894) = 1
1874 ELSEIF(KEY3(1:5)=='sigeq')THEN ! element equivalent stress criteria
1875 ANIM_SE(4895) = 1
1876 ELSEIF(KEY3(1:4)=='bulk.OR.' KEY3(1:4)=='qvis')THEN
1877 ANIM_SE(4896) = 1
1878 IF (LAW51_PHASE1) ANIM_SE(4931) = 1
1879 IF (LAW51_PHASE2) ANIM_SE(4932) = 1
1880 IF (LAW51_PHASE3) ANIM_SE(4933) = 1
1881 IF (LAW51_PHASE4) ANIM_SE(4934) = 1
1882
1883CTentation: numeros reserves anim_save for sub-materiaux law 51
1884C ANIM_SE(4897) ==> !DENS_law51-phase1
1885C ANIM_SE(4898) ==> ! phase2
1886C ANIM_SE(4899) ==> ! phase3
1887C ANIM_SE(4900) ==> ! phase4
1888C ANIM_SE(4901) ==> !ENER_law51-phase1
1889C ANIM_SE(4902) ==> ! phase2
1890C ANIM_SE(4903) ==> ! phase3
1891C ANIM_SE(4904) ==> ! phase4
1892C ANIM_SE(4905) ==> !TEMP_law51-phase1
1893C ANIM_SE(4906) ==> ! phase2
1894C ANIM_SE(4907) ==> ! phase3
1895C ANIM_SE(4908) ==> ! phase4
1896C ANIM_SE(4909) ==> !PRES_law51-phase1
1897C ANIM_SE(4910) ==> ! phase2
1898C ANIM_SE(4911) ==> ! phase3
1899C ANIM_SE(4912) ==> ! phase4
1900C ANIM_SE(4913) ==> !EPSP_law51-phase1
1901C ANIM_SE(4914) ==> ! phase2
1902C ANIM_SE(4915) ==> ! phase3
1903C ANIM_SE(4916) ==> ! phase4
1904C ANIM_SE(4917) ==> ! SSP_law51-phase1
1905C ANIM_SE(4918) ==> ! phase2
1906C ANIM_SE(4919) ==> ! phase3
1907C ANIM_SE(4920) ==> ! phase4
1908
1909 ELSEIF(KEY3(1:3)=='vol.OR.' KEY3(1:4)=='volu')THEN
1910 ANIM_SE(4921) = 1
1911 ANIM_SE(4922) = 1
1912 ANIM_SE(4923) = 1
1913 ANIM_SE(4924) = 1
1914 ANIM_SE(4925) = 1
1915
1916CTentation: numeros reserves anim_save for sub-materiaux law 51
1917C ANIM_SE(4926) ==> !MASS_law51-phase1
1918C ANIM_SE(4927) ==> ! phase2
1919C ANIM_SE(4928) ==> ! phase3
1920C ANIM_SE(4929) ==> ! phase4
1921
1922 ELSEIF(KEY3(1:4)=='tdet')THEN
1923 ANIM_SE(4930) = 1
1924
1925CTentation: numeros reserves anim_save for sub-materiaux law 51
1926C ANIM_SE(4931) ==> !QVIS_law51-phase1
1927C ANIM_SE(4932) ==> ! phase2
1928C ANIM_SE(4933) ==> ! phase3
1929C ANIM_SE(4934) ==> ! phase4
1930C Anim_se (4935) ==>!Law37 - Liquid
1931C ANIM_SE(4936) ==> ! law37 - gas
1932
1933 ELSEIF(KEY3(1:5)=='dt ')THEN
1934 ANIM_SE(4937) = 1
1935
1936 ELSEIF(KEY3(1:3) == 'mom.OR.' KEY3(2:4) == 'mom')THEN
1937 IF(KEY3(1:5) == 'mom ')THEN
1938 ANIM_SE(4938)=1
1939 ANIM_SE(4939)=1
1940 ANIM_SE(4940)=1
1941 ANIM_SE(4941)=1
1942 ANIM_SE(4942)=1
1943 ANIM_SE(4943)=1
1944 ANIM_SE(4944)=1
1945 ELSEIF(KEY3(1:5) == 'momx ')THEN
1946 ANIM_SE(4938)=1
1947 ELSEIF(KEY3(1:5) == 'momy ')THEN
1948 ANIM_SE(4939)=1
1949 ELSEIF(KEY3(1:5) == 'momz ')THEN
1950 ANIM_SE(4940)=1
1951 ELSEIF(KEY3(1:5) == 'momxy')THEN
1952 ANIM_SE(4941)=1
1953 ELSEIF(KEY3(1:5) == 'momyz')THEN
1954 ANIM_SE(4942)=1
1955 ELSEIF(KEY3(1:5) == 'momxz')THEN
1956 ANIM_SE(4943)=1
1957 ELSEIF(KEY3(1:5) == '|mom|')THEN
1958 ANIM_SE(4944)=1
1959 ENDIF
1960 ELSEIF(KEY3(1:3) == 'vel.OR.' KEY3(2:4) == 'vel' )THEN
1961 IF(KEY3(1:5) == 'vel ')THEN
1962 ANIM_SE(4945)=1
1963 ANIM_SE(4946)=1
1964 ANIM_SE(4947)=1
1965 ANIM_SE(4948)=1
1966 ANIM_SE(4949)=1
1967 ANIM_SE(4950)=1
1968 ANIM_SE(4951)=1
1969 ELSEIF(KEY3(1:5) == 'velx ')THEN
1970 ANIM_SE(4945)=1
1971 ELSEIF(KEY3(1:5) == 'vely ')THEN
1972 ANIM_SE(4946)=1
1973 ELSEIF(KEY3(1:5) == 'velz ')THEN
1974 ANIM_SE(4947)=1
1975 ELSEIF(KEY3(1:5) == 'velxy')THEN
1976 ANIM_SE(4948)=1
1977 ELSEIF(KEY3(1:5) == 'velyz')THEN
1978 ANIM_SE(4949)=1
1979 ELSEIF(KEY3(1:5) == 'velxz')THEN
1980 ANIM_SE(4950)=1
1981 ELSEIF(KEY3(1:5) == '|vel|')THEN
1982 ANIM_SE(4951)=1
1983 ENDIF
1984 ELSEIF(KEY3(1:3) == 'fin.OR.' KEY3(2:4) == 'fin')THEN
1985 IF(KEY3(1:5) == 'fint ')THEN
1986 ANIM_SE(4952)=1
1987 ANIM_SE(4953)=1
1988 ANIM_SE(4954)=1
1989 ANIM_SE(4955)=1
1990 ANIM_SE(4956)=1
1991 ANIM_SE(4957)=1
1992 ANIM_SE(4958)=1
1993 ELSEIF(KEY3(1:5) == 'finx ')THEN
1994 ANIM_SE(4952)=1
1995 ELSEIF(KEY3(1:5) == 'finy ')THEN
1996 ANIM_SE(4953)=1
1997 ELSEIF(KEY3(1:5) == 'finz ')THEN
1998 ANIM_SE(4954)=1
1999 ELSEIF(KEY3(1:5) == 'finxy')THEN
2000 ANIM_SE(4955)=1
2001 ELSEIF(KEY3(1:5) == 'finyz')THEN
2002 ANIM_SE(4956)=1
2003 ELSEIF(KEY3(1:5) == 'finxz')THEN
2004 ANIM_SE(4957)=1
2005 ELSEIF(KEY3(1:5) == '|fin|')THEN
2006 ANIM_SE(4958)=1
2007 ENDIF
2008 ELSEIF(KEY3(1:3)=='ams')THEN
2009 ANIM_SE(4959) = 1
2010 ELSEIF(KEY3(1:5)=='vorty')THEN
2011 ANIM_E(4960) = 1
2012 ANIM_SE(4960) = 1
2013 ELSEIF(KEY3(1:5)=='vortz')THEN
2014 ANIM_E(4961) = 1
2015 ANIM_SE(4961) = 1
2016 ELSEIF(KEY3(1:4)=='vort')THEN
2017 ANIM_E(4962) = 1
2018 ANIM_SE(4962) = 1
2019 ELSEIF(KEY3(1:4)=='eint')THEN
2020 ANIM_E(4963) = 1
2021 ANIM_SE(4963) = 1
2022 ELSEIF(KEY3(1:3) == 'off')THEN
2023 ANIM_SE(4965) = 1
2024 ELSEIF(KEY3(1:4)=='mach')THEN
2025 ANIM_SE(4966) = 1
2026 ELSEIF(KEY3(1:5)=='color')THEN
2027 ANIM_SE(4967) = 1
2028 ELSEIF(KEY3(1:4)=='damg')THEN
2029 ANIM_SE(4968) = 1
2030 ELSEIF(KEY3(1:7)=='nl_epsp')THEN
2031 ANIM_SE(4969) = 1
2032 ELSEIF(KEY3(1:7)=='nl_epsd')THEN
2033 ANIM_SE(4970) = 1
2034 ELSEIF(KEY3(1:6)=='tsaiwu')THEN
2035 READ(KEY4(1:5),'(i5)')N1
2036 IF (N1 /= 0 ) THEN
2037 ANIM_SE(4971+N1) = 1
2038 ELSE
2039 ANIM_SE(4971) = 1
2040 ENDIF
2041 ELSEIF(KEY3(1:9)=='tillotson')THEN
2042 ANIM_SE(5172) = 1
2043 ELSEIF(KEY3(1:7)=='vstrain')THEN
2044 read(KEY4, '(i2)',ERR=2206) ILAY
20452206 CONTINUE
2046.and. if(ILAY >= 1 ILAY <= 10)then
2047 ANIM_SE(5173+ILAY) = 1
2048 else
2049 ANIM_SE(5173) = 1
2050 endif
2051!D YOUR NEW ANIMATION CONTOUR HERE
2052 !ELSEIF(...)THEN
2053 ! ANIM_SE(5173+22+1) = 1
2054
2055 ENDIF
2056C
2057C--------------------------
2058C FIL
2059C--------------------------
2060 ELSEIF(KEY2(1:4)=='trus.OR.'KEY2(1:4)=='beam.OR.'
2061 . KEY2(1:5)=='sprin')THEN
2062C--------------------------
2063C FIL:TORSEUR
2064C--------------------------
2065 IF(KEY3(1:4)=='forc')THEN
2066 ANIM_FT(1) = 1
2067C--------------------------
2068C FILS:SCALAIRES
2069C--------------------------
2070 ELSEIF(KEY3(1:4)=='ener')THEN
2071 ANIM_FE(3) = 1
2072 ELSEIF(KEY3(1:4)=='vonm')THEN
2073 ANIM_FE(7) = 1
2074 ELSEIF(KEY3(1:5)=='dam1')THEN
2075 ANIM_FE(11) = 1
2076 ELSEIF(KEY3(1:5)=='dam2')THEN
2077 ANIM_FE(12) = 1
2078 ELSEIF(KEY3(1:5)=='dam3')THEN
2079 ANIM_FE(13) = 1
2080 ELSEIF(KEY3(1:5)=='sigx ')THEN
2081 ANIM_FE(14) = 1
2082 ELSEIF(KEY3(1:5)=='sigy ')THEN
2083 ANIM_FE(15) = 1
2084 ELSEIF(KEY3(1:5)=='sigz ')THEN
2085 ANIM_FE(16) = 1
2086 ELSEIF(KEY3(1:5)=='sigxy.OR.'KEY3(1:5)=='sigyx')THEN
2087 ANIM_FE(17) = 1
2088 ELSEIF(KEY3(1:5)=='sigyz.OR.'KEY3(1:5)=='sigzy')THEN
2089 ANIM_FE(18) = 1
2090 ELSEIF(KEY3(1:5)=='sigzx.OR.'KEY3(1:5)=='sigxz')THEN
2091 ANIM_FE(19) = 1
2092 ELSEIF(KEY3(1:5)=='dt ')THEN
2093 ANIM_FE(20) = 1
2094 ELSEIF(KEY3(1:3)=='ams')THEN
2095 ANIM_FE(21) = 1
2096 ELSEIF(KEY3(1:3)=='off')THEN
2097 ANIM_FE(22) = 1
2098 ELSEIF(KEY3(1:4)=='epsp')THEN
2099 ANIM_FE(1) = 1 ! global element plastic strain
2100 IF (KEY2(1:4)=='beam') THEN
2101 N1 = 0
2102 IF (KEY4(1:3) == 'all') THEN
2103 IBEAM_EPSPALL = 1
2104 ELSE
2105 READ(KEY4(1:5),'(i5)')N1 ! ! global element plastic strain
2106 IF (N1 > 0) ANIM_FE(22 + N1) = 1
2107 ENDIF
2108 ENDIF ! IF (KEY2(1:4)=='beam')
2109 ELSEIF (KEY3(1:5)=='epsx ') THEN
2110 IF (KEY2(1:5)=='truss') THEN
2111 ANIM_FE(123) = 1
2112 ENDIF
2113 ELSEIF(KEY3(1:4)=='epsd')THEN
2114 ANIM_FE(124) = 1
2115 ELSEIF(KEY3(1:4)=='dama')THEN
2116 IF (KEY2(1:4)=='beam') THEN
2117 ANIM_FE(125) = 1 ! damage factor
2118 ENDIF
2119 ELSE
2120 IXITKEY=IXITKEY+1
2121 ENDIF
2122C--------------------------------------------
2123 ELSEIF(KEY2(1:4)=='elem')THEN
2124 IF (KEY3(1:5) == 'law51') THEN
2125 IF (KEY4(1:3) == 'all') THEN
2126 LAW51_PHASE1 = .TRUE.
2127 LAW51_PHASE2 = .TRUE.
2128 LAW51_PHASE3 = .TRUE.
2129 LAW51_PHASE4 = .TRUE.
2130 ELSE IF (KEY4(1:1) == '1') THEN
2131 LAW51_PHASE1 = .TRUE.
2132 ELSE IF (KEY4(1:1) == '2') THEN
2133 LAW51_PHASE2 = .TRUE.
2134 ELSE IF (KEY4(1:1) == '3') THEN
2135 LAW51_PHASE3 = .TRUE.
2136 ELSE IF (KEY4(1:1) == '4') THEN
2137 LAW51_PHASE4 = .TRUE.
2138 ENDIF
2139C--------------------------
2140C GENERIQUE SCALAIRES
2141C--------------------------
2142 ELSE IF (KEY3(1:4)=='epsp')THEN
2143 ANIM_E(1) = 1
2144 ANIM_SE(1) = 1
2145 ANIM_CE(1) = 1
2146 ANIM_FE(1) = 1
2147 IF (LAW51_PHASE1) ANIM_SE(4913) = 1 !EPSP_law51-phase1
2148 IF (LAW51_PHASE2) ANIM_SE(4914) = 1 !EPSP_law51-phase2
2149 IF (LAW51_PHASE3) ANIM_SE(4915) = 1 !EPSP_law51-phase3
2150 IF (LAW51_PHASE4) ANIM_SE(4916) = 1 !EPSP_law51-phase4
2151 ELSEIF(KEY3(1:5)=='dens3')THEN
2152 !/ANIM/ELEM/DENS37
2153 ANIM_SE(4935) = 1 !law37 - rho_liquid
2154 ANIM_SE(4936) = 1 !law37 - rho_gas
2155 ELSEIF(KEY3(1:4)=='dens')THEN
2156 ANIM_E(2) = 1
2157 ANIM_SE(2) = 1
2158 ANIM_CE(2) = 1
2159 IF (LAW51_PHASE1) THEN
2160 ANIM_SE(4897) = 1 !DENS_law51-phase1
2161 ANIM_CE(11890) = 1
2162 ENDIF
2163 IF (LAW51_PHASE2) THEN
2164 ANIM_SE(4898) = 1 !DENS_law51-phase2
2165 ANIM_CE(11891) = 1
2166 ENDIF
2167 IF (LAW51_PHASE3) THEN
2168 ANIM_SE(4899) = 1 !DENS_law51-phase3
2169 ANIM_CE(11892) = 1
2170 ENDIF
2171 IF (LAW51_PHASE4) THEN
2172 ANIM_SE(4900) = 1 !DENS_law51-phase4
2173 ANIM_CE(11893) = 1
2174 ENDIF
2175 !Anim_se (4935) = 1! Law37 -rho_liquid
2176 !ANIM_SE(4936) = 1 !law37 - rho_gas
2177 ELSEIF(KEY3(1:4)=='ener')THEN
2178 ANIM_E(3) = 1
2179 ANIM_SE(3) = 1
2180 ANIM_CE(3) = 1
2181 ANIM_FE(3) = 1
2182 IF (LAW51_PHASE1) THEN
2183 ANIM_SE(4901) = 1 !ENER_law51-phase1
2184 ANIM_CE(11894) = 1
2185 ENDIF
2186 IF (LAW51_PHASE2) THEN
2187 ANIM_SE(4902) = 1 !ENER_law51-phase2
2188 ANIM_CE(11895) = 1
2189 ENDIF
2190 IF (LAW51_PHASE3) THEN
2191 ANIM_SE(4903) = 1 !ENER_law51-phase3
2192 ANIM_CE(11896) = 1
2193 ENDIF
2194 IF (LAW51_PHASE4) THEN
2195 ANIM_SE(4904) = 1 !ENER_law51-phase4
2196 ANIM_CE(11897) = 1
2197 ENDIF
2198 ELSEIF(KEY3(1:4)=='temp')THEN
2199 ANIM_E(4) = 1
2200 ANIM_SE(4) = 1
2201 ANIM_CE(4) = 1
2202 IF (LAW51_PHASE1) THEN
2203 ANIM_SE(4905) = 1 !TEMP_law51-phase1
2204 ANIM_CE(11898) = 1
2205 ENDIF
2206 IF (LAW51_PHASE2) THEN
2207 ANIM_SE(4906) = 1 !TEMP_law51-phase2
2208 ANIM_CE(11899) = 1
2209 ENDIF
2210 IF (LAW51_PHASE3) THEN
2211 ANIM_SE(4907) = 1 !TEMP_law51-phase3
2212 ANIM_CE(11900) = 1
2213 ENDIF
2214 IF (LAW51_PHASE4) THEN
2215 ANIM_SE(4908) = 1 !TEMP_law51-phase4
2216 ANIM_CE(11901) = 1
2217 ENDIF
2218 ELSEIF(KEY3(1:4)=='thic')THEN
2219 ANIM_E(5) = 1
2220 ANIM_CE(5) = 1
2221 ELSEIF(KEY3(1:1)=='p.AND.'KEY3(1:5)/='pla51')THEN
2222 ANIM_E(6) = 1
2223 ANIM_SE(6) = 1
2224 ANIM_CE(6) = 1
2225 IF (LAW51_PHASE1) THEN
2226 ANIM_SE(4909) = 1 !PRES_law51-phase1
2227 ANIM_CE(11902) = 1
2228 ENDIF
2229 IF (LAW51_PHASE2) THEN
2230 ANIM_SE(4910) = 1 !PRES_law51-phase2
2231 ANIM_CE(11903) = 1
2232 ENDIF
2233 IF (LAW51_PHASE3) THEN
2234 ANIM_SE(4911) = 1 !PRES_law51-phase3
2235 ANIM_CE(11904) = 1
2236 ENDIF
2237 IF (LAW51_PHASE4) THEN
2238 ANIM_SE(4912) = 1 !PRES_law51-phase4
2239 ANIM_CE(11905) = 1
2240 ENDIF
2241 ELSEIF(KEY3(1:4)=='vonm')THEN
2242 ANIM_E(7) = 1
2243 ANIM_SE(7) = 1
2244 ANIM_CE(7) = 1
2245 ANIM_FE(7) = 1
2246 ELSEIF(KEY3(1:1)=='k')THEN
2247 ANIM_E(8) = 1
2248 ANIM_SE(8) = 1
2249 ANIM_CE(8) = 1
2250 ELSEIF(KEY3(1:4)=='tvis')THEN
2251 ANIM_E(9) = 1
2252 ANIM_SE(9) = 1
2253 ANIM_CE(9) = 1
2254 ELSEIF(KEY3(1:5)=='vortx')THEN
2255 ANIM_E(10) = 1
2256 ANIM_SE(10) = 1
2257 ANIM_CE(10) = 1
2258 ELSEIF(KEY3(1:5)=='dam1')THEN
2259 ANIM_E(11) = 1
2260 ANIM_SE(11) = 1
2261 ANIM_CE(11) = 1
2262 ANIM_FE(11) = 1
2263 ELSEIF(KEY3(1:5)=='dam2')THEN
2264 ANIM_E(12) = 1
2265 ANIM_SE(12) = 1
2266 ANIM_CE(12) = 1
2267 ANIM_FE(12) = 1
2268 ELSEIF(KEY3(1:5)=='dam3')THEN
2269 ANIM_E(13) = 1
2270 ANIM_SE(13) = 1
2271 ANIM_CE(13) = 1
2272 ANIM_FE(13) = 1
2273 ELSEIF(KEY3(1:5)=='sigx ')THEN
2274 ANIM_E(14) = 1
2275 ANIM_SE(14) = 1
2276 ANIM_CE(14) = 1
2277 ANIM_FE(14) = 1
2278 ELSEIF(KEY3(1:5)=='sigy ')THEN
2279 ANIM_E(15) = 1
2280 ANIM_SE(15) = 1
2281 ANIM_CE(15) = 1
2282 ANIM_FE(15) = 1
2283 ELSEIF(KEY3(1:5)=='sigz ')THEN
2284 ANIM_E(16) = 1
2285 ANIM_SE(16) = 1
2286 ANIM_CE(16) = 1
2287 ANIM_FE(16) = 1
2288 ELSEIF(KEY3(1:5)=='sigxy.OR.'KEY3(1:5)=='sigyx')THEN
2289 ANIM_E(17) = 1
2290 ANIM_SE(17) = 1
2291 ANIM_CE(17) = 1
2292 ANIM_FE(17) = 1
2293 ELSEIF(KEY3(1:5)=='sigyz.OR.'KEY3(1:5)=='sigzy')THEN
2294 ANIM_E(18) = 1
2295 ANIM_SE(18) = 1
2296 ANIM_CE(18) = 1
2297 ANIM_FE(18) = 1
2298 ELSEIF(KEY3(1:5)=='sigzx.OR.'KEY3(1:5)=='sigxz')THEN
2299 ANIM_E(19) = 1
2300 ANIM_SE(19) = 1
2301 ANIM_CE(19) = 1
2302 ANIM_FE(19) = 1
2303 ELSEIF(KEY3(1:5)=='usr1 .OR.'KEY3(1:5)=='user1.OR.'
2304 . KEY3(1:5)=='pla51')THEN
2305 ANIM_E(20) = 1
2306 ANIM_SE(20) = 1
2307 ANIM_CE(20) = 1
2308 ELSEIF(KEY3(1:5)=='usr2 .OR.'KEY3(1:5)=='user2.OR.'
2309 . KEY3(1:5)=='tem51')THEN
2310 ANIM_E(21) = 1
2311 ANIM_SE(21) = 1
2312 ANIM_CE(21) = 1
2313 ELSEIF(KEY3(1:5)=='usr3 .OR.'KEY3(1:5)=='user3.OR.'
2314 . KEY3(1:5)=='bf51 ')THEN
2315 ANIM_E(22) = 1
2316 ANIM_SE(22) = 1
2317 ANIM_CE(22) = 1
2318 ELSEIF(KEY3(1:5)=='usr4 .OR.'KEY3(1:5)=='user4')THEN
2319 ANIM_E(23) = 1
2320 ANIM_SE(23) = 1
2321 ANIM_CE(23) = 1
2322 ELSEIF(KEY3(1:5)=='usr5 .OR.'KEY3(1:5)=='user5')THEN
2323 ANIM_E(24) = 1
2324 ANIM_SE(24) = 1
2325 ANIM_CE(24) = 1
2326 ELSEIF(KEY3(1:4)=='hour')THEN
2327 ANIM_E(25) = 1
2328 ANIM_SE(25) = 1
2329 ANIM_CE(25) = 1
2330 ELSEIF(KEY3(1:4)=='epsd')THEN
2331 ANIM_CE(26) = 1
2332 ANIM_SE(26) = 1
2333 ELSEIF(KEY3(1:5)=='usr6 .OR.'KEY3(1:5)=='user6')THEN
2334 ANIM_SE(27) = 1
2335 ANIM_CE(27) = 1
2336 ELSEIF(KEY3(1:5)=='usr7 .OR.'KEY3(1:5)=='user7')THEN
2337 ANIM_SE(28) = 1
2338 ANIM_CE(28) = 1
2339 ELSEIF(KEY3(1:5)=='usr8 .OR.'KEY3(1:5)=='user8')THEN
2340 ANIM_SE(29) = 1
2341 ANIM_CE(29) = 1
2342 ELSEIF(KEY3(1:5)=='usr9 .OR.'KEY3(1:5)=='user9')THEN
2343 ANIM_SE(30) = 1
2344 ANIM_CE(30) = 1
2345 ELSEIF(KEY3(1:5)=='usr10')THEN
2346 ANIM_SE(31) = 1
2347 ANIM_CE(31) = 1
2348 ELSEIF(KEY3(1:5)=='usr11')THEN
2349 ANIM_SE(32) = 1
2350 ANIM_CE(32) = 1
2351 ELSEIF(KEY3(1:5)=='usr12')THEN
2352 ANIM_SE(33) = 1
2353 ANIM_CE(33) = 1
2354 ELSEIF(KEY3(1:5)=='usr13')THEN
2355 ANIM_SE(34) = 1
2356 ANIM_CE(34) = 1
2357 ELSEIF(KEY3(1:5)=='usr14')THEN
2358 ANIM_SE(35) = 1
2359 ANIM_CE(35) = 1
2360 ELSEIF(KEY3(1:5)=='usr15')THEN
2361 ANIM_SE(36) = 1
2362 ANIM_CE(36) = 1
2363 ELSEIF(KEY3(1:5)=='usr16')THEN
2364 ANIM_SE(37) = 1
2365 ANIM_CE(37) = 1
2366 ELSEIF(KEY3(1:5)=='usr17')THEN
2367 ANIM_SE(38) = 1
2368 ANIM_CE(38) = 1
2369 ELSEIF(KEY3(1:5)=='usr18')THEN
2370 ANIM_SE(39) = 1
2371 ANIM_CE(39) = 1
2372 ELSEIF(KEY3(1:5)=='usr19')THEN
2373 ANIM_SE(40) = 1
2374 ANIM_CE(40) = 1
2375 ELSEIF(KEY3(1:5)=='usr20')THEN
2376 ANIM_SE(41) = 1
2377 ANIM_CE(41) = 1
2378 ELSEIF(KEY3(1:5)=='usr21')THEN
2379 ANIM_SE(42) = 1
2380 ANIM_CE(42) = 1
2381 ELSEIF(KEY3(1:5)=='usr22')THEN
2382 ANIM_SE(43) = 1
2383 ANIM_CE(43) = 1
2384 ELSEIF(KEY3(1:5)=='usr23')THEN
2385 ANIM_SE(44) = 1
2386 ANIM_CE(44) = 1
2387 ELSEIF(KEY3(1:5)=='usr24')THEN
2388 ANIM_SE(45) = 1
2389 ANIM_CE(45) = 1
2390 ELSEIF(KEY3(1:5)=='usr25')THEN
2391 ANIM_SE(46) = 1
2392 ANIM_CE(46) = 1
2393 ELSEIF(KEY3(1:5)=='usr26')THEN
2394 ANIM_SE(47) = 1
2395 ANIM_CE(47) = 1
2396 ELSEIF(KEY3(1:5)=='usr27')THEN
2397 ANIM_SE(48) = 1
2398 ANIM_CE(48) = 1
2399 ELSEIF(KEY3(1:5)=='usr28')THEN
2400 ANIM_SE(49) = 1
2401 ANIM_CE(49) = 1
2402 ELSEIF(KEY3(1:5)=='usr29')THEN
2403 ANIM_SE(50) = 1
2404 ANIM_CE(50) = 1
2405 ELSEIF(KEY3(1:5)=='usr30')THEN
2406 ANIM_SE(51) = 1
2407 ANIM_CE(51) = 1
2408 ELSEIF(KEY3(1:5)=='usr31')THEN
2409 ANIM_SE(52) = 1
2410 ANIM_CE(52) = 1
2411 ELSEIF(KEY3(1:5)=='usr32')THEN
2412 ANIM_SE(53) = 1
2413 ANIM_CE(53) = 1
2414 ELSEIF(KEY3(1:5)=='usr33')THEN
2415 ANIM_SE(54) = 1
2416 ANIM_CE(54) = 1
2417 ELSEIF(KEY3(1:5)=='usr34')THEN
2418 ANIM_SE(55) = 1
2419 ANIM_CE(55) = 1
2420 ELSEIF(KEY3(1:5)=='usr35')THEN
2421 ANIM_SE(56) = 1
2422 ANIM_CE(56) = 1
2423 ELSEIF(KEY3(1:5)=='usr36')THEN
2424 ANIM_SE(57) = 1
2425 ANIM_CE(57) = 1
2426 ELSEIF(KEY3(1:5)=='usr37')THEN
2427 ANIM_SE(58) = 1
2428 ANIM_CE(58) = 1
2429 ELSEIF(KEY3(1:5)=='usr38')THEN
2430 ANIM_SE(59) = 1
2431 ANIM_CE(59) = 1
2432 ELSEIF(KEY3(1:5)=='usr39')THEN
2433 ANIM_SE(60) = 1
2434 ANIM_CE(60) = 1
2435 ELSEIF(KEY3(1:5)=='usr40')THEN
2436 ANIM_SE(61) = 1
2437 ANIM_CE(61) = 1
2438 ELSEIF(KEY3(1:5)=='usr41')THEN
2439 ANIM_SE(62) = 1
2440 ANIM_CE(62) = 1
2441 ELSEIF(KEY3(1:5)=='usr42')THEN
2442 ANIM_SE(63) = 1
2443 ANIM_CE(63) = 1
2444 ELSEIF(KEY3(1:5)=='usr43')THEN
2445 ANIM_SE(64) = 1
2446 ANIM_CE(64) = 1
2447 ELSEIF(KEY3(1:5)=='usr44')THEN
2448 ANIM_SE(65) = 1
2449 ANIM_CE(65) = 1
2450 ELSEIF(KEY3(1:5)=='usr45')THEN
2451 ANIM_SE(66) = 1
2452 ANIM_CE(66) = 1
2453 ELSEIF(KEY3(1:5)=='usr46')THEN
2454 ANIM_SE(67) = 1
2455 ANIM_CE(67) = 1
2456 ELSEIF(KEY3(1:5)=='usr47')THEN
2457 ANIM_SE(68) = 1
2458 ANIM_CE(68) = 1
2459 ELSEIF(KEY3(1:5)=='usr48')THEN
2460 ANIM_SE(69) = 1
2461 ANIM_CE(69) = 1
2462 ELSEIF(KEY3(1:5)=='usr49')THEN
2463 ANIM_SE(70) = 1
2464 ANIM_CE(70) = 1
2465 ELSEIF(KEY3(1:5)=='usr50')THEN
2466 ANIM_SE(71) = 1
2467 ANIM_CE(71) = 1
2468 ELSEIF(KEY3(1:5)=='vfrac')THEN
2469 !BIMAT(law20)
2470 !BIPHAS(law37) : TRIMAT = -2
2471 !MULTIMAT(law51) : TRIMAT = 3|4
2472 ! law51 & law37 : TRIMAT = 3|4
2473 !consistency with material law checked in radioss2.F since TRIMAT common var not yet read
2474 ANIM_SE(283) = 1 !vfrac phase1 USR4 (law51 & 37)
2475 ANIM_SE(284) = 1 !vfrac phase2 USR4 (law51 & 37)
2476 ANIM_SE(285) = 1 !vfrac phase3 USR22 (law51)
2477 ANIM_SE(286) = 1 !vfrac phase4 USR48 (law51)
2478 ANIM_CE(10248) = 1 !VFRAC Phase1 (Law20/51)
2479 ANIM_CE(10249) = 1 !VFRAC Phase2 (Law20/51)
2480 ANIM_CE(10250) = 1 !VFRAC Phase3 (Law20/51)
2481 ANIM_CE(10251) = 1 !VFRAC Phase4 (Law20/51)
2482 ELSEIF(KEY3(1:5)=='bfrac')THEN
2483 ANIM_SE(887) = 1
2484 ANIM_CE(10252) = 1
2485 ELSEIF(KEY3(1:3)=='ssp')THEN
2486 ANIM_SE(4891) = 1
2487 ANIM_CE(10671) = 1
2488 IF (LAW51_PHASE1) THEN
2489 ANIM_SE(4917) = 1 ! SSP_law51-phase1
2490 ANIM_CE(11906) = 1
2491 ENDIF
2492 IF (LAW51_PHASE2) THEN
2493 ANIM_SE(4918) = 1 ! phase2
2494 ANIM_CE(11907) = 1
2495 ENDIF
2496 IF (LAW51_PHASE3) THEN
2497 ANIM_SE(4919) = 1 ! phase3
2498 ANIM_CE(11908) = 1
2499 ENDIF
2500 IF (LAW51_PHASE4) THEN
2501 ANIM_SE(4920) = 1 ! phase4
2502 ANIM_CE(11909) = 1
2503 ENDIF
2504 ELSEIF(KEY3(1:5)=='schli')THEN
2505 ANIM_SE(4892) = 1
2506 ANIM_CE(10672) = 1
2507 ELSEIF(KEY3(1:5)=='domai')THEN
2508 ANIM_SE(4893) = 1
2509 ANIM_CE(10676) = 1
2510 ELSEIF(KEY3(1:3)=='vol.OR.' KEY3(1:4)=='volu')THEN
2511 ANIM_SE(4921) = 1
2512 ANIM_CE(4*MX_PLY_ANIM + 14566) = 1
2513 IF (LAW51_PHASE1) THEN
2514 ANIM_SE(4922) = 1 !VOLU_law51-phase1
2515 ANIM_CE(11914) = 1
2516 ENDIF
2517 IF (LAW51_PHASE2) THEN
2518 ANIM_SE(4923) = 1 ! phase2
2519 ANIM_CE(11915) = 1
2520 ENDIF
2521 IF (LAW51_PHASE3) THEN
2522 ANIM_SE(4924) = 1 ! phase3
2523 ANIM_CE(11916) = 1
2524 ENDIF
2525 IF (LAW51_PHASE4) THEN
2526 ANIM_SE(4925) = 1 ! phase4
2527 ANIM_CE(11917) = 1
2528 ENDIF
2529 ELSEIF(KEY3(1:4)=='tdet')THEN
2530 ANIM_SE(4930) = 1
2531 ANIM_CE(11889) = 1
2532 ELSEIF(KEY3(1:5)=='law20')THEN
2533 ANIM_CE(11890) = 1 !dens
2534 ANIM_CE(11891) = 1
2535 ANIM_CE(11892) = 0
2536 ANIM_CE(11893) = 0
2537 ANIM_CE(11894) = 1 !ener
2538 ANIM_CE(11895) = 1
2539 ANIM_CE(11896) = 0
2540 ANIM_CE(11897) = 0
2541 ANIM_CE(11898) = 1 !temp
2542 ANIM_CE(11899) = 1
2543 ANIM_CE(11900) = 0
2544 ANIM_CE(11901) = 0
2545 ANIM_CE(11902) = 1 !pres
2546 ANIM_CE(11903) = 1
2547 ANIM_CE(11904) = 0
2548 ANIM_CE(11905) = 0
2549 ANIM_CE(11906) = 1 !plas
2550 ANIM_CE(11907) = 1
2551 ANIM_CE(11908) = 0
2552 ANIM_CE(11909) = 0
2553 ANIM_CE(11910) = 1 !ssp
2554 ANIM_CE(11911) = 1
2555 ANIM_CE(11912) = 0
2556 ANIM_CE(11913) = 0
2557 ANIM_CE(11914) = 1 !volume
2558 ANIM_CE(11915) = 1
2559 ANIM_CE(11916) = 0
2560 ANIM_CE(11917) = 0
2561 ANIM_CE(11918) = 1 !mass
2562 ANIM_CE(11919) = 1
2563 ANIM_CE(11920) = 0
2564 ANIM_CE(11921) = 0
2565 ANIM_CE(11922) = 1 !Qvis
2566 ANIM_CE(11923) = 1
2567 ANIM_CE(11924) = 0
2568 ANIM_CE(11925) = 0
2569 ELSEIF(KEY3(1:5)=='sigeq')THEN ! element equivalent stress criteria
2570 ANIM_SE(4895) = 1
2571 ANIM_CE(10677) = 1
2572 ELSEIF(KEY3(1:4)=='bulk.OR.' KEY3(1:4)=='qvis')THEN
2573 ANIM_SE(4896) = 1
2574 IF (LAW51_PHASE1) ANIM_SE(4931) = 1
2575 IF (LAW51_PHASE2) ANIM_SE(4932) = 1
2576 IF (LAW51_PHASE3) ANIM_SE(4933) = 1
2577 IF (LAW51_PHASE4) ANIM_SE(4934) = 1
2578 ANIM_CE(11888) = 1 !PID15 UNIM_CE (10678: 11887)
2579 ELSEIF(KEY3(1:5)=='dt ')THEN
2580 IDX = 13242 + 4*MX_PLY_ANIM
2581 ANIM_SE(4937) = 1
2582 ANIM_CE(IDX) = 1
2583 ANIM_FE(20) = 1
2584 ELSEIF(KEY3(1:3) == 'mom.OR.' KEY3(2:4) == 'mom' )THEN
2585 IF(KEY3(1:5) == 'mom ')THEN
2586 ANIM_SE(4938) =1
2587 ANIM_SE(4939) =1
2588 ANIM_SE(4940) =1
2589 ANIM_SE(4941) =1
2590 ANIM_SE(4942) =1
2591 ANIM_SE(4943) =1
2592 ANIM_SE(4944) =1
2593 ELSEIF(KEY3(1:5) == 'momx ')THEN
2594 ANIM_SE(4938) =1
2595 ELSEIF(KEY3(1:5) == 'momy ')THEN
2596 ANIM_SE(4939) =1
2597 ELSEIF(KEY3(1:5) == 'momz ')THEN
2598 ANIM_SE(4940) =1
2599 ELSEIF(KEY3(1:5) == 'momxy')THEN
2600 ANIM_SE(4941) =1
2601 ELSEIF(KEY3(1:5) == 'momyz')THEN
2602 ANIM_SE(4942) =1
2603 ELSEIF(KEY3(1:5) == 'momxz')THEN
2604 ANIM_SE(4943) =1
2605 ELSEIF(KEY3(1:5) == '|mom|')THEN
2606 ANIM_SE(4944) =1
2607 ENDIF
2608 ELSEIF(KEY3(1:3) == 'vel.OR.' KEY3(2:4) == 'vel' )THEN
2609 IF(KEY3(1:5) == 'vel ')THEN
2610 ! 3d
2611 ANIM_SE(4945) =1
2612 ANIM_SE(4946) =1
2613 ANIM_SE(4947) =1
2614 ANIM_SE(4948) =1
2615 ANIM_SE(4949) =1
2616 ANIM_SE(4950) =1
2617 ANIM_SE(4951) =1
2618 ! 2d
2619 IDX0 = 15921 + 4*MX_PLY_ANIM ! VEL-Y 0:10
2620 IDX1 = IDX0+10+1 ! VEL-Z 0:10
2621 IDX2 = IDX1+10+1 ! VEL 0:10
2622 ANIM_CE(IDX0) = 1
2623 ANIM_CE(IDX1) = 1
2624 ANIM_CE(IDX2) = 1
2625 IF(KEY4(1:2) == '10')THEN
2626 ANIM_CE(IDX2+10) = 1
2627 ANIM_CE(IDX2+10) = 1
2628 ANIM_CE(IDX2+10) = 1
2629 ELSEIF(KEY4(1:1) >= '1.AND.' KEY4(1:1) <= '9') THEN
2630 ILAY = ICHAR(KEY4(1:1)) - ICHAR('0')
2631 ANIM_CE(IDX2+ILAY) = 1
2632 ANIM_CE(IDX2+ILAY) = 1
2633 ANIM_CE(IDX2+ILAY) = 1
2634 ENDIF
2635 ELSEIF(KEY3(1:5) == 'velx ')THEN
2636 ANIM_SE(4945) =1
2637 ELSEIF(KEY3(1:5) == 'vely ')THEN
2638 ANIM_SE(4946) =1
2639 IDX0 = 15921 + 4*MX_PLY_ANIM ! VEL-Y 0:10
2640 IDX1 = IDX0+10+1 ! VEL-Z 0:10
2641 IDX2 = IDX1+10+1 ! VEL 0:10
2642 IF(KEY4(1:2) == '10')THEN
2643 ANIM_CE(IDX0+10) = 1
2644 ANIM_CE(IDX0+10) = 1
2645 ANIM_CE(IDX0+10) = 1
2646 ELSEIF(KEY4(1:1) >= '1.AND.' KEY4(1:1) <= '9') THEN
2647 ILAY = ICHAR(KEY4(1:1)) - ICHAR('0')
2648 ANIM_CE(IDX0+ILAY) = 1
2649 ANIM_CE(IDX0+ILAY) = 1
2650 ANIM_CE(IDX0+ILAY) = 1
2651 ENDIF
2652 ELSEIF(KEY3(1:5) == 'velz ')THEN
2653 ANIM_SE(4947) =1
2654 IDX0 = 15921 + 4*MX_PLY_ANIM ! VEL-Y 0:10
2655 IDX1 = IDX0+10+1 ! VEL-Z 0:10
2656 IDX2 = IDX1+10+1 ! VEL 0:10
2657 IF(KEY4(1:2) == '10')THEN
2658 ANIM_CE(IDX1+10) = 1
2659 ANIM_CE(IDX1+10) = 1
2660 ANIM_CE(IDX1+10) = 1
2661 ELSEIF(KEY4(1:1) >= '1.AND.' KEY4(1:1) <= '9') THEN
2662 ILAY = ICHAR(KEY4(1:1)) - ICHAR('0')
2663 ANIM_CE(IDX1+ILAY) = 1
2664 ANIM_CE(IDX1+ILAY) = 1
2665 ANIM_CE(IDX1+ILAY) = 1
2666 ENDIF
2667 ELSEIF(KEY3(1:5) == 'velxy')THEN
2668 ANIM_SE(4948) =1
2669 ELSEIF(KEY3(1:5) == 'velyz')THEN
2670 ANIM_SE(4949) =1
2671 ELSEIF(KEY3(1:5) == 'velxz')THEN
2672 ANIM_SE(4950) =1
2673 ELSEIF(KEY3(1:5) == '|vel|')THEN
2674 ANIM_SE(4951) =1
2675 ENDIF
2676 ELSEIF(KEY3(1:3) == 'fin.OR.' KEY3(2:4) == 'fin' )THEN
2677 IF(KEY3(1:5) == 'fint ')THEN
2678 ANIM_SE(4952)=1
2679 ANIM_SE(4953)=1
2680 ANIM_SE(4954)=1
2681 ANIM_SE(4955)=1
2682 ANIM_SE(4956)=1
2683 ANIM_SE(4957)=1
2684 ANIM_SE(4958)=1
2685 ELSEIF(KEY3(1:5) == 'finx ')THEN
2686 ANIM_SE(4952)=1
2687 ELSEIF(KEY3(1:5) == 'finy ')THEN
2688 ANIM_SE(4953)=1
2689 ELSEIF(KEY3(1:5) == 'finz ')THEN
2690 ANIM_SE(4954)=1
2691 ELSEIF(KEY3(1:5) == 'finxy')THEN
2692 ANIM_SE(4955)=1
2693 ELSEIF(KEY3(1:5) == 'finyz')THEN
2694 ANIM_SE(4956)=1
2695 ELSEIF(KEY3(1:5) == 'finxz')THEN
2696 ANIM_SE(4957)=1
2697 ELSEIF(KEY3(1:5) == '|fin|')THEN
2698 ANIM_SE(4958)=1
2699 ENDIF
2700 ELSEIF(KEY3(1:3)=='ams')THEN
2701 IDX = 13242 + 4*MX_PLY_ANIM + 1
2702 ANIM_SE(4959) = 1
2703 ANIM_CE(IDX) = 1
2704 ANIM_FE(21) = 1
2705 ELSEIF(KEY3(1:5)=='vorty')THEN
2706 ANIM_E(4960) = 1
2707 ANIM_SE(4960) = 1
2708 ELSEIF(KEY3(1:5)=='vortz')THEN
2709 ANIM_E(4961) = 1
2710 ANIM_SE(4961) = 1
2711 ELSEIF(KEY3(1:4)=='vort')THEN
2712 ANIM_E(4962) = 1
2713 ANIM_SE(4962) = 1
2714 ELSEIF(KEY3(1:4)=='eint')THEN
2715 ANIM_E(4963) = 1
2716 ANIM_SE(4963) = 1
2717 IDX = 13242 + 4*MX_PLY_ANIM + 2
2718 ANIM_CE(IDX) = 1
2719 ELSEIF(KEY3(1:4)=='wpla')THEN
2720 IDX = 13242 + 4*MX_PLY_ANIM + 3
2721 ANIM_CE(IDX) = 1
2722 ANIM_SE(4964) = 1
2723 ELSEIF(KEY3(1:3)=='off')THEN
2724 IDX = 13547 + 4*MX_PLY_ANIM + 1000 + 1
2725 ANIM_CE(IDX) = 1
2726 ANIM_SE(4965) = 1
2727 ANIM_FE(22) = 1
2728 ELSEIF(KEY3(1:4)=='mach')THEN
2729 IDX = 13547 + 4*MX_PLY_ANIM + 1000 + 2
2730 ANIM_CE(IDX) = 1
2731 ANIM_SE(4966) = 1
2732 ELSEIF(KEY3(1:5)=='color')THEN
2733 IDX = 13547 + 4*MX_PLY_ANIM + 1000 + 3
2734 ANIM_CE(IDX) = 1
2735 ANIM_SE(4967) = 1
2736 ELSEIF(KEY3(1:4) == 'damg') THEN
2737 IDX = 14551 + 4*MX_PLY_ANIM
2738 ANIM_CE(IDX) = 1
2739 ANIM_SE(4968) = 1
2740 ELSE IF(KEY3(1:7)=='nl_epsp')THEN
2741 IDX = 14567 + 4*MX_PLY_ANIM
2742 ANIM_CE(IDX) = 1
2743 ANIM_SE(4969) = 1
2744 ELSE IF(KEY3(1:7)=='nl_epsd')THEN
2745 IDX = 14581 + 4*MX_PLY_ANIM
2746 ANIM_CE(IDX) = 1
2747 ANIM_SE(4970) = 1
2748 ELSE IF(KEY3(1:6)=='tsaiwu')THEN
2749 IDX = 14595 + 4*MX_PLY_ANIM
2750 ANIM_CE(IDX) = 1
2751 ANIM_SE(4971) = 1
2752 ELSEIF(KEY3(1:9)=='tillotson')THEN
2753 IDX = 15898 + 4*MX_PLY_ANIM
2754 ANIM_CE(IDX) = 1
2755 ANIM_SE(5172) = 1
2756 ELSEIF(KEY3(1:7)=='vstrain')THEN
2757 IDX = 15899 + 4*MX_PLY_ANIM
2758 ILAY = 0
2759 read(KEY4, '(i2)', ERR=2205) ILAY
27602205 CONTINUE
2761.and. if(ILAY >= 1 ILAY <= 10)then
2762 ANIM_CE(IDX+ILAY) = 1
2763 ANIM_SE(5173+ILAY) = 1
2764 else
2765 ANIM_CE(IDX) = 1
2766 ANIM_SE(5173) = 1
2767 endif
2768 !ELSEIF()
2769 ! last used index in ANIM_CE : IDX = 4*MX_PLY_ANIM + 15943
2770 ! next availability in ANIM_CE : IDX = 4*MX_PLY_ANIM + 15944
2771 ! last introduction : VEL,VELy and VELZ for ANIM_CE (2d)
2772 ELSE
2773 IXITKEY=IXITKEY+1
2774 ENDIF
2775C--------------------------
2776C GPS : GRID POINT STRESS 1,2, use ANIM_N(101-301)
2777C--------------------------
2778 ELSEIF(KEY2(1:4)=='gps1')THEN
2779 IADG = IAD_GPS
2780 IF(KEY3(1:1)=='p')THEN
2781 ANIM_N(IADG+1) = 1
2782 ELSEIF(KEY3(1:4)=='vonm')THEN
2783 ANIM_N(IADG+2) = 1
2784 ELSEIF(KEY3(1:4)=='tens')THEN
2785 ANIM_N(IADG+3) = 1
2786 ANIM_N(IADG+4) = 1
2787 ANIM_N(IADG+5) = 1
2788 ANIM_N(IADG+6) = 1
2789 ANIM_N(IADG+7) = 1
2790 ANIM_N(IADG+8) = 1
2791 ELSEIF(KEY3(1:5)=='sigx ')THEN
2792 ANIM_N(IADG+3) = 1
2793 ELSEIF(KEY3(1:5)=='sigy ')THEN
2794 ANIM_N(IADG+4) = 1
2795 ELSEIF(KEY3(1:5)=='sigz ')THEN
2796 ANIM_N(IADG+5) = 1
2797 ELSEIF(KEY3(1:5)=='sigxy.OR.'KEY3(1:5)=='sigyx')THEN
2798 ANIM_N(IADG+6) = 1
2799 ELSEIF(KEY3(1:5)=='sigyz.OR.'KEY3(1:5)=='sigzy')THEN
2800 ANIM_N(IADG+7) = 1
2801 ELSEIF(KEY3(1:5)=='sigzx.OR.'KEY3(1:5)=='sigxz')THEN
2802 ANIM_N(IADG+8) = 1
2803 ELSEIF(KEY3(1:5)=='shell')THEN
2804 IF(KEY4(1:5)=='upper')THEN
2805 ANIM_N(IADG+9) = 1
2806 ANIM_N(IADG+10) = 1
2807 ANIM_N(IADG+11) = 1
2808 ANIM_N(IADG+12) = 1
2809 ANIM_N(IADG+13) = 1
2810 ANIM_N(IADG+14) = 1
2811 ELSEIF(KEY4(1:5)=='lower')THEN
2812 ANIM_N(IADG+15) = 1
2813 ANIM_N(IADG+16) = 1
2814 ANIM_N(IADG+17) = 1
2815 ANIM_N(IADG+18) = 1
2816 ANIM_N(IADG+19) = 1
2817 ANIM_N(IADG+20) = 1
2818 ELSE
2819 IXITKEY=IXITKEY+1
2820 ENDIF
2821 ELSE
2822 IXITKEY=IXITKEY+1
2823 ENDIF
2824 ELSEIF(KEY2(1:4)=='gps2')THEN
2825 IADG = IAD_GPS+100
2826 IF(KEY3(1:1)=='p')THEN
2827 ANIM_N(IADG+1) = 1
2828 ELSEIF(KEY3(1:4)=='vonm')THEN
2829 ANIM_N(IADG+2) = 1
2830 ELSEIF(KEY3(1:4)=='tens')THEN
2831 ANIM_N(IADG+3) = 1
2832 ANIM_N(IADG+4) = 1
2833 ANIM_N(IADG+5) = 1
2834 ANIM_N(IADG+6) = 1
2835 ANIM_N(IADG+7) = 1
2836 ANIM_N(IADG+8) = 1
2837 ELSEIF(KEY3(1:5)=='sigx ')THEN
2838 ANIM_N(IADG+3) = 1
2839 ELSEIF(KEY3(1:5)=='sigy ')THEN
2840 ANIM_N(IADG+4) = 1
2841 ELSEIF(KEY3(1:5)=='sigz ')THEN
2842 ANIM_N(IADG+5) = 1
2843 ELSEIF(KEY3(1:5)=='sigxy.OR.'KEY3(1:5)=='sigyx')THEN
2844 ANIM_N(IADG+6) = 1
2845 ELSEIF(KEY3(1:5)=='sigyz.OR.'KEY3(1:5)=='sigzy')THEN
2846 ANIM_N(IADG+7) = 1
2847 ELSEIF(KEY3(1:5)=='sigzx.OR.'KEY3(1:5)=='sigxz')THEN
2848 ANIM_N(IADG+8) = 1
2849 ELSEIF(KEY3(1:5)=='shell')THEN
2850 IF(KEY4(1:5)=='upper')THEN
2851 ANIM_N(IADG+9) = 1
2852 ANIM_N(IADG+10) = 1
2853 ANIM_N(IADG+11) = 1
2854 ANIM_N(IADG+12) = 1
2855 ANIM_N(IADG+13) = 1
2856 ANIM_N(IADG+14) = 1
2857 ELSEIF(KEY4(1:5)=='lower')THEN
2858 ANIM_N(IADG+15) = 1
2859 ANIM_N(IADG+16) = 1
2860 ANIM_N(IADG+17) = 1
2861 ANIM_N(IADG+18) = 1
2862 ANIM_N(IADG+19) = 1
2863 ANIM_N(IADG+20) = 1
2864 ELSE
2865 IXITKEY=IXITKEY+1
2866 ENDIF
2867 ELSE
2868 IXITKEY=IXITKEY+1
2869 ENDIF
2870 ELSEIF(KEY2(1:3)=='gps')THEN
2871 IADG = IAD_GPS+200
2872 IF(KEY3(1:4)=='tens')THEN
2873 ANIM_N(IADG+1) = 1
2874 ANIM_N(IADG+2) = 1
2875 ANIM_N(IADG+3) = 1
2876 ANIM_N(IADG+4) = 1
2877 ANIM_N(IADG+5) = 1
2878 ANIM_N(IADG+6) = 1
2879 ELSEIF(KEY3(1:5)=='sigx ')THEN
2880 ANIM_N(IADG+1) = 1
2881 ELSEIF(KEY3(1:5)=='sigy ')THEN
2882 ANIM_N(IADG+2) = 1
2883 ELSEIF(KEY3(1:5)=='sigz ')THEN
2884 ANIM_N(IADG+3) = 1
2885 ELSEIF(KEY3(1:5)=='sigxy.OR.'KEY3(1:5)=='sigyx')THEN
2886 ANIM_N(IADG+4) = 1
2887 ELSEIF(KEY3(1:5)=='sigyz.OR.'KEY3(1:5)=='sigzy')THEN
2888 ANIM_N(IADG+5) = 1
2889 ELSEIF(KEY3(1:5)=='sigzx.OR.'KEY3(1:5)=='sigxz')THEN
2890 ANIM_N(IADG+6) = 1
2891c
2892 ELSEIF(KEY3(1:5)=='stres')THEN
2893 IADG = IAD_GPS+200
2894 IF(KEY4(1:4)=='tens')THEN
2895 ANIM_N(IADG+1) = 1
2896 ANIM_N(IADG+2) = 1
2897 ANIM_N(IADG+3) = 1
2898 ANIM_N(IADG+4) = 1
2899 ANIM_N(IADG+5) = 1
2900 ANIM_N(IADG+6) = 1
2901 ELSEIF(KEY4(1:5)=='sigx ')THEN
2902 ANIM_N(IADG+1) = 1
2903 ELSEIF(KEY4(1:5)=='sigy ')THEN
2904 ANIM_N(IADG+2) = 1
2905 ELSEIF(KEY4(1:5)=='sigz ')THEN
2906 ANIM_N(IADG+3) = 1
2907 ELSEIF(KEY4(1:5)=='sigxy.OR.'KEY3(1:5)=='sigyx')THEN
2908 ANIM_N(IADG+4) = 1
2909 ELSEIF(KEY4(1:5)=='sigyz.OR.'KEY3(1:5)=='sigzy')THEN
2910 ANIM_N(IADG+5) = 1
2911 ELSEIF(KEY4(1:5)=='sigzx'.OR.key3(1:5)=='SIGXZ')THEN
2912 anim_n(iadg+6) = 1
2913 ELSE
2914 ixitkey=ixitkey+1
2915 ENDIF
2916c
2917 ELSEIF(key3(1:5)=='STRAI')THEN
2918 iadg = iad_gps+400
2919 IF(key4(1:4)=='TENS')THEN
2920 anim_n(iadg+1) = 1
2921 anim_n(iadg+2) = 1
2922 anim_n(iadg+3) = 1
2923 anim_n(iadg+4) = 1
2924 anim_n(iadg+5) = 1
2925 anim_n(iadg+6) = 1
2926 ELSEIF(key4(1:5)=='SIGX ')THEN
2927 anim_n(iadg+1) = 1
2928 ELSEIF(key4(1:5)=='SIGY ')THEN
2929 anim_n(iadg+2) = 1
2930 ELSEIF(key4(1:5)=='SIGZ ')THEN
2931 anim_n(iadg+3) = 1
2932 ELSEIF(key4(1:5)=='SIGXY'.OR.key3(1:5)=='SIGYX')THEN
2933 anim_n(iadg+4) = 1
2934 ELSEIF(key4(1:5)=='SIGYZ'.OR.key3(1:5)=='SIGZY')THEN
2935 anim_n(iadg+5) = 1
2936 ELSEIF(key4(1:5)=='SIGZX'.OR.key3(1:5)=='SIGXZ')THEN
2937 anim_n(iadg+6) = 1
2938 ELSE
2939 ixitkey=ixitkey+1
2940 ENDIF
2941 ELSE
2942 ixitkey=ixitkey+1
2943 ENDIF
2944 ELSEIF(key2(1:5)=='ISOST')THEN
2945 iadg = iad_gps+300
2946 IF(key3(1:4)=='TENS')THEN
2947 anim_n(iadg+1) = 1
2948 anim_n(iadg+2) = 1
2949 anim_n(iadg+3) = 1
2950 anim_n(iadg+4) = 1
2951 anim_n(iadg+5) = 1
2952 anim_n(iadg+6) = 1
2953 ELSEIF(key3(1:5)=='SIGX ')THEN
2954 anim_n(iadg+1) = 1
2955 ELSEIF(key3(1:5)=='SIGY ')THEN
2956 anim_n(iadg+2) = 1
2957 ELSEIF(key3(1:5)=='SIGZ ')THEN
2958 anim_n(iadg+3) = 1
2959 ELSEIF(key3(1:5)=='SIGXY'.OR.key3(1:5)=='SIGYX')THEN
2960 anim_n(iadg+4) = 1
2961 ELSEIF(key3(1:5)=='SIGYZ'.OR.key3(1:5)=='SIGZY')THEN
2962 anim_n(iadg+5) = 1
2963 ELSEIF(key3(1:5)=='SIGZX'.OR.key3(1:5)=='SIGXZ')THEN
2964 anim_n(iadg+6) = 1
2965 ELSE
2966 ixitkey=ixitkey+1
2967 ENDIF
2968C-----------------------------
2969C NODAL SCALAIRE
2970C-----------------------------
2971 ELSEIF(key2(1:4)=='NODA')THEN
2972 IF (key3(1:2)=='DT')THEN
2973 anim_n(1) = 1
2974 ELSEIF(key3(1:4)=='DMAS')THEN
2975 anim_n(2) = 1
2976 ELSEIF(key3(1:2)=='P ')THEN
2977 anim_n(3) = 1
2978 ELSEIF(key3(1:4)=='DENS')THEN
2979 anim_n(4) = 1
2980 ELSEIF(key3(1:4)=='ENER')THEN
2981 anim_n(5) = 1
2982 ELSEIF(key3(1:4)=='TEMP')THEN
2983 anim_n(6) = 1
2984 ELSEIF(key3(1:5)=='NVAR1') THEN
2985 anim_n(7) = 1
2986 ELSEIF(key3(1:5)=='NVAR2') THEN
2987 anim_n(8) = 1
2988 ELSEIF(key3(1:5)=='NVAR3') THEN
2989 anim_n(9) = 1
2990 ELSEIF(key3(1:5)=='NVAR4') THEN
2991 anim_n(10) = 1
2992 ELSEIF(key3(1:5)=='NVAR5') THEN
2993 anim_n(11) = 1
2994 ELSEIF(key3(1:5)=='DINER')THEN
2995 anim_n(12) = 1
2996 ELSEIF(key3(1:4)=='FPOT')THEN !potential
2997 anim_n(13) = 1
2998 ELSEIF(key3(1:5)=='NDMAS')THEN
2999 anim_n(14) = 1
3000 ELSEIF(key3(1:5)=='DAMA2')THEN
3001 anim_n(15) = 1
3002 anim_n(16) = 1
3003 ELSEIF(key3(1:5)=='SCHLI')THEN
3004 anim_n(17) = 0 !currently unplugged
3005 ELSEIF(key3(1:5)=='STIFR')THEN
3006 anim_n(18) = 1
3007 ELSEIF(key3(1:4)=='STIF')THEN
3008 anim_n(19) = 1
3009 ELSEIF(key3(1:5)=='VFRAC')THEN
3010 anim_n(20) = 1
3011 anim_n(21) = 1
3012 anim_n(22) = 1
3013 anim_n(23) = 1
3014 ELSEIF(key3(1:5)=='ZVFRA')THEN !inter22
3015 anim_n(24) = 1
3016 anim_n(25) = 1
3017 anim_n(26) = 1
3018 anim_n(27) = 1
3019 ELSEIF(key3(1:4)=='ZVOL')THEN !inter22
3020 anim_n(28) = 1
3021 anim_n(29) = 1
3022 ELSEIF(key3(1:3)=='SSP')THEN
3023 anim_n(30) = 1
3024 ELSEIF(key3(1:4)=='PEXT ')THEN
3025 anim_n(31) = 1
3026 output%DATA%ANIM_HAS_NODA_PEXT = 1
3027 ELSE
3028 ixitkey=ixitkey+1
3029 ENDIF
3030 ELSEIF(key2(1:4)=='MASS')THEN
3031 anim_m = 1
3032 IF (law51_phase1) THEN
3033 anim_se(4926) = 1 !MASS_law51-phase1
3034 anim_ce(11918) = 1
3035 ENDIF
3036 IF (law51_phase2) THEN
3037 anim_se(4927) = 1 !MASS_law51-phase2
3038 anim_ce(11919) = 1
3039 ENDIF
3040 IF (law51_phase3) THEN
3041 anim_se(4928) = 1 !MASS_law51-phase3
3042 anim_ce(11920) = 1
3043 ENDIF
3044 IF (law51_phase4) THEN
3045 anim_se(4929) = 1 !MASS_law51-phase4
3046 anim_ce(11921) = 1
3047 ENDIF
3048 ELSEIF(key2(1:4)=='NUMB')THEN
3049 anim_u = 1
3050 ELSEIF(key2(1:4)=='KEEP')THEN
3051 anim_k = 1
3052 ELSEIF(key2(1:5)=='ALEBA')THEN
3053 ifvani=1
3054 IF(key3(1:4)=='MASS')THEN
3055 anim_ce(2145) = 1
3056 ELSEIF(key3(1:4)=='PRES')THEN
3057 anim_ce(2146) = 1
3058 ELSEIF(key3(1:4)=='MOMX')THEN
3059 anim_ce(2147) = 1
3060 ELSEIF(key3(1:4)=='MOMY')THEN
3061 anim_ce(2148) = 1
3062 ELSEIF(key3(1:4)=='MOMZ')THEN
3063 anim_ce(2149) = 1
3064 ELSEIF(key3(1:4)=='DENS')THEN
3065 anim_ce(2150) = 1
3066 ELSEIF(key3(1:4)=='ENER')THEN
3067 anim_ce(2151) = 1
3068 ELSEIF(key3(1:4)=='SSP')THEN
3069 anim_ce(2152) = 1
3070 ELSEIF(key3(1:4)=='GAMA')THEN
3071 anim_ce(2153) = 1
3072 ELSEIF(key3(1:4)=='VISU')THEN
3073 anim_ce(2154) = 1
3074 ENDIF
3075 ELSE
3076 ixitkey=ixitkey+1
3077 ENDIF
3078C
3079 RETURN
3080 END
anim_dcod_key_0
subroutine anim_dcod_key_0(output, key2, key3, key4, key5, key6, jerr, ixitkey)
Definition anim_dcod_key_0.F:36
my_real
#define my_real
Definition cppsort.cpp:32
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::ncharkey
integer, parameter ncharkey
Definition names_and_titles_mod.F:1005
stack_mod
Definition stack_mod.F:75
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895