477 1 NB ,LEN ,BUFR ,NSV ,A ,
478 2 STIFN ,VISCN ,IBC ,ISECIN ,NOINT ,
479 3 IBAG ,ICODT ,SECFCUM,NSTRF ,ICONTACT,
480 4 FCONT ,INACTI ,IADM ,INTTH ,FTHE ,
481 5 CONDN ,H3D_DATA, MULTI_FVM,NIN ,TAGNCONT,
482 6 KLOADPINTER,LOADPINTER,LOADP_HYD_INTER ,
483 . INTCAREA,FSAV ,PARAMETERS,NODADT_THERM)
494#include "implicit_f.inc"
498#include "scr07_c.inc"
499#include "scr14_c.inc"
500#include "scr16_c.inc"
501#include "scr18_c.inc"
502#include "com01_c.inc"
503#include "com04_c.inc"
504#include "com06_c.inc"
505#include "com08_c.inc"
506#include "impl1_c.inc"
510 TYPE(output_),
intent(inout) :: OUTPUT
511 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT, INACTI,
512 . NSV(*), ICODT(*), NSTRF(*),ICONTACT(*),
513 . TAGNCONT(NLOADP_HYD_INTER,*),KLOADPINTER(*),LOADPINTER(*),
514 . LOADP_HYD_INTER(*),
516 INTEGER,
INTENT(in) :: NIN,
517 INTEGER,
INTENT(IN) :: NODADT_THERM
519 . BUFR(,*), A(3,*), STIFN(*), VISCN(*),
520 . SECFCUM(7,NUMNOD,NSECT),
521 . FCONT(3,*),FTHE(*),CONDN(*)
522 my_real,
INTENT(INOUT) :: FSAV(*)
523 TYPE(H3D_DATABASE) :: H3D_DATA
524 TYPE(MULTI_FVM_STRUCT),
INTENT(INOUT) :: MULTI_FVM
525 TYPE (PARAMETERS_) ,
INTENT(IN):: PARAMETERS
529 my_real :: mass, fsav29
531 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER,PP,PPL,NP
536 IF (multi_fvm%INT18_GLOBAL_LIST(nin))
THEN
540 multi_fvm%FORCE_INT(1, ibric) = multi_fvm%FORCE_INT(1, ibric) + bufr(2,i)
541 multi_fvm%FORCE_INT(2, ibric) = multi_fvm%FORCE_INT(2, ibric) + bufr(3,i)
542 multi_fvm%FORCE_INT(3, ibric) = multi_fvm%FORCE_INT(3, ibric) + bufr(4,i)
549 a(1,nod) = a(1,nod) + bufr(2,i)
550 a(2,nod) = a(2,nod) + bufr(3,i)
551 a(3,nod) = a(3,nod) + bufr(4,i)
552 stifn(nod) = stifn(nod) + bufr(5,i)
553 IF(kdtint /= 0) viscn(nod) = viscn(nod) + bufr(6,i)
560 a(1,nod) = a(1,nod) + bufr(2,i)
561 a(2,nod) = a(2,nod) + bufr(3,i)
562 a(3,nod) = a(3,nod) + bufr(4,i)
563 stifn(nod) = stifn(nod) + bufr(5,i)
564 fthe(nod) = fthe(nod) + bufr(6,i)
565 IF(nodadt_therm == 1) condn(nod) = condn(nod) + bufr(7,i)
571 a(1,nod) = a(1,nod) + bufr(2,i)
572 a(2,nod) = a(2,nod) + bufr(3,i)
573 a(3,nod) = a(3,nod) + bufr(4,i)
574 stifn(nod) = stifn(nod) + bufr(5,i)
575 viscn(nod) = viscn(nod) + bufr(6,i)
576 fthe(nod) = fthe(nod) + bufr(7,i)
577 IF(nodadt_therm == 1) condn(nod) = condn(nod) + bufr(8,i)
584 IF (inconv == 1)
THEN
585 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
586 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
587 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
588 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
593 fcont(1,nod)=fcont(1,nod)+bufr(2,i)
594 fcont(2,nod)=fcont(2,nod)+bufr(3,i)
595 fcont(3,nod)=fcont(3,nod)+bufr(4,i)
600 IF(nintloadp > 0)
THEN
604 DO np = kloadpinter(nin)+1, kloadpinter(nin+1)
606 ppl = loadp_hyd_inter(pp)
607 tagncont(ppl,nod) = 1
615 IF(nstrf(1)+nstrf(2)/=0)
THEN
620 IF(nstrf(k1s)==noint)
THEN
625 IF(secfcum(4,nod,i)==1.)
THEN
626 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
627 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
628 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
640 IF((ibag/=0.AND.inacti/=7).OR.
641 . (iadm/=0).OR.(idamp_rdof/=0))
THEN
644 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
645 + bufr(4,i)/=zero)
THEN
655 ibcs = ibc - 8 * ibcm
661 CALL ibcoff(ibcs,icodt(nod))
679 1 NB ,LEN ,BUFR ,IRECTS ,A ,
680 2 STIFN ,VISCN ,IBC ,ISECIN ,NOINT ,
681 3 IBAG ,ICODT ,SECFCUM,NSTRF ,ICONTACT,
682 4 FCONT ,INTTH ,FTHE ,CONDN ,H3D_DATA,
683 5 TAGNCONT,KLOADPINTER,LOADPINTER,LOADP_HYD_INTER,NODADT_THERM)
692#include "implicit_f.inc"
696#include "scr07_c.inc"
697#include "scr14_c.inc"
698#include "scr16_c.inc"
699#include "scr18_c.inc"
700#include "com01_c.inc"
701#include "com04_c.inc"
702#include "com06_c.inc"
703#include "com08_c.inc"
704#include "impl1_c.inc"
708 TYPE(output_),
intent(inout) :: OUTPUT
709 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT,INTTH,
710 . irects(2,*), icodt(*), nstrf(*),icontact(*),
711 . tagncont(nloadp_hyd_inter,*),
712 . kloadpinter(*),loadpinter(*),loadp_hyd_inter(*)
713 INTEGER,
INTENT(IN) :: NODADT_THERM
715 . BUFR(LEN,*), A(3,*), STIFN(*), VISCN(*),
716 . SECFCUM(7,NUMNOD,NSECT),
717 . FCONT(3,*),FTHE(*),CONDN(*)
718 TYPE(H3D_DATABASE) :: H3D_DATA
722 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, ,PP,PPL,NOD1,NOD2,NP
731 a(1,nod) = a(1,nod) + bufr(2,i)
732 a(2,nod) = a(2,nod) + bufr(3,i)
733 a(3,nod) = a(3,nod) + bufr(4,i)
734 stifn(nod) = stifn(nod) + bufr(5,i)
736 a(1,nod) = a(1,nod) + bufr(6,i)
737 a(2,nod) = a(2,nod) + bufr(7,i)
738 a(3,nod) = a(3,nod) + bufr(8,i)
739 stifn(nod) = stifn(nod) + bufr(9,i)
745 a(1,nod) = a(1,nod) + bufr(2,i)
746 a(2,nod) = a(2,nod) + bufr(3,i)
747 a(3,nod) = a(3,nod) + bufr(4,i)
748 stifn(nod) = stifn(nod) + bufr(5,i)
749 viscn(nod) = viscn(nod) + bufr(6,i)
751 a(1,nod) = a(1,nod) + bufr(7,i)
752 a(2,nod) = a(2,nod) + bufr(8,i)
753 a(3,nod) = a(3,nod) + bufr(9,i)
754 stifn(nod) = stifn(nod) + bufr(10,i)
755 viscn(nod) = viscn(nod) + bufr(11,i)
759 IF(nodadt_therm ==1 )
THEN
764 a(1,nod) = a(1,nod) + bufr(2,i)
765 a(2,nod) = a(2,nod) + bufr(3,i)
766 a(3,nod) = a(3,nod) + bufr(4,i)
767 stifn(nod) = stifn(nod) + bufr(5,i)
768 fthe(nod) = fthe(nod) + bufr(6,i)
769 condn(nod) = condn(nod) + bufr(7,i)
771 a(1,nod) = a(1,nod) + bufr(8,i)
772 a(2,nod) = a(2,nod) + bufr(9,i)
773 a(3,nod) = a(3,nod) + bufr(10,i)
774 stifn(nod) = stifn(nod) + bufr(11,i)
775 fthe(nod) = fthe(nod) + bufr(12,i)
776 condn(nod) = condn(nod) + bufr(13,i)
782 a(1,nod) = a(1,nod) + bufr(2,i)
783 a(2,nod) = a(2,nod) + bufr(3,i)
784 a(3,nod) = a(3,nod) + bufr(4,i)
785 stifn(nod) = stifn(nod) + bufr(5,i)
786 viscn(nod) = viscn(nod) + bufr(6,i)
787 fthe(nod) = fthe(nod) + bufr(7,i)
788 condn(nod) = condn(nod) + bufr(8,i)
790 a(1,nod) = a(1,nod) + bufr(9,i)
791 a(2,nod) = a(2,nod) + bufr(10,i)
792 a(3,nod) = a(3,nod) + bufr(11,i)
793 stifn(nod) = stifn(nod) + bufr(12,i)
794 viscn(nod) = viscn(nod) + bufr(13,i)
795 fthe(nod) = fthe(nod) + bufr(14,i)
796 condn(nod) = condn(nod) + bufr(15,i)
804 a(1,nod) = a(1,nod) + bufr(2,i)
805 a(2,nod) = a(2,nod) + bufr(3,i)
806 a(3,nod) = a(3,nod) + bufr(4,i)
807 stifn(nod) = stifn(nod) + bufr(5,i)
808 fthe(nod) = fthe(nod) + bufr(6,i)
810 a(1,nod) = a(1,nod) + bufr(7,i)
811 a(2,nod) = a(2,nod) + bufr(8,i)
812 a(3,nod) = a(3,nod) + bufr(9,i)
813 stifn(nod) = stifn(nod) + bufr(10,i)
814 fthe(nod) = fthe(nod) + bufr(11,i)
820 a(1,nod) = a(1,nod) + bufr(2,i)
821 a(2,nod) = a(2,nod) + bufr(3,i)
822 a(3,nod) = a(3,nod) + bufr(4,i)
823 stifn(nod) = stifn(nod) + bufr(5,i)
824 viscn(nod) = viscn(nod) + bufr(6,i)
825 fthe(nod) = fthe(nod) + bufr(7,i)
827 a(1,nod) = a(1,nod) + bufr(8,i)
828 a(2,nod) = a(2,nod) + bufr(9,i)
829 a(3,nod) = a(3,nod) + bufr(10,i)
830 stifn(nod) = stifn(nod) + bufr(11,i)
831 viscn(nod) = viscn(nod) + bufr(12,i)
832 fthe(nod) = fthe(nod) + bufr(13,i)
841 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
842 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
843 . (manim>=4.AND.manim
844 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
850 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
851 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
852 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
854 fcont(1,nod)=fcont(1,nod)+ bufr(6,i)
855 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
856 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
862 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
863 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
864 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
866 fcont(1,nod)=fcont(1,nod)+ bufr(7,i)
867 fcont(2,nod)=fcont(2,nod)+ bufr(8,i)
868 fcont(3,nod)=fcont(3,nod)+ bufr(9,i)
874 IF(nintloadp > 0)
THEN
879 DO np = kloadpinter(noint)+1, kloadpinter(noint+1)
881 ppl = loadp_hyd_inter(pp)
882 tagncont(ppl,nod1) = 1
883 tagncont(ppl,nod2) = 1
891 IF(nstrf(1)+nstrf(2)/=0)
THEN
896 IF(nstrf(k1s)==noint)
THEN
902 IF(secfcum(4,nod,i)==1.)
THEN
903 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
904 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
905 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
908 IF(secfcum(4,nod,i)==1.)
THEN
909 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(6,ii)
910 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(7,ii)
911 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(8,ii)
918 IF(secfcum(4,nod,i)==1.)
THEN
919 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
920 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
921 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
924 IF(secfcum(4,nod,i)==1.)
THEN
925 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(7,ii)
926 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(8,ii)
927 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(9,ii)
940 IF ((ibag/=0).OR.(idamp_rdof/=0))
THEN
944 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
945 + bufr(4,i)/=zero)
THEN
950 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
951 + bufr(8,i)/=zero)
THEN
958 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
959 + bufr(4,i)/=zero)
THEN
964 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
965 + bufr(9,i)/=zero)
THEN
986 1 NB ,LEN ,BUFR ,NELEMS ,A ,
987 2 STIFN ,FCONT ,IXS ,IXS16 ,FROTS ,
994 use element_mod ,
only :nixs
998#include "implicit_f.inc"
1002#include "scr07_c.inc"
1003#include "scr14_c.inc"
1004#include "scr16_c.inc"
1005#include "com06_c.inc"
1006#include "com08_c.inc"
1010 TYPE(output_),
intent(inout) :: output
1011 INTEGER nb, len, nn,
1012 . nelems(*), IXS(NIXS,*) ,IXS16(8,*)
1014 . BUFR(LEN,*), A(3,*), STIFN(*),
1015 . FCONT(3,*), FROTS(7,*)
1016 TYPE(H3D_DATABASE) :: H3D_DATA
1025 nn = nint(bufr(1,i))
1029 a(1,nod) = a(1,nod) + bufr(2,i)
1030 a(2,nod) = a(2,nod) + bufr(3,i)
1031 a(3,nod) = a(3,nod) + bufr(4,i)
1032 stifn(nod) = stifn(nod) + bufr(5,i)
1034 a(1,nod) = a(1,nod) + bufr(6,i)
1035 a(2,nod) = a(2,nod) + bufr(7,i)
1036 a(3,nod) = a(3,nod) + bufr(8,i)
1037 stifn(nod) = stifn(nod) + bufr(9,i)
1039 a(1,nod) = a(1,nod) + bufr(10,i)
1040 a(2,nod) = a(2,nod) + bufr(11,i)
1041 a(3,nod) = a(3,nod) + bufr(12,i)
1042 stifn(nod) = stifn(nod) + bufr(13,i)
1044 a(1,nod) = a(1,nod) + bufr(14,i)
1045 a(2,nod) = a(2,nod) + bufr(15,i)
1046 a(3,nod) = a(3,nod) + bufr(16,i)
1047 stifn(nod) = stifn(nod) + bufr(17,i)
1049 a(1,nod) = a(1,nod) + bufr(18,i)
1050 a(2,nod) = a(2,nod) + bufr(19,i)
1051 a(3,nod) = a(3,nod) + bufr(20,i)
1052 stifn(nod) = stifn(nod) + bufr(21,i)
1054 a(1,nod) = a(1,nod) + bufr(22,i)
1055 a(2,nod) = a(2,nod) + bufr(23,i)
1056 a(3,nod) = a(3,nod) + bufr(24,i)
1057 stifn(nod) = stifn(nod) + bufr(25,i)
1059 a(1,nod) = a(1,nod) + bufr(26,i)
1060 a(2,nod) = a(2,nod) + bufr(27,i)
1061 a(3,nod) = a(3,nod) + bufr(28,i)
1062 stifn(nod) = stifn(nod) + bufr(29,i)
1064 a(1,nod) = a(1,nod) + bufr(30,i)
1065 a(2,nod) = a(2,nod) + bufr(31,i)
1066 a(3,nod) = a(3,nod) + bufr(32,i)
1067 stifn(nod) = stifn(nod) + bufr(33,i)
1070 a(1,nod) = a(1,nod) + bufr(34,i)
1071 a(2,nod) = a(2,nod) + bufr(35,i)
1072 a(3,nod) = a(3,nod) + bufr(36,i)
1073 stifn(nod) = stifn(nod) + bufr(37,i)
1075 a(1,nod) = a(1,nod) + bufr(38,i)
1076 a(2,nod) = a(2,nod) + bufr(39,i)
1077 a(3,nod) = a(3,nod) + bufr(40,i)
1078 stifn(nod) = stifn(nod) + bufr(41,i)
1080 a(1,nod) = a(1,nod) + bufr(42,i)
1081 a(2,nod) = a(2,nod) + bufr(43,i)
1082 a(3,nod) = a(3,nod) + bufr(44,i)
1083 stifn(nod) = stifn(nod) + bufr(45,i)
1085 a(1,nod) = a(1,nod) + bufr(46,i)
1086 a(2,nod) = a(2,nod) + bufr(47,i)
1087 a(3,nod) = a(3,nod) + bufr(48,i)
1088 stifn(nod) = stifn(nod) + bufr(49,i)
1090 a(1,nod) = a(1,nod) + bufr(50,i)
1091 a(2,nod) = a(2,nod) + bufr(51,i)
1092 a(3,nod) = a(3,nod) + bufr(52,i)
1093 stifn(nod) = stifn(nod) + bufr(53,i)
1095 a(1,nod) = a(1,nod) + bufr(54,i)
1096 a(2,nod) = a(2,nod) + bufr(55,i)
1097 a(3,nod) = a(3,nod) + bufr(56,i)
1098 stifn(nod) = stifn(nod) + bufr(57,i)
1100 a(1,nod) = a(1,nod) + bufr(58,i)
1101 a(2,nod) = a(2,nod) + bufr(59,i)
1102 a(3,nod) = a(3,nod) + bufr(60,i)
1103 stifn(nod) = stifn(nod) + bufr(61,i)
1105 a(1,nod) = a(1,nod) + bufr(62,i)
1106 a(2,nod) = a(2,nod) + bufr(63,i)
1107 a(3,nod) = a(3,nod) + bufr(64,i)
1108 stifn(nod) = stifn(nod) + bufr(65,i)
1110 frots(1,nn) = frots(1,nn) + bufr(66,i)
1111 frots(2,nn) = frots(2,nn) + bufr(67,i)
1112 frots(3,nn) = frots(3,nn) + bufr(68,i)
1113 frots(4,nn) = frots(4,nn) + bufr(69,i)
1118 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
1119 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP
1120 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
1121 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
1124 ne = nint(bufr(1,i))
1127 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
1128 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
1129 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
1131 fcont(1,nod)=fcont(1,nod)+ bufr(6,i)
1132 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
1133 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
1135 fcont(1,nod)=fcont(1,nod)+ bufr(10,i)
1136 fcont(2,nod)=fcont(2,nod)+ bufr(11,i)
1137 fcont(3,nod)=fcont(3,nod)+ bufr(12,i)
1139 fcont(1,nod)=fcont(1,nod)+ bufr(14,i)
1140 fcont(2,nod)=fcont(2,nod)+ bufr(15,i)
1141 fcont(3,nod)=fcont(3,nod)+ bufr(16,i)
1143 fcont(1,nod)=fcont(1,nod)+ bufr(18,i)
1144 fcont(2,nod)=fcont(2,nod)+ bufr(19,i)
1145 fcont(3,nod)=fcont(3,nod)+ bufr(20,i)
1147 fcont(1,nod)=fcont(1,nod)+ bufr(22,i)
1148 fcont(2,nod)=fcont(2,nod)+ bufr(23,i)
1149 fcont(3,nod)=fcont(3,nod)+ bufr(24,i)
1151 fcont(1,nod)=fcont(1,nod)+ bufr(26,i)
1152 fcont(2,nod)=fcont(2,nod)+ bufr(27,i)
1153 fcont(3,nod)=fcont(3,nod)+ bufr(28,i)
1155 fcont(1,nod)=fcont(1,nod)+ bufr(30,i)
1156 fcont(2,nod)=fcont(2,nod)+ bufr(31,i)
1157 fcont(3,nod)=fcont(3,nod)+ bufr(32,i)
1160 fcont(1,nod)=fcont(1,nod)+ bufr(34,i)
1161 fcont(2,nod)=fcont(2,nod)+ bufr(35,i)
1162 fcont(3,nod)=fcont(3,nod)+ bufr(36,i)
1164 fcont(1,nod)=fcont(1,nod)+ bufr(38,i)
1165 fcont(2,nod)=fcont(2,nod)+ bufr(39,i)
1166 fcont(3,nod)=fcont(3,nod)+ bufr(40,i)
1168 fcont(1,nod)=fcont(1,nod)+ bufr(42,i)
1169 fcont(2,nod)=fcont(2,nod)+ bufr(43,i)
1170 fcont(3,nod)=fcont(3,nod)+ bufr(44,i)
1172 fcont(1,nod)=fcont(1,nod)+ bufr(46,i)
1173 fcont(2,nod)=fcont(2,nod)+ bufr(47,i)
1174 fcont(3,nod)=fcont(3,nod)+ bufr(48,i)
1176 fcont(1,nod)=fcont(1,nod)+ bufr(50,i)
1177 fcont(2,nod)=fcont(2,nod)+ bufr(51,i)
1178 fcont(3,nod)=fcont(3,nod)+ bufr(52,i)
1180 fcont(1,nod)=fcont(1,nod)+ bufr(54,i)
1181 fcont(2,nod)=fcont(2,nod)+ bufr(55,i)
1182 fcont(3,nod)=fcont(3,nod)+ bufr(56,i)
1184 fcont(1,nod)=fcont(1,nod)+ bufr(58,i)
1185 fcont(2,nod)=fcont(2,nod)+ bufr(59,i)
1186 fcont(3,nod)=fcont(3,nod)+ bufr(60,i)
1188 fcont(1,nod)=fcont(1,nod)+ bufr(62,i)
1189 fcont(2,nod)=fcont(2,nod)+ bufr(63,i)
1190 fcont(3,nod)=fcont(3,nod)+ bufr(64,i)
1209 1 NB ,LEN ,BUFR ,NSV ,A ,
1210 2 STIFN ,VISCN ,IBC ,ISECIN ,NOINT ,
1211 3 IBAG ,ICODT ,SECFCUM,NSTRF ,ICONTACT,
1212 4 FCONT ,INACTI ,IADM ,INTTH ,DAANC6 ,
1213 5 FTHE ,NLG ,ALPHAK ,H3D_DATA)
1222#include "implicit_f.inc"
1226#include "scr05_c.inc"
1227#include "scr07_c.inc"
1228#include "scr14_c.inc"
1229#include "scr16_c.inc"
1230#include "scr18_c.inc"
1231#include "com01_c.inc"
1232#include "com04_c.inc"
1233#include "com06_c.inc"
1234#include "com08_c.inc"
1238 TYPE(output_),
intent(inout) :: OUTPUT
1239 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT, INACTI,
1240 . NSV(*), ICODT(*), NSTRF(*), NLG(*),
1241 . ICONTACT(*), IADM,INTTH
1243 . BUFR(LEN,*), A(3,*), STIFN(*), VISCN(*),
1244 . SECFCUM(7,NUMNOD,NSECT),
1245 . FCONT(3,*),FTHE(*), ALPHAK(3,*)
1246 DOUBLE PRECISION DAANC6(3,6,*)
1251 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, , NBINTER, INC,ISIGN
1252 DOUBLE PRECISION DAANC6L(3,6)
1257 IF(INTTH == 0 ) then
1260 n = nsv(nint(bufr(1,i)))
1262 a(1,nod) = a(1,nod) + bufr(2,i)
1263 a(2,nod) = a(2,nod) + bufr(3,i)
1264 a(3,nod) = a(3,nod) + bufr(4,i)
1265 stifn(nod) = stifn(nod) + bufr(5,i)
1268 IF(alphak(2,n) < zero .or. bufr(6,i) < zero)isign = -1
1269 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(6,i)))
1271 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1272 daanc6(1,2,n) = daanc6(1,2,n)+daanc6l(1,2)
1273 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1274 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1275 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1276 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1277 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1278 daanc6(2,2,n) = daanc6(2,2,n)+daanc6l(2,2)
1279 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1280 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1281 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1282 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1283 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1284 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1285 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1286 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1287 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1288 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1292 n = nsv(nint(bufr(1,i)))
1294 a(1,nod) = a(1,nod) + bufr(2,i)
1295 a(2,nod) = a(2,nod) + bufr(3,i)
1296 a(3,nod) = a(3,nod) + bufr(4,i)
1297 stifn(nod) = stifn(nod) + bufr(5,i)
1298 viscn(nod) = viscn(nod) + bufr
1300 IF(alphak(2,n) < zero .or. bufr(7,i) < zero)isign = -1
1301 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(7,i)))
1304 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1305 daanc6(1,2,n) = daanc6(1,2,n)+daanc6l(1,2)
1306 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1307 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1308 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1309 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1310 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1311 daanc6(2,2,n) = daanc6(2,2,n
1312 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1313 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1314 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1315 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1316 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1317 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1318 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1319 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1320 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1321 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1328 n = nsv(nint(bufr(1,i)))
1330 a(1,nod) = a(1,nod) + bufr(2,i)
1331 a(2,nod) = a(2,nod) + bufr(3,i)
1332 a(3,nod) = a(3,nod) + bufr(4,i)
1333 stifn(nod) = stifn(nod) + bufr(5,i)
1334 fthe(nod) = fthe(nod) + bufr(6,i)
1337 IF(alphak(2,n) < zero .or. bufr(7,i) < zero)isign = -1
1338 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(7,i)))
1340 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1341 daanc6(1,2,n) = daanc6(1,2,n)+daanc6l(1,2)
1342 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1343 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1344 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1345 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1346 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1347 daanc6(2,2,n) = daanc6(2,2,n)+daanc6l(2,2)
1348 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1349 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1350 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1351 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1352 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1353 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1354 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1355 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1356 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1357 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1361 n = nsv(nint(bufr(1,i)))
1363 a(1,nod) = a(1,nod) + bufr(2,i)
1364 a(2,nod) = a(2,nod) + bufr(3,i)
1365 a(3,nod) = a(3,nod) + bufr(4,i)
1366 stifn(nod) = stifn(nod) + bufr(5,i)
1367 viscn(nod) = viscn(nod) + bufr(6,i)
1368 fthe(nod) = fthe(nod) + bufr(7,i)
1371 IF(alphak(2,n) < zero .or. bufr(8,i) < zero)isign = -1
1372 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(8,i)))
1374 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1375 daanc6(1,2,n) = daanc6(1,2,n)+daanc6l(1,2)
1376 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1377 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1378 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1379 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1380 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1381 daanc6(2,2,n) = daanc6(2,2,n)+daanc6l(2,2)
1382 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1383 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1384 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1385 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1386 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1387 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1388 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1389 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1390 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1391 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1398 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
1399 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
1400 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
1401 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
1406 fcont(1,nod)=fcont(1,nod)+bufr(2,i)
1407 fcont(2,nod)=fcont(2,nod)+bufr(3,i)
1408 fcont(3,nod)=fcont(3,nod)+bufr(4,i)
1415 IF(nstrf(1)+nstrf(2)/=0)
THEN
1417 nbinter=nstrf(k0+14)
1420 IF(nstrf(k1s)==noint)
THEN
1423 n = nint(bufr(1,ii))
1425 IF(secfcum(4,nod,i)==1.)
THEN
1426 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
1427 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
1428 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
1440 IF((ibag/=0.AND.inacti/=7).OR.
1444 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
1445 + bufr(4,i)/=zero)
THEN
1455 ibcs = ibc - 8 * ibcm
1461 CALL ibcoff(ibcs,icodt(nod))
1480 1 NB ,LEN ,BUFR ,IXLINS ,A ,
1481 2 STIFN ,VISCN ,IBC ,ISECIN ,NOINT ,
1482 3 IBAG ,ICODT ,SECFCUM,NSTRF ,ICONTACT,
1483 4 FCONT ,DAANC6 ,NLG ,ALPHAK ,H3D_DATA)
1492#include "implicit_f.inc"
1496#include "scr05_c.inc"
1497#include "scr07_c.inc"
1498#include "scr14_c.inc"
1499#include "scr16_c.inc"
1500#include "scr18_c.inc"
1501#include "com01_c.inc"
1502#include "com04_c.inc"
1503#include "com06_c.inc"
1504#include "com08_c.inc"
1508 TYPE(output_),
intent(inout) :: output
1509 INTEGER nb, len, ibc ,isecin ,ibag , noint,
1510 . ixlins(2,*), ICODT(*), NSTRF(*),
1511 . ICONTACT(*), NLG(*)
1513 . bufr(len,*), a(3,*), stifn(*), viscn(*),
1514 . secfcum(7,numnod,nsect),
1515 . fcont(3,*), alphak(3,*)
1516 DOUBLE PRECISION DAANC6(3,6,*)
1517 TYPE(H3D_DATABASE) :: H3D_DATA
1521 INTEGER I, J, II, N, N1, N2, N1G, N2G, K0, K1S, NBINTER, INC,ISIGN
1522 DOUBLE PRECISION DAANC6L(3,6)
1534 a(1,n1g) = a(1,n1g) + bufr(2,i)
1535 a(2,n1g) = a(2,n1g) + bufr(3,i)
1536 a(3,n1g) = a(3,n1g) + bufr(4,i)
1537 stifn(n1g) = stifn(n1g) + bufr(5,i)
1539 a(1,n2g) = a(1,n2g) + bufr(6,i)
1540 a(2,n2g) = a(2,n2g) + bufr(7,i)
1541 a(3,n2g) = a(3,n2g) + bufr(8,i)
1542 stifn(n2g) = stifn(n2g) + bufr(9,i)
1545 IF(alphak(2,n1) < zero .or. bufr(10,i) < zero)isign = -1
1546 alphak(2,n1) = isign*
min(abs(alphak(2,n1)),abs(bufr(10,i)))
1548 IF(alphak(2,n2) < zero .or. bufr(11,i) < zero)isign = -1
1549 alphak(2,n2) = isign*
min(abs(alphak(2,n2)),abs(bufr(11,i)))
1550 CALL getdpdaanc(daanc6l,bufr(12,i),iresp,inc)
1551 daanc6(1,1,n1) = daanc6(1,1,n1)+daanc6l(1,1)
1552 daanc6(1,2,n1) = daanc6(1,2,n1)+daanc6l(1,2)
1553 daanc6(1,3,n1) = daanc6(1,3,n1)+daanc6l(1,3)
1554 daanc6(1,4,n1) = daanc6(1,4,n1)+daanc6l(1,4)
1555 daanc6(1,5,n1) = daanc6(1,5,n1)+daanc6l(1,5)
1556 daanc6(1,6,n1) = daanc6(1,6,n1)+daanc6l(1,6)
1557 daanc6(2,1,n1) = daanc6(2,1,n1)+daanc6l(2,1)
1558 daanc6(2,2,n1) = daanc6(2,2,n1)+daanc6l(2,2)
1559 daanc6(2,3,n1) = daanc6(2,3,n1)+daanc6l(2,3)
1560 daanc6(2,4,n1) = daanc6(2,4,n1)+daanc6l(2,4)
1561 daanc6(2,5,n1) = daanc6(2,5,n1)+daanc6l(2,5)
1562 daanc6(2,6,n1) = daanc6(2,6,n1)+daanc6l(2,6)
1563 daanc6(3,1,n1) = daanc6(3,1,n1)+daanc6l(3,1)
1564 daanc6(3,2,n1) = daanc6(3,2,n1)+daanc6l(3,2)
1565 daanc6(3,3,n1) = daanc6(3,3,n1)+daanc6l(3,3)
1566 daanc6(3,4,n1) = daanc6(3,4,n1)+daanc6l(3,4)
1567 daanc6(3,5,n1) = daanc6(3,5,n1)+daanc6l(3,5)
1568 daanc6(3,6,n1) = daanc6(3,6,n1)+daanc6l(3,6)
1569 CALL getdpdaanc(daanc6l,bufr(12+inc,i),iresp,inc)
1570 daanc6(1,1,n2) = daanc6(1,1,n2)+daanc6l(1,1)
1571 daanc6(1,2,n2) = daanc6(1,2,n2)+daanc6l(1,2)
1572 daanc6(1,3,n2) = daanc6(1,3,n2)+daanc6l(1,3)
1573 daanc6(1,4,n2) = daanc6(1,4,n2)+daanc6l(1,4)
1574 daanc6(1,5,n2) = daanc6(1,5,n2)+daanc6l(1,5)
1575 daanc6(1,6,n2) = daanc6(1,6,n2)+daanc6l(1,6)
1576 daanc6(2,1,n2) = daanc6(2,1,n2)+daanc6l(2,1)
1577 daanc6(2,2,n2) = daanc6(2,2,n2)+daanc6l(2,2)
1578 daanc6(2,3,n2) = daanc6(2,3,n2)+daanc6l(2,3)
1579 daanc6(2,4,n2) = daanc6(2,4,n2)+daanc6l(2,4)
1580 daanc6(2,5,n2) = daanc6(2,5,n2)+daanc6l(2,5)
1581 daanc6(2,6,n2) = daanc6(2,6,n2)+daanc6l(2,6)
1582 daanc6(3,1,n2) = daanc6(3,1,n2)+daanc6l(3,1)
1583 daanc6(3,2,n2) = daanc6(3,2,n2)+daanc6l(3,2)
1584 daanc6(3,3,n2) = daanc6(3,3,n2)+daanc6l(3,3)
1585 daanc6(3,4,n2) = daanc6(3,4,n2)+daanc6l(3,4)
1586 daanc6(3,5,n2) = daanc6(3,5,n2)+daanc6l(3,5)
1587 daanc6(3,6,n2) = daanc6(3,6,n2)+daanc6l(3,6)
1596 a(1,n1g) = a(1,n1g) + bufr(2,i)
1597 a(2,n1g) = a(2,n1g) + bufr(3,i)
1598 a(3,n1g) = a(3,n1g) + bufr(4,i)
1599 stifn(n1g) = stifn(n1g) + bufr(5,i)
1600 viscn(n1g) = viscn(n1g) + bufr(6,i)
1602 a(1,n2g) = a(1,n2g) + bufr(7,i)
1603 a(2,n2g) = a(2,n2g) + bufr(8,i)
1604 a(3,n2g) = a(3,n2g) + bufr(9,i)
1605 stifn(n2g) = stifn(n2g) + bufr(10,i)
1606 viscn(n2g) = viscn(n2g) + bufr(11,i)
1609 IF(alphak(2,n1) < zero .or. bufr(12,i) < zero)isign = -1
1610 alphak(2,n1) = isign*
min(abs(alphak(2,n1)),abs(bufr(12,i)))
1612 IF(alphak(2,n2) < zero .or. bufr(13,i) < zero)isign = -1
1613 alphak(2,n2) = isign*
min(abs(alphak(2,n2)),abs(bufr(13,i)))
1614 CALL getdpdaanc(daanc6l,bufr(14,i),iresp,inc)
1615 daanc6(1,1,n1) = daanc6(1,1,n1)+daanc6l(1,1)
1616 daanc6(1,2,n1) = daanc6(1,2,n1)+daanc6l(1,2)
1617 daanc6(1,3,n1) = daanc6(1,3,n1)+daanc6l(1,3)
1618 daanc6(1,4,n1) = daanc6(1,4,n1)+daanc6l(1,4)
1619 daanc6(1,5,n1) = daanc6(1,5,n1)+daanc6l(1,5)
1620 daanc6(1,6,n1) = daanc6(1,6,n1)+daanc6l(1,6)
1621 daanc6(2,1,n1) = daanc6(2,1,n1)+daanc6l(2,1)
1622 daanc6(2,2,n1) = daanc6(2,2,n1)+daanc6l(2,2)
1623 daanc6(2,3,n1) = daanc6(2,3,n1)+daanc6l(2,3)
1624 daanc6(2,4,n1) = daanc6(2,4,n1)+daanc6l(2,4)
1625 daanc6(2,5,n1) = daanc6(2,5,n1)+daanc6l(2,5)
1626 daanc6(2,6,n1) = daanc6(2,6,n1)+daanc6l(2,6)
1627 daanc6(3,1,n1) = daanc6(3,1,n1)+daanc6l(3,1)
1628 daanc6(3,2,n1) = daanc6(3,2,n1)+daanc6l(3,2)
1629 daanc6(3,3,n1) = daanc6(3,3,n1)+daanc6l(3,3)
1630 daanc6(3,4,n1) = daanc6(3,4,n1)+daanc6l(3,4)
1631 daanc6(3,5,n1) = daanc6(3,5,n1)+daanc6l(3,5)
1632 daanc6(3,6,n1) = daanc6(3,6,n1)+daanc6l(3,6)
1633 CALL getdpdaanc(daanc6l,bufr(14+inc,i),iresp,inc)
1634 daanc6(1,1,n2) = daanc6(1,1,n2)+daanc6l(1,1)
1635 daanc6(1,2,n2) = daanc6(1,2,n2)+daanc6l(1,2)
1636 daanc6(1,3,n2) = daanc6(1,3,n2)+daanc6l(1,3)
1637 daanc6(1,4,n2) = daanc6(1,4,n2)+daanc6l(1,4)
1638 daanc6(1,5,n2) = daanc6(1,5,n2)+daanc6l(1,5)
1639 daanc6(1,6,n2) = daanc6(1,6,n2)+daanc6l(1,6)
1640 daanc6(2,1,n2) = daanc6(2,1,n2)+daanc6l(2,1)
1641 daanc6(2,2,n2) = daanc6(2,2,n2)+daanc6l(2,2)
1642 daanc6(2,3,n2) = daanc6(2,3,n2)+daanc6l(2,3)
1643 daanc6(2,4,n2) = daanc6(2,4,n2)+daanc6l(2,4)
1644 daanc6(2,5,n2) = daanc6(2,5,n2)+daanc6l(2,5)
1645 daanc6(2,6,n2) = daanc6(2,6,n2)+daanc6l(2,6)
1646 daanc6(3,1,n2) = daanc6(3,1,n2)+daanc6l(3,1)
1647 daanc6(3,2,n2) = daanc6(3,2,n2)+daanc6l(3,2)
1648 daanc6(3,3,n2) = daanc6(3,3,n2)+daanc6l(3,3)
1649 daanc6(3,4,n2) = daanc6(3,4,n2)+daanc6l(3,4)
1650 daanc6(3,5,n2) = daanc6(3,5,n2)+daanc6l(3,5)
1651 daanc6(3,6,n2) = daanc6(3,6,n2)+daanc6l(3,6)
1658 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
1659 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP
1660 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
1661 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
1670 fcont(1,n1g)=fcont(1,n1g)+bufr(2,i)
1671 fcont(2,n1g)=fcont(2,n1g)+bufr(3,i)
1672 fcont(3,n1g)=fcont(3,n1g)+bufr(4,i)
1674 fcont(1,n2g)=fcont(1,n2g)+bufr(5,i)
1675 fcont(2,n2g)=fcont(2,n2g)+bufr(6,i)
1676 fcont(3,n2g)=fcont(3,n2g)+bufr(7,i)
1683 IF(nstrf(1)+nstrf(2)/=0)
THEN
1685 nbinter=nstrf(k0+14)
1688 IF(nstrf(k1s)==noint)
THEN
1692 n = nint(bufr(1,ii))
1697 IF(secfcum(4,n1g,i)==1.)
THEN
1698 secfcum(1,n1g,i)=secfcum(1,n1g,i)+bufr(2,ii)
1699 secfcum(2,n1g,i)=secfcum(2,n1g,i)+bufr(3,ii)
1700 secfcum(3,n1g,i)=secfcum(3,n1g,i)+bufr(4,ii)
1702 IF(secfcum(4,n2g,i)==1.)
THEN
1703 secfcum(1,n2g,i)=secfcum(1,n2g,i)+bufr(6,ii)
1704 secfcum(2,n2g,i)=secfcum(2,n2g,i)+bufr(7,ii)
1705 secfcum(3,n2g,i)=secfcum(3,n2g,i)+bufr(8,ii)
1710 n = nint(bufr(1,ii))
1715 IF(secfcum(4,n1g,i)==1.)
THEN
1716 secfcum(1,n1g,i)=secfcum(1,n1g,i)+bufr(2,ii)
1717 secfcum(2,n1g,i)=secfcum(2,n1g,i)+bufr(3,ii)
1718 secfcum(3,n1g,i)=secfcum(3,n1g,i)+bufr(4,ii)
1720 IF(secfcum(4,n2g,i)==1.)
THEN
1721 secfcum(1,n2g,i)=secfcum(1,n2g,i)+bufr(7,ii)
1722 secfcum(2,n2g,i)=secfcum(2,n2g,i)+bufr(8,ii)
1723 secfcum(3,n2g,i)=secfcum(3,n2g,i)+bufr(9,ii)
1740 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
1741 + bufr(4,i)/=zero)
THEN
1747 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
1748 + bufr(8,i)/=zero)
THEN
1757 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
1758 + bufr(4,i)/=zero)
THEN
1764 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
1765 + bufr(9,i)/=zero)
THEN
2201 1 NB ,LEN ,BUFR ,IRECTS ,FSKYI,
2202 2 ISKY ,IBC ,ISECIN,NOINT ,IBAG ,
2203 3 ICODT ,SECFCUM,NSTRF ,ICONTACT,FCONT,
2204 4 INTTH ,FTHESKYI,CONDNSKYI,H3D_DATA,
2205 5 NIN ,TAGNCONT,KLOADPINTER,LOADPINTER,
2206 7 LOADP_HYD_INTER,NODADT_THERM)
2216#include "implicit_f.inc"
2220#include "parit_c.inc"
2221#include "scr07_c.inc"
2222#include "scr14_c.inc"
2223#include "scr16_c.inc"
2224#include "scr18_c.inc"
2225#include "com01_c.inc"
2226#include "com04_c.inc"
2227#include "com06_c.inc"
2228#include "com08_c.inc"
2232 TYPE(output_),
intent(inout) :: OUTPUT
2233 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT,INTTH,NIN,
2234 . irects(2,*), isky(*), icodt(*), nstrf(*),
2236 . tagncont(nloadp_hyd_inter,*),kloadpinter(*),loadpinter(*),
2237 . loadp_hyd_inter(*)
2238 INTEGER,
INTENT(IN) :: NODADT_THERM
2241 . FSKYI(LSKYI,NFSKYI), SECFCUM(7,NUMNOD,NSECT),
2242 . FCONT(3,*),FTHESKYI(*),CONDNSKYI(*)
2243 TYPE(H3D_DATABASE) :: H3D_DATA
2247 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER,PP,PPL,NOD1,NOD2,NP
2251 IF ((nisky+nb)> lskyi)
THEN
2252 CALL ancmsg(msgid=26,anmode=aninfo)
2256 IF(intth == 0 )
THEN
2262 fskyi(nisky,1)=bufr(2,i)
2263 fskyi(nisky,2)=bufr(3,i)
2264 fskyi(nisky,3)=bufr(4,i)
2265 fskyi(nisky,4)=bufr(5,i)
2269 fskyi(nisky,1)=bufr(6,i)
2270 fskyi(nisky,2)=bufr(7,i)
2271 fskyi(nisky,3)=bufr(8,i)
2272 fskyi(nisky,4)=bufr(9,i)
2280 fskyi(nisky,1)=bufr(2,i)
2281 fskyi(nisky,2)=bufr(3,i)
2282 fskyi(nisky,3)=bufr(4,i)
2283 fskyi(nisky,4)=bufr(5,i)
2284 fskyi(nisky,5)=bufr(6,i)
2288 fskyi(nisky,1)=bufr(7,i)
2289 fskyi(nisky,2)=bufr(8,i)
2290 fskyi(nisky,3)=bufr(9,i)
2291 fskyi(nisky,4)=bufr(10,i)
2292 fskyi(nisky,5)=bufr(11,i)
2297 IF(nodadt_therm == 1)
THEN
2303 fskyi(nisky,1)=bufr(2,i)
2304 fskyi(nisky,2)=bufr(3,i)
2305 fskyi(nisky,3)=bufr(4,i)
2306 fskyi(nisky,4)=bufr(5,i)
2307 ftheskyi(nisky) =bufr(10,i)
2308 condnskyi(nisky)=bufr(12,i)
2312 fskyi(nisky,1)=bufr(6,i)
2313 fskyi(nisky,2)=bufr(7,i)
2314 fskyi(nisky,3)=bufr(8,i)
2315 fskyi(nisky,4)=bufr(9,i)
2316 ftheskyi(nisky) =bufr(11,i)
2317 condnskyi(nisky)=bufr(13,i)
2325 fskyi(nisky,1)=bufr(2,i)
2326 fskyi(nisky,2)=bufr(3,i)
2327 fskyi(nisky,3)=bufr(4,i)
2328 fskyi(nisky,4)=bufr(5,i)
2329 fskyi(nisky,5)=bufr(6,i)
2330 ftheskyi(nisky) =bufr(12,i)
2331 condnskyi(nisky)=bufr(13,i)
2335 fskyi(nisky,1)=bufr(7,i)
2336 fskyi(nisky,2)=bufr(8,i)
2337 fskyi(nisky,3)=bufr(9,i)
2338 fskyi(nisky,4)=bufr(10,i)
2339 fskyi(nisky,5)=bufr(11,i)
2340 ftheskyi(nisky) =bufr(13,i)
2341 condnskyi(nisky)=bufr(14,i)
2351 fskyi(nisky,1)=bufr(2,i)
2352 fskyi(nisky,2)=bufr(3,i)
2353 fskyi(nisky,3)=bufr(4,i)
2354 fskyi(nisky,4)=bufr(5,i)
2355 ftheskyi(nisky) =bufr(10,i)
2359 fskyi(nisky,1)=bufr(6,i)
2360 fskyi(nisky,2)=bufr(7,i)
2361 fskyi(nisky,3)=bufr(8,i)
2362 fskyi(nisky,4)=bufr(9,i)
2363 ftheskyi(nisky) =bufr(11,i)
2371 fskyi(nisky,1)=bufr(2,i)
2372 fskyi(nisky,2)=bufr(3,i)
2373 fskyi(nisky,3)=bufr(4,i)
2374 fskyi(nisky,4)=bufr(5,i)
2375 fskyi(nisky,5)=bufr(6,i)
2376 ftheskyi(nisky) =bufr(12,i)
2380 fskyi(nisky,1)=bufr(7,i)
2381 fskyi(nisky,2)=bufr(8,i)
2382 fskyi(nisky,3)=bufr(9,i)
2383 fskyi(nisky,4)=bufr(10,i)
2384 fskyi(nisky,5)=bufr(11,i)
2385 ftheskyi(nisky) =bufr(13,i)
2394 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
2395 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
2396 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
2397 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
2403 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2404 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2405 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2407 fcont(1,nod)=fcont(1,nod)+ bufr(6,i)
2408 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
2409 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
2415 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2416 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2417 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2419 fcont(1,nod)=fcont(1,nod)+ bufr(7,i)
2420 fcont(2,nod)=fcont(2,nod)+ bufr(8,i)
2421 fcont(3,nod)=fcont(3,nod)+ bufr(9,i)
2427 IF(nintloadp > 0)
THEN
2432 DO np = kloadpinter(nin)+1, kloadpinter(nin+1)
2434 ppl = loadp_hyd_inter(pp)
2435 tagncont(ppl,nod1) = 1
2436 tagncont(ppl,nod2) = 1
2446 nbinter=nstrf(k0+14)
2449 IF(nstrf(k1s)==noint)
THEN
2453 n = nint(bufr(1,ii))
2455 IF(secfcum(4,nod,i)==1.)
THEN
2456 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2458 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2461 IF(secfcum(4,nod,i)==1.)
THEN
2462 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(6,ii)
2463 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(7,ii)
2464 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(8,ii)
2469 n = nint(bufr(1,ii))
2471 IF(secfcum(4,nod,i)==1.)
THEN
2472 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2473 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
2474 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2477 IF(secfcum(4,nod,i)==1.)
THEN
2478 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(7,ii)
2479 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(8,ii)
2480 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(9,ii)
2493 IF((ibag/=0).OR.(idamp_rdof/=0))
THEN
2497 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2498 + bufr(4,i)/=zero)
THEN
2503 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
2504 + bufr(8,i)/=zero)
THEN
2511 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2512 + bufr(4,i)/=zero)
THEN
2517 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
2518 + bufr(9,i)/=zero)
THEN
2542 1 NB ,LEN ,BUFR ,IRECTS ,FSKYI,
2543 2 ISKY ,IBC ,ISECIN,NOINT ,IBAG ,
2544 3 ICODT ,SECFCUM,NSTRF ,ICONTACT,FCONT,
2555#include "implicit_f.inc"
2559#include "parit_c.inc"
2560#include "scr07_c.inc"
2561#include "scr14_c.inc"
2562#include "scr16_c.inc"
2563#include "scr18_c.inc"
2564#include "com01_c.inc"
2565#include "com04_c.inc"
2566#include "com06_c.inc"
2567#include "com08_c.inc"
2571 TYPE(output_),
INTENT(INOUT) :: OUTPUT
2572 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT,
2573 . IRECTS(2,*), ISKY(*), ICODT(*), NSTRF(*),
2574 . icontact(*),nlg(*)
2577 . fskyi(lskyi,nfskyi), secfcum(7,numnod,nsect),
2583 INTEGER I, J, , N, NOD, K0, K1S, IBCS, IBCM, NBINTER
2587 IF ((NISKY+NB)> LSKYI)THEN
2588 CALL (MSGID=26,anmode=aninfo)
2595 nod = nlg(irects(1,n))
2597 fskyi(nisky,1)=bufr(2,i)
2598 fskyi(nisky,2)=bufr(3,i)
2599 fskyi(nisky,3)=bufr(4,i)
2600 fskyi(nisky,4)=bufr(5,i)
2602 nod = nlg(irects(2,n))
2604 fskyi(nisky,1)=bufr(6,i)
2605 fskyi(nisky,2)=bufr(7,i)
2606 fskyi(nisky,3)=bufr(8,i)
2607 fskyi(nisky,4)=bufr(9,i)
2613 nod = nlg(irects(1,n))
2615 fskyi(nisky,1)=bufr(2,i)
2616 fskyi(nisky,2)=bufr(3,i)
2617 fskyi(nisky,3)=bufr(4,i)
2618 fskyi(nisky,4)=bufr(5,i)
2619 fskyi(nisky,5)=bufr(6,i)
2621 nod = nlg(irects(2,n))
2623 fskyi(nisky,1)=bufr(7,i)
2624 fskyi(nisky,2)=bufr(8,i)
2625 fskyi(nisky,3)=bufr(9,i)
2626 fskyi(nisky,4)=bufr(10,i)
2627 fskyi(nisky,5)=bufr(11,i)
2634 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
2635 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
2636 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
2637 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
2642 nod = nlg(irects(1,n))
2643 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2644 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2645 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2646 nod = nlg(irects(2,n))
2647 fcont(1,nod)=fcont(1,nod)+ bufr(6,i)
2648 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
2649 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
2654 nod = nlg(irects(1,n))
2655 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2656 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2657 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2658 nod = nlg(irects(2,n))
2659 fcont(1,nod)=fcont(1,nod)+ bufr(7,i)
2660 fcont(2,nod)=fcont(2,nod)+ bufr(8,i)
2661 fcont(3,nod)=fcont(3,nod)+ bufr(9,i)
2669 IF(nstrf(1)+nstrf(2)/=0)
THEN
2671 nbinter=nstrf(k0+14)
2674 IF(nstrf(k1s)==noint)
THEN
2678 n = nint(bufr(1,ii))
2679 nod = nlg(irects(1,n))
2680 IF(secfcum(4,nod,i)==1.)
THEN
2681 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2682 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
2683 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2685 nod = nlg(irects(2,n))
2686 IF(secfcum(4,nod,i)==1.)
THEN
2687 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(6,ii)
2688 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(7,ii)
2689 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(8,ii)
2694 n = nint(bufr(1,ii))
2695 nod = nlg(irects(1,n))
2696 IF(secfcum(4,nod,i)==1.)
THEN
2697 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2698 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
2699 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2701 nod = nlg(irects(2,n))
2702 IF(secfcum(4,nod,i)==1.)
THEN
2703 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(7,ii)
2704 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(8,ii)
2705 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(9,ii)
2722 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2723 + bufr(4,i)/=zero)
THEN
2725 nod = nlg(irects(1,n))
2728 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
2729 + bufr(8,i)/=zero)
THEN
2730 nod = nlg(irects(2,n))
2736 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2737 + bufr(4,i)/=zero)
THEN
2739 nod = nlg(irects(1,n))
2742 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
2743 + bufr(9,i)/=zero)
THEN
2744 nod = nlg(irects(2,n))
2872 1 IPARI ,ISENDTO ,IRCVFROM,INTLIST ,NBINTC ,
2873 2 ISIZXV ,ILENXV ,IAD_ELEM,I2SIZE ,ITASK ,
2874 3 ISLEN7 ,IRLEN7 ,ISLEN11 ,IRLEN11 ,IGRBRIC ,
2875 4 NME17 ,ISLEN17 ,IRLEN17 ,IRLEN7T ,ISLEN7T ,
2876 5 LINDIDEL,LBUFIDEL,IRLEN20 ,ISLEN20 ,IRLEN20T,
2877 6 ISLEN20T,NBINT20 ,IRLEN20E,ISLEN20E,FR_RBY ,
2878 7 FR_RBY6 ,NPBY ,IRBKIN_L,NRBYKIN_L,KINDRBY,
2879 8 NSENSOR ,SENSOR_TAB,LBUFIDEL24,INTBUF_TAB,
2880 9 SORT_COMM,NEED_COMM_INT25_SOLID_EROSION,COMM_INT25_SOLID_EROSION )
2896#include "implicit_f.inc"
2900#include "com01_c.inc"
2901#include "com04_c.inc"
2902#include "com08_c.inc"
2903#include "scr17_c.inc"
2904#include "scr18_c.inc"
2905#include "task_c.inc"
2906#include "param_c.inc"
2907#include "units_c.inc"
2908#include "warn_c.inc"
2909#include "tabsiz_c.inc"
2911 COMMON /ring/irecvf,isendt,iring
2912 INTEGER IRECVF,ISENDT,IRING
2916 INTEGER ,
INTENT(IN) :: NSENSOR
2917 INTEGER ISENDTO(NINTER+1,NSPMD+1),
2918 . IRCVFROM(NINTER+1,NSPMD+1), IPARI(NPARI,*),
2919 . INTLIST(NINTER), IAD_ELEM(2,*),
2920 . NBINTC, ISIZXV, ILENXV, I2SIZE,ITASK,
2921 . ISLEN7, IRLEN7, ISLEN11, IRLEN11, NME17, ISLEN17, IRLEN17,
2922 . IRLEN7T, ISLEN7T, LINDIDEL, LBUFIDEL,
2923 . IRLEN20,ISLEN20,IRLEN20T,ISLEN20T,NBINT20,
2924 . IRLEN20E, ISLEN20E, NRBYKIN_L, IRBKIN_L(*),KINDRBY(*),
2925 . FR_RBY(*), FR_RBY6(*), NPBY(NNPBY,*),LBUFIDEL24
2926 TYPE(INTBUF_STRUCT_) INTBUF_TAB(*)
2927 TYPE(sorting_comm_type),
DIMENSION(NINTER),
INTENT(inout) :: SORT_COMM
2928 TYPE (SENSOR_STR_) ,
DIMENSION(NSENSOR) :: SENSOR_TAB
2929 LOGICAL,
DIMENSION(NSPMD),
INTENT(inout) :: NEED_COMM_INT25_SOLID_EROSION
2930 INTEGER,
INTENT(inout) :: COMM_INT25_SOLID_EROSION
2932 TYPE (GROUP_) ,
DIMENSION(NGRBRIC) :: IGRBRIC
2936 INTEGER I, J, K, NSN, NMN, ITYP, IERROR, IERROR1, KK, NIN, P,
2937 . scpmax1,isd,irc,my_rank,last, ige, ign, nme, nmes,intth,
2938 . nrtm, nrts, nlinm, nlins,nsne, nmne,n,m,
2939 . work(70000),ity,noint,inacti,isens,interact
2940 INTEGER,
DIMENSION(:),
ALLOCATABLE:: INDX
2941 INTEGER,
DIMENSION(:),
ALLOCATABLE:: PPAR
2942 INTEGER,
DIMENSION(:),
ALLOCATABLE:: INDP
2943 INTEGER,
DIMENSION(:),
ALLOCATABLE:: NNP
2944 INTEGER,
DIMENSION(:),
ALLOCATABLE:: RANK
2945 INTEGER,
DIMENSION(:),
ALLOCATABLE:: IRBODY
2946 INTEGER,
DIMENSION(:),
ALLOCATABLE:: INDB
2947 INTEGER,
DIMENSION(:),
ALLOCATABLE:: NNSN
2948 my_real STARTT, STOPT,TS
2952 CALL my_alloc(indx,2*ninter)
2953 CALL my_alloc(ppar,ninter)
2954 CALL my_alloc(indp,2*nspmd)
2955 CALL my_alloc(nnp,nspmd)
2956 CALL my_alloc(rank,nspmd)
2957 CALL my_alloc(irbody,numnod)
2958 CALL my_alloc(indb,2*nrbykin)
2959 CALL my_alloc(nnsn,nrbykin)
2985 IF(ity==7.OR.ity==10.OR.ity==11.OR.
2986 . ity==17.OR.ity==20.OR.
2987 . ity==22.OR.ity==23.OR.ity==24.OR.ity==25)
THEN
2989 IF (inacti==5.OR.inacti==6.OR.inacti==7.OR.ity==23.OR.inacti==-1)
THEN
2993 IF(ity == 7.OR.ity == 11.OR.ity == 24.OR.ity == 25)
2994 . isens = ipari(64,i)
2996 ts = sensor_tab(isens)%TSTART
2997 IF (ts>tt) interact = 0
2999 startt = intbuf_tab(i)%VARIABLES(3)
3000 stopt = intbuf_tab(i)%VARIABLES(11)
3001 IF (startt>tt) interact = 0
3004 IF (interact==0)
THEN
3008 CALL ancmsg(msgid=220,anmode=aninfo,
3009 * i1=noint,i2=inacti,r1=startt,r2=startt)
3012 IF (interact==0 .AND. ity /=25)
THEN
3016 nsnfi(i)%P(1:nspmd)=0
3018 nsnsi(i)%P(1:nspmd)=0
3031 IF(ipari(17,i)>0)
THEN
3036 IF(ityp==7.OR.ityp==22.OR.
3037 . ityp==10.OR.ityp==5.OR.
3038 . ityp==23.OR.ityp==24.OR.
3040 lindidel = lindidel+nrtm
3041 lbufidel = lbufidel+4*nrtm+4
3042 IF(ityp==24.OR.ityp==25) lbufidel24 = lbufidel24 + nrtm
3043 ELSEIF(ityp==11)
THEN
3044 lindidel = lindidel+nrtm+nrts
3045 lbufidel = lbufidel+2*(nrtm+nrts)+4
3047 lindidel = lindidel+nrtm+nrts
3048 lbufidel = lbufidel+4*(nrtm+nrts)+4
3050 lindidel = lindidel+nsn
3051 lbufidel = lbufidel+4*nsn+4
3052 ELSEIF(ityp==20)
THEN
3053 lindidel = lindidel+nrtm
3054 lbufidel = lbufidel+4*nrtm+4
3058 lindidel = lindidel+nlinm+nlins
3059 lbufidel = lbufidel+2*(nlinm+nlins)+4
3069 IF (ityp==7.OR.ityp==10.OR.ityp==22.OR.
3070 . ityp==11.OR.ityp==23.OR.ityp==24.OR.
3075 isendto(i,ispmd+1)=nsn
3076 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3077 isendto(ninter+1,ispmd+1)=isendto(ninter+1,ispmd+1)+nsn
3080 ircvfrom(i,ispmd+1) = nmn
3081 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3082 ircvfrom(ninter+1,ispmd+1)=ircvfrom(ninter+1,ispmd+1)
3087 ELSEIF(ityp==17)
THEN
3088 IF(ipari(33,i)==0)
THEN
3091 nmes= igrbric(ign)%NENTITY
3092 nme = igrbric(ige)%NENTITY
3093 nme17 = nme17+nme+nmes
3095 isendto(i,ispmd+1)=nmes
3096 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3097 isendto(ninter+1,ispmd+1)=isendto(ninter+1,ispmd+1)
3101 ircvfrom(i,ispmd+1) = nme
3102 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3103 ircvfrom(ninter+1,ispmd+1)=ircvfrom(ninter+1,ispmd+1)
3107 ELSEIF(ityp==20)
THEN
3113 IF(nsn+nsne/=0)
THEN
3114 isendto(i,ispmd+1)=nsn+nsne
3115 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3116 isendto(ninter+1,ispmd+1)=isendto(ninter+1,ispmd+1)
3119 IF(nmn+nmne/=0)
THEN
3120 ircvfrom(i,ispmd+1) = nmn+nmne
3121 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3122 ircvfrom(ninter+1,ispmd+1)=ircvfrom(ninter+1,ispmd+1)
3130 CALL spmd_ibcast(isendto(1,k),isendto(1,k),ninter+1,1,
3132 CALL spmd_ibcast(ircvfrom(1,k),ircvfrom(1,k),ninter+1,1,
3137 isendto(i,nspmd+1)=0
3138 ircvfrom(i,nspmd+1)=0
3140 IF(ircvfrom(i,k)/=0)
THEN
3141 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3143 IF(isendto(i,k)/=0)
THEN
3144 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3154 IF (ityp==7.OR.ityp==10.OR.ityp==11.OR.
3155 . (ityp==17.AND.ipari(33,i) == 0).OR.ityp==20.OR.
3156 . ityp==22.OR.ityp==23.OR.ityp==24.OR.ityp==25)
THEN
3161 scpmax1=
max(scpmax1,ircvfrom(i,j))
3163 ppar(nbintc)=scpmax1
3167 CALL my_orders(0,work,ppar,indx,nbintc,1)
3169 indx(i+nbintc)=intlist(i)
3172 intlist(i)=indx(indx(i)+nbintc)
3175 IF(debug(3)>=1.AND.ispmd==0.AND.nspmd>1)
THEN
3176 WRITE(istdo,*)'** interfaces nodal decomposition
'
3177 WRITE(ISTDO,*)'#PROC NSN TOT NMN TOT TOTAL'
3178 WRITE(iout,*)
'** INTERFACES NODAL DECOMPOSITION '
3179 WRITE(iout,*)
'#PROC NSN TOT NMN TOT TOTAL'
3184 isd = isd + isendto(i,j)
3185 irc = irc + ircvfrom(i,j)
3187 WRITE(istdo,
'(I4,3X,I8,2X,I8,2X,I8)')j,isd,irc,isd+irc
3188 WRITE(iout,
'(I4,3X,I8,2X,I8,2X,I8)')j,isd,irc,isd+irc
3195 nnp(k) = ircvfrom(ninter+1,k) + isendto(ninter+1,k)
3198 IF(nspmd > 1)
CALL my_orders(0,work,nnp,indp,nspmd,1)
3205 IF(nnp(ispmd+1)>0)
THEN
3206 my_rank = rank(ispmd+1)
3208 last = indp(my_rank-1)
3213 IF(my_rank==nspmd)
THEN
3214 IF(irecvf/=0) isendt = -1
3216 isendt = indp(my_rank+1)
3219 IF(nnp(indp(nspmd-1))>0)iring= indp(nspmd)
3227 isizxv = isizxv + iad_elem(1,i+1)-iad_elem(2,i)
3230 ilenxv = 10 + 3*iroddl
3232 IF(iale/=0.AND.ninter>0) ilenxv = ilenxv + 6
3234 IF(idtmin(11)==3.OR.idtmin(10)==3) ilenxv = ilenxv + 2
3269 intth = ipari(47,nin)
3271 IF(ityp==7 .OR.ityp==10.OR.
3272 . ityp==23.OR.ityp==22.OR.ityp==24.OR.
3276 islen7 = islen7 +
nsnsi(nin)%P(p)
3277 irlen7 = irlen7 +
nsnfi(nin)%P(p)
3281 islen7t = islen7t +
nsnsi(nin)%P(p)
3282 irlen7t = irlen7t +
nsnfi(nin)%P(p)
3286 IF(ityp == 25 .AND. ipari(58,nin) >0)
THEN
3291 ELSEIF(ityp==11)
THEN
3294 islen11 = islen11 +
nsnsi(nin)%P(p)
3295 irlen11 = irlen11 +
nsnfi(nin)%P(p)
3297 ELSEIF(ityp==17.AND.ipari(33,nin)==0)
THEN
3300 islen17 = islen17 +
nsnsi(nin)%P(p)
3301 irlen17 = irlen17 +
nsnfi(nin)%P(p)
3303 ELSEIF(ityp==20)
THEN
3305 nbint20 = nbint20 + 1
3308 islen20 = islen20 +
nsnsi(nin)%P(p)
3309 irlen20 = irlen20 +
nsnfi(nin)%P(p)
3310 islen20e= islen20e+
nsnsie(nin)%P(p)
3311 irlen20e= irlen20e+
nsnfie(nin)%P(p)
3315 islen20t = islen20t +
nsnsi(nin)%P(p)
3316 irlen20t = irlen20t +
nsnfi(nin)%P(p)
3317 islen20e = islen20e +
nsnsie(nin)%P(p)
3318 irlen20e = irlen20e +
nsnfie(nin)%P(p)
3324 IF(nrbykin > 0)
THEN
3338 nrbykin_l=nrbykin_l+1
3339 nnsn(nrbykin_l)=-npby(2,n)
3340 irbkin_l(nrbykin_l)=n
3349 fr_rby6(n)=irbody(fr_rby(n))
3352 CALL my_orders(0,work,nnsn,indb,nrbykin_l,1)
3355 indb(nrbykin_l+n)=irbkin_l(n)
3358 irbkin_l(n)=indb(nrbykin_l+indb(n))
3364 indx(i+nbintc)=intlist(i)
3367 intlist(i)=indx(indx(i)+nbintc)
3373 . ircvfrom(ninter+1,ispmd+1)+isendto(ninter+1,ispmd+1),comm_cont)
3375 IF(ninter/=0)
CALL spmd_split_comm_inter( intbuf_tab, nbintc,intlist,ipari,isendto,ircvfrom,sort_comm,
3376 . need_comm_int25_solid_erosion,comm_int25_solid_erosion )
3411 SUBROUTINE spmd_initfi(IPARI,IFLAG,H3D_DATA,PARAMETERS,IDT_THERM,INTHEAT)
3425#include "implicit_f.inc"
3426#include "i25edge_c.inc"
3430#include "com01_c.inc"
3431#include "com04_c.inc"
3432#include "scr14_c.inc"
3433#include "scr16_c.inc"
3434#include "scr18_c.inc"
3435#include "task_c.inc"
3436#include "param_c.inc"
3437#include "parit_c.inc"
3438#include "spmd_c.inc"
3443 INTEGER IPARI(NPARI,*), IFLAG
3444 INTEGER,
INTENT(IN) :: IDT_THERM
3445 INTEGER,
INTENT(IN) :: INTHEAT
3446 TYPE(H3D_DATABASE) :: H3D_DATA
3447 TYPE (PARAMETERS_) ,
INTENT(IN):: PARAMETERS
3451 INTEGER I, P, NSN, NMN, ITYP, IGAP, IERROR, IERROR1, LENS, LENR,
3452 . lskyfi, inacti, nbintc, leni, j,intth, k, l,iedge4 ,intfric ,
3453 . flagremn ,sizremnorfi, ivis2 ,intnitsche,itied,ipstif
3454 INTEGER :: LENR_EDGE,LENS_EDGE
3462 ALLOCATE(
nsvfi(ninter),stat=ierror1)
3463 ierror = ierror + ierror1
3464 ALLOCATE(
nsnfi(ninter),stat=ierror1)
3465 ierror = ierror + ierror1
3466 ALLOCATE(
edge_fi(ninter),stat=ierror1)
3467 ierror = ierror + ierror1
3468 ALLOCATE(
nsnfi_sav(ninter),stat=ierror1)
3469 ierror = ierror + ierror1
3470 ALLOCATE(
nsnsi_sav(ninter),stat=ierror1)
3471 ierror = ierror + ierror1
3472 ALLOCATE(
itafi(ninter),stat=ierror1)
3473 ierror = ierror + ierror1
3474 ALLOCATE(
pmainfi(ninter),stat=ierror1)
3475 ierror = ierror + ierror1
3476 ALLOCATE(
kinfi(ninter),stat=ierror1)
3477 ierror = ierror + ierror1
3478 ALLOCATE(
msfi(ninter),stat=ierror1)
3479 ierror = ierror + ierror1
3480 ALLOCATE(
stifi(ninter),stat=ierror1)
3481 ierror = ierror + ierror1
3482 ALLOCATE(
gapfi(ninter),stat=ierror1)
3483 ierror = ierror + ierror1
3484 ALLOCATE(
gap_lfi(ninter),stat=ierror1)
3485 ierror = ierror + ierror1
3486 ALLOCATE(
xfi(ninter),stat=ierror1)
3487 ierror = ierror + ierror1
3488 ALLOCATE(
vfi(ninter),stat=ierror1)
3489 ierror = ierror + ierror1
3490 ALLOCATE(
nsvsi(ninter),stat=ierror1)
3491 ierror = ierror + ierror1
3492 ALLOCATE(
nsnsi(ninter),stat=ierror1)
3493 ierror = ierror + ierror1
3494 ALLOCATE(
afi(ninter),stat=ierror1)
3495 ierror = ierror + ierror1
3496 ALLOCATE(
stnfi(ninter),stat=ierror1)
3497 ierror = ierror + ierror1
3498 ALLOCATE(
vscfi(ninter),stat=ierror1)
3499 ierror = ierror + ierror1
3500 ALLOCATE(
penfi(ninter),stat=ierror1)
3501 ierror = ierror + ierror1
3502 ALLOCATE(
fskyfi(ninter),stat=ierror1)
3503 ierror = ierror + ierror1
3504 ALLOCATE(
iskyfi(ninter),stat=ierror1)
3505 ierror = ierror + ierror1
3507 ALLOCATE(
nisubsfi(ninter),stat=ierror1)
3508 ierror = ierror + ierror1
3509 ALLOCATE(
lisubsfi(ninter),stat=ierror1)
3510 ierror = ierror + ierror1
3512 ierror = ierror + ierror1
3513 ALLOCATE(
addsubsfi(ninter),stat=ierror1)
3514 ierror = ierror + ierror1
3515 IF(parameters%INTCAREA > 0 )
THEN
3517 ierror = ierror + ierror1
3520 ALLOCATE(
nisubsfie(ninter),stat=ierror1)
3521 ierror = ierror + ierror1
3522 ALLOCATE(
lisubsfie(ninter),stat=ierror1)
3523 ierror = ierror + ierror1
3525 ierror = ierror + ierror1
3527 ierror = ierror + ierror1
3530 ALLOCATE(
fnconti(ninter),stat=ierror1)
3531 ierror = ierror + ierror1
3532 ALLOCATE(
ftconti(ninter),stat=ierror1)
3533 ierror = ierror + ierror1
3535 ALLOCATE (
nlskyfi(ninter),stat=ierror1)
3536 ierror = ierror + ierror1
3539 ALLOCATE (
eminxfi(ninter),stat=ierror1)
3540 ierror = ierror + ierror1
3541 ALLOCATE (
ksfi(ninter),stat=ierror1)
3542 ierror = ierror + ierror1
3543 ALLOCATE (
frotsfi(ninter),stat=ierror1)
3544 ierror = ierror + ierror1
3545 ALLOCATE (
stnfi17(ninter),stat=ierror1)
3546 ierror = ierror + ierror1
3547 ALLOCATE(
xfi17(ninter),stat=ierror1)
3548 ierror = ierror + ierror1
3549 ALLOCATE(
vfi17(ninter),stat=ierror1)
3550 ierror = ierror + ierror1
3551 ALLOCATE(
afi17(ninter),stat=ierror1)
3552 ierror = ierror + ierror1
3554 ALLOCATE(
nbinflfi(ninter),stat=ierror1)
3555 ierror = ierror + ierror1
3556 ALLOCATE(
daanc6fi(ninter),stat=ierror1)
3557 ierror = ierror + ierror1
3558 ALLOCATE(
dxancfi(ninter),stat=ierror1)
3559 ierror = ierror + ierror1
3560 ALLOCATE(
dvancfi(ninter),stat=ierror1)
3561 ierror = ierror + ierror1
3562 ALLOCATE(
penfia(ninter),stat=ierror1)
3563 ierror = ierror + ierror1
3564 ALLOCATE(
alphakfi(ninter),stat=ierror1)
3565 ierror = ierror + ierror1
3566 ALLOCATE(
daancfi(ninter),stat=ierror1)
3567 ierror = ierror + ierror1
3569 ierror = ierror + ierror1
3571 ALLOCATE(
nsvfie(ninter),stat=ierror1)
3572 ierror = ierror + ierror1
3573 ALLOCATE(
nsnfie(ninter),stat=ierror1)
3574 ierror = ierror + ierror1
3575 ALLOCATE(
itafie(ninter),stat=ierror1)
3576 ierror = ierror + ierror1
3577 ALLOCATE(
msfie(ninter),stat=ierror1)
3578 ierror = ierror + ierror1
3579 ALLOCATE(
stifie(ninter),stat=ierror1)
3580 ierror = ierror + ierror1
3581 ALLOCATE(
gapfie(ninter),stat=ierror1)
3582 ierror = ierror + ierror1
3583 ALLOCATE(
xfie(ninter),stat=ierror1)
3584 ierror = ierror + ierror1
3585 ALLOCATE(
vfie(ninter),stat=ierror1)
3586 ierror = ierror + ierror1
3587 ALLOCATE(
nsvsie(ninter),stat=ierror1)
3588 ierror = ierror + ierror1
3589 ALLOCATE(
nsnsie(ninter),stat=ierror1)
3590 ierror = ierror + ierror1
3591 ALLOCATE(
afie(ninter),stat=ierror1)
3592 ierror = ierror + ierror1
3593 ALLOCATE(
stnfie(ninter),stat=ierror1)
3594 ierror = ierror + ierror1
3595 ALLOCATE(
vscfie(ninter),stat=ierror1)
3596 ierror = ierror + ierror1
3597 ALLOCATE(
penfie(ninter),stat=ierror1)
3598 ierror = ierror + ierror1
3599 ALLOCATE(
fskyfie(ninter),stat=ierror1)
3600 ierror = ierror + ierror1
3601 ALLOCATE(
iskyfie(ninter),stat=ierror1)
3602 ierror = ierror + ierror1
3604 ALLOCATE (
nlskyfie(ninter),stat=ierror1)
3605 ierror = ierror + ierror1
3608 ALLOCATE(
daanc6fie(ninter),stat=ierror1)
3609 ierror = ierror + ierror1
3610 ALLOCATE(
dxancfie(ninter),stat=ierror1)
3611 ierror = ierror + ierror1
3613 ierror = ierror + ierror1
3614 ALLOCATE(
penfiae(ninter),stat=ierror1)
3615 ierror = ierror + ierror1
3616 ALLOCATE(
alphakfie(ninter),stat=ierror1)
3617 ierror = ierror + ierror1
3618 ALLOCATE(
daancfie(ninter),stat=ierror1)
3619 ierror = ierror + ierror1
3621 ierror = ierror + ierror1
3623 ALLOCATE (
nodnxfie(ninter),stat=ierror1)
3624 ierror = ierror + ierror1
3625 ALLOCATE (
nodamsfie(ninter),stat=ierror1)
3626 ierror = ierror + ierror1
3628 ierror = ierror + ierror1
3632 ALLOCATE (
mtfi_pene(ninter),stat=ierror1)
3633 ierror = ierror + ierror1
3635 ierror = ierror + ierror1
3636 ALLOCATE (
mtfi_v(ninter),stat=ierror1)
3637 ierror = ierror + ierror1
3638 ALLOCATE (
mtfi_a(ninter),stat=ierror1)
3639 ierror = ierror + ierror1
3640 ALLOCATE (
i18kafi(ninter),stat=ierror1)
3641 ierror = ierror + ierror1
3642 ALLOCATE (
mtfi_n(ninter),stat=ierror1)
3643 ierror = ierror + ierror1
3645 ALLOCATE (
time_sfi(ninter),stat=ierror1)
3646 ierror = ierror + ierror1
3647 ALLOCATE (
irtlm_fi(ninter),stat=ierror1)
3648 ierror = ierror + ierror1
3650 ierror = ierror + ierror1
3652 ierror = ierror + ierror1
3654 ierror = ierror + ierror1
3656 ierror = ierror + ierror1
3657 ALLOCATE (
isedge_fi(ninter),stat=ierror1)
3658 ierror = ierror + ierror1
3659 ALLOCATE (
irtse_fi(ninter),stat=ierror1)
3660 ierror = ierror + ierror1
3661 ALLOCATE (
is2pt_fi(ninter),stat=ierror1)
3662 ierror = ierror + ierror1
3663 ALLOCATE (
isegpt_fi(ninter),stat=ierror1)
3664 ierror = ierror + ierror1
3665 ALLOCATE (
is2se_fi(ninter),stat=ierror1)
3666 ierror = ierror + ierror1
3667 ALLOCATE (
ispt2_fi(ninter),stat=ierror1)
3668 ierror = ierror + ierror1
3669 ALLOCATE (
remnor_fi(ninter),stat=ierror1)
3670 ierror = ierror + ierror1
3672 ierror = ierror + ierror1
3674 ALLOCATE (
islide_fi(ninter),stat=ierror1)
3675 ierror = ierror + ierror1
3676 ALLOCATE (
icodt_fi(ninter),stat=ierror1)
3677 ierror = ierror + ierror1
3678 ALLOCATE (
iskew_fi(ninter),stat=ierror1)
3679 ierror = ierror + ierror1
3682 IF(ipari(97,i) > 0) ipstif = ipstif + ipari(97,i)
3684 IF(parameters%ISTIF_DT > 0 .OR. ipstif>0 )
THEN
3689 IF(intheat > 0 )
THEN
3690 ALLOCATE (
ftheskyfi(ninter),stat=ierror1)
3691 ierror = ierror + ierror1
3692 ALLOCATE(
fthefi(ninter),stat=ierror1)
3693 ierror = ierror + ierror1
3694 ALLOCATE(
tempfi(ninter),stat=ierror1)
3695 ierror = ierror + ierror1
3696 ALLOCATE(
matsfi(ninter),stat=ierror1)
3697 ierror = ierror + ierror1
3698 ALLOCATE(
nmtemp(ninter),stat=ierror1)
3699 ierror = ierror + ierror1
3700 IF(idt_therm == 1)
THEN
3701 ALLOCATE(
condnfi(ninter),stat=ierror1)
3702 ierror = ierror + ierror1
3704 ierror = ierror + ierror1
3707 IF(intheat > 0 .OR.nintloadp21 > 0)
THEN
3708 ALLOCATE(
nmnfi(ninter),stat=ierror1)
3709 ierror = ierror + ierror1
3710 ALLOCATE(
nmvfi(ninter),stat=ierror1)
3711 ierror = ierror + ierror1
3712 ALLOCATE(
nmnsi(ninter),stat=ierror1)
3713 ierror = ierror + ierror1
3714 ALLOCATE(
nmvsi(ninter),stat=ierror1)
3715 ierror = ierror + ierror1
3716 ALLOCATE(
tempnod(ninter),stat=ierror1)
3717 ierror = ierror + ierror1
3719 IF(nintloadp21 > 0)
THEN
3721 ierror = ierror + ierror1
3723 IF(ninterfric > 0 )
THEN
3725 ierror = ierror + ierror1
3727 ierror = ierror + ierror1
3730 ALLOCATE (
nodnxfi(ninter),stat=ierror1)
3731 ierror = ierror + ierror1
3732 ALLOCATE (
nodamsfi(ninter),stat=ierror1)
3733 ierror = ierror + ierror1
3734 ALLOCATE (
procamsfi(ninter),stat=ierror1)
3735 ierror = ierror + ierror1
3738 ierror = ierror + ierror1
3740 ierror = ierror + ierror1
3742 IF(intheat > 0.OR.interadhesion > 0 )
THEN
3743 ALLOCATE(
areasfi(ninter),stat=ierror1)
3747 ierror = ierror + ierror1
3750 IF(interadhesion > 0)
THEN
3751 ALLOCATE(
if_adhfi(ninter),stat=ierror1)
3755 ierror = ierror + ierror1
3758 ALLOCATE(
candf_si(ninter),stat=ierror1)
3759 ierror = ierror + ierror1
3761 IF(nitsche > 0 )
THEN
3762 ALLOCATE(
forneqsfi(ninter),stat=ierror1)
3763 ierror = ierror + ierror1
3767 ALLOCATE(
efricfi(ninter),stat=ierror1)
3768 ierror = ierror + ierror1
3769 ALLOCATE(
efricgfi(ninter),stat=ierror1)
3770 ierror = ierror + ierror1
3809 IF(parameters%INTCAREA > 0)
NULLIFY(
intareanfi(i)%p)
3897 IF(parameters%ISTIF_DT > 0 .OR. ipari(97,i)>0 )
THEN
3913 IF(intheat /= 0)
THEN
3922 IF(interadhesion /= 0)
THEN
3927 IF(idt_therm == 1)
THEN
3934 IF(intheat /= 0.OR.nintloadp21 > 0)
THEN
3943 IF(nintloadp21 > 0)
THEN
3949 IF(ninterfric > 0 )
THEN
3956 IF(nitsche > 0 )
THEN
3963 CALL ancmsg(msgid=20,anmode=aninfo)
3974 inacti = ipari(22,i)
3976 intfric = ipari(72,i)
3977 flagremn = ipari(63,i)
3978 intnitsche = ipari(86,i)
3983 iedge4 = ipari(59,i)
3987 IF(ityp==7.OR.ityp==10.OR.ityp==11.OR.
3988 + (ityp==17.AND.ipari(33,i)==0).OR.ityp==20.OR.
3989 + ityp==22.OR.ityp==23.OR.ityp==24.OR.
3993 ALLOCATE(
nsnsi(i)%P(nspmd),stat=ierror)
3994 ALLOCATE(
nsnfi(i)%P(nspmd),stat=ierror1)
3997 IF(ierror+ierror1/=0)
THEN
3998 CALL ancmsg(msgid=20,anmode=aninfo)
4006 lens = lens +
nsnsi(i)%P(p)
4007 lenr = lenr +
nsnfi(i)%P(p)
4014 ALLOCATE(
nsvsi(i)%P(lens),stat=ierror1)
4016 ierror = ierror + ierror1
4018 ALLOCATE(
candf_si(i)%P(nsn),stat=ierror1)
4022 ALLOCATE(
nsvfi(i)%P(lenr),stat=ierror1)
4023 ierror = ierror + ierror1
4025 IF(ityp==7.OR.ityp==10.OR.
4026 + ityp==20.OR.ityp==22.OR.
4027 + ityp==23.OR.ityp==24.OR.
4029 ALLOCATE(
itafi(i)%P(lenr),stat=ierror1)
4030 ierror = ierror + ierror1
4032 IF(ityp==7.OR.ityp==20.OR.ityp==22.OR.
4033 + ityp==23.OR.ityp==24.OR.ityp==25)
THEN
4034 ALLOCATE(
kinfi(i)%P(lenr),stat=ierror1)
4035 ierror = ierror + ierror1
4039 ALLOCATE(
nbinflfi(i)%P(lenr),stat=ierror1)
4040 ierror = ierror + ierror1
4045 ALLOCATE(
matsfi(i)%P(lenr),stat=ierror1)
4046 ierror = ierror + ierror
4049 IF((ityp==7.OR.ityp==24.OR.ityp==25).AND.intfric > 0)
THEN
4051 ierror = ierror + ierror
4057 *
ALLOCATE(
irtlm_fi(i)%P(2,lenr),stat=ierror1)
4061 *
ALLOCATE(
icont_i_fi(i)%P(lenr),stat=ierror1)
4065 *
ALLOCATE(
isedge_fi(i)%P(lenr),stat=ierror1)
4070 *
ALLOCATE(
irtse_fi(i)%P(5,lenr),stat=ierror1)
4074 *
ALLOCATE(
is2pt_fi(i)%P(lenr),stat=ierror1)
4078 *
ALLOCATE(
ispt2_fi(i)%P(lenr),stat=ierror1)
4082 *
ALLOCATE(
isegpt_fi(i)%P(lenr),stat=ierror1)
4086 *
ALLOCATE(
is2se_fi(i)%P(2,lenr),stat=ierror1)
4090 IF(intnitsche >0)
THEN
4092 *
ALLOCATE(
forneqsfi(i)%P(3,lenr),stat=ierror1)
4099 IF(.NOT.
ASSOCIATED(
pmainfi(i)%P))
4100 *
ALLOCATE(
pmainfi(i)%P(lenr),stat=ierror1)
4103 *
ALLOCATE(
irtlm_fi(i)%P(4,lenr),stat=ierror1)
4106 *
ALLOCATE(
icont_i_fi(i)%P(lenr),stat=ierror1)
4110 *
ALLOCATE(
icodt_fi(i)%P(lenr),stat=ierror1)
4115 *
ALLOCATE(
iskew_fi(i)%P(lenr),stat=ierror1)
4121 *
ALLOCATE(
islide_fi(i)%P(4,lenr),stat=ierror1)
4125 IF(flagremn==2)
THEN
4130 IF(sizremnorfi /= 0)
THEN
4132 *
ALLOCATE(
remnor_fi(i)%P(sizremnorfi),stat=ierror1)
4134 ELSE IF(sizremnorfi == 0)
THEN
4136 *
ALLOCATE(
remnor_fi(i)%P(sizremnorfi),stat=ierror1)
4142 IF(idtmins_old == 2)
THEN
4143 ALLOCATE (
nodnxfi(i)%P(lenr),stat=ierror1)
4144 ierror = ierror + ierror1
4146 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4147 ierror = ierror + ierror1
4149 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4150 ierror = ierror + ierror1
4154 ierror = ierror + ierror1
4157 ELSEIF(idtmins_int_old /= 0)
THEN
4158 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4159 ierror = ierror + ierror1
4161 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4162 ierror = ierror + ierror1
4166 ierror = ierror + ierror1
4172 IF(.NOT.
ASSOCIATED(
nodnxfi(i)%P))
THEN
4173 ALLOCATE (
nodnxfi(i)%P(lenr),stat=ierror1)
4174 ierror = ierror + ierror1
4177 IF(.NOT.
ASSOCIATED(
nodamsfi(i)%P))
THEN
4178 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4179 ierror = ierror + ierror1
4182 IF(.NOT.
ASSOCIATED(
procamsfi(i)%P))
THEN
4183 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4184 ierror = ierror + ierror1
4190 ierror = ierror + ierror1
4194 ELSEIF(idtmins_int /= 0)
THEN
4195 IF(.NOT.
ASSOCIATED(
nodamsfi(i)%P))
THEN
4196 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4197 ierror = ierror + ierror1
4200 IF(.NOT.
ASSOCIATED(
procamsfi(i)%P))
THEN
4201 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4202 ierror = ierror + ierror1
4208 ierror = ierror + ierror1
4214 ALLOCATE(
msfi(i)%P(lenr),stat=ierror1)
4215 ierror = ierror + ierror1
4216 ALLOCATE(
stifi(i)%P(lenr),stat=ierror1)
4217 ierror = ierror + ierror1
4219 ALLOCATE(
gapfi(i)%P(lenr),stat=ierror1)
4220 ierror = ierror + ierror1
4222 ALLOCATE(
gap_lfi(i)%P(lenr),stat=ierror1)
4223 ierror = ierror + ierror1
4226 ALLOCATE(
xfi(i)%P(3,lenr),stat=ierror1)
4227 ierror = ierror + ierror1
4228 ALLOCATE(
vfi(i)%P(3,lenr),stat=ierror1)
4229 ierror = ierror + ierror1
4232 ALLOCATE(
afi(i)%P(3,lenr*nthread),stat=ierror1)
4233 ierror = ierror + ierror1
4234 ALLOCATE(
stnfi(i)%P(lenr*nthread),stat=ierror1)
4235 ierror = ierror + ierror1
4237 ALLOCATE(
vscfi(i)%P(lenr*nthread),stat=ierror1)
4238 ierror = ierror + ierror1
4249 vscfi(i)%P(1:lenr*nthread)=zero
4256 ALLOCATE(
fthefi(i)%P(lenr*nthread),stat=ierror1)
4257 ierror = ierror + ierror1
4259 fthefi(i)%P(1:lenr*nthread)=zero
4261 ALLOCATE(
tempfi(i)%P(lenr),stat=ierror1)
4262 ierror = ierror + ierror1
4266 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4267 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4268 ierror = ierror + ierror1
4271 IF(ityp == 25.AND.ivis2==-1)
THEN
4272 ALLOCATE(
if_adhfi(i)%P(lenr),stat=ierror1)
4273 ierror = ierror + ierror1
4276 IF(idt_therm ==1.AND.intth > 0)
THEN
4277 ALLOCATE(
condnfi(i)%P(lenr*nthread),stat=ierror1)
4278 ierror = ierror + ierror1
4279 condnfi(i)%P(1:lenr*nthread)=zero
4289 ALLOCATE(
tempfi(i)%P(lenr),stat=ierror1)
4290 ierror = ierror + ierror1
4293 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4294 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4295 ierror = ierror + ierror1
4297 IF(ityp == 25.AND.ivis2==-1)
THEN
4298 ALLOCATE(
if_adhfi(i)%P(lenr),stat=ierror1)
4299 ierror = ierror + ierror1
4304 ALLOCATE(
time_sfi(i)%P(lenr),stat=ierror1)
4306 ierror = ierror + ierror1
4307 ELSEIF(ityp == 25)
THEN
4308 ALLOCATE(
time_sfi(i)%P(2*lenr),stat=ierror1)
4310 ierror = ierror + ierror1
4312 IF(ityp == 24 .OR. ityp == 25)
THEN
4313 ALLOCATE(
secnd_frfi(i)%P(6,lenr),stat=ierror1)
4315 ierror = ierror + ierror1
4316 ALLOCATE(
pene_oldfi(i)%P(5,lenr),stat=ierror1)
4318 ierror = ierror + ierror1
4319 ALLOCATE(
stif_oldfi(i)%P(2,lenr),stat=ierror1)
4321 ierror = ierror + ierror1
4323 IF(ipari(97,i) > 0)
THEN
4328 IF(ityp == 25.AND.parameters%INTCAREA > 0.AND.ipari(36,i)>0)
THEN
4330 ierror = ierror + ierror1
4333 IF(ityp == 24.AND.parameters%INTCAREA > 0)
THEN
4335 ierror = ierror + ierror1
4339 ELSEIF(ityp==11)
THEN
4341 ALLOCATE(
itafi(i)%P(2*lenr),stat=ierror1)
4342 ierror = ierror + ierror1
4345 IF(idtmins_old == 2)
THEN
4346 ALLOCATE (
nodnxfi(i)%P(2*lenr),stat=ierror1)
4347 ierror = ierror + ierror1
4349 ALLOCATE (
nodamsfi(i)%P(2*lenr),stat=ierror1)
4350 ierror = ierror + ierror1
4352 ALLOCATE (
procamsfi(i)%P(2*lenr),stat=ierror1)
4353 ierror = ierror + ierror1
4355 ELSEIF(idtmins_int_old /= 0)
THEN
4356 ALLOCATE (
nodamsfi(i)%P(2*lenr),stat=ierror1)
4357 ierror = ierror + ierror1
4360 ierror = ierror + ierror1
4365 ALLOCATE(
matsfi(i)%P(lenr),stat=ierror1)
4366 ierror = ierror + ierror1
4370 IF(intfric > 0)
THEN
4372 ierror = ierror + ierror
4376 ALLOCATE(
msfi(i)%P(2*lenr),stat=ierror1)
4377 ierror = ierror + ierror1
4378 ALLOCATE(
stifi(i)%P(lenr),stat=ierror1)
4379 ierror = ierror + ierror1
4381 ALLOCATE(
gapfi(i)%P(lenr),stat=ierror1)
4382 ierror = ierror + ierror1
4384 ALLOCATE(
gap_lfi(i)%P(lenr),stat=ierror1)
4385 ierror = ierror + ierror1
4388 ALLOCATE(
xfi(i)%P(3,2*lenr),stat=ierror1)
4389 ierror = ierror + ierror1
4390 ALLOCATE(
vfi(i)%P(3,2*lenr),stat=ierror1)
4391 ierror = ierror + ierror1
4392 IF(inacti==5.OR.inacti==6)
THEN
4393 ALLOCATE(
penfi(i)%P(2,lenr),stat=ierror1)
4396 ALLOCATE(
afi(i)%P(3,2*lenr*nthread),stat=ierror1)
4397 ierror = ierror + ierror1
4398 ALLOCATE(
stnfi(i)%P(2*lenr*nthread),stat=ierror1)
4399 ierror = ierror + ierror1
4401 ALLOCATE(
vscfi(i)%P(2*lenr*nthread),stat=ierror1)
4402 ierror = ierror + ierror1
4405 ALLOCATE(
fthefi(i)%P(2*lenr*nthread),stat=ierror1)
4406 ierror = ierror + ierror1
4407 ALLOCATE(
tempfi(i)%P(2*lenr),stat=ierror1)
4408 ierror = ierror + ierror1
4411 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4412 ierror = ierror + ierror1
4414 IF(idt_therm ==1.AND.intth > 0)
THEN
4415 ALLOCATE(
condnfi(i)%P(2*lenr*nthread),stat=ierror1)
4416 ierror = ierror + ierror1
4425 ALLOCATE(
tempfi(i)%P(2*lenr),stat=ierror1)
4426 ierror = ierror + ierror1
4431 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4432 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4433 ierror = ierror + ierror1
4437 ELSEIF(ityp==17)
THEN
4438 ALLOCATE(
xfi17(i)%P(3,16,lenr),stat=ierror1)
4439 ierror = ierror + ierror1
4440 ALLOCATE(
vfi17(i)%P(3,16,lenr),stat=ierror1)
4441 ierror = ierror + ierror1
4442 ALLOCATE(
frotsfi(i)%P(7,lenr),stat=ierror1)
4443 ierror = ierror + ierror1
4444 ALLOCATE(
ksfi(i)%P(2,lenr),stat=ierror1)
4445 ierror = ierror + ierror1
4446 ALLOCATE(
eminxfi(i)%P(6,lenr),stat=ierror1)
4447 ierror = ierror + ierror1
4449 ALLOCATE(
afi17(i)%P(3,16,lenr),stat=ierror1)
4450 ierror = ierror + ierror1
4451 ALLOCATE(
stnfi17(i)%P(16,lenr),stat=ierror1)
4452 ierror = ierror + ierror1
4455 lskyfi = lenr * multimax
4458 ALLOCATE(
iskyfi(i)%P(lskyfi),stat=ierror1)
4459 ierror = ierror + ierror1
4460 ALLOCATE(
fskyfi(i)%P(40,lskyfi),stat=ierror1)
4461 ierror = ierror + ierror1
4471 CALL ancmsg(msgid=20,anmode=aninfo)
4476 ALLOCATE(
daanc6fi(i)%P(3,6,lenr),stat=ierror1)
4477 ierror = ierror + ierror1
4478 ALLOCATE(
dxancfi(i)%P(3,lenr),stat=ierror1)
4479 ierror = ierror + ierror1
4480 ALLOCATE(
dvancfi(i)%P(3,lenr),stat=ierror1)
4481 ierror = ierror + ierror1
4482 IF((inacti==5.OR.inacti==6) .AND. lenr > 0)
THEN
4483 ALLOCATE(
penfi(i)%P(2,lenr),stat=ierror1)
4484 ierror = ierror + ierror1
4485 ALLOCATE(
penfia(i)%P(5,lenr),stat=ierror1)
4486 ierror = ierror + ierror1
4488 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
4489 ALLOCATE (
diag_smsfi(i)%P(lenr),stat=ierror1)
4490 ierror = ierror + ierror1
4492 ALLOCATE(
alphakfi(i)%P(lenr),stat=ierror1)
4493 ierror = ierror + ierror1
4494 ALLOCATE(
daancfi(i)%P(3,lenr),stat=ierror1)
4496 ierror = ierror + ierror1
4498 CALL ancmsg(msgid=20,anmode=aninfo)
4502 ALLOCATE(
nsnsie(i)%P(nspmd),stat=ierror1)
4503 ierror = ierror + ierror1
4504 ALLOCATE(
nsnfie(i)%P(nspmd),stat=ierror1)
4505 ierror = ierror + ierror1
4511 lens = lens +
nsnsie(i)%P(p)
4512 lenr = lenr +
nsnfie(i)%P(p)
4516 ALLOCATE(
nsvsie(i)%P(lens),stat=ierror1)
4518 ierror = ierror + ierror1
4522 ALLOCATE(
nsvfie(i)%P(lenr),stat=ierror1)
4523 ierror = ierror + ierror1
4525 ALLOCATE(
itafie(i)%P(2*lenr),stat=ierror1)
4526 ierror = ierror + ierror1
4528 ALLOCATE(
msfie(i)%P(2*lenr),stat=ierror1)
4529 ierror = ierror + ierror1
4530 ALLOCATE(
stifie(i)%P(lenr),stat=ierror1)
4531 ierror = ierror + ierror1
4533 ALLOCATE(
gapfie(i)%P(lenr),stat=ierror1)
4534 ierror = ierror + ierror1
4536 ALLOCATE(
xfie(i)%P(3,2*lenr),stat=ierror1)
4537 ierror = ierror + ierror1
4538 ALLOCATE(
vfie(i)%P(3,2*lenr),stat=ierror1)
4539 ierror = ierror + ierror1
4540 IF(inacti==5.OR.inacti==6)
THEN
4541 ALLOCATE(
penfie(i)%P(2,lenr),stat=ierror1)
4542 ierror = ierror + ierror1
4543 ALLOCATE(
penfiae(i)%P(5,2*lenr),stat=ierror1)
4544 ierror = ierror + ierror1
4547 ALLOCATE(
afie(i)%P(3,2*lenr*nthread),stat=ierror1)
4548 ierror = ierror + ierror1
4549 ALLOCATE(
stnfie(i)%P(2*lenr*nthread),stat=ierror1)
4550 ierror = ierror + ierror1
4552 ALLOCATE(
vscfie(i)%P(2*lenr*nthread),stat=ierror1)
4553 ierror = ierror + ierror1
4557 lskyfi = lenr * multimax
4559 ALLOCATE(
iskyfie(i)%P(lskyfi),stat=ierror1)
4560 ierror = ierror + ierror1
4563 ALLOCATE(
fskyfie(i)%P(10,lskyfi),stat=ierror1)
4566 ALLOCATE(
fskyfie(i)%P(10,lskyfi),stat=ierror1)
4568 ierror = ierror + ierror1
4570 ALLOCATE(
daanc6fie(i)%P(3,6,2*lenr),stat=ierror1)
4571 ierror = ierror + ierror1
4572 ALLOCATE(
dxancfie(i)%P(3,2*lenr),stat=ierror1)
4573 ierror = ierror + ierror1
4574 ALLOCATE(
dvancfie(i)%P(3,2*lenr),stat=ierror1)
4575 ierror = ierror + ierror1
4576 ALLOCATE(
alphakfie(i)%P(2*lenr),stat=ierror1)
4577 ierror = ierror + ierror1
4578 ALLOCATE(
daancfie(i)%P(3,lenr),stat=ierror1)
4579 ierror = ierror + ierror1
4580 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
4582 ierror = ierror + ierror1
4585 IF(idtmins_old == 2)
THEN
4586 ALLOCATE (
nodnxfie(i)%P(lenr),stat=ierror1)
4587 ierror = ierror + ierror1
4589 ALLOCATE (
nodamsfie(i)%P(lenr),stat=ierror1)
4590 ierror = ierror + ierror1
4592 ALLOCATE (
procamsfie(i)%P(lenr),stat=ierror1)
4593 ierror = ierror + ierror1
4595 ELSEIF(idtmins_int_old /= 0)
THEN
4596 ALLOCATE (
nodamsfie(i)%P(lenr),stat=ierror1)
4597 ierror = ierror + ierror1
4599 ALLOCATE (
procamsfie(i)%P(lenr),stat=ierror1)
4600 ierror = ierror + ierror1
4605 CALL ancmsg(msgid=20,anmode=aninfo)
4614 IF(ityp == 25 )
THEN
4616 ALLOCATE(
nsnsie(i)%P(nspmd),stat=ierror1)
4617 ALLOCATE(
nsnfie(i)%P(nspmd),stat=ierror1)
4622 IF(ipari(58,i) /=0)
THEN
4630 lens_edge = lens_edge +
nsnsie(i)%P(p)
4631 lenr_edge = lenr_edge +
nsnfie(i)%P(p)
4634 IF(lens_edge>0)
THEN
4635 ALLOCATE(
nsvsie(i)%P(lens_edge),stat=ierror1)
4637 ierror = ierror + ierror1
4640 IF(iparit == 0)
THEN
4642 ALLOCATE(
afie(i)%P(3,2*lenr_edge*nthread),stat=ierror1)
4643 ierror = ierror + ierror1
4644 ALLOCATE(
stnfie(i)%P(2*lenr_edge*nthread),stat=ierror1)
4645 ierror = ierror + ierror1
4647 ALLOCATE(
vscfie(i)%P(2*lenr_edge*nthread),stat=ierror1)
4648 ierror = ierror + ierror1
4651 lskyfi = lenr_edge * multimax
4653 ALLOCATE(
iskyfie(i)%P(lskyfi),stat=ierror1)
4654 ierror = ierror + ierror1
4656 ALLOCATE(
fskyfie(i)%P(8,lskyfi),stat=ierror1)
4658 ALLOCATE(
fskyfie(i)%P(8,lskyfi),stat=ierror1)
4660 ierror = ierror + ierror1
4662 IF(lenr_edge>0)
THEN
4663 ALLOCATE(
nsvfie(i)%P(lenr_edge))
4664 ALLOCATE(
itafie(i)%P(lenr_edge*2))
4665 ALLOCATE(
ledge_fie(i)%P(e_ledge_size,lenr_edge))
4666 ALLOCATE(
xfie(i)%P(3,lenr_edge*2))
4667 ALLOCATE(
vfie(i)%P(3,lenr_edge*2))
4668 ALLOCATE(
msfie(i)%P(lenr_edge*2))
4669 ALLOCATE(
gapfie(i)%P(lenr_edge))
4673 ALLOCATE(
stifie(i)%P(lenr_edge))
4684 ipari(69,i) = lenr_edge
4685 IF(intfric > 0)
THEN
4689 IF(ipari(97,i) > 0)
THEN
4700 IF(ipari(36,i)>0.AND.ipari(7,i)/=17)
THEN
4703 ALLOCATE(
nisubsfi(i)%P(nspmd),stat=ierror1)
4704 ierror = ierror + ierror1
4710 lenr = lenr +
nsnfi(i)%P(p)
4715 ALLOCATE(
addsubsfi(i)%P(lenr+1),stat=ierror1)
4716 ierror = ierror + ierror1
4719 ALLOCATE(
lisubsfi(i)%P(leni),stat=ierror1)
4720 ierror = ierror + ierror1
4722 IF(ipari(7,i)==25.OR.ipari(7,i)==7.OR.ipari(7,i)==24.OR.ipari(7,i)==11)
THEN
4724 ierror = ierror + ierror1
4729 ALLOCATE(
addsubsfi(i)%P(lenr+1),stat=ierror1)
4736 CALL ancmsg(msgid=20,anmode=aninfo)
4741 IF(ipari(7,i) == 25 .AND. ipari(58,i) > 0)
THEN
4743 ALLOCATE(
nisubsfie(i)%P(nspmd),stat=ierror1)
4744 ierror = ierror + ierror1
4751 IF(lenr_edge>0)
THEN
4752 ALLOCATE(
addsubsfie(i)%P(lenr_edge+1),stat=ierror1)
4753 ierror = ierror + ierror1
4756 ALLOCATE(
lisubsfie(i)%P(leni),stat=ierror1)
4757 ierror = ierror + ierror1
4760 ierror = ierror + ierror1
4763 IF(lenr_edge>0)
THEN
4764 ALLOCATE(
addsubsfie(i)%P(lenr_edge+1),stat=ierror1)
4774 IF(anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT >0)
THEN
4776 ALLOCATE(
fnconti(i)%P(3,lenr),stat=ierror1)
4777 ierror = ierror + ierror1
4778 ALLOCATE(
ftconti(i)%P(3,lenr),stat=ierror1)
4779 ierror = ierror + ierror1
4781 CALL ancmsg(msgid=20,anmode=aninfo)
4801 IF(intth==2.OR.ipari(95,i) > 0)
THEN
4807 ALLOCATE(
nmnsi(i)%P(nspmd),stat=ierror1)
4808 ierror = ierror + ierror1
4810 ALLOCATE(
nmnfi(i)%P(nspmd),stat=ierror1)
4811 ierror = ierror + ierror1
4813 CALL ancmsg(msgid=20,anmode=aninfo)
4821 lens = lens +
nmnsi(i)%P(p)
4822 lenr = lenr +
nmnfi(i)%P(p)
4829 ALLOCATE(
nmvsi(i)%P(lens),stat=ierror1)
4831 ierror = ierror + ierror1
4833 ALLOCATE(
tempnod(i)%P(lens),stat=ierror1)
4835 ierror = ierror + ierror1
4840 ALLOCATE(
nmvfi(i)%P(lenr),stat=ierror1)
4841 ierror = ierror + ierror1
4844 CALL ancmsg(msgid=20,anmode=aninfo)
4850 IF(ipari(95,i) > 0)
THEN
4857 IF(ipari(48,i) /= 0 )
THEN
4859 ALLOCATE(
fthefi(i)%P(lenr*nthread),stat=ierror1)
4860 ierror = ierror + ierror1
4862 fthefi(i)%P(1:lenr*nthread)=zero
4864 lskyfi = lenr * multimax
4865 ALLOCATE(
iskyfi(i)%P(lskyfi),stat=ierror1)
4866 ierror = ierror + ierror1
4867 ALLOCATE(
ftheskyfi(i)%P(lskyfi),stat=ierror1)
4868 ierror = ierror + ierror1
4891 inacti = ipari(22,i)
4893 IF(ityp==7.OR.ityp==10.OR.ityp==11.OR.
4894 + (ityp==17.AND.ipari(33,i)==0).OR.ityp==20.OR.
4895 + ityp==22.OR.ityp==23.OR.ityp==24.OR.
4899 lenr = lenr +
nsnfi(i)%P(p)
4902 IF(ityp==7.OR.ityp==22.OR.ityp==23.OR.
4903 + ityp==20.OR.ityp==10.OR.ityp==24.OR.
4916 IF((ityp == 7.OR.ityp == 20.OR.
4917 + ityp==22.OR.ityp==23.OR.ityp==24.OR.
4918 + ityp==25).AND.ipari(47,i) > 0)
THEN
4921 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4924 IF(ityp == 25.AND.ivis2==-1)
THEN
4930 ELSEIF(ityp==11)
THEN
4941 IF(inacti==5.OR.inacti==6)
THEN
4944 IF(ipari(47,i)>0)
THEN
4948 ELSEIF(ityp==17.AND.ipari(33,i)==0)
THEN
4964 IF(inacti==5.OR.inacti==6)
THEN
4970 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
4977 lenr = lenr +
nsnfie(i)%P(p)
4987 IF(inacti==5.OR.inacti==6)
THEN
4994 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
5001 ELSEIF(ityp==25)
THEN
5004 IF (ityp==24.OR.ityp==25)
THEN
5008 IF(ipari(97,i) > 0)
THEN
5011 IF(ityp==25.AND.parameters%INTCAREA > 0.AND.ipari(36,i)>0)
THEN
5014 IF(ityp==24.AND.parameters%INTCAREA > 0)
THEN
5020 IF(idtmins_old == 2)
THEN
5022 ierror = ierror + ierror1
5024 ELSEIF(idtmins_int_old /= 0)
THEN
5026 ierror = ierror + ierror1
5033 ierror = ierror + ierror1
5036 ELSEIF(idtmins_int /= 0)
THEN
5039 ierror = ierror + ierror1
5045 IF(ityp == 25 )
THEN
5046 IF(ipari(58,i) /=0)
THEN
5049 lenr = lenr +
nsnfie(i)%P(p)
5063 IF(idtmins_old == 2)
THEN
5070 ELSEIF(idtmins_int_old /= 0)
THEN
5077 IF(ipari(97,i) > 0)
THEN
5090 lenr = lenr +
nmnfi(i)%P(p)
5098 ALLOCATE(
nmtemp(i)%P(lenr),stat=ierror1)
5099 ierror = ierror + ierror1
5103 CALL ancmsg(msgid=20,anmode=aninfo)