33!||--- uses -----------------------------------------------------
37 . INTHEAT,IDT_THERM,NODADT_THERM)
45#include "implicit_f.inc"
49 INTEGER ,
INTENT(IN) :: INTHEAT
50 INTEGER ,
INTENT(IN) :: IDT_THERM
51 INTEGER ,
INTENT(IN) :: NODADT_THERM
52 INTEGER :: NSNR, NSNROLD
53 INTEGER :: OLDNUM(*), RENUM(*)
69 + oldnum(renum(i)) = i
83 1 CAND_N ,RENUM ,II_STOK,NIN,NSN,
92#include "implicit_f.inc"
100 INTEGER CAND_N(*), NSNFIOLD(*), RENUM(*),
101 . NIN, NSN, NSNROLD,II_STOK
105 INTEGER I, J, P, I_STOK, IDEB, JDEB, NI
118 DO WHILE (j<=
nsnfi(nin)%P(p).AND.i<=nsnfiold(p))
119 IF(nint(xrem(4,j+jdeb))==
120 +
nsvfi(nin)%P(i+ideb))
THEN
121 renum(i+ideb) = j+jdeb
124 ELSEIF(nint(xrem(4,j+jdeb))<
125 +
nsvfi(nin)%P(i+ideb))
THEN
127 ELSEIF(nint(xrem(4,j+jdeb))>
128 +
nsvfi(nin)%P(i+ideb))
THEN
133 jdeb = jdeb +
nsnfi(nin)%P(p)
134 ideb = ideb + nsnfiold(p)
143 cand_n(i) = renum(ni) + nsn
160 1 CAND_N ,RENUM ,II_STOK,NIN,NSN,
169#include "implicit_f.inc"
173#include "com01_c.inc"
177 INTEGER CAND_N(*), NSNFIOLD(*), RENUM(*),
178 . NIN, NSN, NSNROLD,II_STOK
182 INTEGER I, J, P, I_STOK, IDEB, JDEB, NI
195 DO WHILE (j<=
nsnfi(nin)%P(p).AND.i<=nsnfiold(p))
197 +
nsvfi(nin)%P(i+ideb))
THEN
198 renum(i+ideb) = j+jdeb
201 ELSEIF(
irem(1,j+jdeb)<
202 +
nsvfi(nin)%P(i+ideb))
THEN
204 ELSEIF(
irem(1,j+jdeb)>
205 +
nsvfi(nin)%P(i+ideb))
THEN
210 jdeb = jdeb +
nsnfi(nin)%P(p)
211 ideb = ideb + nsnfiold(p)
220 cand_n(i) = renum(ni) + nsn
237 1 CAND_N ,RENUM ,II_STOK,NIN,NSN,
246#include "implicit_f.inc"
250#include "com01_c.inc"
254 INTEGER CAND_N(*), NSNFIOLD(*), RENUM(*),
255 . NIN, NSN, NSNROLD,II_STOK
259 INTEGER I, J, P, I_STOK, IDEB, JDEB, NI
274 DO WHILE (j<=
nsnfi(nin)%P(p).AND.i<=nsnfiold(p))
277 +
nsvfi(nin)%P(i+ideb))
THEN
278 renum(i+ideb) = j+jdeb
281 ELSEIF(
irem(1,j+jdeb)<
282 +
nsvfi(nin)%P(i+ideb))
THEN
284 ELSEIF(
irem(1,j+jdeb)>
285 +
nsvfi(nin)%P(i+ideb))
THEN
291 jdeb = jdeb +
nsnfi(nin)%P(p)
292 ideb = ideb + nsnfiold(p)
301 IF(renum(ni) > 0)
THEN
302 cand_n(i) = renum(ni) + nsn
305 WRITE(6,*)
"possible p/on issue for",
306 .
"/INTER/TYPE10 bouncing node",
itafi(nin)%P(ni)
315!||--- called by -------------------------------------
317!||--- uses -----------------------------------------------------
321 1 RENUM ,NIN,NSN,NSNFIOLD,NSNROLD)
329#include "implicit_f.inc"
333#include "com01_c.inc"
337 INTEGER NSNFIOLD(*), RENUM(*),
342 INTEGER I, J, P, I_STOK, IDEB, JDEB, NI
355 DO WHILE (j<=
nsnfi(nin)%P(p).AND.i<=nsnfiold(p))
357 +
nsvfi(nin)%P(i+ideb))
THEN
358 renum(i+ideb) = j+jdeb
361 ELSEIF(
irem(1,j+jdeb)<
362 +
nsvfi(nin)%P(i+ideb))
THEN
364 ELSEIF(
irem(1,j+jdeb)>
365 +
nsvfi(nin)%P(i+ideb))
THEN
370 jdeb = jdeb +
nsnfi(nin)%P(p)
371 ideb = ideb + nsnfiold(p)
385 1 CAND_N ,RENUM ,II_STOK,NIN,NRTS,
386 2 NSNFIOLD,NSNROLD,ADDCM,CHAINE,CAND_M,NSN4,NRTM)
394#include "implicit_f.inc"
398#include "com01_c.inc"
402 INTEGER CAND_N(*), NSNFIOLD(*), RENUM(*),
403 . NIN, NRTS, NSNROLD,ADDCM(*),CHAINE(2,*),
404 . CAND_M(*),NSN4,NRTM,II_STOK
408 INTEGER I, J, P, I_STOK, IDEB, JDEB, NI, IAD,IADM1
421 DO WHILE (j<=
nsnfi(nin)%P(p).AND.i<=nsnfiold(p))
423 +
nsvfi(nin)%P(i+ideb))
THEN
424 renum(i+ideb) = j+jdeb
427 ELSEIF(
irem(1,j+jdeb)<
428 +
nsvfi(nin)%P(i+ideb))
THEN
430 ELSEIF(
irem(1,j+jdeb)>
431 +
nsvfi(nin)%P(i+ideb))
THEN
436 jdeb = jdeb +
nsnfi(nin)%P(p)
437 ideb = ideb + nsnfiold(p)
446 cand_n(i) = renum(ni) + nrts
452 DO WHILE(iad>0.AND.j<nsn4)
454 IF(chaine(1,iad) > nrts )
THEN
455 chaine(1,iad) = renum(chaine(1,iad)-nrts)+nrts
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 INTEGER , LEN, IBC ,ISECIN ,IBAG , NOINT, INACTI,
511 . NSV(*), ICODT(*), NSTRF(*),ICONTACT(*),
512 . TAGNCONT(NLOADP_HYD_INTER,*),KLOADPINTER(*),LOADPINTER(*),
513 . LOADP_HYD_INTER(*),
515 INTEGER,
INTENT(in) :: NIN, INTCAREA
516 INTEGER,
INTENT(IN) :: NODADT_THERM
518 . BUFR(LEN,*), A(3,*), STIFN(*), VISCN(*),
519 . SECFCUM(7,NUMNOD,NSECT),
520 . FCONT(3,*),FTHE(*),CONDN(*)
521 my_real,
INTENT(INOUT) :: FSAV(*)
522 TYPE(H3D_DATABASE) :: H3D_DATA
523 TYPE(MULTI_FVM_STRUCT),
INTENT(INOUT) :: MULTI_FVM
524 TYPE (PARAMETERS_) ,
INTENT(IN):: PARAMETERS
528 my_real :: mass, fsav29
530 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER,PP,PPL,NP
535 IF (multi_fvm%INT18_GLOBAL_LIST(nin))
THEN
539 multi_fvm%FORCE_INT(1, ibric) = multi_fvm%FORCE_INT(1, ibric) + bufr(2,i)
540 multi_fvm%FORCE_INT(2, ibric) = multi_fvm%FORCE_INT(2, ibric) + bufr(3,i)
541 multi_fvm%FORCE_INT(3, ibric) = multi_fvm%FORCE_INT(3, ibric) + bufr(4,i)
548 a(1,nod) = a(1,nod) + bufr(2,i)
549 a(2,nod) = a(2,nod) + bufr(3,i)
550 a(3,nod) = a(3,nod) + bufr(4,i)
551 stifn(nod) = stifn(nod) + bufr(5,i)
552 IF(kdtint /= 0) viscn(nod) = viscn(nod) + bufr(6,i)
559 a(1,nod) = a(1,nod) + bufr(2,i)
560 a(2,nod) = a(2,nod) + bufr(3,i)
561 a(3,nod) = a(3,nod) + bufr(4,i)
562 stifn(nod) = stifn(nod) + bufr(5,i)
563 fthe(nod) = fthe(nod) + bufr(6,i)
564 IF(nodadt_therm == 1) condn(nod) = condn(nod
570 a(1,nod) = a(1,nod) + bufr(2,i)
571 a(2,nod) = a(2,nod) + bufr(3,i)
572 a(3,nod) = a(3,nod) + bufr(4,i)
573 stifn(nod) = stifn(nod) + bufr(5,i)
574 viscn(nod) = viscn(nod) + bufr(6,i)
575 fthe(nod) = fthe(nod) + bufr(7,i)
576 IF(nodadt_therm == 1) condn(nod) = condn(nod) + bufr(8,i)
583 IF (inconv == 1)
THEN
584 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
585 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
586 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
587 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
592 fcont(1,nod)=fcont(1,nod)+bufr(2,i)
593 fcont(2,nod)=fcont(2,nod)+bufr(3,i)
594 fcont(3,nod)=fcont(3,nod)+bufr(4,i)
599 IF(nintloadp > 0)
THEN
603 DO np = kloadpinter(nin)+1, kloadpinter(nin+1)
605 ppl = loadp_hyd_inter(pp)
606 tagncont(ppl,nod) = 1
614 IF(nstrf(1)+nstrf(2)/=0)
THEN
619 IF(nstrf(k1s)==noint)
THEN
624 IF(secfcum(4,nod,i)==1.)
THEN
625 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
626 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
627 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
639 IF((ibag/=0.AND.inacti/=7).OR.
640 . (iadm/=0).OR.(idamp_rdof/=0))
THEN
643 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
644 + bufr(4,i)/=zero)
THEN
654 ibcs = ibc - 8 * ibcm
660 CALL ibcoff(ibcs,icodt(nod))
678 1 NB ,LEN ,BUFR ,IRECTS ,A ,
679 2 STIFN ,VISCN ,IBC ,ISECIN ,NOINT ,
680 3 IBAG ,ICODT ,SECFCUM,NSTRF ,ICONTACT,
681 4 FCONT ,INTTH ,FTHE ,CONDN ,H3D_DATA,
682 5 TAGNCONT,KLOADPINTER,LOADPINTER,LOADP_HYD_INTER,NODADT_THERM)
691#include "implicit_f.inc"
695#include "scr07_c.inc"
696#include "scr14_c.inc"
697#include "scr16_c.inc"
698#include "scr18_c.inc"
699#include "com01_c.inc"
700#include "com04_c.inc"
701#include "com06_c.inc"
702#include "com08_c.inc"
703#include "impl1_c.inc"
707 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT,INTTH,
708 . irects(2,*), icodt(*), nstrf(*),icontact(*),
709 . tagncont(nloadp_hyd_inter,*),
710 . kloadpinter(*),loadpinter
711 INTEGER,
INTENT(IN) :: NODADT_THERM
713 . BUFR(LEN,*), A(3,*), STIFN(*), VISCN(*),
714 . SECFCUM(7,NUMNOD,NSECT),
715 . FCONT(3,*),FTHE(*),CONDN(*)
716 TYPE(H3D_DATABASE) :: H3D_DATA
720 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER,PP,PPL,NOD1,NOD2,NP
729 a(1,nod) = a(1,nod) + bufr(2,i)
730 a(2,nod) = a(2,nod) +
731 a(3,nod) = a(3,nod) + bufr(4,i)
734 a(1,nod) = a(1,nod) + bufr(6,i)
735 a(2,nod) = a(2,nod) + bufr(7,i)
736 a(3,nod) = a(3,nod) + bufr(8,i)
737 stifn(nod) = stifn(nod) + bufr(9,i)
743 a(1,nod) = a(1,nod) + bufr(2,i)
744 a(2,nod) = a(2,nod) + bufr(3,i)
745 a(3,nod) = a(3,nod) + bufr(4,i)
746 stifn(nod) = stifn(nod) + bufr(5,i)
747 viscn(nod) = viscn(nod) + bufr(6,i)
749 a(1,nod) = a(1,nod) + bufr(7,i)
750 a(2,nod) = a(2,nod) + bufr(8,i)
751 a(3,nod) = a(3,nod) + bufr(9,i)
752 stifn(nod) = stifn(nod) + bufr(10,i)
753 viscn(nod) = viscn(nod) + bufr(11,i)
757 IF(nodadt_therm ==1 )
THEN
762 a(1,nod) = a(1,nod) + bufr(2,i)
763 a(2,nod) = a(2,nod) + bufr(3,i)
764 a(3,nod) = a(3,nod) + bufr(4,i)
765 stifn(nod) = stifn(nod) + bufr(5,i)
766 fthe(nod) = fthe(nod) + bufr(6,i)
767 condn(nod) = condn(nod) + bufr(7,i)
769 a(1,nod) = a(1,nod) + bufr(8,i)
770 a(2,nod) = a(2,nod) + bufr(9,i)
771 a(3,nod) = a(3,nod) + bufr(10,i)
772 stifn(nod) = stifn(nod) + bufr(11,i)
773 fthe(nod) = fthe(nod) + bufr(12,i)
774 condn(nod) = condn(nod) + bufr(13,i)
780 a(1,nod) = a(1,nod) + bufr(2,i)
781 a(2,nod) = a(2,nod) + bufr(3,i)
782 a(3,nod) = a(3,nod) + bufr(4,i)
783 stifn(nod) = stifn(nod) + bufr(5,i)
784 viscn(nod) = viscn(nod) + bufr(6,i)
785 fthe(nod) = fthe(nod) + bufr(7,i)
786 condn(nod) = condn(nod) + bufr(8,i)
788 a(1,nod) = a(1,nod) + bufr(9,i)
789 a(2,nod) = a(2,nod) + bufr(10,i)
790 a(3,nod) = a(3,nod) + bufr(11,i)
791 stifn(nod) = stifn(nod) + bufr(12,i)
792 viscn(nod) = viscn(nod) + bufr(13,i)
793 fthe(nod) = fthe(nod) + bufr(14,i)
794 condn(nod) = condn(nod) + bufr(15,i)
802 a(1,nod) = a(1,nod) + bufr(
803 a(2,nod) = a(2,nod) + bufr(3,i)
804 a(3,nod) = a(3,nod) + bufr(4,i)
805 stifn(nod) = stifn(nod) + bufr(5,i)
806 fthe(nod) = fthe(nod) + bufr(6,i)
808 a(1,nod) = a(1,nod) + bufr(7,i)
809 a(2,nod) = a(2,nod) + bufr(8,i)
810 a(3,nod) = a(3,nod) + bufr(9,i)
811 stifn(nod) = stifn(nod) + bufr(10,i)
812 fthe(nod) = fthe(nod) + bufr(11,i)
818 a(1,nod) = a(1,nod) + bufr(2,i)
819 a(2,nod) = a(2,nod) + bufr(3,i)
820 a(3,nod) = a(3,nod) + bufr(4,i)
821 stifn(nod) = stifn(nod) + bufr(5,i)
822 viscn(nod) = viscn(nod) + bufr(6,i)
823 fthe(nod) = fthe(nod) + bufr(7,i)
825 a(1,nod) = a(1,nod) + bufr(8,i)
826 a(2,nod) = a(2,nod) + bufr(9,i)
827 a(3,nod) = a(3,nod) + bufr(10,i)
828 stifn(nod) = stifn(nod) + bufr(11,i)
829 viscn(nod) = viscn(nod) + bufr(12,i)
830 fthe(nod) = fthe(nod) + bufr(13,i)
839 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
840 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
841 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
842 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
848 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
849 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
850 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
852 fcont(1,nod)=fcont(1,nod)+ bufr
853 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
854 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
860 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
861 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
862 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
864 fcont(1,nod)=fcont(1,nod)+ bufr(7,i)
865 fcont(2,nod)=fcont(2,nod)+ bufr(8,i)
866 fcont(3,nod)=fcont(3,nod)+ bufr(9,i)
872 IF(nintloadp > 0)
THEN
877 DO np = kloadpinter(noint)+1, kloadpinter(noint+1)
879 ppl = loadp_hyd_inter(pp)
880 tagncont(ppl,nod1) = 1
881 tagncont(ppl,nod2) = 1
889 IF(nstrf(1)+nstrf(2)/=0)
THEN
894 IF(nstrf(k1s)==noint)
THEN
900 IF(secfcum(4,nod,i)==1.)
THEN
901 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
902 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
903 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
906 IF(secfcum(4,nod,i)==1.)
THEN
907 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(6,ii)
908 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(7,ii)
909 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(8,ii)
916 IF(secfcum(4,nod,i)==1.)
THEN
917 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
918 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
919 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
922 IF(secfcum(4,nod,i)==1.)
THEN
923 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(7,ii)
924 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(8,ii)
925 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(9,ii)
938 IF ((ibag/=0).OR.(idamp_rdof/=0))
THEN
942 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
943 + bufr(4,i)/=zero)
THEN
948 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
949 + bufr(8,i)/=zero)
THEN
956 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
957 + bufr(4,i)/=zero)
THEN
962 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
963 + bufr(9,i)/=zero)
THEN
983 1 NB ,LEN ,BUFR ,NELEMS ,A ,
984 2 STIFN ,FCONT ,IXS ,IXS16 ,FROTS ,
994#include "implicit_f.inc"
998#include "scr07_c.inc"
999#include "scr14_c.inc"
1000#include "scr16_c.inc"
1001#include "com06_c.inc"
1002#include "com08_c.inc"
1006 INTEGER NB, LEN, NN,
1007 . nelems(*), ixs(nixs,*) ,ixs16(8,*)
1009 . bufr(len,*), a(3,*), stifn(*),
1010 . fcont(3,*), frots(7,*)
1020 nn = nint(bufr(1,i))
1024 a(1,nod) = a(1,nod) + bufr(2,i)
1025 a(2,nod) = a(2,nod) + bufr(3,i)
1026 a(3,nod) = a(3,nod) + bufr(4,i)
1027 stifn(nod) = stifn(nod) + bufr(5,i)
1029 a(1,nod) = a(1,nod) + bufr(6,i)
1030 a(2,nod) = a(2,nod) + bufr(7,i)
1031 a(3,nod) = a(3,nod) + bufr(8,i)
1032 stifn(nod) = stifn(nod) + bufr(9,i)
1034 a(1,nod) = a(1,nod) + bufr(10,i)
1035 a(2,nod) = a(2,nod) + bufr(11,i)
1036 a(3,nod) = a(3,nod) + bufr(12,i)
1037 stifn(nod) = stifn(nod) + bufr(13,i)
1039 a(1,nod) = a(1,nod) + bufr(14,i)
1040 a(2,nod) = a(2,nod) + bufr(15,i)
1041 a(3,nod) = a(3,nod) + bufr(16,i)
1042 stifn(nod) = stifn(nod) + bufr(17,i)
1044 a(1,nod) = a(1,nod) + bufr(18,i)
1045 a(2,nod) = a(2,nod) + bufr(19,i)
1046 a(3,nod) = a(3,nod) + bufr(20,i)
1047 stifn(nod) = stifn(nod) + bufr(21,i)
1049 a(1,nod) = a(1,nod) + bufr(22,i)
1051 a(3,nod) = a(3,nod) + bufr(24,i)
1052 stifn(nod) = stifn(nod) + bufr(25,i)
1054 a(1,nod) = a(1,nod) + bufr(26,i)
1055 a(2,nod) = a(2,nod) + bufr(27,i)
1056 a(3,nod) = a(3,nod) + bufr(28,i)
1057 stifn(nod) = stifn(nod) + bufr(29,i)
1059 a(1,nod) = a(1,nod) + bufr(30,i)
1060 a(2,nod) = a(2,nod) + bufr(31,i)
1061 a(3,nod) = a(3,nod) + bufr(32,i)
1062 stifn(nod) = stifn(nod) + bufr(33,i)
1065 a(1,nod) = a(1,nod) + bufr(34,i)
1066 a(2,nod) = a(2,nod) + bufr(35,i)
1067 a(3,nod) = a(3,nod) + bufr(36,i)
1068 stifn(nod) = stifn(nod) + bufr(37,i)
1070 a(1,nod) = a(1,nod) + bufr(38,i)
1071 a(2,nod) = a(2,nod) + bufr(39,i)
1072 a(3,nod) = a(3,nod) + bufr(40,i)
1073 stifn(nod) = stifn(nod) + bufr(41,i)
1075 a(1,nod) = a(1,nod) + bufr(42,i)
1076 a(2,nod) = a(2,nod) + bufr(43,i)
1078 stifn(nod) = stifn(nod) + bufr(45,i)
1080 a(1,nod) = a(1,nod) + bufr(46,i)
1081 a(2,nod) = a(2,nod) + bufr(47,i)
1082 a(3,nod) = a(3,nod) + bufr(48,i)
1083 stifn(nod) = stifn(nod) + bufr(49,i)
1085 a(1,nod) = a(1,nod) + bufr(50,i)
1086 a(2,nod) = a(2,nod) + bufr(51,i)
1087 a(3,nod) = a(3,nod) + bufr(52,i)
1088 stifn(nod) = stifn(nod) + bufr(53,i)
1090 a(1,nod) = a(1,nod) + bufr(54,i)
1091 a(2,nod) = a(2,nod) + bufr(55,i)
1092 a(3,nod) = a(3,nod) + bufr(56,i)
1093 stifn(nod) = stifn(nod) + bufr(57,i)
1095 a(1,nod) = a(1,nod) + bufr(58,i)
1096 a(2,nod) = a(2,nod) + bufr(59,i)
1097 a(3,nod) = a(3,nod) + bufr(60,i)
1098 stifn(nod) = stifn(nod) + bufr(61,i)
1100 a(1,nod) = a(1,nod) + bufr(62,i)
1101 a(2,nod) = a(2,nod) + bufr(63,i)
1102 a(3,nod) = a(3,nod) + bufr(64,i)
1103 stifn(nod) = stifn(nod) + bufr(65,i)
1105 frots(1,nn) = frots(1,nn) + bufr(66,i)
1106 frots(2,nn) = frots(2,nn) + bufr(67,i)
1107 frots(3,nn) = frots(3,nn) + bufr(68,i)
1108 frots(4,nn) = frots(4,nn) + bufr(69,i)
1113 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
1114 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
1115 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
1116 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
1119 ne = nint(bufr(1,i))
1122 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
1123 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
1124 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
1126 fcont(1,nod)=fcont(1,nod)+ bufr(6,i)
1127 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
1128 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
1130 fcont(1,nod)=fcont(1,nod)+ bufr(10,i)
1131 fcont(2,nod)=fcont(2,nod)+ bufr(11,i)
1132 fcont(3,nod)=fcont(3,nod)+ bufr(12,i)
1134 fcont(1,nod)=fcont(1,nod)+ bufr(14,i)
1135 fcont(2,nod)=fcont(2,nod)+ bufr(15,i)
1136 fcont(3,nod)=fcont(3,nod)+ bufr(16,i)
1138 fcont(1,nod)=fcont(1,nod)+ bufr(18,i)
1139 fcont(2,nod)=fcont(2,nod)+ bufr(19,i)
1140 fcont(3,nod)=fcont(3,nod)+ bufr(20,i)
1142 fcont(1,nod)=fcont(1,nod)+ bufr(22,i)
1143 fcont(2,nod)=fcont(2,nod)+ bufr(23,i)
1144 fcont(3,nod)=fcont(3,nod)+ bufr(24,i)
1146 fcont(1,nod)=fcont(1,nod)+ bufr(26,i)
1147 fcont(2,nod)=fcont(2,nod)+ bufr(27,i)
1148 fcont(3,nod)=fcont(3,nod)+ bufr(28,i)
1150 fcont(1,nod)=fcont(1,nod)+ bufr(30,i)
1151 fcont(2,nod)=fcont(2,nod)+ bufr(31,i)
1152 fcont(3,nod)=fcont(3,nod)+ bufr(32,i)
1155 fcont(1,nod)=fcont(1,nod)+ bufr(34,i)
1156 fcont(2,nod)=fcont(2,nod)+ bufr(35,i)
1157 fcont(3,nod)=fcont(3,nod)+ bufr(36,i)
1159 fcont(1,nod)=fcont(1,nod)+ bufr(38,i)
1160 fcont(2,nod)=fcont(2,nod)+ bufr(39,i)
1161 fcont(3,nod)=fcont(3,nod)+ bufr(40,i)
1163 fcont(1,nod)=fcont(1,nod)+ bufr(42,i)
1164 fcont(2,nod)=fcont(2,nod)+ bufr(43,i)
1165 fcont(3,nod)=fcont(3,nod)+ bufr(44,i)
1167 fcont(1,nod)=fcont(1,nod)+ bufr(46,i)
1168 fcont(2,nod)=fcont(2,nod)+ bufr(47,i)
1169 fcont(3,nod)=fcont(3,nod)+ bufr(48,i)
1171 fcont(1,nod)=fcont(1,nod)+ bufr(50,i)
1172 fcont(2,nod)=fcont(2,nod)+ bufr(51,i)
1173 fcont(3,nod)=fcont(3,nod)+ bufr(52,i)
1175 fcont(1,nod)=fcont(1,nod)+ bufr(54,i)
1176 fcont(2,nod)=fcont(2,nod)+ bufr(55,i)
1177 fcont(3,nod)=fcont(3,nod)+ bufr(56,i)
1179 fcont(1,nod)=fcont(1,nod)+ bufr(58,i)
1180 fcont(2,nod)=fcont(2,nod)+ bufr
1181 fcont(3,nod)=fcont(3,nod)+ bufr(60,i)
1183 fcont(1,nod)=fcont(1,nod)+ bufr(62,i)
1185 fcont(3,nod)=fcont(3,nod)+ bufr(64,i)
1197!||
getdpdaanc ../engine/source/mpi/interfaces/spmd_i7tool.f
1204 1 NB ,LEN ,BUFR ,NSV ,A ,
1205 2 STIFN ,VISCN ,IBC ,ISECIN ,NOINT ,
1206 3 IBAG ,ICODT ,SECFCUM,NSTRF ,ICONTACT,
1207 4 FCONT ,INACTI ,IADM ,INTTH ,DAANC6 ,
1208 5 FTHE ,NLG ,ALPHAK ,H3D_DATA)
1217#include "implicit_f.inc"
1221#include "scr05_c.inc"
1222#include "scr07_c.inc"
1223#include "scr14_c.inc"
1224#include "scr16_c.inc"
1225#include "scr18_c.inc"
1226#include "com01_c.inc"
1227#include "com04_c.inc"
1228#include "com06_c.inc"
1229#include "com08_c.inc"
1233 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT, INACTI,
1234 . NSV(*), ICODT(*), NSTRF(*), NLG(*),
1235 . ICONTACT(*), IADM,INTTH
1237 . BUFR(LEN,*), A(3,*), STIFN(*), VISCN(*),
1238 . SECFCUM(7,NUMNOD,NSECT),
1239 . FCONT(3,*),FTHE(*), ALPHAK(3,*)
1240 DOUBLE PRECISION DAANC6(3,6,*)
1241 TYPE(H3D_DATABASE) :: H3D_DATA
1245 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER, INC,ISIGN
1246 DOUBLE PRECISION DAANC6L(3,6)
1251 IF(INTTH == 0 ) then
1254 n = nsv(nint(bufr(1,i)))
1256 a(1,nod) = a(1,nod) + bufr(2,i)
1257 a(2,nod) = a(2,nod) + bufr(3,i)
1258 a(3,nod) = a(3,nod) + bufr(4,i)
1259 stifn(nod) = stifn(nod) + bufr(5,i)
1262 IF(alphak(2,n) < zero .or. bufr(6,i) < zero)isign = -1
1263 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(6,i)))
1265 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1266 daanc6(1,2,n) = daanc6(1,2,n)+daanc6l(1,2)
1267 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1268 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1269 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1270 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1271 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1272 daanc6(2,2,n) = daanc6(2,2,n)+daanc6l(2,2)
1273 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1274 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1275 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1276 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1277 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1278 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1279 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1280 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1281 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1282 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1286 n = nsv(nint(bufr(1,i)))
1288 a(1,nod) = a(1,nod) + bufr(2,i)
1289 a(2,nod) = a(2,nod) + bufr(3,i)
1290 a(3,nod) = a(3,nod) + bufr(4,i)
1291 stifn(nod) = stifn(nod) + bufr(5,i)
1292 viscn(nod) = viscn(nod) + bufr(6,i)
1294 IF(alphak(2,n) < zero .or. bufr(7,i) < zero)isign = -1
1295 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(7,i)))
1298 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1300 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1301 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1302 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1303 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1304 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1306 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1307 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1308 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1309 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1310 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1311 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1312 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1313 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1314 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1315 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1322 n = nsv(nint(bufr(1,i)))
1324 a(1,nod) = a(1,nod) + bufr(2,i)
1325 a(2,nod) = a(2,nod) + bufr(3,i)
1326 a(3,nod) = a(3,nod) + bufr(4,i)
1327 stifn(nod) = stifn(nod) + bufr(5,i)
1328 fthe(nod) = fthe(nod) + bufr(6,i)
1331 IF(alphak(2,n) < zero .or. bufr(7,i) < zero)isign = -1
1332 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(7,i)))
1334 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1335 daanc6(1,2,n) = daanc6(1,2,n)+daanc6l(1,2)
1336 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1337 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1338 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1339 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1340 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1341 daanc6(2,2,n) = daanc6(2,2,n)+daanc6l(2,2)
1342 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1343 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1344 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1345 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1346 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1347 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1348 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1349 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1350 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1351 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1355 n = nsv(nint(bufr(1,i)))
1357 a(1,nod) = a(1,nod) + bufr(2,i)
1358 a(2,nod) = a(2,nod) + bufr(3,i)
1359 a(3,nod) = a(3,nod) + bufr(4,i)
1360 stifn(nod) = stifn(nod) + bufr(5,i)
1361 viscn(nod) = viscn(nod) + bufr(6,i)
1362 fthe(nod) = fthe(nod) + bufr(7,i)
1365 IF(alphak(2,n) < zero .or. bufr(8,i) < zero)isign = -1
1366 alphak(2,n) = isign*
min(abs(alphak(2,n)),abs(bufr(8,i)))
1368 daanc6(1,1,n) = daanc6(1,1,n)+daanc6l(1,1)
1369 daanc6(1,2,n) = daanc6(1,2,n)+daanc6l(1,2)
1370 daanc6(1,3,n) = daanc6(1,3,n)+daanc6l(1,3)
1371 daanc6(1,4,n) = daanc6(1,4,n)+daanc6l(1,4)
1372 daanc6(1,5,n) = daanc6(1,5,n)+daanc6l(1,5)
1373 daanc6(1,6,n) = daanc6(1,6,n)+daanc6l(1,6)
1374 daanc6(2,1,n) = daanc6(2,1,n)+daanc6l(2,1)
1375 daanc6(2,2,n) = daanc6(2,2,n)+daanc6l(2,2)
1376 daanc6(2,3,n) = daanc6(2,3,n)+daanc6l(2,3)
1377 daanc6(2,4,n) = daanc6(2,4,n)+daanc6l(2,4)
1378 daanc6(2,5,n) = daanc6(2,5,n)+daanc6l(2,5)
1379 daanc6(2,6,n) = daanc6(2,6,n)+daanc6l(2,6)
1380 daanc6(3,1,n) = daanc6(3,1,n)+daanc6l(3,1)
1381 daanc6(3,2,n) = daanc6(3,2,n)+daanc6l(3,2)
1382 daanc6(3,3,n) = daanc6(3,3,n)+daanc6l(3,3)
1383 daanc6(3,4,n) = daanc6(3,4,n)+daanc6l(3,4)
1384 daanc6(3,5,n) = daanc6(3,5,n)+daanc6l(3,5)
1385 daanc6(3,6,n) = daanc6(3,6,n)+daanc6l(3,6)
1392 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
1393 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
1394 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
1395 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
1401 fcont(2,nod)=fcont(2,nod)+bufr(3,i)
1402 fcont(3,nod)=fcont(3,nod)+bufr(4,i)
1409 IF(nstrf(1)+nstrf(2)/=0)
THEN
1411 nbinter=nstrf(k0+14)
1414 IF(nstrf(k1s)==noint)
THEN
1417 n = nint(bufr(1,ii))
1419 IF(secfcum(4,nod,i)==1.)
THEN
1420 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
1421 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
1422 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
1434 IF((ibag/=0.AND.inacti/=7).OR.
1438 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
1439 + bufr(4,i)/=zero)
THEN
1449 ibcs = ibc - 8 * ibcm
1455 CALL ibcoff(ibcs,icodt(nod))
1467!||--- calls -----------------------------------------------------
1474 1 NB ,LEN ,BUFR ,IXLINS ,A ,
1475 2 STIFN ,VISCN ,IBC ,ISECIN ,NOINT ,
1476 3 IBAG ,ICODT ,SECFCUM,NSTRF ,ICONTACT,
1477 4 FCONT ,DAANC6 ,NLG ,ALPHAK ,H3D_DATA)
1486#include "implicit_f.inc"
1490#include "scr05_c.inc"
1491#include "scr07_c.inc"
1492#include "scr14_c.inc"
1493#include "scr16_c.inc"
1494#include "scr18_c.inc"
1495#include "com01_c.inc"
1496#include "com04_c.inc"
1497#include "com06_c.inc"
1498#include "com08_c.inc"
1502 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT,
1503 . (2,*), ICODT(*), NSTRF(*),
1504 . ICONTACT(*), NLG(*)
1506 . BUFR(LEN,*), A(3,*), STIFN(*), VISCN(*),
1507 . secfcum(7,numnod,nsect),
1508 . fcont(3,*), alphak(3,*)
1509 DOUBLE PRECISION DAANC6(3,6,*)
1510 TYPE(H3D_DATABASE) :: H3D_DATA
1514 INTEGER I, J, II, N, N1, N2, N1G, N2G, K0, K1S, NBINTER, INC,ISIGN
1515 DOUBLE PRECISION DAANC6L(3,6)
1527 a(1,n1g) = a(1,n1g) + bufr(2,i)
1528 a(2,n1g) = a(2,n1g) + bufr(3,i)
1529 a(3,n1g) = a(3,n1g) + bufr(4,i)
1530 stifn(n1g) = stifn(n1g) + bufr(5,i)
1532 a(1,n2g) = a(1,n2g) + bufr(6,i)
1533 a(2,n2g) = a(2,n2g) + bufr(7,i)
1534 a(3,n2g) = a(3,n2g) + bufr(8,i)
1535 stifn(n2g) = stifn(n2g) + bufr(9,i)
1538 IF(alphak(2,n1) < zero .or. bufr(10,i) < zero)isign = -1
1539 alphak(2,n1) = isign*
min(abs(alphak(2,n1)),abs(bufr(10,i)))
1541 IF(alphak(2,n2) < zero .or. bufr(11,i) <
1542 alphak(2,n2) = isign*
min(abs(alphak(2,n2)),abs(bufr(11,i)))
1543 CALL getdpdaanc(daanc6l,bufr(12,i),iresp,inc)
1544 daanc6(1,1,n1) = daanc6(1,1,n1)+daanc6l(1,1)
1545 daanc6(1,2,n1) = daanc6(1,2,n1)+daanc6l(1,2)
1546 daanc6(1,3,n1) = daanc6(1,3,n1)+daanc6l(1,3)
1547 daanc6(1,4,n1) = daanc6(1,4,n1)+daanc6l(1,4)
1548 daanc6(1,5,n1) = daanc6(1,5,n1)+daanc6l(1,5)
1549 daanc6(1,6,n1) = daanc6(1,6,n1)+daanc6l(1,6)
1550 daanc6(2,1,n1) = daanc6(2,1,n1)+daanc6l(2,1)
1551 daanc6(2,2,n1) = daanc6(2,2,n1)+daanc6l(2,2)
1552 daanc6(2,3,n1) = daanc6(2,3,n1)+daanc6l(2,3)
1553 daanc6(2,4,n1) = daanc6(2,4,n1)+daanc6l(2,4)
1554 daanc6(2,5,n1) = daanc6(2,5,n1)+daanc6l(2,5)
1555 daanc6(2,6,n1) = daanc6(2,6,n1)+daanc6l(2,6)
1556 daanc6(3,1,n1) = daanc6(3,1,n1)+daanc6l(3,1)
1557 daanc6(3,2,n1) = daanc6(3,2,n1)+daanc6l(3,2)
1558 daanc6(3,3,n1) = daanc6(3,3,n1)+daanc6l(3,3)
1559 daanc6(3,4,n1) = daanc6(3,4,n1)+daanc6l(3,4)
1560 daanc6(3,5,n1) = daanc6(3,5,n1)+daanc6l(3,5)
1561 daanc6(3,6,n1) = daanc6(3,6,n1)+daanc6l(3,6)
1562 CALL getdpdaanc(daanc6l,bufr(12+inc,i),iresp,inc)
1563 daanc6(1,1,n2) = daanc6(1,1,n2)+daanc6l(1,1)
1564 daanc6(1,2,n2) = daanc6(1,2,n2)+daanc6l(1,2)
1565 daanc6(1,3,n2) = daanc6(1,3,n2)+daanc6l(1,3)
1566 daanc6(1,4,n2) = daanc6(1,4,n2)+daanc6l(1,4)
1567 daanc6(1,5,n2) = daanc6(1,5,n2)+daanc6l(1,5)
1568 daanc6(1,6,n2) = daanc6(1,6,n2)+daanc6l(1,6)
1569 daanc6(2,1,n2) = daanc6(2,1,n2)+daanc6l(2,1)
1570 daanc6(2,2,n2) = daanc6(2,2,n2)+daanc6l(2,2)
1571 daanc6(2,3,n2) = daanc6(2,3,n2)+daanc6l(2,3)
1572 daanc6(2,4,n2) = daanc6(2,4,n2)+daanc6l(2,4)
1573 daanc6(2,5,n2) = daanc6(2,5,n2)+daanc6l(2,5)
1574 daanc6(2,6,n2) = daanc6(2,6,n2)+daanc6l(2,6)
1575 daanc6(3,1,n2) = daanc6(3,1,n2)+daanc6l(3,1)
1576 daanc6(3,2,n2) = daanc6(3,2,n2)+daanc6l(3,2)
1577 daanc6(3,3,n2) = daanc6(3,3,n2)+daanc6l(3,3)
1578 daanc6(3,4,n2) = daanc6(3,4,n2)+daanc6l(3,4)
1579 daanc6(3,5,n2) = daanc6(3,5,n2)+daanc6l(3,5)
1580 daanc6(3,6,n2) = daanc6(3,6,n2)+daanc6l(3,6)
1589 a(1,n1g) = a(1,n1g) + bufr(2,i)
1590 a(2,n1g) = a(2,n1g) + bufr(3,i)
1591 a(3,n1g) = a(3,n1g) + bufr(4,i)
1592 stifn(n1g) = stifn(n1g) + bufr(5,i)
1593 viscn(n1g) = viscn(n1g) + bufr(6,i)
1595 a(1,n2g) = a(1,n2g) + bufr(7,i)
1597 a(3,n2g) = a(3,n2g) + bufr(9,i)
1598 stifn(n2g) = stifn(n2g) + bufr(10,i)
1599 viscn(n2g) = viscn(n2g) + bufr(11,i)
1602 IF(alphak(2,n1) < zero .or. bufr(12,i) < zero)isign = -1
1603 alphak(2,n1) = isign*
min(abs(alphak(2,n1)),abs(bufr(12,i)))
1605 IF(alphak(2,n2) < zero .or. bufr(13,i) < zero)isign = -1
1606 alphak(2,n2) = isign*
min(abs(alphak(2,n2)),abs(bufr(13,i)))
1607 CALL getdpdaanc(daanc6l,bufr(14,i),iresp,inc)
1608 daanc6(1,1,n1) = daanc6(1,1,n1)+daanc6l(1,1)
1609 daanc6(1,2,n1) = daanc6(1,2,n1)+daanc6l(1,2)
1610 daanc6(1,3,n1) = daanc6(1,3,n1)+daanc6l(1,3)
1611 daanc6(1,4,n1) = daanc6(1,4,n1)+daanc6l(1,4)
1612 daanc6(1,5,n1) = daanc6(1,5,n1)+daanc6l(1,5)
1613 daanc6(1,6,n1) = daanc6(1,6,n1)+daanc6l(1,6)
1614 daanc6(2,1,n1) = daanc6(2,1,n1)+daanc6l(2,1)
1615 daanc6(2,2,n1) = daanc6(2,2,n1)+daanc6l(2,2)
1616 daanc6(2,3,n1) = daanc6(2,3,n1)+daanc6l(2,3)
1617 daanc6(2,4,n1) = daanc6(2,4,n1)+daanc6l(2,4)
1618 daanc6(2,5,n1) = daanc6(2,5,n1)+daanc6l(2,5)
1619 daanc6(2,6,n1) = daanc6(2,6,n1)+daanc6l(2,6)
1620 daanc6(3,1,n1) = daanc6(3,1,n1)+daanc6l(3,1)
1621 daanc6(3,2,n1) = daanc6(3,2,n1)+daanc6l(3,2)
1622 daanc6(3,3,n1) = daanc6(3,3,n1)+daanc6l(3,3)
1623 daanc6(3,4,n1) = daanc6(3,4,n1)+daanc6l(3,4)
1624 daanc6(3,5,n1) = daanc6(3,5,n1)+daanc6l(3,5)
1625 daanc6(3,6,n1) = daanc6(3,6,n1)+daanc6l(3,6)
1626 CALL getdpdaanc(daanc6l,bufr(14+inc,i),iresp,inc)
1627 daanc6(1,1,n2) = daanc6(1,1,n2)+daanc6l(1,1)
1628 daanc6(1,2,n2) = daanc6(1,2,n2)+daanc6l(1,2)
1629 daanc6(1,3,n2) = daanc6(1,3,n2)+daanc6l(1,3)
1630 daanc6(1,4,n2) = daanc6(1,4,n2)+daanc6l(1,4)
1631 daanc6(1,5,n2) = daanc6(1,5,n2)+daanc6l(1,5)
1632 daanc6(1,6,n2) = daanc6(1,6,n2)+daanc6l(1,6)
1633 daanc6(2,1,n2) = daanc6(2,1,n2)+daanc6l(2,1)
1634 daanc6(2,2,n2) = daanc6(2,2,n2)+daanc6l(2,2)
1635 daanc6(2,3,n2) = daanc6(2,3,n2)+daanc6l(2,3)
1636 daanc6(2,4,n2) = daanc6(2,4,n2)+daanc6l(2,4)
1637 daanc6(2,5,n2) = daanc6(2,5,n2)+daanc6l(2,5)
1638 daanc6(2,6,n2) = daanc6(2,6,n2)+daanc6l(2,6)
1639 daanc6(3,1,n2) = daanc6(3,1,n2)+daanc6l(3,1)
1640 daanc6(3,2,n2) = daanc6(3,2,n2)+daanc6l(3,2)
1641 daanc6(3,3,n2) = daanc6(3,3,n2)+daanc6l(3,3)
1642 daanc6(3,4,n2) = daanc6(3,4,n2)+daanc6l(3,4)
1643 daanc6(3,5,n2) = daanc6(3,5,n2)+daanc6l(3,5)
1644 daanc6(3,6,n2) = daanc6(3,6,n2)+daanc6l(3,6)
1651 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
1652 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
1653 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
1654 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
1663 fcont(1,n1g)=fcont(1,n1g)+bufr(2,i)
1664 fcont(2,n1g)=fcont(2,n1g)+bufr(3,i)
1665 fcont(3,n1g)=fcont(3,n1g)+bufr(4,i)
1667 fcont(1,n2g)=fcont(1,n2g)+bufr(5,i)
1668 fcont(2,n2g)=fcont(2,n2g)+bufr(6,i)
1669 fcont(3,n2g)=fcont(3,n2g)+bufr(7,i)
1676 IF(nstrf(1)+nstrf(2)/=0)
THEN
1678 nbinter=nstrf(k0+14)
1681 IF(nstrf(k1s)==noint)
THEN
1685 n = nint(bufr(1,ii))
1690 IF(secfcum(4,n1g,i)==1.)
THEN
1691 secfcum(1,n1g,i)=secfcum(1,n1g,i)+bufr(2,ii)
1692 secfcum(2,n1g,i)=secfcum(2,n1g,i)+bufr(3,ii)
1693 secfcum(3,n1g,i)=secfcum(3,n1g,i)+bufr(4,ii)
1695 IF(secfcum(4,n2g,i)==1.)
THEN
1696 secfcum(1,n2g,i)=secfcum(1,n2g,i)+bufr(6,ii)
1697 secfcum(2,n2g,i)=secfcum(2,n2g,i)+bufr(7,ii)
1698 secfcum(3,n2g,i)=secfcum(3,n2g,i)+bufr(8,ii)
1703 n = nint(bufr(1,ii))
1708 IF(secfcum(4,n1g,i)==1.)
THEN
1709 secfcum(1,n1g,i)=secfcum(1,n1g,i)+bufr(2,ii)
1710 secfcum(2,n1g,i)=secfcum(2,n1g,i)+bufr(3,ii)
1711 secfcum(3,n1g,i)=secfcum(3,n1g,i)+bufr(4,ii)
1713 IF(secfcum(4,n2g,i)==1.)
THEN
1714 secfcum(1,n2g,i)=secfcum(1,n2g,i)+bufr(7,ii)
1715 secfcum(2,n2g,i)=secfcum(2,n2g,i)+bufr(8,ii)
1716 secfcum(3,n2g,i)=secfcum(3,n2g,i)+bufr(9,ii)
1733 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
1734 + bufr(4,i)/=zero)
THEN
1740 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
1741 + bufr(8,i)/=zero)
THEN
1750 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
1751 + bufr(4,i)/=zero)
THEN
1757 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
1758 + bufr(9,i)/=zero)
THEN
1777!||
ancmsg ../engine/source/output/message/message.f
1787 1 NB ,LEN ,BUFR ,NSV ,FSKYI,
1788 2 ISKY ,IBC ,ISECIN ,NOINT ,IBAG ,
1789 3 ICODT ,SECFCUM,NSTRF ,ICONTACT,FCONT,
1790 4 INACTI ,IADM ,INTTH ,FTHESKYI,CONDNSKYI,
1791 5 H3D_DATA,NIN ,TAGNCONT,KLOADPINTER,LOADPINTER ,
1792 6 LOADP_HYD_INTER,INTCAREA,FSAV ,PARAMETERS,NODADT_THERM)
1803#include "implicit_f.inc"
1807#include "parit_c.inc"
1808#include "scr07_c.inc"
1809#include "scr14_c.inc"
1810#include "scr16_c.inc"
1811#include "scr18_c.inc"
1812#include "com01_c.inc"
1813#include "com04_c.inc"
1814#include "com06_c.inc"
1815#include "com08_c.inc"
1819 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT, INACTI,NIN,
1820 . NSV(*), ISKY(*), ICODT(*), NSTRF(*),ICONTACT(*),
1821 . TAGNCONT(NLOADP_HYD_INTER,*),KLOADPINTER(*),LOADPINTER(*),
1822 . LOADP_HYD_INTER(*),
1824 INTEGER ,
INTENT(IN):: INTCAREA
1825 INTEGER,
INTENT(IN) :: NODADT_THERM
1828 . fskyi(lskyi,nfskyi), secfcum(7,numnod,nsect),
1829 . fcont(3,*),ftheskyi(lskyi),condnskyi(lskyi)
1830 my_real,
INTENT(INOUT) :: fsav(*)
1832 TYPE (PARAMETERS_) ,
INTENT(IN):: PARAMETERS
1836 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER,
1837 . nisky_sav,temp_siz,ierror,pp,ppl,np
1842 IF ((nisky+nb)> lskyi)
THEN
1843 CALL ancmsg(msgid=26,anmode=aninfo)
1852 fskyi(nisky,1)=bufr(2,i)
1853 fskyi(nisky,2)=bufr(3,i)
1854 fskyi(nisky,3)=bufr(4,i)
1855 fskyi(nisky,4)=bufr(5,i)
1860 IF(kdtint /= 0 )
THEN
1864 fskyi(nisky,5)=bufr(temp_siz,i)
1869 IF(intth /= 0 )
THEN
1873 ftheskyi(nisky)=bufr(temp_siz,i)
1877 IF(nodadt_therm ==1)
THEN
1881 condnskyi(nisky)=bufr(temp_siz,i)
1891 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
1892 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
1893 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
1894 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
1899 fcont(1,nod)=fcont(1,nod)+bufr(2,i)
1900 fcont(2,nod)=fcont(2,nod)+bufr(3,i)
1901 fcont(3,nod)=fcont(3,nod)+bufr(4,i)
1906 IF(nintloadp > 0)
THEN
1910 DO np = kloadpinter(nin)+1, kloadpinter(nin+1)
1912 ppl = loadp_hyd_inter(pp)
1913 tagncont(ppl,nod) = 1
1932 IF(nstrf(1)+nstrf(2)/=0)
THEN
1934 nbinter=nstrf(k0+14)
1937 IF(nstrf(k1s)==noint)
THEN
1940 n = nint(bufr(1,ii))
1942 IF(secfcum(4,nod,i)==1.)
THEN
1943 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
1944 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
1945 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
1957 IF((ibag/=0.AND.inacti/=7).OR.
1958 . (iadm/=0).OR.(idamp_rdof/=0))
THEN ! warning conflict inacti=7 and ibag=3
1961 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
1962 + bufr(4,i)/=zero)
THEN
1972 ibcs = ibc - 8 * ibcm
1978 CALL ibcoff(ibcs,icodt(nod))
1999 1 NB ,LEN ,BUFR ,NSV ,FSKYI,
2000 2 ISKY ,IBC ,ISECIN ,NOINT ,IBAG ,
2001 3 ICODT ,SECFCUM,NSTRF ,ICONTACT,FCONT,
2002 4 INACTI ,IADM ,INTTH ,FTHESKYI,NLG ,
2013#include "implicit_f.inc"
2017#include "parit_c.inc"
2018#include "scr07_c.inc"
2019#include "scr14_c.inc"
2020#include "scr16_c.inc"
2021#include "scr18_c.inc"
2022#include "com01_c.inc"
2023#include "com04_c.inc"
2024#include "com06_c.inc"
2025#include "com08_c.inc"
2029 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT, INACTI,
2030 . NSV(*), ISKY(*), ICODT(*), NSTRF(*), NLG(*),
2031 . ICONTACT(*), IADM,INTTH
2034 . FSKYI(LSKYI,NFSKYI), SECFCUM(7,NUMNOD,NSECT),
2035 . FCONT(3,*),FTHESKYI(LSKYI)
2036 TYPE(H3D_DATABASE) :: H3D_DATA
2040 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER
2044 IF ((nisky+nb)> lskyi)
THEN
2045 CALL ancmsg(msgid=26,anmode=aninfo)
2049 IF(intth == 0 )
THEN
2055 fskyi(nisky,1)=bufr(2,i)
2056 fskyi(nisky,2)=bufr(3,i)
2057 fskyi(nisky,3)=bufr(4,i)
2058 fskyi(nisky,4)=bufr(5,i)
2066 fskyi(nisky,1)=bufr(2,i)
2067 fskyi(nisky,2)=bufr(3,i)
2068 fskyi(nisky,3)=bufr(4,i)
2069 fskyi(nisky,4)=bufr(5,i)
2070 fskyi(nisky,5)=bufr(6,i)
2083 fskyi(nisky,1)=bufr(2,i)
2084 fskyi(nisky,2)=bufr(3,i)
2085 fskyi(nisky,3)=bufr(4,i)
2086 fskyi(nisky,4)=bufr(5,i)
2087 ftheskyi(nisky) =bufr(6,i)
2095 fskyi(nisky,1)=bufr(2,i)
2096 fskyi(nisky,2)=bufr(3,i)
2097 fskyi(nisky,3)=bufr(4,i)
2098 fskyi(nisky,4)=bufr(5,i)
2099 fskyi(nisky,5)=bufr(6,i)
2100 ftheskyi(nisky) =bufr(7,i)
2108 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT
2109 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
2110 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
2111 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
2116 fcont(1,nod)=fcont(1,nod)+bufr(2,i)
2117 fcont(2,nod)=fcont(2,nod)+bufr(3,i)
2118 fcont(3,nod)=fcont(3,nod)+bufr(4,i)
2125 IF(nstrf(1)+nstrf(2)/=0)
THEN
2127 nbinter=nstrf(k0+14)
2130 IF(nstrf(k1s)==noint)
THEN
2133 n = nint(bufr(1,ii))
2135 IF(secfcum(4,nod,i)==1.)
THEN
2136 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2137 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
2138 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2150 IF((ibag/=0.AND.inacti/=7).OR.
2154 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2155 + bufr(4,i)/=zero)
THEN
2165 ibcs = ibc - 8 * ibcm
2171 CALL ibcoff(ibcs,icodt(nod))
2192 1 NB ,LEN ,BUFR ,IRECTS ,FSKYI,
2193 2 ISKY ,IBC ,ISECIN,NOINT ,IBAG ,
2194 3 ICODT ,SECFCUM,NSTRF ,ICONTACT,FCONT,
2195 4 INTTH ,FTHESKYI,CONDNSKYI,H3D_DATA,
2196 5 NIN ,TAGNCONT,KLOADPINTER,LOADPINTER,
2197 7 LOADP_HYD_INTER,NODADT_THERM)
2207#include "implicit_f.inc"
2211#include "parit_c.inc"
2212#include "scr07_c.inc"
2213#include "scr14_c.inc"
2214#include "scr16_c.inc"
2215#include "scr18_c.inc"
2216#include "com01_c.inc"
2217#include
"com04_c.inc"
2218#include "com06_c.inc"
2219#include "com08_c.inc"
2223 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT,INTTH,NIN,
2224 . irects(2,*), isky(*), icodt(*), nstrf(*),
2226 . tagncont(nloadp_hyd_inter,*),kloadpinter(*),loadpinter(*),
2227 . loadp_hyd_inter(*)
2228 INTEGER,
INTENT(IN) :: NODADT_THERM
2231 . FSKYI(LSKYI,NFSKYI), SECFCUM(7,NUMNOD,NSECT),
2232 . FCONT(3,*),FTHESKYI(*),CONDNSKYI(*)
2233 TYPE(H3D_DATABASE) :: H3D_DATA
2237 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, , NBINTER,PP,PPL,NOD1,NOD2,NP
2241 IF ((NISKY+NB)> LSKYI)THEN
2242 CALL ancmsg(msgid=26,anmode=aninfo)
2246 IF(intth == 0 )
THEN
2252 fskyi(nisky,1)=bufr(2,i)
2253 fskyi(nisky,2)=bufr(3,i)
2254 fskyi(nisky,3)=bufr(4,i)
2255 fskyi(nisky,4)=bufr(5,i)
2259 fskyi(nisky,1)=bufr(6,i)
2260 fskyi(nisky,2)=bufr(7,i)
2261 fskyi(nisky,3)=bufr(8,i)
2262 fskyi(nisky,4)=bufr(9,i)
2270 fskyi(nisky,1)=bufr(2,i)
2271 fskyi(nisky,2)=bufr(3,i)
2272 fskyi(nisky,3)=bufr(4,i)
2273 fskyi(nisky,4)=bufr(5,i)
2274 fskyi(nisky,5)=bufr(6,i)
2278 fskyi(nisky,1)=bufr(7,i)
2279 fskyi(nisky,2)=bufr(8,i)
2280 fskyi(nisky,3)=bufr(9,i)
2281 fskyi(nisky,4)=bufr(10,i)
2282 fskyi(nisky,5)=bufr(11,i)
2287 IF(nodadt_therm == 1)
THEN
2293 fskyi(nisky,1)=bufr(2,i)
2294 fskyi(nisky,2)=bufr(3,i)
2295 fskyi(nisky,3)=bufr(4,i)
2296 fskyi(nisky,4)=bufr(5,i)
2297 ftheskyi(nisky) =bufr(10,i)
2298 condnskyi(nisky)=bufr(12,i)
2302 fskyi(nisky,1)=bufr(6,i)
2303 fskyi(nisky,2)=bufr(7,i)
2304 fskyi(nisky,3)=bufr(8,i)
2305 fskyi(nisky,4)=bufr(9,i)
2306 ftheskyi(nisky) =bufr(11,i)
2307 condnskyi(nisky)=bufr(13,i)
2315 fskyi(nisky,1)=bufr(2,i)
2316 fskyi(nisky,2)=bufr(3,i)
2317 fskyi(nisky,3)=bufr(4,i)
2318 fskyi(nisky,4)=bufr(5,i)
2319 fskyi(nisky,5)=bufr(6,i)
2320 ftheskyi(nisky) =bufr(12,i)
2321 condnskyi(nisky)=bufr(13,i)
2325 fskyi(nisky,1)=bufr(7,i)
2326 fskyi(nisky,2)=bufr(8,i)
2327 fskyi(nisky,3)=bufr(9,i)
2328 fskyi(nisky,4)=bufr(10,i)
2329 fskyi(nisky,5)=bufr(11,i)
2330 ftheskyi(nisky) =bufr(13,i)
2331 condnskyi(nisky)=bufr(14,i)
2341 fskyi(nisky,1)=bufr(2,i)
2342 fskyi(nisky,2)=bufr(3,i)
2343 fskyi(nisky,3)=bufr(4,i)
2344 fskyi(nisky,4)=bufr(5,i)
2345 ftheskyi(nisky) =bufr(10,i)
2349 fskyi(nisky,1)=bufr(6,i)
2350 fskyi(nisky,2)=bufr(7,i)
2351 fskyi(nisky,3)=bufr(8,i)
2352 fskyi(nisky,4)=bufr(9,i)
2353 ftheskyi(nisky) =bufr(11,i)
2361 fskyi(nisky,1)=bufr(2,i)
2362 fskyi(nisky,2)=bufr(3,i)
2363 fskyi(nisky,3)=bufr(4,i)
2364 fskyi(nisky,4)=bufr(5,i)
2365 fskyi(nisky,5)=bufr(6,i)
2366 ftheskyi(nisky) =bufr(12,i)
2370 fskyi(nisky,1)=bufr(7,i)
2371 fskyi(nisky,2)=bufr(8,i)
2372 fskyi(nisky,3)=bufr(9,i)
2373 fskyi(nisky,4)=bufr(10,i)
2374 fskyi(nisky,5)=bufr(11,i)
2375 ftheskyi(nisky) =bufr(13,i)
2384 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
2385 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
2386 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
2387 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
2393 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2394 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2395 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2397 fcont(1,nod)=fcont(1,nod)+ bufr(6,i)
2398 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
2399 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
2405 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2406 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2407 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2409 fcont(1,nod)=fcont(1,nod)+ bufr(7,i)
2410 fcont(2,nod)=fcont(2,nod)+ bufr(8,i)
2411 fcont(3,nod)=fcont(3,nod)+ bufr(9,i)
2417 IF(nintloadp > 0)
THEN
2422 DO np = kloadpinter(nin)+1, kloadpinter(nin+1)
2424 ppl = loadp_hyd_inter(pp)
2425 tagncont(ppl,nod1) = 1
2426 tagncont(ppl,nod2) = 1
2434 IF(nstrf(1)+nstrf(2)/=0)
THEN
2436 nbinter=nstrf(k0+14)
2439 IF(nstrf(k1s)==noint)
THEN
2443 n = nint(bufr(1,ii))
2445 IF(secfcum(4,nod,i)==1.)
THEN
2446 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2447 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
2448 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2451 IF(secfcum(4,nod,i)==1.)
THEN
2452 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(6,ii)
2453 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(7,ii)
2454 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(8,ii)
2461 IF(secfcum(4,nod,i)==1.)
THEN
2462 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2464 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2467 IF(secfcum(4,nod,i)==1.)
THEN
2468 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(7,ii)
2469 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(8,ii)
2470 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(9,ii)
2483 IF((ibag/=0).OR.(idamp_rdof/=0))
THEN
2487 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2488 + bufr(4,i)/=zero)
THEN
2493 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
2494 + bufr(8,i)/=zero)
THEN
2501 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2502 + bufr(4,i)/=zero)
THEN
2507 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
2508 + bufr(9,i)/=zero)
THEN
2532 1 NB ,LEN ,BUFR ,IRECTS ,FSKYI,
2533 2 ISKY ,IBC ,ISECIN,NOINT ,IBAG ,
2534 3 ICODT ,SECFCUM,NSTRF ,ICONTACT,FCONT,
2545#include "implicit_f.inc"
2549#include "parit_c.inc"
2550#include "scr07_c.inc"
2551#include "scr14_c.inc"
2552#include "scr16_c.inc"
2553#include "scr18_c.inc"
2554#include "com01_c.inc"
2555#include "com04_c.inc"
2556#include "com06_c.inc"
2557#include "com08_c.inc"
2561 INTEGER NB, LEN, IBC ,ISECIN ,IBAG , NOINT,
2562 . irects(2,*), isky(*), icodt(*), nstrf(*),
2563 . icontact(*),nlg(*)
2566 . fskyi(lskyi,nfskyi), secfcum(7,numnod,nsect),
2572 INTEGER I, J, II, N, NOD, K0, K1S, IBCS, IBCM, NBINTER
2576 IF ((+NB)> LSKYI)THEN
2577 CALL ANCMSG(MSGID=26,anmode=aninfo)
2584 nod = nlg(irects(1,n))
2586 fskyi(nisky,1)=bufr(2,i)
2587 fskyi(nisky,2)=bufr(3,i)
2588 fskyi(nisky,3)=bufr(4,i)
2589 fskyi(nisky,4)=bufr(5,i)
2591 nod = nlg(irects(2,n))
2593 fskyi(nisky,1)=bufr(6,i)
2594 fskyi(nisky,2)=bufr(7,i)
2595 fskyi(nisky,3)=bufr(8,i)
2596 fskyi(nisky,4)=bufr(9,i)
2602 nod = nlg(irects(1,n))
2604 fskyi(nisky,1)=bufr(2,i)
2605 fskyi(nisky,2)=bufr(3,i)
2606 fskyi(nisky,3)=bufr(4,i)
2607 fskyi(nisky,4)=bufr(5,i)
2608 fskyi(nisky,5)=bufr(6,i)
2610 nod = nlg(irects(2,n))
2612 fskyi(nisky,1)=bufr(7,i)
2613 fskyi(nisky,2)=bufr(8,i)
2614 fskyi(nisky,3)=bufr(9,i)
2615 fskyi(nisky,4)=bufr(10,i)
2616 fskyi(nisky,5)=bufr(11,i)
2623 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
2624 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
2625 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
2626 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
2631 nod = nlg(irects(1,n))
2632 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2633 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2634 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2635 nod = nlg(irects(2,n))
2636 fcont(1,nod)=fcont(1,nod)+ bufr(6,i)
2637 fcont(2,nod)=fcont(2,nod)+ bufr(7,i)
2638 fcont(3,nod)=fcont(3,nod)+ bufr(8,i)
2643 nod = nlg(irects(1,n))
2644 fcont(1,nod)=fcont(1,nod)+ bufr(2,i)
2645 fcont(2,nod)=fcont(2,nod)+ bufr(3,i)
2646 fcont(3,nod)=fcont(3,nod)+ bufr(4,i)
2647 nod = nlg(irects(2,n))
2648 fcont(1,nod)=fcont(1,nod)+ bufr(7,i)
2649 fcont(2,nod)=fcont(2,nod)+ bufr(8,i)
2650 fcont(3,nod)=fcont(3,nod)+ bufr(9,i)
2658 IF(nstrf(1)+nstrf(2)/=0)
THEN
2660 nbinter=nstrf(k0+14)
2663 IF(nstrf(k1s)==noint)
THEN
2667 n = nint(bufr(1,ii))
2668 nod = nlg(irects(1,n))
2669 IF(secfcum(4,nod,i)==1.)
THEN
2670 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2671 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
2672 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2674 nod = nlg(irects(2,n))
2675 IF(secfcum(4,nod,i)==1.)
THEN
2676 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(6,ii)
2677 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(7,ii)
2678 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(8,ii)
2683 n = nint(bufr(1,ii))
2684 nod = nlg(irects(1,n))
2685 IF(secfcum(4,nod,i)==1.)
THEN
2686 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(2,ii)
2687 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(3,ii)
2688 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(4,ii)
2690 nod = nlg(irects(2,n))
2691 IF(secfcum(4,nod,i)==1.)
THEN
2692 secfcum(1,nod,i)=secfcum(1,nod,i)+bufr(7,ii)
2693 secfcum(2,nod,i)=secfcum(2,nod,i)+bufr(8,ii)
2694 secfcum(3,nod,i)=secfcum(3,nod,i)+bufr(9,ii)
2711 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2712 + bufr(4,i)/=zero)
THEN
2714 nod = nlg(irects(1,n))
2717 IF(bufr(6,i)/=zero.OR.bufr(7,i)/=zero.OR.
2718 + bufr(8,i)/=zero)
THEN
2719 nod = nlg(irects(2,n))
2725 IF(bufr(2,i)/=zero.OR.bufr(3,i)/=zero.OR.
2726 + bufr(4,i)/=zero)
THEN
2728 nod = nlg(irects(1,n))
2731 IF(bufr(7,i)/=zero.OR.bufr(8,i)/=zero.OR.
2732 + bufr(9,i)/=zero)
THEN
2733 nod = nlg(irects(2,n))
2743!||====================================================================
2752 1 NB ,LEN ,BUFR ,NELEMS ,FSKYI ,
2753 2 ISKY ,FCONT ,IXS ,IXS16 ,H3D_DATA )
2762#include "implicit_f.inc"
2766#include "parit_c.inc"
2767#include "scr07_c.inc"
2768#include "scr14_c.inc"
2769#include "scr16_c.inc"
2770#include "com04_c.inc"
2771#include "com06_c.inc"
2772#include "com08_c.inc"
2777 . nelems(*), isky(*), ixs(nixs,*), ixs16(8,*)
2779 . fskyi(lskyi,nfskyi),
2785 INTEGER I, J, II, IIIS, NOD, IES, NN
2790 nn = nint(bufr(1,i))
2793 iiis = nint(bufr(6+(ii-1)*5,i))
2795 nod = ixs(iiis+1,ies)
2797 nod = ixs16(iiis-8,ies-numels8-numels10-numels20)
2800 fskyi(nisky,1)=bufr(2+(ii-1)*5,i)
2801 fskyi(nisky,2)=bufr(3+(ii-1)*5,i)
2802 fskyi(nisky,3)=bufr(4+(ii-1)*5,i)
2803 fskyi(nisky,4)=bufr(5+(ii-1)*5,i)
2810 IF((anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0.AND.
2811 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
2812 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))
2813 . .OR.anim_v(26)+h3d_data%N_VECT_CONT_MAX>0)
THEN
2816 ies = nelems(nint(bufr(1,i)))
2818 iiis = nint(bufr(6+(ii-1)*5,i))
2821 nod = ixs(iiis+1,ies)
2823 nod = ixs16(iiis-8,ies-numels8-numels10-numels20)
2825 fcont(1,nod)=fcont(1,nod)+ bufr(2+(ii-1)*5,i)
2826 fcont(2,nod)=fcont(2,nod)+ bufr(3+(ii-1)*5,i)
2827 fcont(3,nod)=fcont(3,nod)+ bufr(4+(ii-1)*5,i)
2858 1 IPARI ,ISENDTO ,IRCVFROM,INTLIST ,NBINTC ,
2859 2 ISIZXV ,ILENXV ,IAD_ELEM,I2SIZE ,ITASK ,
2860 3 ISLEN7 ,IRLEN7 ,ISLEN11 ,IRLEN11 ,IGRBRIC ,
2861 4 NME17 ,ISLEN17 ,IRLEN17 ,IRLEN7T ,ISLEN7T ,
2862 5 LINDIDEL,LBUFIDEL,IRLEN20 ,ISLEN20 ,IRLEN20T,
2863 6 ISLEN20T,NBINT20 ,IRLEN20E,ISLEN20E,FR_RBY ,
2864 7 FR_RBY6 ,NPBY ,IRBKIN_L,NRBYKIN_L,KINDRBY,
2865 8 NSENSOR ,SENSOR_TAB,LBUFIDEL24,INTBUF_TAB,
2866 9 SORT_COMM,NEED_COMM_INT25_SOLID_EROSION,COMM_INT25_SOLID_EROSION )
2882#include "implicit_f.inc"
2886#include "com01_c.inc"
2887#include "com04_c.inc"
2888#include "com08_c.inc"
2889#include "scr17_c.inc"
2890#include "scr18_c.inc"
2891#include "task_c.inc"
2892#include "param_c.inc"
2893#include "units_c.inc"
2894#include "warn_c.inc"
2895#include "tabsiz_c.inc"
2897 COMMON /ring/irecvf,isendt,iring
2898 INTEGER IRECVF,ISENDT,IRING
2902 INTEGER ,
INTENT(IN) :: NSENSOR
2903 INTEGER ISENDTO(NINTER+1,NSPMD+1),
2904 . IRCVFROM(NINTER+1,NSPMD+1), IPARI(NPARI,*),
2905 . INTLIST(NINTER), IAD_ELEM(2,*),
2906 . NBINTC, ISIZXV, ILENXV, I2SIZE,ITASK,
2907 . ISLEN7, IRLEN7, ISLEN11, IRLEN11, NME17, ISLEN17, IRLEN17,
2908 . IRLEN7T, ISLEN7T, LINDIDEL, LBUFIDEL,
2909 . IRLEN20,ISLEN20,IRLEN20T,ISLEN20T,NBINT20,
2910 . IRLEN20E, ISLEN20E, NRBYKIN_L, IRBKIN_L(*),KINDRBY(*),
2911 . FR_RBY(*), FR_RBY6(*), NPBY(NNPBY,*),LBUFIDEL24
2912 TYPE(INTBUF_STRUCT_) INTBUF_TAB(*)
2913 TYPE(sorting_comm_type),
DIMENSION(NINTER),
INTENT(inout) :: SORT_COMM
2914 TYPE (SENSOR_STR_) ,
DIMENSION(NSENSOR) :: SENSOR_TAB
2915 LOGICAL,
INTENT(inout) :: NEED_COMM_INT25_SOLID_EROSION
2916 INTEGER,
INTENT(inout) :: COMM_INT25_SOLID_EROSION
2918 TYPE (GROUP_) ,
DIMENSION(NGRBRIC) :: IGRBRIC
2922 INTEGER I, J, K, NSN, NMN, ITYP, IERROR, IERROR1, KK, NIN, P,
2923 . scpmax1,isd,irc,my_rank,last, ige, ign, nme, nmes,intth,
2924 . nrtm, nrts, nlinm, nlins,nsne, nmne,n,m,
2925 . work(70000),ity,noint,inacti,isens,interact
2926 INTEGER,
DIMENSION(:),
ALLOCATABLE:: INDX
2927 INTEGER,
DIMENSION(:),
ALLOCATABLE:: PPAR
2928 INTEGER,
DIMENSION(:),
ALLOCATABLE:: INDP
2929 INTEGER,
DIMENSION(:),
ALLOCATABLE:: NNP
2930 INTEGER,
DIMENSION(:),
ALLOCATABLE:: RANK
2931 INTEGER,
DIMENSION(:),
ALLOCATABLE:: IRBODY
2932 INTEGER,
DIMENSION(:),
ALLOCATABLE:: INDB
2933 INTEGER,
DIMENSION(:),
ALLOCATABLE:: NNSN
2934 my_real STARTT, STOPT,TS
2938 CALL my_alloc(indx,2*ninter)
2939 CALL my_alloc(ppar,ninter)
2940 CALL my_alloc(indp,2*nspmd)
2941 CALL my_alloc(nnp,nspmd)
2942 CALL my_alloc(rank,nspmd)
2943 CALL my_alloc(irbody,numnod)
2944 CALL my_alloc(indb,2*nrbykin)
2945 CALL my_alloc(nnsn,nrbykin)
2971 IF(ity==7.OR.ity==10.OR.ity==11.OR.
2972 . ity==17.OR.ity==20.OR.
2973 . ity==22.OR.ity==23.OR.ity==24.OR.ity==25)
THEN
2975 IF (inacti==5.OR.inacti==6.OR.inacti==7.OR.ity==23.OR.inacti==-1)
THEN
2979 IF(ity == 7.OR.ity == 11.OR.ity == 24.OR.ity == 25)
2980 . isens = ipari(64,i)
2982 ts = sensor_tab(isens)%TSTART
2983 IF (ts>tt) interact = 0
2985 startt = intbuf_tab(i)%VARIABLES(3)
2986 stopt = intbuf_tab(i)%VARIABLES(11)
2987 IF (startt>tt) interact = 0
2990 IF (interact==0)
THEN
2994 CALL ancmsg(msgid=220,anmode=aninfo,
2995 * i1=noint,i2=inacti,r1=startt,r2=startt)
2998 IF (interact==0 .AND. ity /=25)
THEN
3002 nsnfi(i)%P(1:nspmd)=0
3004 nsnsi(i)%P(1:nspmd)=0
3017 IF(ipari(17,i)>0)
THEN
3022 IF(ityp==7.OR.ityp==22.OR.
3023 . ityp==10.OR.ityp==5.OR.
3024 . ityp==23.OR.ityp==24.OR.
3026 lindidel = lindidel+nrtm
3027 lbufidel = lbufidel+4*nrtm+4
3028 IF(ityp==24.OR.ityp==25) lbufidel24 = lbufidel24
3029 ELSEIF(ityp==11)
THEN
3030 lindidel = lindidel+nrtm+nrts
3031 lbufidel = lbufidel+2*(nrtm+nrts)+4
3033 lindidel = lindidel+nrtm+nrts
3034 lbufidel = lbufidel+4*(nrtm+nrts)+4
3036 lindidel = lindidel+nsn
3037 lbufidel = lbufidel+4*nsn+4
3038 ELSEIF(ityp==20)
THEN
3039 lindidel = lindidel+nrtm
3040 lbufidel = lbufidel+4*nrtm+4
3044 lindidel = lindidel+nlinm+nlins
3045 lbufidel = lbufidel+2*(nlinm+nlins)+4
3055 IF (ityp==7.OR.ityp==10.OR.ityp==22.OR.
3056 . ityp==11.OR.ityp==23.OR.ityp==24.OR.
3061 isendto(i,ispmd+1)=nsn
3062 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3063 isendto(ninter+1,ispmd+1)=isendto(ninter+1,ispmd+1)+nsn
3066 ircvfrom(i,ispmd+1) = nmn
3067 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3068 ircvfrom(ninter+1,ispmd+1)=ircvfrom(ninter+1,ispmd+1)
3073 ELSEIF(ityp==17)
THEN
3074 IF(ipari(33,i)==0)
THEN
3077 nmes= igrbric(ign)%NENTITY
3078 nme = igrbric(ige)%NENTITY
3079 nme17 = nme17+nme+nmes
3081 isendto(i,ispmd+1)=nmes
3082 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3083 isendto(ninter+1,ispmd+1)=isendto(ninter+1,ispmd+1)
3087 ircvfrom(i,ispmd+1) = nme
3088 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3089 ircvfrom(ninter+1,ispmd+1)=ircvfrom(ninter+1,ispmd+1)
3093 ELSEIF(ityp==20)
THEN
3099 IF(nsn+nsne/=0)
THEN
3100 isendto(i,ispmd+1)=nsn+nsne
3101 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3102 isendto(ninter+1,ispmd+1)=isendto(ninter+1,ispmd+1)
3105 IF(nmn+nmne/=0)
THEN
3106 ircvfrom(i,ispmd+1) = nmn+nmne
3107 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3108 ircvfrom(ninter+1,ispmd+1)=ircvfrom(ninter+1,ispmd+1)
3116 CALL spmd_ibcast(isendto(1,k),isendto(1,k),ninter+1,1,
3118 CALL spmd_ibcast(ircvfrom(1,k),ircvfrom(1,k),ninter+1,1,
3123 isendto(i,nspmd+1)=0
3124 ircvfrom(i,nspmd+1)=0
3126 IF(ircvfrom(i,k)/=0)
THEN
3127 ircvfrom(i,nspmd+1)=ircvfrom(i,nspmd+1)+1
3129 IF(isendto(i,k)/=0)
THEN
3130 isendto(i,nspmd+1)=isendto(i,nspmd+1)+1
3140 IF (ityp==7.OR.ityp==10.OR.ityp==11.OR.
3141 . (ityp==17.AND.ipari(33,i) == 0).OR.ityp==20.OR.
3142 . ityp==22.OR.ityp==23.OR.ityp==24.OR.ityp==25)
THEN
3147 scpmax1=
max(scpmax1,ircvfrom(i,j))
3149 ppar(nbintc)=scpmax1
3153 CALL my_orders(0,work,ppar,indx,nbintc,1)
3155 indx(i+nbintc)=intlist(i)
3158 intlist(i)=indx(indx(i)+nbintc)
3161 IF(debug(3)>=1.AND.ispmd==0.AND.nspmd>1)
THEN
3162 WRITE(istdo,*)
'** INTERFACES NODAL DECOMPOSITION '
3163 WRITE(istdo,*)
'#PROC NSN TOT NMN TOT TOTAL'
3164 WRITE(iout,*)
'** INTERFACES NODAL DECOMPOSITION '
3165 WRITE(iout,*)
'#PROC NSN TOT NMN TOT TOTAL'
3170 isd = isd + isendto(i,j)
3171 irc = irc + ircvfrom(i,j)
3173 WRITE(istdo,
'(I4,3X,I8,2X,I8,2X,I8)')j,isd,irc,isd+irc
3174 WRITE(iout,
'(I4,3X,I8,2X,I8,2X,I8)')j,isd,irc,isd+irc
3181 nnp(k) = ircvfrom(ninter+1,k) + isendto(ninter+1,k)
3184 IF(nspmd > 1)
CALL my_orders(0,work,nnp,indp,nspmd,1)
3191 IF(nnp(ispmd+1)>0)
THEN
3192 my_rank = rank(ispmd+1)
3194 last = indp(my_rank-1)
3199 IF(my_rank==nspmd)
THEN
3200 IF(irecvf/=0) isendt = -1
3202 isendt = indp(my_rank+1)
3205 IF(nnp(indp(nspmd-1))>0)iring= indp(nspmd)
3213 isizxv = isizxv + iad_elem(1,i+1)-iad_elem(2,i)
3216 ilenxv = 10 + 3*iroddl
3218 IF(iale/=0.AND.ninter>0) ilenxv = ilenxv + 6
3220 IF(idtmin(11)==3.OR.idtmin(10)==3) ilenxv = ilenxv + 2
3255 intth = ipari(47,nin)
3257 IF(ityp==7 .OR.ityp==10.OR.
3258 . ityp==23.OR.ityp==22.OR.ityp==24.OR.
3262 islen7 = islen7 +
nsnsi(nin)%P(p)
3263 irlen7 = irlen7 +
nsnfi(nin)%P(p)
3267 islen7t = islen7t +
nsnsi(nin)%P(p)
3268 irlen7t = irlen7t +
nsnfi(nin)%P(p)
3272 IF(ityp == 25 .AND. ipari(58,nin) >0)
THEN
3277 ELSEIF(ityp==11)
THEN
3280 islen11 = islen11 +
nsnsi(nin)%P(p)
3281 irlen11 = irlen11 +
nsnfi(nin)%P(p)
3283 ELSEIF(ityp==17.AND.ipari(33,nin)==0)
THEN
3286 islen17 = islen17 +
nsnsi(nin)%P(p)
3287 irlen17 = irlen17 +
nsnfi(nin)%P(p)
3289 ELSEIF(ityp==20)
THEN
3291 nbint20 = nbint20 + 1
3294 islen20 = islen20 +
nsnsi(nin)%P(p)
3295 irlen20 = irlen20 +
nsnfi(nin)%P(p)
3296 islen20e= islen20e+
nsnsie(nin)%P(p)
3297 irlen20e= irlen20e+
nsnfie(nin)%P(p)
3301 islen20t = islen20t +
nsnsi(nin)%P(p)
3302 irlen20t = irlen20t +
nsnfi(nin)%P(p)
3303 islen20e = islen20e +
nsnsie(nin)%P(p)
3304 irlen20e = irlen20e +
nsnfie(nin)%P(p)
3310 IF(nrbykin > 0)
THEN
3324 nrbykin_l=nrbykin_l+1
3325 nnsn(nrbykin_l)=-npby(2,n)
3326 irbkin_l(nrbykin_l)=n
3335 fr_rby6(n)=irbody(fr_rby(n))
3338 CALL my_orders(0,work,nnsn,indb,nrbykin_l,1)
3341 indb(nrbykin_l+n)=irbkin_l(n)
3344 irbkin_l(n)=indb(nrbykin_l+indb(n))
3350 indx(i+nbintc)=intlist(i)
3353 intlist(i)=indx(indx(i)+nbintc)
3359 . ircvfrom(ninter+1,ispmd+1)+isendto(ninter+1,ispmd+1),comm_cont)
3362 . need_comm_int25_solid_erosion,comm_int25_solid_erosion )
3379!||
spmd_initfi ../engine/source/mpi/interfaces/spmd_i7tool.f
3382!||--- calls -----------------------------------------------------
3397 SUBROUTINE spmd_initfi(IPARI,IFLAG,H3D_DATA,PARAMETERS,IDT_THERM,INTHEAT)
3411#include "implicit_f.inc"
3412#include "i25edge_c.inc"
3416#include "com01_c.inc"
3417#include "com04_c.inc"
3418#include "scr14_c.inc"
3419#include "scr16_c.inc"
3420#include "scr18_c.inc"
3421#include "task_c.inc"
3422#include "param_c.inc"
3423#include "parit_c.inc"
3424#include "spmd_c.inc"
3429 INTEGER IPARI(NPARI,*), IFLAG
3430 INTEGER,
INTENT(IN) :: IDT_THERM
3431 INTEGER,
INTENT(IN) :: INTHEAT
3432 TYPE(H3D_DATABASE) :: H3D_DATA
3433 TYPE (PARAMETERS_) ,
INTENT(IN):: PARAMETERS
3437 INTEGER I, P, NSN, NMN, ITYP, IGAP, IERROR, IERROR1, LENS, LENR,
3438 . lskyfi, inacti, nbintc, leni, j,intth, k, l,iedge4 ,intfric ,
3439 . flagremn ,sizremnorfi, ivis2 ,intnitsche,itied,ipstif
3440 INTEGER :: LENR_EDGE,LENS_EDGE
3448 ALLOCATE(
nsvfi(ninter),stat=ierror1)
3449 ierror = ierror + ierror1
3450 ALLOCATE(
nsnfi(ninter),stat=ierror1)
3451 ierror = ierror + ierror1
3452 ALLOCATE(
edge_fi(ninter),stat=ierror1)
3453 ierror = ierror + ierror1
3454 ALLOCATE(
nsnfi_sav(ninter),stat=ierror1)
3455 ierror = ierror + ierror1
3456 ALLOCATE(
nsnsi_sav(ninter),stat=ierror1)
3457 ierror = ierror + ierror1
3458 ALLOCATE(
itafi(ninter),stat=ierror1)
3459 ierror = ierror + ierror1
3460 ALLOCATE(
pmainfi(ninter),stat=ierror1)
3461 ierror = ierror + ierror1
3462 ALLOCATE(
kinfi(ninter),stat=ierror1)
3463 ierror = ierror + ierror1
3464 ALLOCATE(
msfi(ninter),stat=ierror1)
3465 ierror = ierror + ierror1
3466 ALLOCATE(
stifi(ninter),stat=ierror1)
3467 ierror = ierror + ierror1
3468 ALLOCATE(
gapfi(ninter),stat=ierror1)
3469 ierror = ierror + ierror1
3470 ALLOCATE(
gap_lfi(ninter),stat=ierror1)
3471 ierror = ierror + ierror1
3472 ALLOCATE(
xfi(ninter),stat=ierror1)
3473 ierror = ierror + ierror1
3474 ALLOCATE(
vfi(ninter),stat=ierror1)
3475 ierror = ierror + ierror1
3476 ALLOCATE(
nsvsi(ninter),stat=ierror1)
3477 ierror = ierror + ierror1
3478 ALLOCATE(
nsnsi(ninter),stat=ierror1)
3479 ierror = ierror + ierror1
3480 ALLOCATE(
afi(ninter),stat=ierror1)
3481 ierror = ierror + ierror1
3482 ALLOCATE(
stnfi(ninter),stat=ierror1)
3483 ierror = ierror + ierror1
3484 ALLOCATE(
vscfi(ninter),stat=ierror1)
3485 ierror = ierror + ierror1
3486 ALLOCATE(
penfi(ninter),stat=ierror1)
3487 ierror = ierror + ierror1
3488 ALLOCATE(
fskyfi(ninter),stat=ierror1)
3489 ierror = ierror + ierror1
3490 ALLOCATE(
iskyfi(ninter),stat=ierror1)
3491 ierror = ierror + ierror1
3493 ALLOCATE(
nisubsfi(ninter),stat=ierror1)
3494 ierror = ierror + ierror1
3495 ALLOCATE(
lisubsfi(ninter),stat=ierror1)
3496 ierror = ierror + ierror1
3498 ierror = ierror + ierror1
3499 ALLOCATE(
addsubsfi(ninter),stat=ierror1)
3500 ierror = ierror + ierror1
3501 IF(parameters%INTCAREA > 0 )
THEN
3503 ierror = ierror + ierror1
3506 ALLOCATE(
nisubsfie(ninter),stat=ierror1)
3507 ierror = ierror + ierror1
3508 ALLOCATE(
lisubsfie(ninter),stat=ierror1)
3509 ierror = ierror + ierror1
3511 ierror = ierror + ierror1
3513 ierror = ierror + ierror1
3516 ALLOCATE(
fnconti(ninter),stat=ierror1)
3517 ierror = ierror + ierror1
3518 ALLOCATE(
ftconti(ninter),stat=ierror1)
3519 ierror = ierror + ierror1
3521 ALLOCATE (
nlskyfi(ninter),stat=ierror1)
3522 ierror = ierror + ierror1
3525 ALLOCATE (
eminxfi(ninter),stat=ierror1)
3526 ierror = ierror + ierror1
3527 ALLOCATE (
ksfi(ninter),stat=ierror1)
3528 ierror = ierror + ierror1
3529 ALLOCATE (
frotsfi(ninter),stat
3530 ierror = ierror + ierror1
3531 ALLOCATE (
stnfi17(ninter),stat=ierror1)
3532 ierror = ierror + ierror1
3533 ALLOCATE(
xfi17(ninter),stat=ierror1)
3534 ierror = ierror + ierror1
3535 ALLOCATE(
vfi17(ninter),stat=ierror1)
3536 ierror = ierror + ierror1
3537 ALLOCATE(
afi17(ninter),stat=ierror1)
3538 ierror = ierror + ierror1
3540 ALLOCATE(
nbinflfi(ninter),stat=ierror1)
3541 ierror = ierror + ierror1
3542 ALLOCATE(
daanc6fi(ninter),stat=ierror1)
3543 ierror = ierror + ierror1
3544 ALLOCATE(
dxancfi(ninter),stat=ierror1)
3545 ierror = ierror + ierror1
3546 ALLOCATE(
dvancfi(ninter),stat=ierror1)
3547 ierror = ierror + ierror1
3548 ALLOCATE(
penfia(ninter),stat=ierror1)
3549 ierror = ierror + ierror1
3550 ALLOCATE(
alphakfi(ninter),stat=ierror1)
3551 ierror = ierror + ierror1
3553 ierror = ierror + ierror1
3555 ierror = ierror + ierror1
3557 ALLOCATE(
nsvfie(ninter),stat=ierror1)
3558 ierror = ierror + ierror1
3559 ALLOCATE(
nsnfie(ninter),stat=ierror1)
3560 ierror = ierror + ierror1
3561 ALLOCATE(
itafie(ninter),stat=ierror1)
3562 ierror = ierror + ierror1
3563 ALLOCATE(
msfie(ninter),stat=ierror1)
3564 ierror = ierror + ierror1
3565 ALLOCATE(
stifie(ninter),stat=ierror1)
3566 ierror = ierror + ierror1
3567 ALLOCATE(
gapfie(ninter),stat=ierror1)
3568 ierror = ierror + ierror1
3569 ALLOCATE(
xfie(ninter),stat=ierror1)
3570 ierror = ierror + ierror1
3571 ALLOCATE(
vfie(ninter),stat=ierror1)
3572 ierror = ierror + ierror1
3573 ALLOCATE(
nsvsie(ninter),stat=ierror1)
3574 ierror = ierror + ierror1
3575 ALLOCATE(
nsnsie(ninter),stat=ierror1)
3576 ierror = ierror + ierror1
3577 ALLOCATE(
afie(ninter),stat=ierror1)
3578 ierror = ierror + ierror1
3579 ALLOCATE(
stnfie(ninter),stat=ierror1)
3580 ierror = ierror + ierror1
3581 ALLOCATE(
vscfie(ninter),stat=ierror1)
3582 ierror = ierror + ierror1
3583 ALLOCATE(
penfie(ninter),stat=ierror1)
3584 ierror = ierror + ierror1
3585 ALLOCATE(
fskyfie(ninter),stat=ierror1)
3586 ierror = ierror + ierror1
3587 ALLOCATE(
iskyfie(ninter),stat=ierror1)
3588 ierror = ierror + ierror1
3590 ALLOCATE (
nlskyfie(ninter),stat=ierror1)
3591 ierror = ierror + ierror1
3594 ALLOCATE(
daanc6fie(ninter),stat=ierror1)
3595 ierror = ierror + ierror1
3596 ALLOCATE(
dxancfie(ninter),stat=ierror1)
3597 ierror = ierror + ierror1
3598 ALLOCATE(
dvancfie(ninter),stat=ierror1)
3600 ALLOCATE(
penfiae(ninter),stat=ierror1)
3601 ierror = ierror + ierror1
3602 ALLOCATE(
alphakfie(ninter),stat=ierror1)
3603 ierror = ierror + ierror1
3604 ALLOCATE(
daancfie(ninter),stat=ierror1)
3605 ierror = ierror + ierror1
3607 ierror = ierror + ierror1
3609 ALLOCATE (
nodnxfie(ninter),stat=ierror1)
3610 ierror = ierror + ierror1
3611 ALLOCATE (
nodamsfie(ninter),stat=ierror1)
3612 ierror = ierror + ierror1
3614 ierror = ierror + ierror1
3618 ALLOCATE (
mtfi_pene(ninter),stat=ierror1)
3619 ierror = ierror + ierror1
3621 ierror = ierror + ierror1
3622 ALLOCATE (
mtfi_v(ninter),stat=ierror1)
3623 ierror = ierror + ierror1
3624 ALLOCATE (
mtfi_a(ninter),stat=ierror1)
3625 ierror = ierror + ierror1
3626 ALLOCATE (
i18kafi(ninter),stat=ierror1)
3627 ierror = ierror + ierror1
3628 ALLOCATE (
mtfi_n(ninter),stat=ierror1)
3629 ierror = ierror + ierror1
3631 ALLOCATE (
time_sfi(ninter),stat=ierror1)
3632 ierror = ierror + ierror1
3633 ALLOCATE (
irtlm_fi(ninter),stat=ierror1)
3634 ierror = ierror + ierror1
3636 ierror = ierror + ierror1
3638 ierror = ierror + ierror1
3640 ierror = ierror + ierror1
3642 ierror = ierror + ierror1
3643 ALLOCATE (
isedge_fi(ninter),stat=ierror1)
3644 ierror = ierror + ierror1
3645 ALLOCATE (
irtse_fi(ninter),stat=ierror1)
3646 ierror = ierror + ierror1
3647 ALLOCATE (
is2pt_fi(ninter),stat=ierror1)
3648 ierror = ierror + ierror1
3649 ALLOCATE (
isegpt_fi(ninter),stat=ierror1)
3650 ierror = ierror + ierror1
3651 ALLOCATE (
is2se_fi(ninter),stat=ierror1)
3652 ierror = ierror + ierror1
3653 ALLOCATE (
ispt2_fi(ninter),stat=ierror1)
3654 ierror = ierror + ierror1
3655 ALLOCATE (
remnor_fi(ninter),stat=ierror1)
3656 ierror = ierror + ierror1
3658 ierror = ierror + ierror1
3660 ALLOCATE (
islide_fi(ninter),stat=ierror1)
3661 ierror = ierror + ierror1
3662 ALLOCATE (
icodt_fi(ninter),stat=ierror1)
3663 ierror = ierror + ierror1
3664 ALLOCATE (
iskew_fi(ninter),stat=ierror1)
3665 ierror = ierror + ierror1
3668 IF(ipari(97,i) > 0) ipstif = ipstif + ipari(97,i)
3670 IF(parameters%ISTIF_DT > 0 .OR. ipstif>0 )
THEN
3675 IF(intheat > 0 )
THEN
3676 ALLOCATE (
ftheskyfi(ninter),stat=ierror1)
3677 ierror = ierror + ierror1
3678 ALLOCATE(
fthefi(ninter),stat=ierror1)
3679 ierror = ierror + ierror1
3680 ALLOCATE(
tempfi(ninter),stat=ierror1)
3681 ierror = ierror + ierror1
3682 ALLOCATE(
matsfi(ninter),stat=ierror1)
3683 ierror = ierror + ierror1
3684 ALLOCATE(
nmtemp(ninter),stat=ierror1)
3685 ierror = ierror + ierror1
3686 IF(idt_therm == 1)
THEN
3687 ALLOCATE(
condnfi(ninter),stat=ierror1)
3688 ierror = ierror + ierror1
3690 ierror = ierror + ierror1
3693 IF(intheat > 0 .OR.nintloadp21 > 0)
THEN
3694 ALLOCATE(
nmnfi(ninter),stat=ierror1)
3695 ierror = ierror + ierror1
3696 ALLOCATE(
nmvfi(ninter),stat=ierror1)
3697 ierror = ierror + ierror1
3698 ALLOCATE(
nmnsi(ninter),stat=ierror1)
3699 ierror = ierror + ierror1
3700 ALLOCATE(
nmvsi(ninter),stat=ierror1)
3701 ierror = ierror + ierror1
3702 ALLOCATE(
tempnod(ninter),stat=ierror1)
3703 ierror = ierror + ierror1
3705 IF(nintloadp21 > 0)
THEN
3707 ierror = ierror + ierror1
3709 IF(ninterfric > 0 )
THEN
3711 ierror = ierror + ierror1
3713 ierror = ierror + ierror1
3716 ALLOCATE (
nodnxfi(ninter),stat=ierror1)
3717 ierror = ierror + ierror1
3718 ALLOCATE (
nodamsfi(ninter),stat=ierror1)
3719 ierror = ierror + ierror1
3720 ALLOCATE (
procamsfi(ninter),stat=ierror1)
3721 ierror = ierror + ierror1
3724 ierror = ierror + ierror1
3726 ierror = ierror + ierror1
3728 IF(intheat > 0.OR.interadhesion > 0 )
THEN
3729 ALLOCATE(
areasfi(ninter),stat=ierror1)
3733 ierror = ierror + ierror1
3736 IF(interadhesion > 0)
THEN
3737 ALLOCATE(
if_adhfi(ninter),stat=ierror1)
3741 ierror = ierror + ierror1
3744 ALLOCATE(
candf_si(ninter),stat=ierror1)
3745 ierror = ierror + ierror1
3747 IF(nitsche > 0 )
THEN
3748 ALLOCATE(
forneqsfi(ninter),stat=ierror1)
3749 ierror = ierror + ierror1
3753 ALLOCATE(
efricfi(ninter),stat=ierror1)
3754 ierror = ierror + ierror1
3755 ALLOCATE(
efricgfi(ninter),stat=ierror1)
3756 ierror = ierror + ierror1
3795 IF(parameters%INTCAREA > 0)
NULLIFY(
intareanfi(i)%p)
3883 IF(parameters%ISTIF_DT > 0 .OR. ipari(97,i)>0 )
THEN
3899 IF(intheat /= 0)
THEN
3908 IF(interadhesion /= 0)
THEN
3913 IF(idt_therm == 1)
THEN
3920 IF(intheat /= 0.OR.nintloadp21 > 0)
THEN
3929 IF(nintloadp21 > 0)
THEN
3935 IF(ninterfric > 0 )
THEN
3942 IF(nitsche > 0 )
THEN
3949 CALL ancmsg(msgid=20,anmode=aninfo)
3960 inacti = ipari(22,i)
3962 intfric = ipari(72,i)
3964 intnitsche = ipari(86,i)
3969 iedge4 = ipari(59,i)
3973 IF(ityp==7.OR.ityp==10.OR.ityp==11.OR.
3974 + (ityp==17.AND.ipari(33,i)==0).OR.ityp==20.OR.
3975 + ityp==22.OR.ityp==23.OR.ityp==24.OR.
3979 ALLOCATE(
nsnsi(i)%P(nspmd),stat=ierror)
3980 ALLOCATE(
nsnfi(i)%P(nspmd),stat=ierror1)
3983 IF(ierror+ierror1/=0)
THEN
3984 CALL ancmsg(msgid=20,anmode=aninfo)
3992 lens = lens +
nsnsi(i)%P(p)
3993 lenr = lenr +
nsnfi(i)%P(p)
4000 ALLOCATE(
nsvsi(i)%P(lens),stat=ierror1)
4002 ierror = ierror + ierror1
4004 ALLOCATE(
candf_si(i)%P(nsn),stat=ierror1)
4008 ALLOCATE(
nsvfi(i)%P(lenr),stat=ierror1)
4009 ierror = ierror + ierror1
4011 IF(ityp==7.OR.ityp==10.OR.
4012 + ityp==20.OR.ityp==22.OR.
4013 + ityp==23.OR.ityp==24.OR.
4015 ALLOCATE(
itafi(i)%P(lenr),stat=ierror1)
4016 ierror = ierror + ierror1
4018 IF(ityp==7.OR.ityp==20.OR.ityp==22.OR.
4019 + ityp==23.OR.ityp==24.OR.ityp==25)
THEN
4020 ALLOCATE(
kinfi(i)%P(lenr),stat=ierror1)
4021 ierror = ierror + ierror1
4025 ALLOCATE(
nbinflfi(i)%P(lenr),stat=ierror1)
4026 ierror = ierror + ierror1
4031 ALLOCATE(
matsfi(i)%P(lenr),stat=ierror1)
4032 ierror = ierror + ierror
4035 IF((ityp==7.OR.ityp==24.OR.ityp==25).AND.intfric > 0)
THEN
4037 ierror = ierror + ierror
4043 *
ALLOCATE(
irtlm_fi(i)%P(2,lenr),stat=ierror1)
4047 *
ALLOCATE(
icont_i_fi(i)%P(lenr),stat=ierror1)
4051 *
ALLOCATE(
isedge_fi(i)%P(lenr),stat=ierror1)
4056 *
ALLOCATE(
irtse_fi(i)%P(5,lenr),stat=ierror1)
4060 *
ALLOCATE(
is2pt_fi(i)%P(lenr),stat=ierror1)
4064 *
ALLOCATE(
ispt2_fi(i)%P(lenr),stat=ierror1)
4068 *
ALLOCATE(
isegpt_fi(i)%P(lenr),stat=ierror1)
4072 *
ALLOCATE(
is2se_fi(i)%P(2,lenr),stat=ierror1)
4076 IF(intnitsche >0)
THEN
4078 *
ALLOCATE(
forneqsfi(i)%P(3,lenr),stat=ierror1)
4085 IF(.NOT.
ASSOCIATED(
pmainfi(i)%P))
4086 *
ALLOCATE(
pmainfi(i)%P(lenr),stat=ierror1)
4089 *
ALLOCATE(
irtlm_fi(i)%P(4,lenr),stat=ierror1)
4092 *
ALLOCATE(
icont_i_fi(i)%P(lenr),stat=ierror1)
4096 *
ALLOCATE(
icodt_fi(i)%P(lenr),stat=ierror1)
4101 *
ALLOCATE(
iskew_fi(i)%P(lenr),stat=ierror1)
4107 *
ALLOCATE(
islide_fi(i)%P(4,lenr),stat=ierror1)
4111 IF(flagremn==2)
THEN
4113 *
ALLOCATE(
kremnor_fi(i)%P(lenr+1),stat=ierror1)
4116 IF(sizremnorfi /= 0)
THEN
4118 *
ALLOCATE(
remnor_fi(i)%P(sizremnorfi),stat=ierror1)
4120 ELSE IF(sizremnorfi == 0)
THEN
4122 *
ALLOCATE(
remnor_fi(i)%P(sizremnorfi),stat=ierror1)
4128 IF(idtmins_old == 2)
THEN
4129 ALLOCATE (
nodnxfi(i)%P(lenr),stat=ierror1)
4130 ierror = ierror + ierror1
4132 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4133 ierror = ierror + ierror1
4135 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4136 ierror = ierror + ierror1
4140 ierror = ierror + ierror1
4143 ELSEIF(idtmins_int_old /= 0)
THEN
4144 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4145 ierror = ierror + ierror1
4147 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4148 ierror = ierror + ierror1
4152 ierror = ierror + ierror1
4158 IF(.NOT.
ASSOCIATED(
nodnxfi(i)%P))
THEN
4159 ALLOCATE (
nodnxfi(i)%P(lenr),stat=ierror1)
4160 ierror = ierror + ierror1
4163 IF(.NOT.
ASSOCIATED(
nodamsfi(i)%P))
THEN
4164 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4165 ierror = ierror + ierror1
4168 IF(.NOT.
ASSOCIATED(
procamsfi(i)%P))
THEN
4169 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4170 ierror = ierror + ierror1
4176 ierror = ierror + ierror1
4180 ELSEIF(idtmins_int /= 0)
THEN
4181 IF(.NOT.
ASSOCIATED(
nodamsfi(i)%P))
THEN
4182 ALLOCATE (
nodamsfi(i)%P(lenr),stat=ierror1)
4183 ierror = ierror + ierror1
4186 IF(.NOT.
ASSOCIATED(
procamsfi(i)%P))
THEN
4187 ALLOCATE (
procamsfi(i)%P(lenr),stat=ierror1)
4188 ierror = ierror + ierror1
4194 ierror = ierror + ierror1
4200 ALLOCATE(
msfi(i)%P(lenr),stat=ierror1)
4201 ierror = ierror + ierror1
4202 ALLOCATE(
stifi(i)%P(lenr),stat=ierror1)
4203 ierror = ierror + ierror1
4205 ALLOCATE(
gapfi(i)%P(lenr),stat=ierror1)
4206 ierror = ierror + ierror1
4208 ALLOCATE(
gap_lfi(i)%P(lenr),stat=ierror1)
4209 ierror = ierror + ierror1
4212 ALLOCATE(
xfi(i)%P(3,lenr),stat=ierror1)
4213 ierror = ierror + ierror1
4214 ALLOCATE(
vfi(i)%P(3,lenr),stat=ierror1)
4215 ierror = ierror + ierror1
4218 ALLOCATE(
afi(i)%P(3,lenr*nthread),stat=ierror1)
4219 ierror = ierror + ierror1
4220 ALLOCATE(
stnfi(i)%P(lenr*nthread),stat=ierror1)
4221 ierror = ierror + ierror1
4223 ALLOCATE(
vscfi(i)%P(lenr*nthread),stat=ierror1)
4224 ierror = ierror + ierror1
4235 vscfi(i)%P(1:lenr*nthread)=zero
4242 ALLOCATE(
fthefi(i)%P(lenr*nthread),stat=ierror1)
4243 ierror = ierror + ierror1
4247 ALLOCATE(
tempfi(i)%P(lenr),stat=ierror1)
4248 ierror = ierror + ierror1
4252 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4253 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4254 ierror = ierror + ierror1
4257 IF(ityp == 25.AND.ivis2==-1)
THEN
4258 ALLOCATE(
if_adhfi(i)%P(lenr),stat=ierror1)
4259 ierror = ierror + ierror1
4262 IF(idt_therm ==1.AND.intth > 0)
THEN
4263 ALLOCATE(
condnfi(i)%P(lenr*nthread),stat=ierror1)
4264 ierror = ierror + ierror1
4265 condnfi(i)%P(1:lenr*nthread)=zero
4275 ALLOCATE(
tempfi(i)%P(lenr),stat=ierror1)
4276 ierror = ierror + ierror1
4279 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4280 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4281 ierror = ierror + ierror1
4283 IF(ityp == 25.AND.ivis2==-1)
THEN
4284 ALLOCATE(
if_adhfi(i)%P(lenr),stat=ierror1)
4285 ierror = ierror + ierror1
4290 ALLOCATE(
time_sfi(i)%P(lenr),stat=ierror1)
4292 ierror = ierror + ierror1
4293 ELSEIF(ityp == 25)
THEN
4294 ALLOCATE(
time_sfi(i)%P(2*lenr),stat=ierror1)
4296 ierror = ierror + ierror1
4298 IF(ityp == 24 .OR. ityp == 25)
THEN
4299 ALLOCATE(
secnd_frfi(i)%P(6,lenr),stat=ierror1)
4301 ierror = ierror + ierror1
4302 ALLOCATE(
pene_oldfi(i)%P(5,lenr),stat=ierror1)
4304 ierror = ierror + ierror1
4305 ALLOCATE(
stif_oldfi(i)%P(2,lenr),stat=ierror1)
4307 ierror = ierror + ierror1
4309 IF(ipari(97,i) > 0)
THEN
4314 IF(ityp == 25.AND.parameters%INTCAREA > 0.AND.ipari(36,i)>0)
THEN
4316 ierror = ierror + ierror1
4319 IF(ityp == 24.AND.parameters%INTCAREA > 0)
THEN
4321 ierror = ierror + ierror1
4325 ELSEIF(ityp==11)
THEN
4327 ALLOCATE(
itafi(i)%P(2*lenr),stat=ierror1)
4328 ierror = ierror + ierror1
4331 IF(idtmins_old == 2)
THEN
4332 ALLOCATE (
nodnxfi(i)%P(2*lenr),stat=ierror1)
4333 ierror = ierror + ierror1
4335 ALLOCATE (
nodamsfi(i)%P(2*lenr),stat=ierror1)
4336 ierror = ierror + ierror1
4338 ALLOCATE (
procamsfi(i)%P(2*lenr),stat=ierror1)
4339 ierror = ierror + ierror1
4341 ELSEIF(idtmins_int_old /= 0)
THEN
4342 ALLOCATE (
nodamsfi(i)%P(2*lenr),stat=ierror1)
4343 ierror = ierror + ierror1
4345 ALLOCATE (
procamsfi(i)%P(2*lenr),stat=ierror1)
4346 ierror = ierror + ierror1
4351 ALLOCATE(
matsfi(i)%P(lenr),stat=ierror1)
4352 ierror = ierror + ierror1
4356 IF(intfric > 0)
THEN
4358 ierror = ierror + ierror
4362 ALLOCATE(
msfi(i)%P(2*lenr),stat=ierror1)
4363 ierror = ierror + ierror1
4364 ALLOCATE(
stifi(i)%P(lenr),stat=ierror1)
4365 ierror = ierror + ierror1
4367 ALLOCATE(
gapfi(i)%P(lenr),stat=ierror1)
4368 ierror = ierror + ierror1
4370 ALLOCATE(
gap_lfi(i)%P(lenr),stat=ierror1)
4371 ierror = ierror + ierror1
4374 ALLOCATE(
xfi(i)%P(3,2*lenr),stat=ierror1)
4375 ierror = ierror + ierror1
4376 ALLOCATE(
vfi(i)%P(3,2*lenr),stat=ierror1)
4377 ierror = ierror + ierror1
4378 IF(inacti==5.OR.inacti==6)
THEN
4379 ALLOCATE(
penfi(i)%P(2,lenr),stat=ierror1)
4382 ALLOCATE(
afi(i)%P(3,2*lenr*nthread),stat=ierror1)
4383 ierror = ierror + ierror1
4384 ALLOCATE(
stnfi(i)%P(2*lenr*nthread),stat=ierror1)
4385 ierror = ierror + ierror1
4387 ALLOCATE(
vscfi(i)%P(2*lenr*nthread),stat=ierror1)
4388 ierror = ierror + ierror1
4391 ALLOCATE(
fthefi(i)%P(2*lenr*nthread),stat=ierror1)
4392 ierror = ierror + ierror1
4393 ALLOCATE(
tempfi(i)%P(2*lenr),stat=ierror1)
4394 ierror = ierror + ierror1
4397 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4398 ierror = ierror + ierror1
4400 IF(idt_therm ==1.AND.intth > 0)
THEN
4401 ALLOCATE(
condnfi(i)%P(2*lenr*nthread),stat=ierror1)
4402 ierror = ierror + ierror1
4411 ALLOCATE(
tempfi(i)%P(2*lenr),stat=ierror1)
4412 ierror = ierror + ierror1
4417 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4418 ALLOCATE(
areasfi(i)%P(lenr),stat=ierror1)
4419 ierror = ierror + ierror1
4423 ELSEIF(ityp==17)
THEN
4424 ALLOCATE(
xfi17(i)%P(3,16,lenr),stat=ierror1)
4425 ierror = ierror + ierror1
4426 ALLOCATE(
vfi17(i)%P(3,16,lenr),stat=ierror1)
4427 ierror = ierror + ierror1
4428 ALLOCATE(
frotsfi(i)%P(7,lenr),stat=ierror1)
4429 ierror = ierror + ierror1
4430 ALLOCATE(
ksfi(i)%P(2,lenr),stat=ierror1)
4431 ierror = ierror + ierror1
4432 ALLOCATE(
eminxfi(i)%P(6,lenr),stat=ierror1)
4433 ierror = ierror + ierror1
4435 ALLOCATE(
afi17(i)%P(3,16,lenr),stat=ierror1)
4436 ierror = ierror + ierror1
4437 ALLOCATE(
stnfi17(i)%P(16,lenr),stat=ierror1)
4438 ierror = ierror + ierror1
4441 lskyfi = lenr * multimax
4444 ALLOCATE(
iskyfi(i)%P(lskyfi),stat=ierror1)
4445 ierror = ierror + ierror1
4446 ALLOCATE(
fskyfi(i)%P(40,lskyfi),stat=ierror1)
4447 ierror = ierror + ierror1
4457 CALL ancmsg(msgid=20,anmode=aninfo)
4462 ALLOCATE(
daanc6fi(i)%P(3,6,lenr),stat=ierror1)
4463 ierror = ierror + ierror1
4464 ALLOCATE(
dxancfi(i)%P(3,lenr),stat=ierror1)
4465 ierror = ierror + ierror1
4466 ALLOCATE(
dvancfi(i)%P(3,lenr),stat=ierror1)
4467 ierror = ierror + ierror1
4468 IF((inacti==5.OR.inacti==6) .AND. lenr > 0)
THEN
4469 ALLOCATE(
penfi(i)%P(2,lenr),stat=ierror1)
4470 ierror = ierror + ierror1
4471 ALLOCATE(
penfia(i)%P(5,lenr),stat=ierror1)
4472 ierror = ierror + ierror1
4474 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
4475 ALLOCATE (
diag_smsfi(i)%P(lenr),stat=ierror1)
4476 ierror = ierror + ierror1
4478 ALLOCATE(
alphakfi(i)%P(lenr),stat=ierror1)
4479 ierror = ierror + ierror1
4480 ALLOCATE(
daancfi(i)%P(3,lenr),stat=ierror1)
4482 ierror = ierror + ierror1
4484 CALL ancmsg(msgid=20,anmode=aninfo)
4488 ALLOCATE(
nsnsie(i)%P(nspmd),stat=ierror1)
4489 ierror = ierror + ierror1
4491 ierror = ierror + ierror1
4497 lens = lens +
nsnsie(i)%P(p)
4498 lenr = lenr +
nsnfie(i)%P(p)
4502 ALLOCATE(
nsvsie(i)%P(lens),stat=ierror1)
4504 ierror = ierror + ierror1
4508 ALLOCATE(
nsvfie(i)%P(lenr),stat=ierror1)
4509 ierror = ierror + ierror1
4511 ALLOCATE(
itafie(i)%P(2*lenr),stat=ierror1)
4512 ierror = ierror + ierror1
4514 ALLOCATE(
msfie(i)%P(2*lenr),stat=ierror1)
4515 ierror = ierror + ierror1
4516 ALLOCATE(
stifie(i)%P(lenr),stat=ierror1)
4517 ierror = ierror + ierror1
4519 ALLOCATE(
gapfie(i)%P(lenr),stat=ierror1)
4520 ierror = ierror + ierror1
4522 ALLOCATE(
xfie(i)%P(3,2*lenr),stat=ierror1)
4523 ierror = ierror + ierror1
4524 ALLOCATE(
vfie(i)%P(3,2*lenr),stat=ierror1)
4525 ierror = ierror + ierror1
4526 IF(inacti==5.OR.inacti==6)
THEN
4527 ALLOCATE(
penfie(i)%P(2,lenr),stat=ierror1)
4528 ierror = ierror + ierror1
4529 ALLOCATE(
penfiae(i)%P(5,2*lenr),stat=ierror1)
4530 ierror = ierror + ierror1
4533 ALLOCATE(
afie(i)%P(3,2*lenr*nthread),stat=ierror1)
4534 ierror = ierror + ierror1
4535 ALLOCATE(
stnfie(i)%P(2*lenr*nthread),stat=ierror1)
4536 ierror = ierror + ierror1
4538 ALLOCATE(
vscfie(i)%P(2*lenr*nthread),stat=ierror1)
4539 ierror = ierror + ierror1
4543 lskyfi = lenr * multimax
4545 ALLOCATE(
iskyfie(i)%P(lskyfi),stat=ierror1)
4546 ierror = ierror + ierror1
4549 ALLOCATE(
fskyfie(i)%P(10,lskyfi),stat=ierror1)
4552 ALLOCATE(
fskyfie(i)%P(10,lskyfi),stat=ierror1)
4554 ierror = ierror + ierror1
4556 ALLOCATE(
daanc6fie(i)%P(3,6,2*lenr),stat=ierror1)
4557 ierror = ierror + ierror1
4558 ALLOCATE(
dxancfie(i)%P(3,2*lenr),stat=ierror1)
4559 ierror = ierror + ierror1
4560 ALLOCATE(
dvancfie(i)%P(3,2*lenr),stat=ierror1)
4561 ierror = ierror + ierror1
4562 ALLOCATE(
alphakfie(i)%P(2*lenr),stat=ierror1)
4563 ierror = ierror + ierror1
4564 ALLOCATE(
daancfie(i)%P(3,lenr),stat=ierror1)
4565 ierror = ierror + ierror1
4566 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
4568 ierror = ierror + ierror1
4571 IF(idtmins_old == 2)
THEN
4572 ALLOCATE (
nodnxfie(i)%P(lenr),stat=ierror1)
4573 ierror = ierror + ierror1
4575 ALLOCATE (
nodamsfie(i)%P(lenr),stat=ierror1)
4576 ierror = ierror + ierror1
4578 ALLOCATE (
procamsfie(i)%P(lenr),stat=ierror1)
4579 ierror = ierror + ierror1
4581 ELSEIF(idtmins_int_old /= 0)
THEN
4582 ALLOCATE (
nodamsfie(i)%P(lenr),stat=ierror1)
4583 ierror = ierror + ierror1
4585 ALLOCATE (
procamsfie(i)%P(lenr),stat=ierror1)
4586 ierror = ierror + ierror1
4591 CALL ancmsg(msgid=20,anmode=aninfo)
4600 IF(ityp == 25 )
THEN
4602 ALLOCATE(
nsnsie(i)%P(nspmd),stat=ierror1)
4603 ALLOCATE(
nsnfie(i)%P(nspmd),stat=ierror1)
4608 IF(ipari(58,i) /=0)
THEN
4616 lens_edge = lens_edge +
nsnsie(i)%P(p)
4617 lenr_edge = lenr_edge +
nsnfie(i)%P(p)
4620 IF(lens_edge>0)
THEN
4621 ALLOCATE(
nsvsie(i)%P(lens_edge),stat=ierror1)
4623 ierror = ierror + ierror1
4626 IF(iparit == 0)
THEN
4628 ALLOCATE(
afie(i)%P(3,2*lenr_edge*nthread),stat=ierror1)
4629 ierror = ierror + ierror1
4630 ALLOCATE(
stnfie(i)%P(2*lenr_edge*nthread),stat=ierror1)
4631 ierror = ierror + ierror1
4633 ALLOCATE(
vscfie(i)%P(2*lenr_edge*nthread),stat=ierror1)
4634 ierror = ierror + ierror1
4637 lskyfi = lenr_edge * multimax
4639 ALLOCATE(
iskyfie(i)%P(lskyfi),stat=ierror1)
4640 ierror = ierror + ierror1
4642 ALLOCATE(
fskyfie(i)%P(8,lskyfi),stat=ierror1)
4644 ALLOCATE(
fskyfie(i)%P(8,lskyfi),stat
4646 ierror = ierror + ierror1
4648 IF(lenr_edge>0)
THEN
4649 ALLOCATE(
nsvfie(i)%P(lenr_edge))
4650 ALLOCATE(
itafie(i)%P(lenr_edge*2))
4651 ALLOCATE(
ledge_fie(i)%P(e_ledge_size,lenr_edge))
4652 ALLOCATE(
xfie(i)%P(3,lenr_edge*2))
4653 ALLOCATE(
vfie(i)%P(3,lenr_edge*2))
4654 ALLOCATE(
msfie(i)%P(lenr_edge*2))
4659 ALLOCATE(
stifie(i)%P(lenr_edge))
4670 ipari(69,i) = lenr_edge
4671 IF(intfric > 0)
THEN
4675 IF(ipari(97,i) > 0)
THEN
4686 IF(ipari(36,i)>0.AND.ipari(7,i)/=17)
THEN
4689 ALLOCATE(
nisubsfi(i)%P(nspmd),stat=ierror1)
4690 ierror = ierror + ierror1
4696 lenr = lenr +
nsnfi(i)%P(p)
4701 ALLOCATE(
addsubsfi(i)%P(lenr+1),stat=ierror1)
4702 ierror = ierror + ierror1
4705 ALLOCATE(
lisubsfi(i)%P(leni),stat=ierror1)
4706 ierror = ierror + ierror1
4708 IF(ipari(7,i)==25.OR.ipari(7,i)==7.OR.ipari(7,i)==24.OR.ipari(7,i)==11)
THEN
4710 ierror = ierror + ierror1
4715 ALLOCATE(
addsubsfi(i)%P(lenr+1),stat=ierror1)
4722 CALL ancmsg(msgid=20,anmode=aninfo)
4727 IF(ipari(7,i) == 25 .AND. ipari(58,i) > 0)
THEN
4729 ALLOCATE(
nisubsfie(i)%P(nspmd),stat=ierror1)
4730 ierror = ierror + ierror1
4737 IF(lenr_edge>0)
THEN
4738 ALLOCATE(
addsubsfie(i)%P(lenr_edge+1),stat=ierror1)
4739 ierror = ierror + ierror1
4742 ALLOCATE(
lisubsfie(i)%P(leni),stat=ierror1)
4743 ierror = ierror + ierror1
4746 ierror = ierror + ierror1
4749 IF(lenr_edge>0)
THEN
4750 ALLOCATE(
addsubsfie(i)%P(lenr_edge+1),stat=ierror1)
4760 IF(anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT >0)
THEN
4762 ALLOCATE(
fnconti(i)%P(3,lenr),stat=ierror1)
4763 ierror = ierror + ierror1
4764 ALLOCATE(
ftconti(i)%P(3,lenr),stat=ierror1)
4765 ierror = ierror + ierror1
4767 CALL ancmsg(msgid=20,anmode=aninfo)
4787 IF(intth==2.OR.ipari(95,i) > 0)
THEN
4793 ALLOCATE(
nmnsi(i)%P(nspmd),stat=ierror1)
4794 ierror = ierror + ierror1
4796 ALLOCATE(
nmnfi(i)%P(nspmd),stat=ierror1)
4797 ierror = ierror + ierror1
4799 CALL ancmsg(msgid=20,anmode=aninfo)
4807 lens = lens +
nmnsi(i)%P(p)
4808 lenr = lenr +
nmnfi(i)%P(p)
4815 ALLOCATE(
nmvsi(i)%P(lens),stat=ierror1)
4817 ierror = ierror + ierror1
4819 ALLOCATE(
tempnod(i)%P(lens),stat=ierror1)
4821 ierror = ierror + ierror1
4826 ALLOCATE(
nmvfi(i)%P(lenr),stat=ierror1)
4827 ierror = ierror + ierror1
4830 CALL ancmsg(msgid=20,anmode=aninfo)
4836 IF(ipari(95,i) > 0)
THEN
4843 IF(ipari(48,i) /= 0 )
THEN
4845 ALLOCATE(
fthefi(i)%P(lenr*nthread),stat=ierror1)
4846 ierror = ierror + ierror1
4848 fthefi(i)%P(1:lenr*nthread)=zero
4850 lskyfi = lenr * multimax
4851 ALLOCATE(
iskyfi(i)%P(lskyfi),stat=ierror1)
4852 ierror = ierror + ierror1
4853 ALLOCATE(
ftheskyfi(i)%P(lskyfi),stat=ierror1)
4854 ierror = ierror + ierror1
4877 inacti = ipari(22,i)
4879 IF(ityp==7.OR.ityp==10.OR.ityp==11.OR.
4880 + (ityp==17.AND.ipari(33,i)==0).OR.ityp==20.OR.
4881 + ityp==22.OR.ityp==23.OR.ityp==24.OR.
4885 lenr = lenr +
nsnfi(i)%P(p)
4888 IF(ityp==7.OR.ityp==22.OR.ityp==23.OR.
4889 + ityp==20.OR.ityp==10.OR.ityp==24.OR.
4902 IF((ityp == 7.OR.ityp == 20.OR.
4903 + ityp==22.OR.ityp==23.OR.ityp
4904 + ityp==25).AND.ipari(47,i) >
THEN
4907 IF(intth>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
4910 IF(ityp == 25.AND.ivis2==-1)
THEN
4916 ELSEIF(ityp==11)
THEN
4927 IF(inacti==5.OR.inacti==6)
THEN
4930 IF(ipari(47,i)>0)
THEN
4934 ELSEIF(ityp==17.AND.ipari(33,i)==0)
THEN
4950 IF(inacti==5.OR.inacti==6)
THEN
4956 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
4963 lenr = lenr +
nsnfie(i)%P(p)
4980 IF(idtmins_old > 0 .OR. idtmins_int_old /= 0)
THEN
4987 ELSEIF(ityp==25)
THEN
4990 IF (ityp==24.OR.ityp==25)
THEN
4994 IF(ipari(97,i) > 0)
THEN
4997 IF(ityp==25.AND.parameters%INTCAREA > 0.AND.ipari(36,i)>0)
THEN
5000 IF(ityp==24.AND.parameters%INTCAREA > 0)
THEN
5006 IF(idtmins_old == 2)
THEN
5008 ierror = ierror + ierror1
5010 ELSEIF(idtmins_int_old /= 0)
THEN
5012 ierror = ierror + ierror1
5019 ierror = ierror + ierror1
5022 ELSEIF(idtmins_int /= 0)
THEN
5025 ierror = ierror + ierror1
5031 IF(ityp == 25 )
THEN
5032 IF(ipari(58,i) /=0)
THEN
5035 lenr = lenr +
nsnfie(i)%P(p)
5049 IF(idtmins_old == 2)
THEN
5056 ELSEIF(idtmins_int_old /= 0)
THEN
5063 IF(ipari(97,i) > 0)
THEN
5076 lenr = lenr +
nmnfi(i)%P(p)
5084 ALLOCATE(
nmtemp(i)%P(lenr),stat=ierror1)
5085 ierror = ierror + ierror1
5089 CALL ancmsg(msgid=20,anmode=aninfo)
5103!||
spmd_savefi ../engine/source/mpi/interfaces/spmd_i7tool.f
5110!||--- uses -----------------------------------------------------
5111!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.f90
5118!||====================================================================
5119 SUBROUTINE spmd_savefi(IPARI, IFLAG,INTBUF_TAB,NSENSOR,SENSOR_TAB,PARAMETERS)
5133#include "implicit_f.inc"
5137#include "com01_c.inc"
5138#include "com04_c.inc"
5139#include "com08_c.inc"
5140#include "param_c.inc"
5145 INTEGER ,
INTENT(IN) :: NSENSOR
5146 INTEGER IPARI(NPARI,*), IFLAG
5147 TYPE(intbuf_struct_) INTBUF_TAB(*)
5148 TYPE (SENSOR_STR_) ,
DIMENSION(NSENSOR) :: SENSOR_TAB
5149 TYPE (PARAMETERS_) ,
INTENT(IN)::
5153 INTEGER I, P, NSN, NMN, ITYP, IGAP, LENS, LENR
5155 INTEGER INTERACT,ISENS
5157 INTEGER,
DIMENSION(NSPMD) :: SAVE_NSNFI,SAVE_NSNSI
5172 + ityp==10.OR.ityp==11.OR.
5173 + (ityp==17.AND.ipari(33,i)==0).OR.ityp==20.OR.
5174 + ityp==22.OR.ityp==23.OR.ityp==24.OR.
5180 save_nsnfi(1:nspmd) =
nsnfi(i)%P(1:nspmd)
5181 save_nsnsi(1:nspmd) =
nsnsi(i)%P(1:nspmd)
5184 IF (inacti==5.OR.inacti==6.OR.inacti==7.OR.ityp==23.OR.inacti==-1)
THEN
5187 IF(ityp == 7.OR.ityp == 11.OR.ityp == 24)
5188 . isens = ipari(64,i)
5190 ts = sensor_tab(isens)%TSTART
5191 IF (ts>tt) interact = 0
5193 startt = intbuf_tab(i)%VARIABLES(3)
5194 stopt = intbuf_tab(i)%VARIABLES(11)
5195 IF (startt>tt) interact = 0
5198 IF((interact == 0 .OR.
nsnfi_flag(i)==1).AND.ityp/=25)
THEN
5202 nsnfi_sav(i)%P(1:nspmd) = save_nsnfi(1:nspmd)
5203 nsnsi_sav(i)%P(1:nspmd) = save_nsnsi(1:nspmd)
5211 intnitsche=ipari(86,i)
5216 lens = lens +
nsnsi(i)%P(p)
5217 lenr = lenr +
nsnfi(i)%P(p)
5228 IF(ityp==7.OR.ityp==20.OR.ityp==22.OR.
5229 + ityp==23.OR.ityp==24.OR.ityp==25)
THEN
5238 IF((ityp == 7.OR.ityp == 24.OR.ityp == 25).AND.ipari(72,i) > 0)
THEN
5268 IF(ipari(63,i)==2)
THEN
5281 ELSEIF(idtmins_int/=0)
THEN
5287 ELSEIF(ityp==10)
THEN
5294 ELSEIF(idtmins_int/=0)
THEN
5299 ELSEIF(ityp==11)
THEN
5306 ELSEIF(idtmins_int/=0)
THEN
5315 ELSEIF(ityp==17)
THEN
5323 lens = lens +
nsnsie(i)%P(p)
5324 lenr = lenr +
nsnfie(i)%P(p)
5341 ELSEIF(idtmins_int/=0)
THEN
5348 IF(ityp == 25 )
THEN
5349 IF(ipari(58,i) /=0)
THEN
5355 lens = lens +
nsnsie(i)%P(p)
5356 lenr = lenr +
nsnfie(i)%P(p)
5370 IF(ipari(36,i)>0.AND.ipari(7,i)/=17)
THEN
5376 lenr = lenr +
nsnfi(i)%P(p)
5383 IF(ipari(7,i)==25.OR.ipari(7,i)==7.OR.ipari(7,i)==24.OR.ipari(7,i)==11)
THEN
5389 IF(ipari(7,i) == 25 .AND. ipari(58,i) > 0)
THEN
5397 lenr = lenr +
nsnfie(i)%P(p)
5399 IF(leni > 0 .AND. lenr > 0)
THEN
5418 IF(intth==2.OR.ipari(95,i) > 0)
THEN
5422 lens = lens +
nmnsi(i)%P(p)
5423 lenr = lenr +
nmnfi(i)%P(p)
5450 inacti = ipari(22,i)
5454 + ityp==10.OR.ityp==11.OR.
5455 + (ityp==17.AND.ipari(33,i)==0).OR.ityp==20.OR.
5456 + ityp==22.OR.ityp==23.OR.ityp==24.OR.
5460 lenr = lenr +
nsnfi(i)%P(p)
5468 IF(ityp == 7.OR.ityp == 11.OR.ityp == 24) isens = ipari(64,i)
5471 ts = sensor_tab(isens)%TSTART
5472 IF (ts>tt) interact = 0
5476 IF (inacti==5.OR.inacti==6.OR.inacti==7.OR.ityp==23.OR.inacti==-1)
THEN
5479 IF(ityp == 7.OR.ityp == 11.OR.ityp == 24) isens = ipari(64,i)
5481 startt = intbuf_tab(i)%VARIABLES(3)
5482 stopt = intbuf_tab(i)%VARIABLES(11)
5483 IF (startt>tt) interact = 0
5487 IF((interact == 0 .OR.
nsnfi_flag(i)==1).AND.ityp/=25)
THEN
5495 IF(ityp==7.OR.ityp==22.OR.ityp==23.OR.
5496 + ityp==10.OR.ityp==20.OR.ityp==24.OR.
5505 IF((ityp == 7.OR.ityp == 20.OR.ityp == 22.OR.
5506 + ityp == 23.OR.ityp == 24.OR.ityp == 25).AND.
5507 + ipari(47,i)>0)
THEN
5510 IF(ipari(47,i)>0.OR.(ityp == 25.AND.ivis2==-1))
THEN
5513 IF(ityp == 25.AND.ivis2==-1)
THEN
5517 ELSEIF(ityp==11)
THEN
5524 IF(inacti==5.OR.inacti==6)
5526 IF(ipari(47,i)>0)
THEN
5531 ELSEIF(ityp==17)
THEN
5547 IF(inacti==5.OR.inacti==6)
THEN
5553 IF(idtmins > 0 .OR. idtmins_int/=0)
THEN
5560 lenr = lenr +
nsnfie(i)%P(p)
5569 IF(inacti==5.OR.inacti==6)
THEN
5576 IF(idtmins > 0 .OR. idtmins_int/=0)
THEN
5585 ELSEIF(ityp==25)
THEN
5588 IF (ityp==24.OR.ityp==25)
THEN
5594 IF(ityp==25.AND.parameters%INTCAREA > 0.AND.ipari(36,i)>0)
CALL write_db(
intareanfi(i)%P(1),lenr)
5600 ELSEIF(idtmins_int/=0)
THEN
5606 IF(ityp == 25 )
THEN
5607 IF(ipari(58,i) /=0)
THEN
5610 lenr = lenr +
nsnfie(i)%P(p)
5629 ELSEIF(idtmins_int/=0)
THEN
5634 IF(ipari(97,i) > 0)
THEN
5648 lenr = lenr +
nmnfi(i)%P(p)
5675#include
"implicit_f.inc"
5679 INTEGER N, ISKY(*),INDEX(*)
5683 INTEGER I, J, , IR, IIA,IIA2
5718 IF(isky(j)<isky(j+1))j=j+1
5739!||====================================================================
5740 SUBROUTINE sorti7(N,ISKY,FSKYI,FTHESKYI,NFSKYI)
5745#include "implicit_f.inc"
5749 INTEGER N, NFSKYI, ISKY(*)
5751 . fskyi(nfskyi,*),ftheskyi(*)
5755 INTEGER I, J, L, IR, IIA
5757 . rra1, rra2, rra3, rra4, rra5,ppa
5783 fskyi(1,ir)=fskyi(1,1)
5784 fskyi(2,ir)=fskyi(2,1)
5785 fskyi(3,ir)=fskyi(3,1)
5786 fskyi(4,ir)=fskyi(4,1)
5787 ftheskyi(ir)=ftheskyi(1)
5804 IF(isky(j)<isky(j+1))j=j+1
5808 fskyi(1,i)=fskyi(1,j)
5809 fskyi(2,i)=fskyi(2,j)
5810 fskyi(3,i)=fskyi(3,j)
5811 fskyi(4,i)=fskyi(4,j)
5812 ftheskyi(i)=ftheskyi(j)
5846 fskyi(1,ir)=fskyi(1,1)
5847 fskyi(2,ir)=fskyi(2,1)
5848 fskyi(3,ir)=fskyi(3,1)
5849 fskyi(4,ir)=fskyi(4,1)
5850 fskyi(5,ir)=fskyi(5,1)
5851 ftheskyi(ir)=ftheskyi(1)
5869 IF(isky(j)<isky(j+1))j=j+1
5873 fskyi(1,i)=fskyi(1,j)
5874 fskyi(2,i)=fskyi(2,j)
5875 fskyi(3,i)=fskyi(3,j)
5876 fskyi(4,i)=fskyi(4,j)
5877 fskyi(5,i)=fskyi(5,j)
5878 ftheskyi(i)=ftheskyi(j)
5901!||====================================================================
5902 SUBROUTINE sorti7t(N,ISKY,FSKYI,FTHESKYI,CONDNSKYI,NFSKYI)
5907#include "implicit_f.inc"
5911 INTEGER N, NFSKYI, ISKY(
5913 . fskyi(nfskyi,*),ftheskyi(*),condnskyi(*)
5917 INTEGER I, J, L, IR, IIA
5919 . rra1, rra2, rra3, rra4, rra5,ppa,kka
5946 fskyi(1,ir)=fskyi(1,1)
5947 fskyi(2,ir)=fskyi(2,1)
5948 fskyi(3,ir)=fskyi(3,1)
5949 fskyi(4,ir)=fskyi(4,1)
5950 ftheskyi(ir)=ftheskyi(1)
5951 condnskyi(ir)=condnskyi(1)
5969 IF(isky(j)<isky(j+1))j=j+1
5973 fskyi(1,i)=fskyi(1,j)
5974 fskyi(2,i)=fskyi(2,j)
5975 fskyi(3,i)=fskyi(3,j)
5976 fskyi(4,i)=fskyi(4,j)
5977 ftheskyi(i)=ftheskyi(j)
5978 condnskyi(i)=condnskyi(j)
6014 fskyi(1,ir)=fskyi(1,1)
6015 fskyi(2,ir)=fskyi(2,1)
6016 fskyi(3,ir)=fskyi(3,1)
6017 fskyi(4,ir)=fskyi(4,1)
6018 fskyi(5,ir)=fskyi(5,1)
6019 ftheskyi(ir)=ftheskyi(1)
6020 condnskyi(ir)=condnskyi(1)
6039 IF(isky(j)<isky(j+1))j=j+1
6043 fskyi(1,i)=fskyi(1,j)
6044 fskyi(2,i)=fskyi(2,j)
6045 fskyi(3,i)=fskyi(3,j)
6046 fskyi(4,i)=fskyi(4,j)
6047 fskyi(5,i)=fskyi(5,j)
6048 ftheskyi(i)=ftheskyi(j)
6049 condnskyi(i)=condnskyi(j)
6082#include "implicit_f.inc"
6086 INTEGER N, NFSKYI, ISKY(*)
6092 INTEGER I, J, L, IR, IIA
6094 . rra1, rra2, rra3, rra4, rra5
6118 fskyi(1,ir)=fskyi(1,1)
6119 fskyi(2,ir)=fskyi(2,1)
6120 fskyi(3,ir)=fskyi(3,1)
6121 fskyi(4,ir)=fskyi(4,1)
6137 IF(isky(j)<isky(j+1))j=j+1
6141 fskyi(1,i)=fskyi(1,j)
6142 fskyi(2,i)=fskyi(2,j)
6144 fskyi(4,i)=fskyi(4,j)
6176 fskyi(1,ir)=fskyi(1,1)
6177 fskyi(2,ir)=fskyi(2,1)
6178 fskyi(3,ir)=fskyi(3,1)
6179 fskyi(4,ir)=fskyi(4,1)
6180 fskyi(5,ir)=fskyi(5,1)
6197 IF(isky(j)<isky(j+1))j=j+1
6201 fskyi(1,i)=fskyi(1,j)
6202 fskyi(2,i)=fskyi(2,j)
6203 fskyi(3,i)=fskyi(3,j)
6204 fskyi(4,i)=fskyi(4,j)
6205 fskyi(5,i)=fskyi(5,j)
6237#include "implicit_f.inc"
6241#include "com01_c.inc"
6246 . nsnfi(*), isky(*), isizenv(2,*), nsnfitot(*)
6250 INTEGER IG, P, I, LASTIG
6262 isizenv(1,p) = isizenv(1,p) + len
6263 isizenv(2,p) = isizenv(2,p) + 1
6267 lastig = lastig+nsnfi(p)
6273 nsnfitot(p) = nsnfitot(p) + nsnfi(p)
6289#include
"implicit_f.inc"
6293#include "com01_c.inc"
6298 . nsnfi(*), nsvfi(*), isizenv(2,*)
6302 INTEGER IDEB, , N, NB
6310 isizenv(1,p) = isizenv(1,p) + nb*leni20
6333#include
"implicit_f.inc"
6337 INTEGER N, NFSKYI, ISKY(*)
6343 INTEGER I, J, L, IR, IIA
6345 . RRA1, RRA2, RRA3, RRA4, ,
6346 . RRA6, RRA7, RRA8, RRA9, RRA10
6378 fskyi(1,ir)=fskyi(1,1)
6379 fskyi(2,ir)=fskyi(2,1)
6380 fskyi(3,ir)=fskyi(3,1)
6381 fskyi(4,ir)=fskyi(4,1)
6382 fskyi(5,ir)=fskyi(5,1)
6383 fskyi(6,ir)=fskyi(6,1)
6384 fskyi(7,ir)=fskyi(7,1)
6385 fskyi(8,ir)=fskyi(8,1)
6405 IF(isky(j)<isky(j+1))j=j+1
6409 fskyi(1,i)=fskyi(1,j)
6410 fskyi(2,i)=fskyi(2,j)
6411 fskyi(3,i)=fskyi(3,j)
6412 fskyi(4,i)=fskyi(4,j)
6413 fskyi(5,i)=fskyi(5,j)
6414 fskyi(6,i)=fskyi(6,j)
6415 fskyi(7,i)=fskyi(7,j)
6416 fskyi(8,i)=fskyi(8,j)
6462 fskyi(1,ir)=fskyi(1,1)
6463 fskyi(2,ir)=fskyi(2,1)
6464 fskyi(3,ir)=fskyi(3,1)
6465 fskyi(4,ir)=fskyi(4,1)
6466 fskyi(5,ir)=fskyi(5,1)
6467 fskyi(6,ir)=fskyi(6,1)
6468 fskyi(7,ir)=fskyi(7,1)
6469 fskyi(8,ir)=fskyi(8,1)
6470 fskyi(9,ir)=fskyi(9,1)
6471 fskyi(10,ir)=fskyi(10,1)
6493 IF(isky(j)<isky(j+1))j=j+1
6497 fskyi(1,i)=fskyi(1,j)
6498 fskyi(2,i)=fskyi(2,j)
6499 fskyi(3,i)=fskyi(3,j)
6500 fskyi(4,i)=fskyi(4,j)
6501 fskyi(5,i)=fskyi(5,j)
6502 fskyi(6,i)=fskyi(6,j)
6503 fskyi(7,i)=fskyi(7,j)
6504 fskyi(8,i)=fskyi(8,j)
6505 fskyi(9,i)=fskyi(9,j)
6506 fskyi(10,i)=fskyi(10,j)
6541#include "implicit_f.inc"
6545 INTEGER , NFSKYI, ISKY(*)
6547 . fskyi(2*nfskyi,*),ftheskyi(*)
6551 INTEGER I, J, L, IR, IIA
6553 . rra1, rra2, rra3, rra4, rra5,
6554 . rra6, rra7, rra8, rra9, rra10,
6589 fskyi(1,ir)=fskyi(1,1)
6590 fskyi(2,ir)=fskyi(2,1)
6591 fskyi(3,ir)=fskyi(3,1)
6592 fskyi(4,ir)=fskyi(4,1)
6593 fskyi(5,ir)=fskyi(5,1)
6594 fskyi(6,ir)=fskyi(6,1)
6595 fskyi(7,ir)=fskyi(7,1)
6596 fskyi(8,ir)=fskyi(8,1)
6597 ftheskyi(ir)=ftheskyi(1)
6618 IF(isky(j)<isky(j+1))j=j+1
6622 fskyi(1,i)=fskyi(1,j)
6623 fskyi(2,i)=fskyi(2,j)
6624 fskyi(3,i)=fskyi(3,j)
6625 fskyi(4,i)=fskyi(4,j)
6626 fskyi(5,i)=fskyi(5,j)
6627 fskyi(6,i)=fskyi(6,j)
6628 fskyi(7,i)=fskyi(7,j)
6629 fskyi(8,i)=fskyi(8,j)
6630 ftheskyi(i)=ftheskyi(j)
6679 fskyi(1,ir)=fskyi(1,1)
6680 fskyi(2,ir)=fskyi(2,1)
6681 fskyi(3,ir)=fskyi(3,1)
6682 fskyi(4,ir)=fskyi(4,1)
6683 fskyi(5,ir)=fskyi(5,1)
6685 fskyi(7,ir)=fskyi(7,1)
6686 fskyi(8,ir)=fskyi(8,1)
6687 fskyi(9,ir)=fskyi(9,1)
6688 fskyi(10,ir)=fskyi(10,1)
6689 ftheskyi(ir)=ftheskyi(1)
6712 IF(isky(j)<isky(j+1))j=j+1
6716 fskyi(1,i)=fskyi(1,j)
6717 fskyi(2,i)=fskyi(2,j)
6718 fskyi(3,i)=fskyi(3,j)
6719 fskyi(4,i)=fskyi(4,j)
6720 fskyi(5,i)=fskyi(5,j)
6721 fskyi(6,i)=fskyi(6,j)
6722 fskyi(7,i)=fskyi(7,j)
6723 fskyi(8,i)=fskyi(8,j)
6724 fskyi(9,i)=fskyi(9,j)
6725 fskyi(10,i)=fskyi(10,j)
6726 ftheskyi(i)=ftheskyi(j)
6757 SUBROUTINE sorti11tt(N,ISKY,FSKYI,FTHESKYI,CONDNSKYI,NFSKYI)
6762#include "implicit_f.inc"
6766 INTEGER N, NFSKYI, ISKY(*)
6768 . fskyi(2*nfskyi,*),ftheskyi(*),condnskyi(*)
6772 INTEGER I, J, L, IR, IIA
6774 . rra1, rra2, rra3, rra4, rra5,
6775 . rra6, rra7, rra8, rra9, rra10,
6812 fskyi(1,ir)=fskyi(1,1)
6813 fskyi(2,ir)=fskyi(2,1)
6814 fskyi(3,ir)=fskyi(3,1)
6815 fskyi(4,ir)=fskyi(4,1)
6816 fskyi(5,ir)=fskyi(5,1)
6817 fskyi(6,ir)=fskyi(6,1)
6818 fskyi(7,ir)=fskyi(7,1)
6819 fskyi(8,ir)=fskyi(8,1)
6820 ftheskyi(ir)=ftheskyi(1)
6821 condnskyi(ir)=condnskyi(1)
6843 IF(isky(j)<isky(j+1))j=j+1
6847 fskyi(1,i)=fskyi(1,j)
6848 fskyi(2,i)=fskyi(2,j)
6849 fskyi(3,i)=fskyi(3,j)
6850 fskyi(4,i)=fskyi(4,j)
6851 fskyi(5,i)=fskyi(5,j)
6852 fskyi(6,i)=fskyi(6,j)
6853 fskyi(7,i)=fskyi(7,j)
6854 fskyi(8,i)=fskyi(8,j)
6855 ftheskyi(i)=ftheskyi(j)
6856 condnskyi(i)=condnskyi(j)
6908 fskyi(1,ir)=fskyi(1,1)
6909 fskyi(2,ir)=fskyi(2,1)
6910 fskyi(3,ir)=fskyi(3,1)
6911 fskyi(4,ir)=fskyi(4,1)
6912 fskyi(5,ir)=fskyi(5,1)
6913 fskyi(6,ir)=fskyi(6,1)
6914 fskyi(7,ir)=fskyi(7,1)
6915 fskyi(8,ir)=fskyi(8,1)
6916 fskyi(9,ir)=fskyi(9,1)
6917 fskyi(10,ir)=fskyi(10,1)
6918 ftheskyi(ir)=ftheskyi(1)
6919 condnskyi(ir)=condnskyi(1)
6943 IF(isky(j)<isky(j+1))j=j+1
6947 fskyi(1,i)=fskyi(1,j)
6948 fskyi(2,i)=fskyi(2,j)
6949 fskyi(3,i)=fskyi(3,j)
6950 fskyi(4,i)=fskyi(4,j)
6951 fskyi(5,i)=fskyi(5,j)
6952 fskyi(6,i)=fskyi(6,j)
6953 fskyi(7,i)=fskyi(7,j)
6954 fskyi(8,i)=fskyi(8,j)
6955 fskyi(9,i)=fskyi(9,j)
6956 fskyi(10,i)=fskyi(10,j)
6957 ftheskyi(i)=ftheskyi(j)
6958 condnskyi(i)=condnskyi(j)
6995#include "implicit_f.inc"
6999 INTEGER N, NFSKYI, ISKY(*)
7005 INTEGER I, J, L, IR, IIA, II
7024 rra(ii)=fskyi(ii,ir)
7028 fskyi(ii,ir)=fskyi(ii,1)
7043 IF(isky(j)<isky(j+1))j=j+1
7048 fskyi(ii,i)=fskyi(ii,j)
7072!||====================================================================
7078#include "implicit_f.inc"
7120#include "implicit_f.inc"
7126 . daanc6(*), daanc6l(*)
7137 datmp(i) = daanc6(i) + daanc6l(i)
7141 daanc6(i) = datmp(i)
7147!||====================================================================
7148!||
conversion7 ../engine/source/mpi/interfaces/spmd_i7tool.f
7157#include "implicit_f.inc"
7161 INTEGER LEN, SIZ, IREM(1,*)
7173 irem(1,i) = nint(xrem_dp(10,i))
7177 xrem(j,i) = xrem_dp(j,i)
7184!||====================================================================
7186!||--- called by ------------------------------------------------------
7195#include "implicit_f.inc"
7199 INTEGER LEN, SIZ, IREM(2,*)
7211 irem(1,i) = nint(xrem_dp(9,i))
7212 irem(2,i) = nint(xrem_dp(17,i))
7216 xrem(j,i) = xrem_dp(j,i)
7241#include "implicit_f.inc"
7245#include "com04_c.inc"
7246#include "scr18_c.inc"
7247#include "param_c.inc"
7251 INTEGER IPARI(NPARI,*),(*),COUNT_REMSLVE(*)
7252 INTEGER,
INTENT(IN) :: NODADT_THERM
7256 INTEGER NI,ITYP,LSKYFI,IERROR1,INTTH,IERROR
7263 IF (count_remslv(ni) >
nlskyfi(ni))
THEN
7267 nlskyfi(ni)=nint(count_remslv(ni)*1.20d0)
7271 IF(ityp==7.OR.ityp==10.OR.ityp==20.OR.
7272 . ityp==22.OR.ityp==23.OR.ityp==24.OR.
7274 intth = ipari(47,ni)
7279 ALLOCATE(
iskyfi(ni)%P(lskyfi),stat=ierror1)
7280 ierror = ierror + ierror1
7282 ALLOCATE(
fskyfi(ni)%P(4,lskyfi),stat=ierror1)
7283 ierror = ierror + ierror1
7285 ALLOCATE(
fskyfi(ni)%P(5,lskyfi),stat=ierror1)
7286 ierror = ierror + ierror1
7288 IF(ityp==7 .OR. ityp==25)
THEN
7291 ALLOCATE(
ftheskyfi(ni)%P(lskyfi),stat=ierror1)
7292 ierror = ierror + ierror1
7294 IF(intth /= 0 .AND. nodadt_therm ==1)
THEN
7296 ALLOCATE(
condnskyfi(ni)%P(lskyfi),stat=ierror1)
7297 ierror = ierror + ierror1
7300 ELSEIF(ityp==11)
THEN
7301 intth = ipari(47,ni)
7305 ALLOCATE(
iskyfi(ni)%P(lskyfi),stat=ierror1)
7306 ierror = ierror + ierror1
7309 ALLOCATE(
ftheskyfi(ni)%P(2*lskyfi),stat=ierror1)
7310 ierror = ierror + ierror1
7313 ALLOCATE(
fskyfi(ni)%P(8,lskyfi),stat=ierror1)
7314 ierror = ierror + ierror1
7316 ALLOCATE(
fskyfi(ni)%P(10,lskyfi),stat=ierror1)
7317 ierror = ierror + ierror1
7319 IF(intth /= 0 .AND. nodadt_therm ==1)
THEN
7321 ALLOCATE(
condnskyfi(ni)%P(2*lskyfi),stat=ierror1)
7322 ierror = ierror + ierror1
7326 IF( count_remslve(ni) >
nlskyfie(ni) )
THEN
7327 nlskyfie(ni) = nint(1.2d0 * count_remslve(ni))
7328 lskyfi = count_remslve(ni)
7331 ALLOCATE(
iskyfie(ni)%P(lskyfi),stat=ierror1)
7334 ALLOCATE(
fskyfie(ni)%P(8,lskyfi),stat=ierror1)
7335 ierror = ierror + ierror1
7337 ALLOCATE(
fskyfie(ni)%P(8,lskyfi),stat=ierror1)
7338 ierror = ierror + ierror1
7343 CALL ancmsg(msgid=20,anmode=aninfo)
7365#include
"implicit_f.inc"
7403#include "implicit_f.inc"
subroutine ibcoff(ibc, icodt)
subroutine inter_sort_07(timers, ipari, x, nin, itask, isendto, ircvfrom, retri, itab, nrtm_t, renum, renum_siz, nsnfiold, eshift, multi_fvm, intbuf_tab, h3d_data, inter_struct, sort_comm, intheat, idt_therm, nodadt_therm)
void my_orders(int *mode, int *iwork, int *data, int *index, int *n, int *irecl)
type(int_pointer), dimension(:), allocatable nsnsi_sav
type(int_pointer), dimension(:), allocatable nmvsi
type(int_pointer), dimension(:), allocatable nmnsi
type(int_pointer), dimension(:), allocatable nmvfi
type(int_pointer), dimension(:), allocatable nmnfi
type(real_pointer), dimension(:), allocatable nmtemp
type(int_pointer), dimension(:), allocatable tempnod
type(int_pointer), dimension(:), allocatable nisubsfie
type(int_pointer), dimension(:), allocatable inflg_subsfie
type(real_pointer), dimension(:), allocatable gape_l_fie
type(int_pointer), dimension(:), allocatable lisubsfie
type(real4_pointer3), dimension(:), allocatable edg_bisector_fie
type(real4_pointer3), dimension(:), allocatable vtx_bisector_fie
type(int_pointer), dimension(:), allocatable addsubsfie
type(real_pointer3), dimension(:), allocatable x_seg_fie
type(int_pointer), dimension(:), allocatable edge_fi
type(int_pointer2), dimension(:), allocatable ledge_fie
type(int_pointer), dimension(:), allocatable ispt2_fi
type(real_pointer2), dimension(:), allocatable stif_oldfi
type(real_pointer2), dimension(:), allocatable dxancfi
type(real_pointer2), dimension(:), allocatable stnfi17
type(real_pointer), dimension(:), allocatable gapfie
type(real_pointer2), dimension(:), allocatable eminxfi
type(real_pointer2), dimension(:), allocatable secnd_frfi
type(real_pointer), dimension(:), allocatable stif_msdt_fi
type(real_pointer), dimension(:), allocatable condnfi
type(int_pointer), dimension(:), allocatable is2pt_fi
type(int_pointer), dimension(:), allocatable inflg_subsfi
type(int_pointer), dimension(:), allocatable iskew_fi
type(real_pointer), dimension(:), allocatable ftheskyfi
type(int_pointer), dimension(:), allocatable isegpt_fi
type(real_pointer2), dimension(:), allocatable forneqsfi
type(real_pointer2), dimension(:), allocatable vfie
type(real_pointer), dimension(:), allocatable time_sfi
type(real_pointer3), dimension(:), allocatable afi17
type(real_pointer2), dimension(:), allocatable daancfi
type(real_pointer2), dimension(:), allocatable fskyfie
type(real_pointer2), dimension(:), allocatable i18kafi
type(real_pointer2), dimension(:), allocatable vfi
type(real_pointer2), dimension(:), allocatable fnconti
type(int_pointer2), dimension(:), allocatable irtlm_fi
type(real_pointer2), dimension(:), allocatable frotsfi
type(int_pointer), dimension(:), allocatable matsfi
type(int_pointer2), dimension(:), allocatable is2se_fi
type(int_pointer), dimension(:), allocatable ipartfric_fie
type(real_pointer), dimension(:), allocatable efricgfi
type(real_pointer2), dimension(:), allocatable dvancfi
type(real_pointer), dimension(:), allocatable tempfi
type(int_pointer), dimension(:), allocatable candf_si
type(real_pointer), dimension(:), allocatable stifi
type(real_pointer3), dimension(:), allocatable xfi17
type(real_pointer3), dimension(:), allocatable vfi17
type(real_pointer2), dimension(:), allocatable ksfi
type(int_pointer), dimension(:), allocatable nsvsi
type(real_pointer2), dimension(:), allocatable xfie
type(int_pointer), dimension(:), allocatable iskyfie
type(int_pointer), dimension(:), allocatable kremnor_fi
type(real_pointer), dimension(:), allocatable gap_lfi
type(int_pointer), dimension(:), allocatable nsnfie
type(real_pointer2), dimension(:), allocatable dvancfie
type(real_pointer2), dimension(:), allocatable mtfi_a
type(r8_pointer3), dimension(:), allocatable daanc6fi
type(real_pointer), dimension(:), allocatable condnskyfi
type(real_pointer2), dimension(:), allocatable penfi
type(int_pointer), dimension(:), allocatable nsnsie
type(real_pointer), dimension(:), allocatable alphakfi
type(real_pointer), dimension(:), allocatable stnfi
type(int_pointer), dimension(:), allocatable lisubsfi
type(int_pointer2), dimension(:), allocatable irtse_fi
type(real_pointer2), dimension(:), allocatable penfia
type(real_pointer2), dimension(:), allocatable afi
type(int_pointer), dimension(:), allocatable nsvsie
type(real_pointer), dimension(:), allocatable mtfi_pene
type(int_pointer), dimension(:), allocatable nodamsfi
type(int_pointer), dimension(:), allocatable nsnsi
type(real_pointer2), dimension(:), allocatable fskyfi
type(int_pointer), dimension(:), allocatable nisubsfi
type(int_pointer), dimension(:), allocatable remnor_fi
type(real_pointer), dimension(:), allocatable stifie
type(int_pointer), dimension(:), allocatable pmainfi
type(real_pointer), dimension(:), allocatable stnfie
type(real_pointer), dimension(:), allocatable gapfi
type(int_pointer), dimension(:), allocatable procamsfie
type(int_pointer), dimension(:), allocatable nodnxfi
type(int_pointer), dimension(:), allocatable nodnxfie
type(int_pointer), dimension(:), allocatable nsvfi
type(real_pointer), dimension(:), allocatable intareanfi
type(real_pointer), dimension(:), allocatable areasfi
type(int_pointer), dimension(:), allocatable icodt_fi
type(int_pointer), dimension(:), allocatable iskyfi
type(int_pointer), dimension(:), allocatable isedge_fi
type(real_pointer2), dimension(:), allocatable afie
type(int_pointer), dimension(:), allocatable nbinflfi
type(int_pointer), dimension(:), allocatable nsnfi_sav
type(real_pointer), dimension(:), allocatable efricfi
type(real_pointer), dimension(:), allocatable msfi
type(real_pointer), dimension(:), allocatable diag_smsfi
type(real_pointer), dimension(:), allocatable vscfi
type(real_pointer2), dimension(:), allocatable mtfi_n
type(int_pointer), dimension(:), allocatable nsvfie
type(int_pointer), dimension(:), allocatable addsubsfi
type(int_pointer), dimension(:), allocatable ipartfricsfi
type(int_pointer), dimension(:), allocatable nodamsfie
type(real_pointer), dimension(:), allocatable t2fac_sms_fi
type(r8_pointer3), dimension(:), allocatable daanc6fie
type(real_pointer2), dimension(:), allocatable xfi
type(int_pointer), dimension(:), allocatable kinfi
type(real_pointer), dimension(:), allocatable vscfie
type(real_pointer), dimension(:), allocatable diag_smsfie
type(real_pointer), dimension(:), allocatable alphakfie
type(int_pointer), dimension(:), allocatable itafie
integer, dimension(:), allocatable nlskyfie
type(real_pointer2), dimension(:), allocatable penfie
integer, dimension(:), allocatable nsnfi_flag
type(real_pointer2), dimension(:), allocatable dxancfie
type(real_pointer), dimension(:), allocatable stife_msdt_fi
type(real_pointer), dimension(:), allocatable mtfi_penemin
integer, dimension(:), allocatable nlskyfi
type(real_pointer), dimension(:), allocatable fthefi
type(real_pointer2), dimension(:), allocatable ftconti
type(int_pointer), dimension(:), allocatable procamsfi
type(int_pointer2), dimension(:), allocatable islide_fi
type(real_pointer2), dimension(:), allocatable pene_oldfi
type(int_pointer), dimension(:), allocatable icont_i_fi
type(real_pointer2), dimension(:), allocatable mtfi_v
type(int_pointer), dimension(:), allocatable tagncontfi
type(int_pointer), dimension(:), allocatable nsnfi
type(int_pointer), dimension(:), allocatable if_adhfi
type(int_pointer2), dimension(:), allocatable t2main_sms_fi
integer, dimension(:,:), allocatable irem
type(real_pointer2), dimension(:), allocatable daancfie
type(real_pointer), dimension(:), allocatable msfie
type(int_pointer), dimension(:), allocatable itafi
type(real_pointer2), dimension(:), allocatable penfiae
subroutine spmd_exch_da20(intbuf_tab, ipari, iad_elem, fr_elem, len20, nbint20, lenr, intlist, nbintc)
subroutine spmd_i7fcom_pon(ipari, intlist, nbintc, niskyfi, icodt, secfcum, nstrf, icontact, fcont, igrbric, ixs, ixs16, niskyfie, nbint20, iflag, intbuf_tab, sfskyi, sisky, h3d_data, multi_fvm, tagncont, kloadpinter, loadpinter, loadp_hyd_inter, fsav, interfaces, glob_therm)
subroutine spmd_ibcast(tabi, tabr, n1, n2, from, add)
subroutine spmd_split_comm(isactive, comm_local)
subroutine spmd_split_comm_inter(nbintc, intlist, ipari, isendto, ircvfrom, sort_comm, need_comm_int25_solid_erosion, comm_int25_solid_erosion)
subroutine spmd_split_comm_joint()
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
subroutine write_db(a, n)
void write_i_c(int *w, int *len)
void write_r_c(float *w, int *len)
void read_i_c(int *w, int *len)
void read_r_c(float *w, int *len)
subroutine wrrestp(elements, nodes, af, iaf, ich, addcne, elbuf_tab, xfem_tab, intbuf_tab, multi_fvm, mat_elem, h3d_data, intbuf_fric_tab, subset, pinch_data, ale_connectivity, t_monvol, sensors, ebcs_tab, dynain_data, user_windows, output, interfaces, loads, python, names_and_titles, eigipm, eigibuf, eigrpm, neipm, leibuf, nerpm, iflow, rflow, liflow, lrflow, impbuf_tab, impl_s, impl_s0, mcp, temp, forneqs, unitab, stack, ndrape, drape_sh3n, drape_sh4n, drapeg, restsize, skews, glob_therm, pblast, rbe3)