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