1 : C ELPDM2 SOURCE PV 09/03/12 21:20:58 6325 2 : SUBROUTINE ELPDM2(SBORD,SMAT 3 : & ,MCHELM) 4 : IMPLICIT INTEGER(I-N) 5 : IMPLICIT REAL*8(A-B,D-H,O-Z) 6 : IMPLICIT COMPLEX*16(C) 7 : ************************************************************************ 8 : * 9 : * PLAQUES PAR EQUATION INTEGRALE: 10 : * 11 : * REMPLISSAGE DES TERMES C.L. DE LA MATRICE CMA1(NS4,NS4) 12 : * ET DU SECOND MEMBRE 13 : * | A1 B1 C1 D1 A2.............. 14 : * | E1 F1 G1 H1 A2............. 15 : * | -------C.L.------------ 16 : * | -------C.L.------------ 17 : * | ... A1 B1 C1 D1 ... 18 : * | .... E1 F1 G1 H1.... 19 : * | -------C.L.-- 20 : * | -------C.L.-- 21 : * 22 : * 23 : * A1 TEL QUE A1 + A2 + .... AN = 0. 24 : * 25 : * 26 : * 27 : * 28 : * 29 : ************************************************************************ 30 : * 31 : -INC SMELEME 32 : -INC SMCHAML 33 : * 34 : SEGMENT SBORD 35 : REAL*8 XBORD(15,NS) 36 : INTEGER IBORD (2 ,NS) 37 : ENDSEGMENT 38 : SEGMENT SMAT 39 : COMPLEX*16 CMA1(NS4,NS4) 40 : COMPLEX*16 CSM (NS4) 41 : COMPLEX*16 CSOM2 (NS) 42 : COMPLEX*16 CRE (NS4) 43 : ENDSEGMENT 44 : * 45 : CHARACTER*8 NOMC0(4) 46 : NS = XBORD(/2) 47 : NS4= CSM (/1) 48 : 49 : * 50 : *--1. MISE A 0 DES LIGNES CONCERNEES 51 : * 52 : DO 100 I=1,NS 53 : DO 110 J=1,NS4 54 : CMA1(4*(I-1) + 3,J) =(0D0,0d0) 55 : CMA1(4*(I-1) + 4,J) =(0D0,0d0) 56 : 110 CONTINUE 57 : 100 CONTINUE 58 : * 59 : *--2. ON BALAIE LE CONTENU DU CHAMELEME 60 : * 61 : NOMC0(1) = 'WW' 62 : NOMC0(2) = 'WN' 63 : NOMC0(3) = 'MN' 64 : NOMC0(4) = 'KN' 65 : SEGACT MCHELM 66 : N1 = IMACHE(/1) 67 : 68 : DO 200 IM0 =1,N1 69 : IPT1 =IMACHE(IM0) 70 : MCHAML = ICHAML(IM0) 71 : SEGACT IPT1 72 : SEGACT MCHAML 73 : 74 : DO 210 I=1,IPT1.NUM(/2) 75 : I1 = IPT1.NUM(1,I) 76 : I2 = IPT1.NUM(3,I) 77 : 78 : DO 211 J=1,NS 79 : J1 = IBORD(1,J) 80 : J2 = IBORD(2,J) 81 : IF ( (J1.EQ.I1).AND.(J2.EQ.I2)) THEN 82 : JELEM = J 83 : GOTO 212 84 : ENDIF 85 : 211 CONTINUE 86 : 212 CONTINUE 87 : 88 : ILL0= 4*(JELEM - 1) 89 : DO 220 JC=1,2 90 : ILL = ILL0 + 2 + JC 91 : MELVAL = IELVAL(JC) 92 : SEGACT MELVAL 93 : IF ( VELCHE(/2) .EQ. 1) THEN 94 : CSM (ILL) =VELCHE(1,1) 95 : ELSE 96 : CSM (ILL) =VELCHE(1,I) 97 : ENDIF 98 : SEGDES MELVAL 99 : DO 221 KC=1,4 100 : IF (NOMCHE(JC).EQ.NOMC0(KC)) THEN 101 : CMA1(ILL,ILL0 + KC ) = (1d0 , 0d0) 102 : ENDIF 103 : 221 CONTINUE 104 : 220 CONTINUE 105 : 210 CONTINUE 106 : SEGDES IPT1 107 : SEGDES MCHAML 108 : 200 CONTINUE 109 : SEGDES MCHELM 110 : 111 : 112 : * 113 : RETURN 114 : END 115 : 116 : 117 :
© Cast3M 2003 - All rights reserved.
Disclaimer