Télécharger bnpor2.eso

Retour à la liste

Numérotation des lignes :

  1. C BNPOR2 SOURCE CHAT 05/01/12 21:41:44 5004
  2. SUBROUTINE BNPOR2 (YGENE,NCOT,IELE)
  3. *______________________________________________________________________
  4. *
  5. * POUR PENALISATION (MILIEUX POREUX) APPELE PAR RIGI2
  6. *______________________________________________________________________
  7. *
  8. IMPLICIT INTEGER(I-N)
  9. IMPLICIT REAL*8(A-H,O-Z)
  10. *
  11. -INC CCOPTIO
  12. *
  13. DIMENSION YGENE(NCOT,*)
  14. *
  15. * CONSTRUCTION DU TABLEAU SELON LES ELEMENTS
  16. * LES SOMMETS SONT RANGES AVANT LES MILIEUX
  17. *
  18. * T6
  19. *
  20. IF(IELE.EQ.6) THEN
  21. YGENE(1,1) = 1.D0
  22. YGENE(1,2) = 1.D0
  23. YGENE(1,4) =-2.D0
  24. YGENE(2,2) = 1.D0
  25. YGENE(2,3) = 1.D0
  26. YGENE(2,5) =-2.D0
  27. YGENE(3,3) = 1.D0
  28. YGENE(3,1) = 1.D0
  29. YGENE(3,6) =-2.D0
  30. *
  31. * Q8
  32. *
  33. ELSE IF(IELE.EQ.10) THEN
  34. YGENE(1,1) = 1.D0
  35. YGENE(1,2) = 1.D0
  36. YGENE(1,5) =-2.D0
  37. YGENE(2,2) = 1.D0
  38. YGENE(2,3) = 1.D0
  39. YGENE(2,6) =-2.D0
  40. YGENE(3,3) = 1.D0
  41. YGENE(3,4) = 1.D0
  42. YGENE(3,7) =-2.D0
  43. YGENE(4,4) = 1.D0
  44. YGENE(4,1) = 1.D0
  45. YGENE(4,8) =-2.D0
  46. *
  47. * C20
  48. *
  49. ELSE IF(IELE.EQ.15) THEN
  50. YGENE(1,1) = 1.D0
  51. YGENE(1,2) = 1.D0
  52. YGENE(1,9) =-2.D0
  53. YGENE(2,2) = 1.D0
  54. YGENE(2,3) = 1.D0
  55. YGENE(2,10) =-2.D0
  56. YGENE(3,3) = 1.D0
  57. YGENE(3,4) = 1.D0
  58. YGENE(3,11) =-2.D0
  59. YGENE(4,4) = 1.D0
  60. YGENE(4,1) = 1.D0
  61. YGENE(4,12) =-2.D0
  62. YGENE(5,5) = 1.D0
  63. YGENE(5,6) = 1.D0
  64. YGENE(5,17) =-2.D0
  65. YGENE(6,6) = 1.D0
  66. YGENE(6,7) = 1.D0
  67. YGENE(6,18) =-2.D0
  68. YGENE(7,7) = 1.D0
  69. YGENE(7,8) = 1.D0
  70. YGENE(7,19) =-2.D0
  71. YGENE(8,8) = 1.D0
  72. YGENE(8,5) = 1.D0
  73. YGENE(8,20) =-2.D0
  74. YGENE(9,1) = 1.D0
  75. YGENE(9,5) = 1.D0
  76. YGENE(9,13) =-2.D0
  77. YGENE(10,2) = 1.D0
  78. YGENE(10,6) = 1.D0
  79. YGENE(10,14)=-2.D0
  80. YGENE(11,3) = 1.D0
  81. YGENE(11,7) = 1.D0
  82. YGENE(11,15)=-2.D0
  83. YGENE(12,4) = 1.D0
  84. YGENE(12,8) = 1.D0
  85. YGENE(12,16)=-2.D0
  86. *
  87. * P15
  88. *
  89. ELSE IF(IELE.EQ.17) THEN
  90. YGENE(1,1) = 1.D0
  91. YGENE(1,2) = 1.D0
  92. YGENE(1,7) =-2.D0
  93. YGENE(2,2) = 1.D0
  94. YGENE(2,3) = 1.D0
  95. YGENE(2,8) =-2.D0
  96. YGENE(3,3) = 1.D0
  97. YGENE(3,1) = 1.D0
  98. YGENE(3,9) =-2.D0
  99. YGENE(4,4) = 1.D0
  100. YGENE(4,5) = 1.D0
  101. YGENE(4,13)=-2.D0
  102. YGENE(5,5) = 1.D0
  103. YGENE(5,6) = 1.D0
  104. YGENE(5,14)=-2.D0
  105. YGENE(6,6) = 1.D0
  106. YGENE(6,4) = 1.D0
  107. YGENE(6,15)=-2.D0
  108. YGENE(7,1) = 1.D0
  109. YGENE(7,4) = 1.D0
  110. YGENE(7,10)=-2.D0
  111. YGENE(8,2) = 1.D0
  112. YGENE(8,5) = 1.D0
  113. YGENE(8,11)=-2.D0
  114. YGENE(9,3) = 1.D0
  115. YGENE(9,6) = 1.D0
  116. YGENE(9,12)=-2.D0
  117. *
  118. * T10
  119. *
  120. ELSE IF(IELE.EQ.24) THEN
  121. YGENE(1,1) = 1.D0
  122. YGENE(1,2) = 1.D0
  123. YGENE(1,5) =-2.D0
  124. YGENE(2,2) = 1.D0
  125. YGENE(2,3) = 1.D0
  126. YGENE(2,6) =-2.D0
  127. YGENE(3,3) = 1.D0
  128. YGENE(3,1) = 1.D0
  129. YGENE(3,7) =-2.D0
  130. YGENE(4,1) = 1.D0
  131. YGENE(4,4) = 1.D0
  132. YGENE(4,8) =-2.D0
  133. YGENE(5,2) = 1.D0
  134. YGENE(5,4) = 1.D0
  135. YGENE(5,9) =-2.D0
  136. YGENE(6,3) = 1.D0
  137. YGENE(6,4) = 1.D0
  138. YGENE(6,10)=-2.D0
  139. *
  140. * RP3
  141. *
  142. ELSE IF(IELE.EQ.29) THEN
  143. YGENE(1,1) = 1.D0
  144. YGENE(1,2) = 1.D0
  145. YGENE(1,7) =-2.D0
  146. YGENE(2,3) = 1.D0
  147. YGENE(2,4) = 1.D0
  148. YGENE(2,8) =-2.D0
  149. *
  150. * LP6
  151. *
  152. ELSE IF(IELE.EQ.30) THEN
  153. YGENE(1,1) = 1.D0
  154. YGENE(1,2) = 1.D0
  155. YGENE(1,10)=-2.D0
  156. YGENE(2,2) = 1.D0
  157. YGENE(2,3) = 1.D0
  158. YGENE(2,11)=-2.D0
  159. YGENE(3,3) = 1.D0
  160. YGENE(3,1) = 1.D0
  161. YGENE(3,12)=-2.D0
  162. YGENE(4,4) = 1.D0
  163. YGENE(4,5) = 1.D0
  164. YGENE(4,13)=-2.D0
  165. YGENE(5,5) = 1.D0
  166. YGENE(5,6) = 1.D0
  167. YGENE(5,14)=-2.D0
  168. YGENE(6,6) = 1.D0
  169. YGENE(6,4) = 1.D0
  170. YGENE(6,15)=-2.D0
  171. *
  172. * LP8
  173. *
  174. ELSE IF(IELE.EQ.31) THEN
  175. YGENE(1,1) = 1.D0
  176. YGENE(1,2) = 1.D0
  177. YGENE(1,13)=-2.D0
  178. YGENE(2,2) = 1.D0
  179. YGENE(2,3) = 1.D0
  180. YGENE(2,14)=-2.D0
  181. YGENE(3,3) = 1.D0
  182. YGENE(3,4) = 1.D0
  183. YGENE(3,15)=-2.D0
  184. YGENE(4,4) = 1.D0
  185. YGENE(4,1) = 1.D0
  186. YGENE(4,16)=-2.D0
  187. YGENE(5,5) = 1.D0
  188. YGENE(5,6) = 1.D0
  189. YGENE(5,17)=-2.D0
  190. YGENE(6,6) = 1.D0
  191. YGENE(6,7) = 1.D0
  192. YGENE(6,18)=-2.D0
  193. YGENE(7,7) = 1.D0
  194. YGENE(7,8) = 1.D0
  195. YGENE(7,19)=-2.D0
  196. YGENE(8,8) = 1.D0
  197. YGENE(8,5) = 1.D0
  198. YGENE(8,20)=-2.D0
  199. *
  200. ELSE
  201. CALL ERREUR(5)
  202. ENDIF
  203. RETURN
  204. END
  205.  
  206.  
  207.  
  208.  

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