Télécharger elpdm2.eso

Retour à la liste

Numérotation des lignes :

  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.  
  118.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales