C ETALCH SOURCE CHAT 05/01/12 23:46:02 5004 SUBROUTINE ETALCH(MCHPOI,IINC,ICPR,ICONTR,MVA,IPB,NPR2) C C======================================================================= C CE SUBROUTINE ETALE LE CHPOINT DANS LE TABLEAU MVA C ***** ENTREES ***** C IINC CONTIENT LA LISTE DES COMPOSANTES DU TABLEAU C ICPR(I)=K : LE POINT I EST LE KIEME DU TABLEAU MVA C ON VERIFIE QUE LES POINTS DE MCHPOI MUNIS DE LEUR COMPOSANTES SONT C CONTENUS DANS MVA, C.A.D. QUE MCONTR(I,J)=1 POUR CES POINTS C ***** SORTIES ***** C LES VALEURS DE MVA DANS LE CHPOINT C EVENTUELLEMENT IPB LE TABLEAU DE CORRESPONDANCE DES POINTS C IPB(I)=II LE IEME POINT DU CHPOINT EST A LA PLACE II DANS MVA C NPR2 = NOMBRE DE POINTS DU MCHPOI C C ATTENTION : SEGACT MVA,IPB ET MISE A ZERO DE CES 2 TABLEAUX AVANT C UTILISATION. C======================================================================= IMPLICIT INTEGER(I-N) -INC CCOPTIO -INC SMCHPOI -INC SMELEME SEGMENT ICPR(XCOOR(/1)/(IDIM+1)) SEGMENT IINC CHARACTER*4 CIINC(0) ENDSEGMENT SEGMENT/MVA/(VA(NNI1,IPR1)*D),MVA1.MVA SEGMENT/ICONTR/(MCONTR(NNI1,IPR1)) SEGMENT IPB(IPR1) C SEGACT ICPR SEGACT IINC SEGACT MVA*MOD SEGACT ICONTR IF(IPB.NE.0) SEGACT IPB*MOD NNI1=VA(/1) IPR1=VA(/2) NPR2=0 C DO 1 J=1,IPR1 DO 1 K=1,NNI1 VA(K,J)=0. 1 CONTINUE C IF(IPB.EQ.0) GO TO 2 DO 3 J=1,IPR1 IPB(J)=0 3 CONTINUE 2 CONTINUE SEGACT MCHPOI NSOUP=IPCHP(/1) DO 60 ISOUP=1,NSOUP MSOUPO=IPCHP(ISOUP) SEGACT MSOUPO MELEME=IGEOC SEGACT MELEME N2=NUM(/2) MPOVAL=IPOVAL SEGACT MPOVAL DO 61 J=1,NOCOMP(/2) DO 62 K=1,NNI1 IF(NOCOMP(J).EQ.CIINC(K)) GO TO 63 62 CONTINUE IJ=1 GOTO 66 63 CONTINUE IF(IPB.NE.0) GO TO 65 DO 64 IJ=1,N2 KI=ICPR(NUM(1,IJ)) IF(MCONTR(K,KI).NE.1) GO TO 66 VA(K,KI)=VPOCHA(IJ,J) 64 CONTINUE GO TO 61 65 CONTINUE DO 75 IJ=1,N2 KI=ICPR(NUM(1,IJ)) IF(MCONTR(K,KI).NE.1) GO TO 66 IF(J.EQ.1) IPB(NPR2+IJ)=KI VA(K,KI)=VPOCHA(IJ,J) 75 CONTINUE GO TO 61 66 CONTINUE MOTERR(1:4)=NOCOMP(J) INTERR(1)=NUM(1,IJ) CALL ERREUR(140) C INCOMPATIBILITE ENTRE LES POINTS ET COMPOSANTES DES 2 CHPOINTS GO TO 5000 61 CONTINUE NPR2=NPR2+N2 SEGDES MELEME SEGDES MSOUPO SEGDES MPOVAL 60 CONTINUE SEGDES MCHPOI C SEGDES ICPR,IINC SEGDES ICONTR SEGDES MVA IF(IPB.NE.0) SEGDES IPB 5000 CONTINUE RETURN END