37!|| bcs_wall_trigger ../engine/source/boundary_conditions/bcs_wall_trigger.f90
40!|| brokmann_elem_renum ../
starter/source/materials/fail/windshield_alter/brokmann_elem_spmd_renum.f90
110!|| get_neighbour_surface_from_remote_proc ../engine/source/interfaces/interf/get_neighbour_surface_from_remote_proc.f90
172!|| law133_upd ../
starter/source/materials/mat/mat133/law133_upd.f90
173!|| law163_upd ../
starter/source/materials/mat/mat163/law163_upd.f90
185!|| mat25_tsaiwu_c ../engine/source/materials/mat/mat025/mat25_tsaiwu_c.f90
296#include "my_real.inc"
297 my_real,
parameter :: infinity = 1e20
298 my_real,
parameter :: zero = 0.
301 my_real,
parameter :: three = 3.
302 my_real,
parameter :: four = 4.
305 my_real,
parameter :: seven = 7.
306 my_real,
parameter :: eight = 8.
307 my_real,
parameter :: nine = 9.
308 my_real,
parameter :: ten = 10.
309 my_real,
parameter :: eleven = 11.
310 my_real,
parameter :: twelve = 12.
311 my_real,
parameter :: fourteen = 14.
312 my_real,
parameter :: fifteen = 15.
313 my_real,
parameter :: sixteen = 16.
314 my_real,
parameter :: eighteen = 18.
315 my_real,
parameter :: twenty = 20.
316 my_real,
parameter :: twenty4 = 24.
317 my_real,
parameter :: twenty5 = 25.
318 my_real,
parameter :: twenty7 = 27.
319 my_real,
parameter :: thirty = 30.
320 my_real,
parameter :: thirty2 = 32.
321 my_real,
parameter :: thirty6 = 36.
322 my_real,
parameter :: fourty = 40.
323 my_real,
parameter :: fourty5 = 45.
324 my_real,
parameter :: fourty8 = 48.
325 my_real,
parameter :: fifty = 50.
326 my_real,
parameter :: fifty4 = 54.
327 my_real,
parameter :: sixty4 = 64.
328 my_real,
parameter :: eighty = 80.
329 my_real,
parameter :: ninty = 90.
330 my_real,
parameter :: eighty16 = 96.
331 my_real,
parameter :: eighty19 = 99.
332 my_real,
parameter :: hundred = 100.
333 my_real,
parameter :: hundred28 = 128.
334 my_real,
parameter :: hundred50 = 150.
335 my_real,
parameter :: hundred80 = 180.
336 my_real,
parameter :: three100 = 300.
337 my_real,
parameter :: twohundred25 = 225.
338 my_real,
parameter :: twohundred95 = 295.
339 my_real,
parameter :: twohundred98 = 298.
340 my_real,
parameter :: twohundred93 = 293.
341 my_real,
parameter :: threehundred24 = 324.
342 my_real,
parameter :: fivehundred12 = 512.
343 my_real,
parameter :: two1000 = 2000.
344 my_real,
parameter :: three1000 = 3000.
345 my_real,
parameter :: ep02 = hundred
346 my_real,
parameter :: ep03 = ep02 * ten
347 my_real,
parameter :: ep04 = ep03 * ten
348 my_real,
parameter :: ep05 = ep04 * ten
349 my_real,
parameter :: ep06 = ep05 * ten
350 my_real,
parameter :: ep07 = ep06 * ten
351 my_real,
parameter :: ep08 = ep07 * ten
352 my_real,
parameter :: ep09 = ep08 * ten
353 my_real,
parameter :: ep10 = ep09 * ten
354 my_real,
parameter :: ep11 = ep10 * ten
355 my_real,
parameter :: ep12 = ep11 * ten
356 my_real,
parameter :: ep13 = ep12 * ten
357 my_real,
parameter :: ep14 = ep13 * ten
358 my_real,
parameter :: ep15 = ep14 * ten
359 my_real,
parameter :: ep16 = ep15 * ten
360 my_real,
parameter :: ep17 = ep16 * ten
361 my_real,
parameter :: ep18 = ep17 * ten
362 my_real,
parameter :: ep19 = ep18 * ten
363 my_real,
parameter :: ep20 = ep19 * ten
364 my_real,
parameter :: ep21 = ep20 * ten
365 my_real,
parameter :: ep29 = ep20 * ep09
366 my_real,
parameter :: ep30 = ep20 * ep10
367 my_real,
parameter :: ep31 = ep30 * ten
368 my_real,
parameter :: ep35 = ep30 * ep05
369 my_real,
parameter :: ep36 = ep30 * ep06
370 my_real,
parameter :: ep37 = ep30 * ep07
371 my_real,
parameter :: em01 = one / ten
372 my_real,
parameter :: em02 = one / ep02
373 my_real,
parameter :: em03 = one / ep03
374 my_real,
parameter :: em04 = one / ep04
375 my_real,
parameter :: em05 = one / ep05
376 my_real,
parameter :: em06 = one / ep06
377 my_real,
parameter :: em07 = one / ep07
378 my_real,
parameter :: em08 = one / ep08
379 my_real,
parameter :: em09 = one / ep09
380 my_real,
parameter :: em10 = one / ep10
381 my_real,
parameter :: em11 = one / ep11
382 my_real,
parameter :: em12 = one / ep12
383 my_real,
parameter :: em13 = one / ep13
384 my_real,
parameter :: em14 = one / ep14
385 my_real,
parameter :: em15 = one / ep15
386 my_real,
parameter :: em16 = one / ep16
387 my_real,
parameter :: em17 = one / ep17
388 my_real,
parameter :: em18 = one / ep18
389 my_real,
parameter :: em19 = one / ep19
390 my_real,
parameter :: em20 = one / ep20
391 my_real,
parameter :: em30 = one / ep30
392 my_real,
parameter :: em35 = one / ep35
393 my_real,
parameter :: zep1 = em01
394 my_real,
parameter :: zep2 = two / ten
395 my_real,
parameter :: zep3 = three / ten
396 my_real,
parameter :: zep4 = four / ten
397 my_real,
parameter :: zep5 = five / ten
398 my_real,
parameter :: zep6 = one - zep4
399 my_real,
parameter :: zep7 = one - zep3
400 my_real,
parameter :: zep8 = one - zep2
401 my_real,
parameter :: zep9 = one - zep1
402 my_real,
parameter :: zep01 = em02
403 my_real,
parameter :: zep02 = two / ep02
404 my_real,
parameter :: zep03 = three / ep02
405 my_real,
parameter :: zep04 = four / ep02
406 my_real,
parameter :: zep05 = five / ep02
407 my_real,
parameter :: zep06 = zep1 - zep04
408 my_real,
parameter :: zep07 = zep1 - zep03
409 my_real,
parameter :: zep08 = zep1 - zep02
410 my_real,
parameter :: zep09 = zep1 - zep01
411 my_real,
parameter :: twoem3 = two / ep03
412 my_real,
parameter :: threeem3 = three / ep03
413 my_real,
parameter :: fourem3 = four / ep03
414 my_real,
parameter :: fiveem3 = five / ep03
415 my_real,
parameter :: sixem3 = zep01 - fourem3
416 my_real,
parameter :: sevenem3 = zep01 - threeem3
417 my_real,
parameter :: eightem3 = zep01 - twoem3
418 my_real,
parameter :: nineem3 = zep01 - em03
419 my_real,
parameter :: twoem4 = two / ep04
420 my_real,
parameter :: threeem4 = three / ep04
421 my_real,
parameter :: fourem4 = four / ep04
422 my_real,
parameter :: fiveem4 = five / ep04
423 my_real,
parameter :: sixem4 = em03 - fourem4
424 my_real,
parameter :: sevenem4 = em03 - threeem4
425 my_real,
parameter :: eightem4 = em03 - twoem4
426 my_real,
parameter :: nineem4 = em03 - em04
427 my_real,
parameter :: twoem5 = two / ep05
428 my_real,
parameter :: threeem5 = three / ep05
429 my_real,
parameter :: fourem5 = four / ep05
430 my_real,
parameter :: fiveem5 = five / ep05
431 my_real,
parameter :: sixem5 = em04 - fourem5
432 my_real,
parameter :: sevenem5 = em04 - threeem5
433 my_real,
parameter :: eightem5 = em04 - twoem5
434 my_real,
parameter :: nineem5 = em04 - em05
435 my_real,
parameter :: twoem6 = two / ep06
436 my_real,
parameter :: threeem6 = three / ep06
437 my_real,
parameter :: fourem6 = four / ep06
438 my_real,
parameter :: fiveem6 = five / ep06
439 my_real,
parameter :: sixem6 = em05 - fourem6
440 my_real,
parameter :: sevenem6 = em05 - threeem6
441 my_real,
parameter :: eightem6 = em05 - twoem6
442 my_real,
parameter :: nineem6 = em05 - em06
443 my_real,
parameter :: twoep30 = two *ep30
444 my_real,
parameter :: threeep30 = three *ep30
445 my_real,
parameter :: fourep30 = four*ep30
446 my_real,
parameter :: zep9ep30 = nine *ep29
447 my_real,
parameter :: onep1ep30 = eleven * ep29
448 my_real,
parameter :: half = one/two
449 my_real,
parameter :: third = one/three
450 my_real,
parameter :: fourth = one/four
451 my_real,
parameter :: one_fifth = one/five
452 my_real,
parameter :: one_over_6 = one/six
453 my_real,
parameter :: one_over_8 = one/eight
454 my_real,
parameter :: one_over_9 = one/nine
455 my_real,
parameter :: one_over_12 = one/twelve
456 my_real,
parameter :: one_over_16 = one/sixteen
457 my_real,
parameter :: one_over_64 = one/sixty4
458 my_real,
parameter :: one_over_512 = one/fivehundred12
459 my_real,
parameter :: two_third = two/three
460 my_real,
parameter :: three_half = three/two
461 my_real,
parameter :: three_over_4 = three/four
462 my_real,
parameter :: trhee_over_14 = three/fourteen
463 my_real,
parameter :: tseize = three/sixteen
464 my_real,
parameter :: four_over_3 = four/three
465 my_real,
parameter :: four_over_5 = four/five
466 my_real,
parameter :: qneuf = four/nine
467 my_real,
parameter :: five_over_3 = five/three
468 my_real,
parameter :: five_over_4 = five/four
469 my_real,
parameter :: five_over_6 = five/six
470 my_real,
parameter :: six_over_5 = six/five
471 my_real,
parameter :: seven_half = seven/two
472 my_real,
parameter :: seven_over_5 = seven/five
473 my_real,
parameter :: seven_over_9 = seven/nine
474 my_real,
parameter :: nine_over_2 = nine/two
475 my_real,
parameter :: nine_over_10 = nine/ten
476 my_real,
parameter :: nseize = nine/sixteen
477 my_real,
parameter :: nine_over_20 = nine/twenty
478 my_real,
parameter :: zep11 = zep1 + zep01
479 my_real,
parameter :: zep12 = zep1 + zep02
480 my_real,
parameter :: zep13 = zep1 + zep03
481 my_real,
parameter :: zep14 = zep1 + zep04
482 my_real,
parameter :: zep15 = zep1 + zep05
483 my_real,
parameter :: zep16 = zep1 + zep06
484 my_real,
parameter :: zep17 = zep1 + zep07
485 my_real,
parameter :: zep18 = zep1 + zep08
486 my_real,
parameter :: zep19 = zep1 + zep09
487 my_real,
parameter :: zep21 = zep2 + zep01
488 my_real,
parameter :: zep23 = zep2 + zep03
489 my_real,
parameter :: zep25 = zep2 + zep05
490 my_real,
parameter :: zep26 = zep2 + zep06
491 my_real,
parameter :: zep213 = zep21 + threeem3
492 my_real,
parameter :: zep29 = zep2 + zep09
493 my_real,
parameter :: zep33 = zep3 + zep03
494 my_real,
parameter :: zep34 = zep3 + zep04
495 my_real,
parameter :: zep35 = zep3 + zep05
496 my_real,
parameter :: zep37 = zep3 + zep07
497 my_real,
parameter :: zep38 = zep3 + zep08
498 my_real,
parameter :: zep333 = zep33 + threeem3
499 my_real,
parameter :: zep373 = zep37 + threeem3
500 my_real,
parameter :: zep41 = zep4 + zep01
501 my_real,
parameter :: zep4187 = zep41 + eightem3 + sevenem4
502 my_real,
parameter :: zep43 = zep4 + zep03
503 my_real,
parameter :: zep45 = zep4 + zep05
504 my_real,
parameter :: zep46 = zep4 + zep06
505 my_real,
parameter :: zep47 = zep4 + zep07
506 my_real,
parameter :: zep499 = zep4 + zep09 + nineem3
507 my_real,
parameter :: zep475 = zep4 + zep07 + fiveem3
508 my_real,
parameter :: zep62 = one - zep38
509 my_real,
parameter :: zep65 = one - zep35
510 my_real,
parameter :: zep66 = one - zep34
511 my_real,
parameter :: zep67 = one - zep33
512 my_real,
parameter :: zep71 = one - zep29
513 my_real,
parameter :: zep74 = one - zep26
514 my_real,
parameter :: zep707 = one - zep3 + sevenem3
515 my_real,
parameter :: zep81 = one - zep19
516 my_real,
parameter :: zep82 = one - zep18
517 my_real,
parameter :: zep84 = one - zep16
518 my_real,
parameter :: zep85 = one - zep15
519 my_real,
parameter :: zep87 = one - zep13
520 my_real,
parameter :: zep89 = one - zep11
521 my_real,
parameter :: zep92 = one - zep08
522 my_real,
parameter :: zep93 = one - zep07
523 my_real,
parameter :: zep94 = one - zep06
524 my_real,
parameter :: zep95 = one - zep05
525 my_real,
parameter :: zep96 = one - zep04
526 my_real,
parameter :: zep97 = one - zep03
527 my_real,
parameter :: zep98 = one - zep02
528 my_real,
parameter :: zep99 = one - zep01
529 my_real,
parameter :: zep995 = one - fiveem3
530 my_real,
parameter :: zep997 = one - threeem3
531 my_real,
parameter :: zep999 = one - em03
532 my_real,
parameter :: zep9999 = one - em04
533 my_real,
parameter :: zep99999 = one - em05
534 my_real,
parameter :: zep931 = zep93 + em03
535 my_real,
parameter :: zep015 = zep15 / ten
536 my_real,
parameter :: zep015625 = one / sixty4
537 my_real,
parameter :: zep025 = one / fourty
538 my_real,
parameter :: zep037 = zep03 + sevenem3
539 my_real,
parameter :: zep072169 = zep07 + twoem3 + em04 + sixem5 + nineem6
540 my_real,
parameter :: zep075 = three / fourty
541 my_real,
parameter :: onep01 = one + zep01
542 my_real,
parameter :: onep001 = one + em03
543 my_real,
parameter :: onep0001 = one + em04
544 my_real,
parameter :: onep00001 = one + em05
545 my_real,
parameter :: onep02 = one + zep02
546 my_real,
parameter :: onep03 = one + zep03
547 my_real,
parameter :: onep04 = one + zep04
548 my_real,
parameter :: onep05 = one + zep05
549 my_real,
parameter :: onep06 = one + zep06
550 my_real,
parameter :: onep08 = one + zep08
551 my_real,
parameter :: onep09 = one + zep09
552 my_real,
parameter :: onep1 = one + zep1
553 my_real,
parameter :: onep11 = onep1 + zep01
554 my_real,
parameter :: onep12 = onep1 + zep02
555 my_real,
parameter :: onep13 = onep1 + zep03
556 my_real,
parameter :: onep14 = onep1 + zep04
557 my_real,
parameter :: onep15 = onep1 + zep05
558 my_real,
parameter :: onep16 = onep1 + zep06
559 my_real,
parameter :: onep17 = onep1 + zep07
560 my_real,
parameter :: onep18 = onep1 + zep08
561 my_real,
parameter :: onep19 = onep1 + zep09
562 my_real,
parameter :: onep136 = onep13 + sixem3
563 my_real,
parameter :: onep138 = onep13 + eightem3
564 my_real,
parameter :: onep2 = one + zep2
565 my_real,
parameter :: onep21 = onep2 + zep01
566 my_real,
parameter :: onep23 = onep2 + zep03
567 my_real,
parameter :: onep24 = onep2 + zep04
568 my_real,
parameter :: onep25 = five / four
569 my_real,
parameter :: onep26 = onep2 + zep06
570 my_real,
parameter :: onep27 = onep2 + zep07
571 my_real,
parameter :: onep28 = onep2 + zep08
572 my_real,
parameter :: onep251 = onep25 + em03
573 my_real,
parameter :: onep252 = onep25 + twoem3
574 my_real,
parameter :: onep3 = one + zep3
575 my_real,
parameter :: onep31 = onep3 + zep01
576 my_real,
parameter :: onep32 = onep3 + zep02
577 my_real,
parameter :: onep33 = onep3 + zep03
578 my_real,
parameter :: onep34 = onep3 + zep04
579 my_real,
parameter :: onep35 = onep3 + zep05
580 my_real,
parameter :: onep36 = onep3 + zep06
581 my_real,
parameter :: onep37 = onep3 + zep07
582 my_real,
parameter :: onep38 = onep3 + zep08
583 my_real,
parameter :: onep39 = onep3 + zep09
584 my_real,
parameter :: onep333 = onep33 + threeem3
585 my_real,
parameter :: onep3348= onep33 + fourem3 + eightem4
586 my_real,
parameter :: onep369 = onep36 + nineem3
587 my_real,
parameter :: onep4 = one + zep4
588 my_real,
parameter :: onep41 = onep4 + zep01
589 my_real,
parameter :: onep42 = onep4 + zep02
590 my_real,
parameter :: onep44 = onep4 + zep04
591 my_real,
parameter :: onep45 = onep4 + zep05
592 my_real,
parameter :: onep46 = onep4 + zep06
593 my_real,
parameter :: onep48 = onep4 + zep08
594 my_real,
parameter :: onep49 = onep4 + zep09
595 my_real,
parameter :: onep414 = onep41 + fourem3
596 my_real,
parameter :: onep484 = onep48 + fourem3
597 my_real,
parameter :: onep5 = three_half
598 my_real,
parameter :: onep51 = onep5 + zep01
599 my_real,
parameter :: onep52 = onep5 + zep02
600 my_real,
parameter :: onep57 = onep5 + zep07
601 my_real,
parameter :: onep59 = onep5 + zep09
602 my_real,
parameter :: onep6 = one + zep6
603 my_real,
parameter :: onep62 = onep6 + zep02
604 my_real,
parameter :: onep65 = onep6 + zep05
605 my_real,
parameter :: onep66 = onep6 + zep06
606 my_real,
parameter :: onep67 = onep6 + zep07
607 my_real,
parameter :: onep68 = onep6 + zep08
608 my_real,
parameter :: onep7 = one + zep7
609 my_real,
parameter :: onep72 = onep7 + zep02
610 my_real,
parameter :: onep73 = onep7 + zep03
611 my_real,
parameter :: onep75 = onep7 + zep05
612 my_real,
parameter :: onep78 = onep7 + zep08
613 my_real,
parameter :: onep732 = onep73 + twoem3
614 my_real,
parameter :: onep8 = one + zep8
615 my_real,
parameter :: onep82 = onep8 + zep02
616 my_real,
parameter :: onep84 = onep8 + zep04
617 my_real,
parameter :: onep88 = onep8 + zep08
618 my_real,
parameter :: onep833 = onep8 + zep03 + threeem3
619 my_real,
parameter :: onep8688 = onep8 + zep06 + eightem3 + eightem4
620 my_real,
parameter :: onep9 = one + zep9
621 my_real,
parameter :: onep92 = onep9 + zep02
622 my_real,
parameter :: onep96 = onep9 + zep06
623 my_real,
parameter :: onep98 = onep9 + zep08
624 my_real,
parameter :: twop03 = two + zep03
625 my_real,
parameter :: twop05 = two + zep05
626 my_real,
parameter :: twop06 = two + zep06
627 my_real,
parameter :: twop07 = two + zep07
628 my_real,
parameter :: twop08 = two + zep08
629 my_real,
parameter :: twop1 = two + zep1
630 my_real,
parameter :: twop2 = two + zep2
631 my_real,
parameter :: twop30 = two + zep3
632 my_real,
parameter :: twop4 = two + zep4
633 my_real,
parameter :: twop5 = five / two
634 my_real,
parameter :: twop60 = two + zep6
635 my_real,
parameter :: twop8 = two + zep8
636 my_real,
parameter :: twop11 = twop1 + zep01
637 my_real,
parameter :: twop12 = twop1 + zep02
638 my_real,
parameter :: twop13 = twop1 + zep03
639 my_real,
parameter :: twop16 = twop1 + zep06
640 my_real,
parameter :: twop19 = twop1 + zep09
641 my_real,
parameter :: twop23 = twop2 + zep03
642 my_real,
parameter :: twop24 = twop2 + zep04
643 my_real,
parameter :: twop25 = twop2 + zep05
644 my_real,
parameter :: twop33 = twop30 + zep03
645 my_real,
parameter :: twop38 = twop30 + zep08
646 my_real,
parameter :: twop45 = twop4 + zep05
647 my_real,
parameter :: twop47 = twop4 + zep07
648 my_real,
parameter :: twop52 = twop5 + zep02
649 my_real,
parameter :: twop83 = twop8 + zep03
650 my_real,
parameter :: twop84 = twop8 + zep04
651 my_real,
parameter :: twop85 = twop8 + zep05
652 my_real,
parameter :: twop185 = twop1 + zep08 + fiveem3
653 my_real,
parameter :: twop261 = twop2 + zep06 + em03
654 my_real,
parameter :: twop645 = twop60 + zep04 + fiveem3
655 my_real,
parameter :: threep02 = three + zep02
656 my_real,
parameter :: threep06 = three + zep06
657 my_real,
parameter :: threep14 = three + zep14
658 my_real,
parameter :: threep15 = three + zep15
659 my_real,
parameter :: threep25 = three + zep2 + zep05
660 my_real,
parameter :: threep27 = three + zep2 + zep07
661 my_real,
parameter :: threep44 = three + zep4 + zep04
662 my_real,
parameter :: threep46 = three + zep4 + zep06
663 my_real,
parameter :: threep52 = three + zep5 + zep02
664 my_real,
parameter :: threep72 = three + zep7 + zep02
665 my_real,
parameter :: threep76 = three + zep7 + zep06
666 my_real,
parameter :: threep8 = three + zep8
667 my_real,
parameter :: threep94 = three + zep9 + zep04
668 my_real,
parameter :: threep266 = three + zep2 + zep06 + sixem3
669 my_real,
parameter :: threep464 = three + zep4 + zep06 + fourem3
670 my_real,
parameter :: threep4641= threep464 + em04
671 my_real,
parameter :: fourp186 = four + zep1 + zep08 + sixem3
672 my_real,
parameter :: fourp40 = four + zep4
673 my_real,
parameter :: fourp76 = four + zep7 + zep06
674 my_real,
parameter :: fivep5 = five + zep5
675 my_real,
parameter :: fivep56 = five + zep5 + zep06
676 my_real,
parameter :: fivep76 = five + zep7 + zep06
677 my_real,
parameter :: sip2832 = six + zep2 + zep08 + threeem3 + twoem4
678 my_real,
parameter :: ninep24 = nine + zep2 + zep04
679 my_real,
parameter :: ninep38 = nine + zep3 + zep08
680 my_real,
parameter :: ninep793 = nine + zep7 + zep09 + threeem3
681 my_real,
parameter :: elevenp5 = eleven + zep5
682 my_real,
parameter :: onzep225 = eleven + zep2 + zep02+fiveem3
683 my_real,
parameter :: twelvep6 = twelve + zep6
684 my_real,
parameter :: fourteenp2 = fourteen + zep2
685 my_real,
parameter :: sixteenp5 = sixteen + zep5
686 my_real,
parameter :: o88p9844 = 88. + zep9 + zep08 + fourem3 + fourem4
687 my_real,
parameter :: eightp314em5 = 8314. * em08
688 my_real,
parameter :: onep2ep30 = onep2 * ep30
689 my_real,
parameter :: zep999ep31 = zep999 * ep31
690 my_real,
parameter :: zep54 = one - zep46
691 my_real,
parameter :: onep001em32 = 1001. * em35
692 my_real,
parameter :: ninep637em5 = 9637. * em08
693 my_real,
parameter :: onep3em5 = onep3 * em05
694 my_real,
parameter :: zep426667 = hundred28 / three100
695 my_real,
parameter :: zep222 = ten / fourty5
696 my_real,
parameter :: zep00666666667 = two / three100
697 my_real,
parameter :: zep01666666667 = five / three100
698 my_real,
parameter :: zep444 = hundred / twohundred25
699 my_real,
parameter :: onep222 = one + zep222
700 my_real,
parameter :: onep8333 = one + twenty5 / thirty
701 my_real,
parameter :: twop444 = two + zep444
702 my_real,
parameter :: twop6666666667= eight / three
703 my_real,
parameter :: twop6667 = twop6666666667
704 my_real,
parameter :: fivep333 = sixteen / three
705 my_real,
parameter :: zep1381966= zep13 + eightem3 + em04 + nineem5 + sixem6 + six *em07
706 my_real,
parameter :: zep5854102= zep5 + zep08 + fiveem3 + fourem4 + em05 + two * em07
707 my_real,
parameter :: r_igc = 8.314472
708 my_real,
parameter :: stefboltz= 5.6704 / ep08
709 my_real,
parameter :: sqr2 = sqrt(two)
710 my_real,
parameter :: sqr3 = sqrt(three)
711 my_real,
parameter :: sqr5 = sqrt(five)
712 my_real,
parameter :: sqr6 = sqrt(six)
713 my_real,
parameter :: sqr7 = sqrt(seven)
714 my_real,
parameter :: sqr3_2 = sqrt(three_half)
715 my_real,
parameter :: sqr4_3 = two * sqrt(third)
716 my_real,
parameter :: sqr16_3 = two * sqr4_3
717 my_real,
parameter :: pi = atan2(zero,-one)
719 my_real,
parameter :: em3 = em03
720 my_real,
parameter :: em4 = em04
721 my_real,
parameter :: em5 = em05
722 my_real,
parameter :: em6 = em06
723 my_real,
parameter :: em7 = em07
724 my_real,
parameter :: em8 = em08
725 my_real,
parameter :: em9 = em09
726 my_real,
parameter :: ep9 = ep09
727 my_real,
parameter :: zep80 = zep8
728 my_real,
parameter :: zep90 = zep9
729 my_real,
parameter :: unp20 = onep2
730 my_real,
parameter :: unp30 = onep3
731 my_real,
parameter :: unp70 = onep7
732 my_real,
parameter :: unp90 = onep9
733 my_real,
parameter :: unem3 = em03
735 my_real,
parameter :: unem5 = em05
736 my_real,
parameter :: unem6 = em06
737 my_real,
parameter :: twoem2 = zep02
738 my_real,
parameter :: fiveem2 = zep05
739 my_real,
parameter :: nineem2 = zep09
subroutine m5in2(pm, mat, m151_id, detonators, tb, x, ix, nix)
subroutine puff(iflag, nel, pm, off, eint, mu, mu2, espe, dvol, df, vnew, mat, pnew, dpdm, dpde)