34
35
36
40 USE multi_fvm_mod
41
42
43
44#include "implicit_f.inc"
45
46
47
48 TYPE (H3D_KEYWORD) H3D_KEYWORD_QUAD_SCALAR(*)
50 INTEGER IGEO(NPROPGI,NUMGEO)
51 TYPE(MULTI_FVM_STRUCT), INTENT(IN) :: MULTI_FVM
52 INTEGER, INTENT(INOUT) :: NKEY
53
54
55
56#include "com04_c.inc"
57#include "param_c.inc"
58
59
60
61 INTEGER I,II
62
63
64
65 i = 0
66
67 i = i + 1
68 h3d_keyword_quad_scalar(i)%KEY3 = 'MASS'
69 h3d_keyword_quad_scalar(i)%TEXT1 = 'Mass'
70
71 i = i + 1
72 h3d_keyword_quad_scalar(i)%KEY3 = 'DENS'
73 h3d_keyword_quad_scalar(i)%TEXT1 = 'Density'
74
75 i = i + 1
76 h3d_keyword_quad_scalar(i)%KEY3 = 'ENER'
77 h3d_keyword_quad_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
78
79 i = i + 1
80 h3d_keyword_quad_scalar(i)%KEY3 = 'TEMP'
81 h3d_keyword_quad_scalar(i)%TEXT1 = 'Temperature'
82
83 i = i + 1
84 h3d_keyword_quad_scalar(i)%KEY3 = 'THICK'
85 h3d_keyword_quad_scalar(i)%TEXT1 = 'Thickness'
86
87 i = i + 1
88 h3d_keyword_quad_scalar(i)%KEY3 = 'P'
89 h3d_keyword_quad_scalar(i)%TEXT1 = 'Pressure'
90
91 i = i + 1
92 h3d_keyword_quad_scalar(i)%KEY3 = 'VONM'
93 h3d_keyword_quad_scalar(i)%TEXT1 = 'Von Mises'
94
95 i = i + 1
96 h3d_keyword_quad_scalar(i)%KEY3 = 'K'
97 h3d_keyword_quad_scalar(i)%TEXT1 = 'Turbulent Energy'
98
99 i = i + 1
100 h3d_keyword_quad_scalar(i)%KEY3 = 'TVIS'
101 h3d_keyword_quad_scalar(i)%TEXT1 = 'Turbulent Viscosity'
102
103 i = i + 1
104 h3d_keyword_quad_scalar(i)%KEY3 = 'VORTX'
105 h3d_keyword_quad_scalar(i)%TEXT1 = 'Vorticity'
106
107 i = i + 1
108 h3d_keyword_quad_scalar(i)%KEY3 = 'DAM1'
109 h3d_keyword_quad_scalar(i)%TEXT1 = 'Damage 1'
110
111 i = i + 1
112 h3d_keyword_quad_scalar(i)%KEY3 = 'DAM2'
113 h3d_keyword_quad_scalar(i)%TEXT1 = 'Damage 2'
114
115 i = i + 1
116 h3d_keyword_quad_scalar(i)%KEY3 = 'DAM3'
117 h3d_keyword_quad_scalar(i)%TEXT1 = 'Damage 3'
118
119 i = i + 1
120 h3d_keyword_quad_scalar(i)%KEY3 = 'SIGX'
121 h3d_keyword_quad_scalar(i)%TEXT1 = 'Stress X'
122
123 i = i + 1
124 h3d_keyword_quad_scalar(i)%KEY3 = 'SIGY'
125 h3d_keyword_quad_scalar(i)%TEXT1 = 'Stress Y'
126
127 i = i + 1
128 h3d_keyword_quad_scalar(i)%KEY3 = 'SIGZ'
129 h3d_keyword_quad_scalar(i)%TEXT1 = 'Stress Z'
130
131 i = i + 1
132 h3d_keyword_quad_scalar(i)%KEY3 = 'SIGXY'
133 h3d_keyword_quad_scalar(i)%TEXT1 = 'Stress XY'
134
135 i = i + 1
136 h3d_keyword_quad_scalar(i)%KEY3 = 'SIGYZ'
137 h3d_keyword_quad_scalar(i)%TEXT1 = 'Stress YZ'
138
139 i = i + 1
140 h3d_keyword_quad_scalar(i)%KEY3 = 'SIGZX'
141 h3d_keyword_quad_scalar(i)%TEXT1 = 'Stress ZX'
142
143 i = i + 1
144 h3d_keyword_quad_scalar(i)%KEY3 = 'USER'
145 h3d_keyword_quad_scalar(i)%IS_UVAR = 1
146 h3d_keyword_quad_scalar(i)%TEXT1 = 'User Var'
147
148 i = i + 1
149 h3d_keyword_quad_scalar(i)%KEY3 = 'HOURGLASS'
150 h3d_keyword_quad_scalar(i)%TEXT1 = 'Hourglass Energy per unit mass'
151
152 IF (multi_fvm%IS_USED) THEN
153
154
155
156
157 DO ii = 1, multi_fvm%NBMAT
158
159 i = i + 1
160 WRITE(h3d_keyword_quad_scalar(i)%KEY3, '(A9,I0)') 'M151VFRAC', ii
161 WRITE(h3d_keyword_quad_scalar(i)%TEXT1, '(A17,I0)') 'Volume fraction '
162 h3d_keyword_quad_scalar(i)%COMMENT = 'Volume fractions (for ALE multi-material Law151)'
163
164 i = i + 1
165 WRITE(h3d_keyword_quad_scalar(i)%KEY3, '(A8,I0)') 'M151DENS', ii
166 WRITE(h3d_keyword_quad_scalar(i)%TEXT1, '(A8,I0)') 'Density ', ii
167 h3d_keyword_quad_scalar(i)%COMMENT = 'Density (for ALE multi-material Law151)'
168
169 i = i + 1
170 WRITE(h3d_keyword_quad_scalar(i)%KEY3, '(A8,I0)') 'M151ENER', ii
171 WRITE(h3d_keyword_quad_scalar(i)%TEXT1, '(A7,I0)') 'Energy ', ii
172 h3d_keyword_quad_scalar(i)%COMMENT = 'Energy (for ALE multi-material Law151)'
173
174 i = i + 1
175 WRITE(h3d_keyword_quad_scalar(i)%KEY3, '(A8,I0)') 'M151PRES', ii
176 WRITE(h3d_keyword_quad_scalar(i)%TEXT1, '(A10,I0)') 'Pressure ', ii
177 h3d_keyword_quad_scalar(i)%COMMENT = 'Pressure (for ALE multi-material Law151)'
178
179 ENDDO
180 ELSE
181
182 i = i + 1
183 h3d_keyword_quad_scalar(i)%KEY3 = 'VFRAC1'
184 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Fraction 1'
185 h3d_keyword_quad_scalar'Volumetric fractions (for ALE multi-material laws: LAW20, LAW37 and LAW51)'
186
187 i = i + 1
188 h3d_keyword_quad_scalar(i)%KEY3 = 'VFRAC2'
189 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Fraction 2'
190 h3d_keyword_quad_scalar(i)%COMMENT = 'Volumetric fractions (for ALE multi-material laws: LAW20, LAW37 and LAW51)'
191
192 i = i + 1
193 h3d_keyword_quad_scalar(i)%KEY3 = 'VFRAC3'
194 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Fraction 3'
195 h3d_keyword_quad_scalar(i)%COMMENT = 'Volumetric fractions (for ALE multi-material laws: LAW20, LAW37 and LAW51)'
196
197 i = i + 1
198 h3d_keyword_quad_scalar(i)%KEY3 = 'VFRAC4'
199 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Fraction 4'
200 h3d_keyword_quad_scalar(i)%COMMENT = 'Volumetric fractions (for ALE multi-material laws: LAW20, LAW37 and LAW51)'
201
202 ENDIF
203
204 i = i + 1
205 h3d_keyword_quad_scalar(i)%KEY3 = 'BFRAC'
206 h3d_keyword_quad_scalar(i)%TEXT1 = 'Burn Fraction'
207
208 i = i + 1
209 h3d_keyword_quad_scalar(i)%KEY3 = 'SSP'
210 h3d_keyword_quad_scalar(i)%TEXT1 = 'Sound Speed'
211
212 i = i + 1
213 h3d_keyword_quad_scalar(i)%KEY3 = 'EPSD'
214 h3d_keyword_quad_scalar(i)%TEXT1 = 'Strain Rate'
215
216 i = i + 1
217 h3d_keyword_quad_scalar(i)%KEY3 = 'EPSP'
218 h3d_keyword_quad_scalar(i)%IS_LAYER = 1
219 h3d_keyword_quad_scalar(i)%IS_LAYER_ALL = 1
220 h3d_keyword_quad_scalar(i)%IS_IPT = 1
221 h3d_keyword_quad_scalar(i)%IS_IPT_ALL = 1
222 h3d_keyword_quad_scalar(i)%IS_PLY = 1
223 h3d_keyword_quad_scalar(i)%IS_PLY_ALL = 1
224 h3d_keyword_quad_scalar(i)%TEXT1 = 'Plastic Strain '
225
226 i = i + 1
227 h3d_keyword_quad_scalar(i)%KEY3 = 'TSAIWU'
228 h3d_keyword_quad_scalar(i)%TEXT1 = 'Tsai-Wu Criterion'
229
230 i = i + 1
231 h3d_keyword_quad_scalar(i)%KEY3 = 'FAIL'
232 h3d_keyword_quad_scalar(i)%IS_PLY = 1
233 h3d_keyword_quad_scalar(i)%IS_IPT = 1
234 h3d_keyword_quad_scalar(i)%TEXT1 = 'Nb of Failed layers'
235
236 i = i + 1
237 h3d_keyword_quad_scalar(i)%KEY3 = 'CRMAS'
238 h3d_keyword_quad_scalar(i)%TEXT1 = 'Airbag crossing mass'
239
240 i = i + 1
241 h3d_keyword_quad_scalar(i)%KEY3 = 'CRVEL'
242 h3d_keyword_quad_scalar(i)%TEXT1 = 'Airbag crossing velocity'
243
244 i = i + 1
245 h3d_keyword_quad_scalar(i)%KEY3 = 'THIN'
246 h3d_keyword_quad_scalar(i)%TEXT1 = '%Thinning'
247
248 i = i + 1
249 h3d_keyword_quad_scalar(i)%KEY3 = 'ERROR/THICK'
250 h3d_keyword_quad_scalar(i)%TEXT1 = 'Estimated Error on Thickness'
251
252 i = i + 1
253 h3d_keyword_quad_scalar(i)%KEY3 = 'PHI'
254 h3d_keyword_quad_scalar(i)%IS_PLY = 1
255 h3d_keyword_quad_scalar(i)%IS_PLY_ALL = 1
256 h3d_keyword_quad_scalar(i)%IS_LAYER = 1
257 h3d_keyword_quad_scalar(i)%IS_LAYER_ALL = 1
258 h3d_keyword_quad_scalar(i)%IS_IPT = 1
259 h3d_keyword_quad_scalar(i)%IS_IPT_ALL = 1
260 h3d_keyword_quad_scalar(i)%TEXT1 = 'PHI'
261
262 i = i + 1
263 h3d_keyword_quad_scalar(i)%KEY3 = 'DAMA'
264 h3d_keyword_quad_scalar(i)%IS_PLY = 1
265 h3d_keyword_quad_scalar(i)%IS_PLY_ALL = 1
266 h3d_keyword_quad_scalar(i)%IS_LAYER = 1
267 h3d_keyword_quad_scalar(i)%IS_LAYER_ALL = 1
268 h3d_keyword_quad_scalar(i)%IS_IPT = 1
269 h3d_keyword_quad_scalar(i)%IS_IPT_ALL = 1
270 h3d_keyword_quad_scalar(i)%TEXT1 = 'DAMAGE'
271
272 i = i + 1
273 h3d_keyword_quad_scalar(i)%KEY3 = 'FAILURE'
274 h3d_keyword_quad_scalar(i)%IS_PLY = 1
275 h3d_keyword_quad_scalar(i)%IS_PLY_ALL = 1
276 h3d_keyword_quad_scalar(i)%IS_LAYER = 1
277 h3d_keyword_quad_scalar(i)%IS_LAYER_ALL = 1
278 h3d_keyword_quad_scalar(i)%IS_IPT = 1
279 h3d_keyword_quad_scalar(i)%IS_IPT_ALL = 1
280 h3d_keyword_quad_scalar(i)%IS_ID = 1
281 h3d_keyword_quad_scalar(i)%IS_MODE = 1
282 h3d_keyword_quad_scalar(i)%TEXT1 = 'Failure'
283 h3d_keyword_quad_scalar(i)%COMMENT = 'Damage of current failure criterion'
284
285 i = i + 1
286 h3d_keyword_quad_scalar(i)%KEY3 = 'DAMA/MEMB'
287 h3d_keyword_quad_scalar(i)%TEXT1 = 'MAX DAMAGE MEMBRANE'
288
289 i = i + 1
290 h3d_keyword_quad_scalar(i)%KEY3 = 'NXTF'
291 h3d_keyword_quad_scalar(i)%IS_LAYER = 1
292 h3d_keyword_quad_scalar(i)%IS_LAYER_ALL = 1
293 h3d_keyword_quad_scalar(i)%IS_IPT = 1
294 h3d_keyword_quad_scalar(i)%IS_IPT_ALL = 1
295 h3d_keyword_quad_scalar(i)%TEXT1 = 'NXT FAILURE FACTOR'
296
297 i = i + 1
298 h3d_keyword_quad_scalar(i)%KEY3 = 'NXTF/MEMB'
299 h3d_keyword_quad_scalar(i)%TEXT1 = 'NXT FAILURE FACTOR MEMBRANE'
300
301 i = i + 1
302 h3d_keyword_quad_scalar(i)%KEY3 = 'TDEL'
303 h3d_keyword_quad_scalar(i)%TEXT1 = 'TIME DELETION ELEMENT'
304
305 i = i + 1
306 h3d_keyword_quad_scalar(i)%KEY3 = 'SCHLIEREN'
307 h3d_keyword_quad_scalar(i)%TEXT1 = 'Schlieren'
308 h3d_keyword_quad_scalar(i)%COMMENT =
309 .'Schlieren image (optical method widely used in CFD field). ALE material laws'
310
311 i = i + 1
312 h3d_keyword_quad_scalar(i)%KEY3 = 'PHI/MEMB'
313 h3d_keyword_quad_scalar(i)%TEXT1 = 'phi MEMBRANE'
314
315 i = i + 1
316 h3d_keyword_quad_scalar(i)%KEY3 = 'DOMAIN'
317 h3d_keyword_quad_scalar(i)%TEXT1 = 'Domain'
318
319 i = i + 1
320 h3d_keyword_quad_scalar(i)%KEY3 = 'SIGEQ'
321 h3d_keyword_quad_scalar(i)%TEXT1 = 'Equiv stress'
322
323 i = i + 1
324 h3d_keyword_quad_scalar(i)%KEY3 = 'BULK'
325 h3d_keyword_quad_scalar(i)%TEXT1 = 'Artificial Viscosity'
326
327 i = i + 1
328 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/DENS1'
329 h3d_keyword_quad_scalar(i)%TEXT1 = 'Density-1'
330
331 i = i + 1
332 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/DENS2'
333 h3d_keyword_quad_scalar(i)%TEXT1 = 'Density-2'
334
335 i = i + 1
336 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/ENER1'
337 h3d_keyword_quad_scalar(i)%TEXT1 = 'Specific Energy-1'
338
339 i = i + 1
340 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/ENER2'
341 h3d_keyword_quad_scalar(i)%TEXT1 = 'Specific Energy-2'
342
343 i = i + 1
344 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/TEMP1'
345 h3d_keyword_quad_scalar(i)%TEXT1 = 'Temperature-1'
346
347 i = i + 1
348 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/TEMP2'
349 h3d_keyword_quad_scalar(i)%TEXT1 = 'Temperature-2'
350
351 i = i + 1
352 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/P1'
353 h3d_keyword_quad_scalar(i)%TEXT1 = 'Pressure-1'
354
355 i = i + 1
356 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/P2'
357 h3d_keyword_quad_scalar(i)%TEXT1 = 'Pressure-2'
358
359 i = i + 1
360 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/EPSP1'
361 h3d_keyword_quad_scalar(i)%TEXT1 = 'Plastic Strain-1'
362
363 i = i + 1
364 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/EPSP1'
365 h3d_keyword_quad_scalar(i)%TEXT1 = 'Plastic Strain-2'
366
367 i = i + 1
368 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/SSP1'
369 h3d_keyword_quad_scalar(i)%TEXT1 = 'Sound Speed-1'
370
371 i = i + 1
372 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/SSP2'
373 h3d_keyword_quad_scalar(i)%TEXT1 = 'Sound Speed-2'
374
375 i = i + 1
376 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/VOLUM1'
377 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volume-1'
378
379 i = i + 1
380 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/VOLUM2'
381 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volume-2'
382
383 i = i + 1
384 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/QVIS1'
385 h3d_keyword_quad_scalar(i)%TEXT1 = 'Artificial Viscosity-1'
386
387 i = i + 1
388 h3d_keyword_quad_scalar(i)%KEY3 = 'LAW20/QVIS2'
389 h3d_keyword_quad_scalar(i)%TEXT1 = 'Artificial Viscosity-2'
390
391 i = i + 1
392 h3d_keyword_quad_scalar(i)%KEY3 = 'TDET'
393 h3d_keyword_quad_scalar(i)%TEXT1 = 'Detonation Time'
394
395 i = i + 1
396 h3d_keyword_quad_scalar(i)%KEY3 = 'DT'
397 h3d_keyword_quad_scalar(i)%TEXT1 = 'Time Step'
398
399 i = i + 1
400 h3d_keyword_quad_scalar(i)%KEY3 = 'AMS'
401 h3d_keyword_quad_scalar(i)%TEXT1 = 'AMS selection'
402
403 i = i + 1
404 h3d_keyword_quad_scalar(i)%KEY3 = 'ENTH'
405 h3d_keyword_quad_scalar(i)%TEXT1 = 'Enthalpy'
406
407 i = i + 1
408 h3d_keyword_quad_scalar(i)%KEY3 = 'ENTHM'
409 h3d_keyword_quad_scalar(i)%TEXT1 = 'Enthalpy per unit mass'
410
411 i = i + 1
412 h3d_keyword_quad_scalar(i)%KEY3 = 'ENTHV'
413 h3d_keyword_quad_scalar(i)%TEXT1 = 'Enthalpy per unit volume'
414
415 i = i + 1
416 h3d_keyword_quad_scalar(i)%KEY3 = 'EINTM'
417 h3d_keyword_quad_scalar(i)%TEXT1 = 'Internal Energy per unit mass'
418
419 i = i + 1
420 h3d_keyword_quad_scalar(i)%KEY3 = 'EINTV'
421 h3d_keyword_quad_scalar(i)%TEXT1 = 'Internal Energy per unit volume'
422
423 i = i + 1
424 h3d_keyword_quad_scalar(i)%KEY3 = 'EINT'
425 h3d_keyword_quad_scalar(i)%TEXT1 = 'Internal Energy'
426
427 i = i + 1
428 h3d_keyword_quad_scalar(i)%KEY3 = 'OFF'
429 h3d_keyword_quad_scalar(i)%TEXT1 = 'Element status'
430 h3d_keyword_quad_scalar(i)%COMMENT = '0.0:deleted >0.&<1.:under failure process 1.0:activated <0:Standby'
431
432 i = i + 1
433 h3d_keyword_quad_scalar(i)%KEY3 = 'MACH'
434 h3d_keyword_quad_scalar(i)%TEXT1 = 'Mach Number'
435 h3d_keyword_quad_scalar(i)%COMMENT = 'Mach Number. Only available with ALE material law 151'
436
437 i = i + 1
438 h3d_keyword_quad_scalar(i)%KEY3 = 'COLOR'
439 h3d_keyword_quad_scalar(i)%TEXT1 = 'Color Function'
440 h3d_keyword_quad_scalar(i)%COMMENT = 'Color function for multimaterial and multifluid ALE laws'
441
442 i = i + 1
443 h3d_keyword_quad_scalar(i)%KEY3 = 'VOLU'
444 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volume'
445 h3d_keyword_quad_scalar(i)%COMMENT = 'Element Volume'
446
447 i = i + 1
448 h3d_keyword_quad_scalar(i)%KEY3 = 'GROUP'
449 h3d_keyword_quad_scalar(i)%TEXT1 = 'Group identifier'
450
451 i = i + 1
452 h3d_keyword_quad_scalar(i)%KEY3 = 'INTERNAL.ID'
453 h3d_keyword_quad_scalar(i)%TEXT1 = 'Internal identifier'
454
455 i = i + 1
456 h3d_keyword_quad_scalar(i)%KEY3 = 'LOCAL.ID'
457 h3d_keyword_quad_scalar(i)%TEXT1 = 'Local identifier within the group'
458
459 i
460 h3d_keyword_quad_scalar(i)%KEY3 = 'TILLOTSON'
461 h3d_keyword_quad_scalar(i)%TEXT1 = 'Region identifier in p,v diagram'
462 h3d_keyword_quad_scalar(i)%COMMENT = '1:compressed, 2:expanded cold, 3:transition, 4:hot expansion'
463
464 i = i + 1
465 h3d_keyword_quad_scalar(i)%KEY3 = 'DIV(U)'
466 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric dilatation rate'
467 h3d_keyword_quad_scalar(i)%COMMENT = 'div(u)'
468
469 i = i + 1
470 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN'
471 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain'
472 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
473
474
475 i = i + 1
476 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/1'
477 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 1'
478 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
479
480 i = i + 1
481 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/2'
482 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 2'
483 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
484
485 i = i + 1
486 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/3'
487 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 3'
488 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
489
490 i = i + 1
491 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/4'
492 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 4'
493 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
494
495 i = i + 1
496 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/5'
497 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 5'
498 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
499
500 i = i + 1
501 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/6'
502 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 6'
503 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
504
505 i = i + 1
506 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/7'
507 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 7'
508 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
509
510 i = i + 1
511 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/8'
512 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 8'
513 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
514
515 i = i + 1
516 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/9'
517 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 9'
518 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
519
520 i = i + 1
521 h3d_keyword_quad_scalar(i)%KEY3 = 'VSTRAIN/10'
522 h3d_keyword_quad_scalar(i)%TEXT1 = 'Volumetric Strain 10'
523 h3d_keyword_quad_scalar(i)%COMMENT = 'mu=rho/rho0-1'
524
525
526 nkey=i
527