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 200 : CALL ERREUR(5) 201 : ENDIF 202 : RETURN 203 : END 204 : 205 : 206 :
© Cast3M 2003 - All rights reserved.
Disclaimer