56!|| damping_funct_ini ../engine/source/assembly/damping_funct_ini.f90
67!|| detach_node ../engine/source/engine/node_spliting/detach_node.f90
85!|| eikonal_remove_first ../
starter/source/initial_conditions/detonation/eikonal_remove_first.f90
99!|| fail_lemaitre_s ../engine/source/materials/fail/lemaitre/fail_lemaitre_s.f90
111!|| fpcont2_min_output ../engine/source/output/h3d/h3d_results/fpcont2_min_output.f90
114!|| fractal_element_neighbor ../
starter/source/materials/fail/fractal/fractal_element_neighbor.f90
284!|| s6zforc3 ../engine/source/elements
326!|| sigeps88 ../engine/source/materials/mat/mat088
337!|| spmd_xv_inter_type1 ../engine/source/mpi/nodes/spmd_sd_xv_inter1.f90
357#include "my_real.inc"
358 my_real,
parameter :: infinity = 1e20
359 my_real,
parameter :: zero = 0.
362 my_real,
parameter :: three = 3.
363 my_real,
parameter :: four = 4.
364 my_real,
parameter :: five = 5.
366 my_real,
parameter :: seven = 7.
367 my_real,
parameter :: eight = 8.
368 my_real,
parameter :: nine = 9.
369 my_real,
parameter :: ten = 10.
370 my_real,
parameter :: eleven = 11.
371 my_real,
parameter :: twelve = 12.
372 my_real,
parameter :: fourteen = 14.
373 my_real,
parameter :: fifteen = 15.
374 my_real,
parameter :: sixteen = 16.
375 my_real,
parameter :: eighteen = 18.
376 my_real,
parameter :: twenty = 20.
377 my_real,
parameter :: twenty4 = 24.
378 my_real,
parameter :: twenty5 = 25.
379 my_real,
parameter :: twenty7 = 27.
380 my_real,
parameter :: thirty = 30.
381 my_real,
parameter :: thirty2 = 32.
382 my_real,
parameter :: thirty6 = 36.
383 my_real,
parameter :: forty = 40.
384 my_real,
parameter :: forty5 = 45.
385 my_real,
parameter :: forty8 = 48.
386 my_real,
parameter :: fifty = 50.
387 my_real,
parameter :: fifty4 = 54.
388 my_real,
parameter :: sixty4 = 64.
389 my_real,
parameter :: eighty = 80.
390 my_real,
parameter :: ninety = 90.
391 my_real,
parameter :: eighty16 = 96.
392 my_real,
parameter :: eighty19 = 99.
393 my_real,
parameter :: hundred = 100.
394 my_real,
parameter :: hundred28 = 128.
395 my_real,
parameter :: hundred50 = 150.
396 my_real,
parameter :: hundred80 = 180.
397 my_real,
parameter :: three100 = 300.
398 my_real,
parameter :: twohundred25 = 225.
399 my_real,
parameter :: twohundred95 = 295.
400 my_real,
parameter :: twohundred98 = 298.
401 my_real,
parameter :: twohundred93 = 293.
402 my_real,
parameter :: threehundred24 = 324.
403 my_real,
parameter :: fivehundred12 = 512.
404 my_real,
parameter :: two1000 = 2000.
405 my_real,
parameter :: three1000 = 3000.
406 my_real,
parameter :: ep02 = hundred
407 my_real,
parameter :: ep03 = ep02 * ten
408 my_real,
parameter :: ep04 = ep03 * ten
409 my_real,
parameter :: ep05 = ep04 * ten
410 my_real,
parameter :: ep06 = ep05 * ten
411 my_real,
parameter :: ep07 = ep06 * ten
412 my_real,
parameter :: ep08 = ep07 * ten
413 my_real,
parameter :: ep09 = ep08 * ten
414 my_real,
parameter :: ep10 = ep09 * ten
415 my_real,
parameter :: ep11 = ep10 * ten
416 my_real,
parameter :: ep12 = ep11 * ten
417 my_real,
parameter :: ep13 = ep12 * ten
418 my_real,
parameter :: ep14 = ep13 * ten
419 my_real,
parameter :: ep15 = ep14 * ten
420 my_real,
parameter :: ep16 = ep15 * ten
421 my_real,
parameter :: ep17 = ep16 * ten
422 my_real,
parameter :: ep18 = ep17 * ten
423 my_real,
parameter :: ep19 = ep18 * ten
424 my_real,
parameter :: ep20 = ep19 * ten
425 my_real,
parameter :: ep21 = ep20 * ten
426 my_real,
parameter :: ep29 = ep20 * ep09
427 my_real,
parameter :: ep30 = ep20 * ep10
428 my_real,
parameter :: ep31 = ep30 * ten
429 my_real,
parameter :: ep35 = ep30 * ep05
430 my_real,
parameter :: ep36 = ep30 * ep06
431 my_real,
parameter :: ep37 = ep30 * ep07
432 my_real,
parameter :: em01 = one / ten
433 my_real,
parameter :: em02 = one / ep02
434 my_real,
parameter :: em03 = one / ep03
435 my_real,
parameter :: em04 = one / ep04
436 my_real,
parameter :: em05 = one / ep05
437 my_real,
parameter :: em06 = one / ep06
438 my_real,
parameter :: em07 = one / ep07
439 my_real,
parameter :: em08 = one / ep08
440 my_real,
parameter :: em09 = one / ep09
441 my_real,
parameter :: em10 = one / ep10
442 my_real,
parameter :: em11 = one / ep11
443 my_real,
parameter :: em12 = one / ep12
444 my_real,
parameter :: em13 = one / ep13
445 my_real,
parameter :: em14 = one / ep14
446 my_real,
parameter :: em15 = one / ep15
447 my_real,
parameter :: em16 = one / ep16
448 my_real,
parameter :: em17 = one / ep17
449 my_real,
parameter :: em18 = one / ep18
450 my_real,
parameter :: em19 = one / ep19
451 my_real,
parameter :: em20 = one / ep20
452 my_real,
parameter :: em30 = one / ep30
453 my_real,
parameter :: em35 = one / ep35
454 my_real,
parameter :: zep1 = em01
455 my_real,
parameter :: zep2 = two / ten
456 my_real,
parameter :: zep3 = three / ten
457 my_real,
parameter :: zep4 = four / ten
458 my_real,
parameter :: zep5 = five / ten
459 my_real,
parameter :: zep6 = one - zep4
460 my_real,
parameter :: zep7 = one - zep3
461 my_real,
parameter :: zep8 = one - zep2
462 my_real,
parameter :: zep9 = one - zep1
463 my_real,
parameter :: zep01 = em02
464 my_real,
parameter :: zep02 = two / ep02
465 my_real,
parameter :: zep03 = three / ep02
466 my_real,
parameter :: zep04 = four / ep02
467 my_real,
parameter :: zep05 = five / ep02
468 my_real,
parameter :: zep06 = zep1 - zep04
469 my_real,
parameter :: zep07 = zep1 - zep03
470 my_real,
parameter :: zep08 = zep1 - zep02
471 my_real,
parameter :: zep09 = zep1 - zep01
472 my_real,
parameter :: twoem3 = two / ep03
473 my_real,
parameter :: threeem3 = three / ep03
474 my_realparameter :: fourem3 = four / ep03
475 my_real,
parameter :: fiveem3 = five / ep03
476 my_real,
parameter :: sixem3 = zep01 - fourem3
477 my_real,
parameter :: sevenem3 = zep01 - threeem3
478 my_real,
parameter :: eightem3 = zep01 - twoem3
479 my_real,
parameter :: nineem3 = zep01 - em03
480 my_real,
parameter :: twoem4 = two / ep04
481 my_real,
parameter :: threeem4 = three / ep04
482 my_real,
parameter :: fourem4 = four / ep04
483 my_real,
parameter :: fiveem4 = five / ep04
484 my_real,
parameter :: sixem4 = em03 - fourem4
485 my_real,
parameter :: sevenem4 = em03 - threeem4
486 my_real,
parameter :: eightem4 = em03 - twoem4
487 my_real,
parameter :: nineem4 = em03 - em04
488 my_real,
parameter :: twoem5 = two / ep05
489 my_real,
parameter :: threeem5 = three / ep05
490 my_real,
parameter :: fourem5 = four / ep05
491 my_real,
parameter :: fiveem5 = five / ep05
492 my_real,
parameter :: sixem5 = em04 - fourem5
493 my_real,
parameter :: sevenem5 = em04 - threeem5
494 my_real,
parameter :: eightem5 = em04 - twoem5
495 my_real,
parameter :: nineem5 = em04 - em05
496 my_real,
parameter :: twoem6 = two / ep06
497 my_real,
parameter :: threeem6 = three / ep06
498 my_real,
parameter :: fourem6 = four / ep06
499 my_real,
parameter :: fiveem6 = five / ep06
500 my_real,
parameter :: sixem6 = em05 - fourem6
501 my_real,
parameter :: sevenem6 = em05 - threeem6
502 my_real,
parameter :: eightem6 = em05 - twoem6
503 my_real,
parameter :: nineem6 = em05 - em06
504 my_real,
parameter :: twoep30 = two *ep30
505 my_real,
parameter :: threeep30 = three *ep30
506 my_real,
parameter :: fourep30 = four*ep30
507 my_real,
parameter :: zep9ep30 = nine *ep29
508 my_real,
parameter :: onep1ep30 = eleven * ep29
509 my_real,
parameter :: half = one/two
510 my_real,
parameter :: third = one/three
511 my_real,
parameter :: fourth = one/four
512 my_real,
parameter :: one_fifth = one/five
513 my_real,
parameter :: one_over_6 = one/six
514 my_real,
parameter :: one_over_7 = one/seven
515 my_real,
parameter :: one_over_8 = one/eight
516 my_real,
parameter :: one_over_9 = one/nine
517 my_real,
parameter :: one_over_12 = one/twelve
518 my_real,
parameter :: one_over_16 = one/sixteen
519 my_real,
parameter :: one_over_64 = one/sixty4
520 my_real,
parameter :: one_over_512 = one/fivehundred12
521 my_real,
parameter :: two_third = two/three
522 my_real,
parameter :: three_half = three/two
523 my_real,
parameter :: three_over_4 = three/four
524 my_real,
parameter :: three_over_14 = three/fourteen
525 my_real,
parameter :: three_over_16 = three/sixteen
526 my_real,
parameter :: four_over_3 = four/three
527 my_real,
parameter :: four_over_5 = four/five
528 my_real,
parameter :: four_over_9 = four/nine
529 my_real,
parameter :: five_over_3 = five/three
530 my_real,
parameter :: five_over_4 = five/four
531 my_real,
parameter :: five_over_6 = five/six
532 my_real,
parameter :: six_over_5 = six/five
533 my_real,
parameter :: seven_half = seven/two
534 my_real,
parameter :: seven_over_5 = seven/five
535 my_real,
parameter :: seven_over_9 = seven/nine
536 my_real,
parameter :: nine_over_2 = nine/two
537 my_real,
parameter :: nine_over_10 = nine/ten
538 my_real,
parameter :: nine_over_16 = nine/sixteen
539 my_real,
parameter :: nine_over_20 = nine/twenty
540 my_real,
parameter :: zep11 = zep1 + zep01
541 my_real,
parameter :: zep12 = zep1 + zep02
542 my_real,
parameter :: zep13 = zep1 + zep03
543 my_real,
parameter :: zep14 = zep1 + zep04
544 my_real,
parameter :: zep15 = zep1 + zep05
545 my_real,
parameter :: zep16 = zep1 + zep06
546 my_real,
parameter :: zep17 = zep1 + zep07
547 my_real,
parameter :: zep18 = zep1 + zep08
548 my_real,
parameter :: zep19 = zep1 + zep09
549 my_real,
parameter :: zep21 = zep2 + zep01
550 my_real,
parameter :: zep23 = zep2 + zep03
551 my_real,
parameter :: zep25 = zep2 + zep05
552 my_real,
parameter :: zep26 = zep2 + zep06
553 my_real,
parameter :: zep213 = zep21 + threeem3
554 my_real,
parameter :: zep29 = zep2 + zep09
555 my_real,
parameter :: zep33 = zep3 + zep03
556 my_real,
parameter :: zep34 = zep3 + zep04
557 my_real,
parameter :: zep35 = zep3 + zep05
558 my_real,
parameter :: zep37 = zep3 + zep07
559 my_real,
parameter :: zep38 = zep3 + zep08
560 my_real,
parameter :: zep333 = zep33 + threeem3
561 my_real,
parameter :: zep373 = zep37 + threeem3
562 my_real,
parameter :: zep41 = zep4 + zep01
563 my_real,
parameter :: zep4187 = zep41 + eightem3 + sevenem4
564 my_real,
parameter :: zep43 = zep4 + zep03
565 my_real,
parameter :: zep45 = zep4 + zep05
566 my_real,
parameter :: zep46 = zep4 + zep06
567 my_real,
parameter :: zep47 = zep4 + zep07
568 my_real,
parameter :: zep499 = zep4 + zep09 + nineem3
569 my_real,
parameter :: zep475 = zep4 + zep07 + fiveem3
570 my_real,
parameter :: zep62 = one - zep38
571 my_real,
parameter :: zep65 = one - zep35
572 my_real,
parameter :: zep66 = one - zep34
573 my_real,
parameter :: zep67 = one - zep33
574 my_real,
parameter :: zep71 = one - zep29
575 my_real,
parameter :: zep74 = one - zep26
576 my_real,
parameter :: zep707 = one - zep3 + sevenem3
577 my_real,
parameter :: zep81 = one - zep19
578 my_real,
parameter :: zep82 = one - zep18
579 my_real,
parameter :: zep84 = one - zep16
580 my_real,
parameter :: zep85 = one - zep15
581 my_real,
parameter :: zep87 = one - zep13
582 my_real,
parameter :: zep89 = one - zep11
583 my_real,
parameter :: zep92 = one - zep08
584 my_real,
parameter :: zep93 = one - zep07
585 my_real,
parameter :: zep94 = one - zep06
586 my_real,
parameter :: zep95 = one - zep05
587 my_real,
parameter :: zep96 = one - zep04
588 my_real,
parameter :: zep97 = one - zep03
589 my_real,
parameter :: zep98 = one - zep02
590 my_real,
parameter :: zep99 = one - zep01
591 my_real,
parameter :: zep995 = one - fiveem3
592 my_real,
parameter :: zep997 = one - threeem3
593 my_real,
parameter :: zep999 = one - em03
594 my_real,
parameter :: zep9999 = one - em04
595 my_real,
parameter :: zep99999 = one - em05
596 my_real,
parameter :: zep931 = zep93 + em03
597 my_real,
parameter :: zep015 = zep15 / ten
598 my_real,
parameter :: zep015625 = one / sixty4
599 my_real,
parameter :: zep025 = one / forty
600 my_real,
parameter :: zep037 = zep03 + sevenem3
601 my_real,
parameter :: zep072169 = zep07 + twoem3 + em04 + sixem5 + nineem6
602 my_real,
parameter :: zep075 = three / forty
603 my_real,
parameter :: onep01 = one + zep01
604 my_real,
parameter :: onep001 = one + em03
605 my_real,
parameter :: onep0001 = one + em04
606 my_real,
parameter :: onep00001 = one + em05
607 my_real,
parameter :: onep02 = one + zep02
608 my_real,
parameter :: onep03 = one + zep03
609 my_real,
parameter :: onep04 = one + zep04
610 my_real,
parameter :: onep05 = one + zep05
611 my_real,
parameter :: onep06 = one + zep06
612 my_real,
parameter :: onep08 = one + zep08
613 my_real,
parameter :: onep09 = one + zep09
614 my_real,
parameter :: onep1 = one + zep1
615 my_real,
parameter :: onep11 = onep1 + zep01
616 my_real,
parameter :: onep12 = onep1 + zep02
617 my_real,
parameter :: onep13 = onep1 + zep03
618 my_real,
parameter :: onep14 = onep1 + zep04
619 my_real,
parameter :: onep15 = onep1 + zep05
620 my_real,
parameter :: onep16 = onep1 + zep06
621 my_real,
parameter :: onep17 = onep1 + zep07
622 my_real,
parameter :: onep18 = onep1 + zep08
623 my_real,
parameter :: onep19 = onep1 + zep09
624 my_real,
parameter :: onep136 = onep13 + sixem3
625 my_real,
parameter :: onep138 = onep13 + eightem3
626 my_real,
parameter :: onep2 = one + zep2
627 my_real,
parameter :: onep21 = onep2 + zep01
628 my_real,
parameter :: onep23 = onep2 + zep03
629 my_real,
parameter :: onep24 = onep2 + zep04
630 my_real,
parameter :: onep25 = five / four
631 my_real,
parameter :: onep26 = onep2 + zep06
632 my_real,
parameter :: onep27 = onep2 + zep07
633 my_real,
parameter :: onep28 = onep2 + zep08
634 my_real,
parameter :: onep251 = onep25 + em03
635 my_real,
parameter :: onep252 = onep25 + twoem3
636 my_real,
parameter :: onep3 = one + zep3
637 my_real,
parameter :: onep31 = onep3 + zep01
638 my_real,
parameter :: onep32 = onep3 + zep02
639 my_real,
parameter :: onep33 = onep3 + zep03
640 my_real,
parameter :: onep34 = onep3 + zep04
641 my_real,
parameter :: onep35 = onep3 + zep05
642 my_real,
parameter :: onep36 = onep3 + zep06
643 my_real,
parameter :: onep37 = onep3 + zep07
644 my_real,
parameter :: onep38 = onep3 + zep08
645 my_real,
parameter :: onep39 = onep3 + zep09
646 my_real,
parameter :: onep333 = onep33 + threeem3
647 my_real,
parameter :: onep3348= onep33 + fourem3 + eightem4
648 my_real,
parameter :: onep369 = onep36 + nineem3
649 my_real,
parameter :: onep4 = one + zep4
650 my_real,
parameter :: onep41 = onep4 + zep01
651 my_real,
parameter :: onep42 = onep4 + zep02
652 my_real,
parameter :: onep44 = onep4 + zep04
653 my_real,
parameter :: onep45 = onep4 + zep05
654 my_real,
parameter :: onep46 = onep4 + zep06
655 my_real,
parameter :: onep48 = onep4 + zep08
656 my_real,
parameter :: onep49 = onep4 + zep09
657 my_real,
parameter :: onep414 = onep41 + fourem3
658 my_real,
parameter :: onep484 = onep48 + fourem3
659 my_real,
parameter :: onep5 = three_half
660 my_real,
parameter :: onep51 = onep5 + zep01
661 my_real,
parameter :: onep52 = onep5 + zep02
662 my_real,
parameter :: onep57 = onep5 + zep07
663 my_real,
parameter :: onep59 = onep5 + zep09
664 my_real,
parameter :: onep6 = one + zep6
665 my_real,
parameter :: onep62 = onep6 + zep02
666 my_real,
parameter :: onep65 = onep6 + zep05
667 my_real,
parameter :: onep66 = onep6 + zep06
668 my_real,
parameter :: onep67 = onep6 + zep07
669 my_real,
parameter :: onep68 = onep6 + zep08
670 my_real,
parameter :: onep7 = one + zep7
671 my_real,
parameter :: onep72 = onep7 + zep02
672 my_real,
parameter :: onep73 = onep7 + zep03
673 my_real,
parameter :: onep75 = onep7 + zep05
674 my_real,
parameter :: onep78 = onep7 + zep08
675 my_real,
parameter :: onep732 = onep73 + twoem3
676 my_real,
parameter :: onep8 = one + zep8
677 my_real,
parameter :: onep82 = onep8 + zep02
678 my_real,
parameter :: onep84 = onep8 + zep04
679 my_real,
parameter :: onep88 = onep8 + zep08
680 my_real,
parameter :: onep833 = onep8 + zep03 + threeem3
681 my_real,
parameter :: onep8688 = onep8 + zep06 + eightem3 + eightem4
682 my_real,
parameter :: onep9 = one + zep9
683 my_real,
parameter :: onep92 = onep9 + zep02
684 my_real,
parameter :: onep96 = onep9 + zep06
685 my_real,
parameter :: onep98 = onep9 + zep08
686 my_real,
parameter :: twop03 = two + zep03
687 my_real,
parameter :: twop05 = two + zep05
688 my_real,
parameter :: twop06 = two + zep06
689 my_real,
parameter :: twop07 = two + zep07
690 my_real,
parameter :: twop08 = two + zep08
691 my_real,
parameter :: twop1 = two + zep1
692 my_real,
parameter :: twop2 = two + zep2
693 my_real,
parameter :: twop30 = two + zep3
694 my_real,
parameter :: twop4 = two + zep4
695 my_real,
parameter :: twop5 = five / two
696 my_real,
parameter :: twop60 = two + zep6
697 my_real,
parameter :: twop8 = two + zep8
698 my_real,
parameter :: twop11 = twop1 + zep01
699 my_real,
parameter :: twop12 = twop1 + zep02
700 my_real,
parameter :: twop13 = twop1 + zep03
701 my_real,
parameter :: twop16 = twop1 + zep06
702 my_real,
parameter :: twop19 = twop1 + zep09
703 my_real,
parameter :: twop23 = twop2 + zep03
704 my_real,
parameter :: twop24 = twop2 + zep04
705 my_real,
parameter :: twop25 = twop2 + zep05
706 my_real,
parameter :: twop33 = twop30 + zep03
707 my_real,
parameter :: twop38 = twop30 + zep08
708 my_real,
parameter :: twop45 = twop4 + zep05
709 my_real,
parameter :: twop47 = twop4 + zep07
710 my_real,
parameter :: twop52 = twop5 + zep02
711 my_real,
parameter :: twop83 = twop8 + zep03
712 my_real,
parameter :: twop84 = twop8 + zep04
713 my_real,
parameter :: twop85 = twop8 + zep05
714 my_real,
parameter :: twop185 = twop1 + zep08 + fiveem3
715 my_real,
parameter :: twop261 = twop2 + zep06 + em03
716 my_real,
parameter :: twop645 = twop60 + zep04 + fiveem3
717 my_real,
parameter :: threep02 = three + zep02
718 my_real,
parameter :: threep06 = three + zep06
719 my_real,
parameter :: threep14 = three + zep14
720 my_real,
parameter :: threep15 = three + zep15
722 my_real,
parameter :: threep27 = three + zep2 + zep07
723 my_realparameter :: threep44 = three + zep4
724 my_real,
parameter :: threep46 = three + zep4 + zep06
725 my_real,
parameter :: threep52 = three + zep5 + zep02
726 my_real,
parameter :: threep72 = three + zep7 + zep02
727 my_real,
parameter :: threep76 = three + zep7 + zep06
730 my_real,
parameter :: threep266
731 my_real,
parameter :: threep464 = three + zep4 + zep06 + fourem3
732 my_real,
parameter :: threep4641= threep464 + em04
733 my_real,
parameter :: fourp186 = four + zep1 + zep08 + sixem3
735 my_real,
parameter :: fourp76 = four + zep7 + zep06
736 my_real,
parameter :: fivep5 = five + zep5
737 my_real,
parameter :: fivep56 = five + zep5 + zep06
739 my_real,
parameter :: sip2832 = six + zep2 + zep08 + threeem3 + twoem4
740 my_real,
parameter :: ninep24 = nine + zep2 + zep04
745 my_real,
parameter :: twelvep6 = twelve + zep6
746 my_real,
parameter :: fourteenp2 = fourteen + zep2
747 my_real,
parameter :: sixteenp5 = sixteen + zep5
748 my_real,
parameter :: o88p9844 = 88. + zep9 + zep08 + fourem3 + fourem4
749 my_real,
parameter :: eightp314em5 = 8314. * em08
751 parameter :: zep999ep31 = zep999 * ep31
753 my_realparameter :: onep001em32
754 my_realparameter :: ninep637em5 = 9637. * em08
756 my_real,
parameter :: zep426667 = hundred28 / three100
757 my_real,
parameter :: zep222 = ten / forty5
758 parameter :: zep00666666667 = two / three100
759 my_real,
parameter :: zep01666666667 = five / three100
760 my_real,
parameter :: zep444 = hundred / twohundred25
761 my_real,
parameter :: onep222 = one
763 my_real,
parameter :: twop444 = two + zep444
764 my_real,
parameter :: twop6666666667= eight / three
765 my_real,
parameter :: twop6667 = twop6666666667
766 my_real,
parameter :: fivep333 = sixteen / three
767 my_real,
parameter :: zep1381966= zep13 + eightem3 + em04 + nineem5 + sixem6 + six *em07
768 my_real,
parameter :: zep5854102= zep5 + zep08 + fiveem3 + fourem4 + em05 + two * em07
769 my_real,
parameter :: r_igc = 8.314472
770 my_real,
parameter :: stefboltz= 5.6704 / ep08
771 my_real,
parameter :: sqr2 = sqrt(two)
772 my_real,
parameter :: sqr3 = sqrt(three)
773 my_real,
parameter :: sqr5 = sqrt(five)
774 my_real,
parameter :: sqr6 = sqrt(six)
775 my_real,
parameter :: sqr7 = sqrt(seven)
776 my_real,
parameter :: sqr3_2 = sqrt(three_half)
777 my_real,
parameter :: sqr4_3 = two * sqrt(third)
778 my_real,
parameter :: sqr16_3 = two * sqr4_3
779 my_real,
parameter :: pi = atan2(zero,-one)
780 my_real,
parameter :: invpi = one / pi
781 my_real,
parameter :: em3 = em03
782 my_real,
parameter :: em4 = em04
783 my_real,
parameter :: em5 = em05
784 my_real,
parameter :: em6 = em06
785 my_real,
parameter :: em7 = em07
786 my_real,
parameter :: em8 = em08
787 my_real,
parameter :: em9 = em09
788 my_real,
parameter :: ep9 = ep09
789 my_real,
parameter :: zep80 = zep8
790 my_real,
parameter :: zep90 = zep9
791 my_real,
parameter :: unp20 = onep2
792 my_real,
parameter :: unp30 = onep3
793 my_real,
parameter :: unp70 = onep7
794 my_real,
parameter :: unp90 = onep9
795 my_real,
parameter :: unem3 = em03
796 my_real,
parameter :: unem4 = em04
797 my_real,
parameter :: unem5 = em05
798 my_real,
parameter :: unem6 = em06
799 my_real,
parameter :: twoem2 = zep02
800 my_real,
parameter :: fiveem2 = zep05
801 my_real,
parameter :: nineem2 = zep09
802 my_real,
parameter :: zep3571429 = half- one_over_7
803 my_real,
parameter :: zep2142857 = half- two*one_over_7
804 my_real,
parameter :: zep0714286 = half- three*one_over_7
805 my_real,
parameter :: zep375 = half- one_over_8
806 my_real,
parameter :: zep3888889 = half- one_over_9
807 my_real,
parameter :: zep2777778 = half- two*one_over_9
808 my_real,
parameter :: zep1666667 = half- three*one_over_9
809 my_real,
parameter :: zep0555555 = half- four*one_over_9