1 : C SHAPEQ SOURCE GOUNAND 08/06/10 21:15:06 6125 2 : SUBROUTINE SHAPEQ(QSI,ETA,DZE,IELE,SHP,IRET) 3 : C======================================================================= 4 : C 5 : C FONCTIONS DE FORME TRIDIMENSIONNELLES QUAF 6 : C EBERSOLT 7 : C======================================================================= 8 : IMPLICIT INTEGER(I-N) 9 : IMPLICIT REAL*8(A-H,O-Z) 10 : PARAMETER (NDIM=3) 11 : DIMENSION SHP(6,*) 12 : * Des constantes 13 : REAL*8 UN 14 : PARAMETER (UN=1.D0) 15 : 16 : * Les matrices de passage : base polynomiale -> base des fonctions de 17 : * forme 18 : REAL*8 PNM27(27,27) 19 : REAL*8 PNM21(21,21) 20 : REAL*8 PNM15(15,15) 21 : * 22 : * REAL*8 PM271(27),PM272(27) 23 : * EQUIVALENCE (PM27(1,1),PM271(1)),(PM27(1,2),PM272(1)) 24 : * 25 : * Travail (on ne peut pas utiliser de segments ici : fortran pur) 26 : * 27 : REAL*8 XCO1(NDIM) 28 : REAL*8 VECTP(27) 29 : * Pour le cube 30 : PARAMETER (NMMO27=1) 31 : INTEGER NMON27(27) 32 : REAL*8 CMON27(NMMO27,27) 33 : INTEGER XMON27(NDIM,NMMO27,27) 34 : * Pour le prisme 35 : PARAMETER (NMMO21=4) 36 : INTEGER NMON21(21) 37 : REAL*8 CMON21(NMMO21,21) 38 : INTEGER XMON21(NDIM,NMMO21,21) 39 : * Pour le tétra 40 : PARAMETER (NMMO15=4) 41 : INTEGER NMON15(15) 42 : REAL*8 CMON15(NMMO15,15) 43 : INTEGER XMON15(NDIM,NMMO15,15) 44 : * 45 : INTEGER ORDDER(NDIM) 46 : * 47 : DATA ((PNM27(I,J),I=1, 27),J= 1, 1)/ 48 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 49 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 50 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 51 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 52 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 53 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 54 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 55 : $ 1.2500000000000000D-01, 0.0000000000000000D+00, 56 : $ 1.2500000000000000D-01, -1.2500000000000000D-01, 57 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 58 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 59 : $ 1.2500000000000000D-01, -1.2500000000000000D-01, 60 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 61 : $ 1.2500000000000000D-01 62 : $ / 63 : DATA ((PNM27(I,J),I=1, 27),J= 2, 2)/ 64 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 65 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 66 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 67 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 68 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 69 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 70 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 71 : $ -2.5000000000000000D-01, -2.5000000000000000D-01, 72 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 73 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 74 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 75 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 76 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 77 : $ -2.5000000000000000D-01 78 : $ / 79 : DATA ((PNM27(I,J),I=1, 27),J= 3, 3)/ 80 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 81 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 82 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 83 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 84 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 85 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 86 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 87 : $ 1.2500000000000000D-01, 0.0000000000000000D+00, 88 : $ -1.2500000000000000D-01, -1.2500000000000000D-01, 89 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 90 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 91 : $ -1.2500000000000000D-01, -1.2500000000000000D-01, 92 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 93 : $ 1.2500000000000000D-01 94 : $ / 95 : DATA ((PNM27(I,J),I=1, 27),J= 4, 4)/ 96 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 97 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 98 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 99 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 100 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 101 : $ -2.5000000000000000D-01, -2.5000000000000000D-01, 102 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 103 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 104 : $ 2.5000000000000000D-01, 2.5000000000000000D-01, 105 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 106 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 107 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 108 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 109 : $ -2.5000000000000000D-01 110 : $ / 111 : DATA ((PNM27(I,J),I=1, 27),J= 5, 5)/ 112 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 113 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 114 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 115 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 116 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 117 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 118 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 119 : $ -1.2500000000000000D-01, 0.0000000000000000D+00, 120 : $ -1.2500000000000000D-01, -1.2500000000000000D-01, 121 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 122 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 123 : $ 1.2500000000000000D-01, 1.2500000000000000D-01, 124 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 125 : $ 1.2500000000000000D-01 126 : $ / 127 : DATA ((PNM27(I,J),I=1, 27),J= 6, 6)/ 128 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 129 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 130 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 131 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 132 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 133 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 134 : $ -2.5000000000000000D-01, 0.0000000000000000D+00, 135 : $ 2.5000000000000000D-01, -2.5000000000000000D-01, 136 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 137 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 138 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 139 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 140 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 141 : $ -2.5000000000000000D-01 142 : $ / 143 : DATA ((PNM27(I,J),I=1, 27),J= 7, 7)/ 144 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 145 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 146 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 147 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 148 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 149 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 150 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 151 : $ -1.2500000000000000D-01, 0.0000000000000000D+00, 152 : $ 1.2500000000000000D-01, -1.2500000000000000D-01, 153 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 154 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 155 : $ -1.2500000000000000D-01, 1.2500000000000000D-01, 156 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 157 : $ 1.2500000000000000D-01 158 : $ / 159 : DATA ((PNM27(I,J),I=1, 27),J= 8, 8)/ 160 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 161 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 162 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 163 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 164 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 165 : $ 2.5000000000000000D-01, -2.5000000000000000D-01, 166 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 167 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 168 : $ -2.5000000000000000D-01, 2.5000000000000000D-01, 169 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 170 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 171 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 172 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 173 : $ -2.5000000000000000D-01 174 : $ / 175 : DATA ((PNM27(I,J),I=1, 27),J= 9, 9)/ 176 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 177 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 178 : $ 2.5000000000000000D-01, -2.5000000000000000D-01, 179 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 180 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 181 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 182 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 183 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 184 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 185 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 186 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 187 : $ -2.5000000000000000D-01, 2.5000000000000000D-01, 188 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 189 : $ -2.5000000000000000D-01 190 : $ / 191 : DATA ((PNM27(I,J),I=1, 27),J= 10, 10)/ 192 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 193 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 194 : $ -2.5000000000000000D-01, -2.5000000000000000D-01, 195 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 196 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 197 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 198 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 199 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 200 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 201 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 202 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 203 : $ 2.5000000000000000D-01, 2.5000000000000000D-01, 204 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 205 : $ -2.5000000000000000D-01 206 : $ / 207 : DATA ((PNM27(I,J),I=1, 27),J= 11, 11)/ 208 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 209 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 210 : $ 2.5000000000000000D-01, 2.5000000000000000D-01, 211 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 212 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 213 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 214 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 215 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 216 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 217 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 218 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 219 : $ -2.5000000000000000D-01, -2.5000000000000000D-01, 220 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 221 : $ -2.5000000000000000D-01 222 : $ / 223 : DATA ((PNM27(I,J),I=1, 27),J= 12, 12)/ 224 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 225 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 226 : $ -2.5000000000000000D-01, 2.5000000000000000D-01, 227 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 228 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 229 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 230 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 231 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 232 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 233 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 234 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 235 : $ 2.5000000000000000D-01, -2.5000000000000000D-01, 236 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 237 : $ -2.5000000000000000D-01 238 : $ / 239 : DATA ((PNM27(I,J),I=1, 27),J= 13, 13)/ 240 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 241 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 242 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 243 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 244 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 245 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 246 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 247 : $ -1.2500000000000000D-01, 0.0000000000000000D+00, 248 : $ -1.2500000000000000D-01, 1.2500000000000000D-01, 249 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 250 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 251 : $ 1.2500000000000000D-01, -1.2500000000000000D-01, 252 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 253 : $ 1.2500000000000000D-01 254 : $ / 255 : DATA ((PNM27(I,J),I=1, 27),J= 14, 14)/ 256 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 257 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 258 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 259 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 260 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 261 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 262 : $ -2.5000000000000000D-01, 0.0000000000000000D+00, 263 : $ 2.5000000000000000D-01, 2.5000000000000000D-01, 264 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 265 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 266 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 267 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 268 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 269 : $ -2.5000000000000000D-01 270 : $ / 271 : DATA ((PNM27(I,J),I=1, 27),J= 15, 15)/ 272 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 273 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 274 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 275 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 276 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 277 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 278 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 279 : $ -1.2500000000000000D-01, 0.0000000000000000D+00, 280 : $ 1.2500000000000000D-01, 1.2500000000000000D-01, 281 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 282 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 283 : $ -1.2500000000000000D-01, -1.2500000000000000D-01, 284 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 285 : $ 1.2500000000000000D-01 286 : $ / 287 : DATA ((PNM27(I,J),I=1, 27),J= 16, 16)/ 288 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 289 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 290 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 291 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 292 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 293 : $ 2.5000000000000000D-01, 2.5000000000000000D-01, 294 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 295 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 296 : $ -2.5000000000000000D-01, -2.5000000000000000D-01, 297 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 298 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 299 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 300 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 301 : $ -2.5000000000000000D-01 302 : $ / 303 : DATA ((PNM27(I,J),I=1, 27),J= 17, 17)/ 304 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 305 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 306 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 307 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 308 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 309 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 310 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 311 : $ 1.2500000000000000D-01, 0.0000000000000000D+00, 312 : $ 1.2500000000000000D-01, 1.2500000000000000D-01, 313 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 314 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 315 : $ 1.2500000000000000D-01, 1.2500000000000000D-01, 316 : $ 0.0000000000000000D+00, 1.2500000000000000D-01, 317 : $ 1.2500000000000000D-01 318 : $ / 319 : DATA ((PNM27(I,J),I=1, 27),J= 18, 18)/ 320 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 321 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 322 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 323 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 324 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 325 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 326 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 327 : $ -2.5000000000000000D-01, 2.5000000000000000D-01, 328 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 329 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 330 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 331 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 332 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 333 : $ -2.5000000000000000D-01 334 : $ / 335 : DATA ((PNM27(I,J),I=1, 27),J= 19, 19)/ 336 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 337 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 338 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 339 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 340 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 341 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 342 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 343 : $ 1.2500000000000000D-01, 0.0000000000000000D+00, 344 : $ -1.2500000000000000D-01, 1.2500000000000000D-01, 345 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 346 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 347 : $ -1.2500000000000000D-01, 1.2500000000000000D-01, 348 : $ 0.0000000000000000D+00, -1.2500000000000000D-01, 349 : $ 1.2500000000000000D-01 350 : $ / 351 : DATA ((PNM27(I,J),I=1, 27),J= 20, 20)/ 352 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 353 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 354 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 355 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 356 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 357 : $ -2.5000000000000000D-01, 2.5000000000000000D-01, 358 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 359 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 360 : $ 2.5000000000000000D-01, -2.5000000000000000D-01, 361 : $ 0.0000000000000000D+00, -2.5000000000000000D-01, 362 : $ 2.5000000000000000D-01, 0.0000000000000000D+00, 363 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 364 : $ 0.0000000000000000D+00, 2.5000000000000000D-01, 365 : $ -2.5000000000000000D-01 366 : $ / 367 : DATA ((PNM27(I,J),I=1, 27),J= 21, 21)/ 368 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 369 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 370 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 371 : $ 5.0000000000000000D-01, 0.0000000000000000D+00, 372 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 373 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 374 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 375 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 376 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 377 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 378 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 379 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 380 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 381 : $ 5.0000000000000000D-01 382 : $ / 383 : DATA ((PNM27(I,J),I=1, 27),J= 22, 22)/ 384 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 385 : $ 5.0000000000000000D-01, 0.0000000000000000D+00, 386 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 387 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 388 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 389 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 390 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 391 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 392 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 393 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 394 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 395 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 396 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 397 : $ 5.0000000000000000D-01 398 : $ / 399 : DATA ((PNM27(I,J),I=1, 27),J= 23, 23)/ 400 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 401 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 402 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 403 : $ 5.0000000000000000D-01, 0.0000000000000000D+00, 404 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 405 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 406 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 407 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 408 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 409 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 410 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 411 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 412 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 413 : $ 5.0000000000000000D-01 414 : $ / 415 : DATA ((PNM27(I,J),I=1, 27),J= 24, 24)/ 416 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 417 : $ 5.0000000000000000D-01, 0.0000000000000000D+00, 418 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 419 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 420 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 421 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 422 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 423 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 424 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 425 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 426 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 427 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 428 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 429 : $ 5.0000000000000000D-01 430 : $ / 431 : DATA ((PNM27(I,J),I=1, 27),J= 25, 25)/ 432 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 433 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 434 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 435 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 436 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 437 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 438 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 439 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 440 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 441 : $ 5.0000000000000000D-01, 0.0000000000000000D+00, 442 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 443 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 444 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 445 : $ 5.0000000000000000D-01 446 : $ / 447 : DATA ((PNM27(I,J),I=1, 27),J= 26, 26)/ 448 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 449 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 450 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 451 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 452 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 453 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 454 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 455 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 456 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 457 : $ 5.0000000000000000D-01, 0.0000000000000000D+00, 458 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 459 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 460 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 461 : $ 5.0000000000000000D-01 462 : $ / 463 : DATA ((PNM27(I,J),I=1, 27),J= 27, 27)/ 464 : $ 1.0000000000000000D+00, 0.0000000000000000D+00, 465 : $ -1.0000000000000000D+00, 0.0000000000000000D+00, 466 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 467 : $ -1.0000000000000000D+00, 0.0000000000000000D+00, 468 : $ 1.0000000000000000D+00, 0.0000000000000000D+00, 469 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 470 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 471 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 472 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 473 : $ -1.0000000000000000D+00, 0.0000000000000000D+00, 474 : $ 1.0000000000000000D+00, 0.0000000000000000D+00, 475 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 476 : $ 1.0000000000000000D+00, 0.0000000000000000D+00, 477 : $ -1.0000000000000000D+00 478 : $ / 479 : DATA ((PNM21(I,J),I=1, 21),J= 1, 1)/ 480 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 481 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 482 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 483 : $ 1.5543122344752192D-15, -5.0000000000000000D-01, 484 : $ 1.5000000000000000D+00, 1.5000000000000000D+00, 485 : $ -1.0000000000000000D+00, -2.0000000000000000D+00, 486 : $ -1.0000000000000000D+00, -1.4999999999999982D+00, 487 : $ 5.0000000000000000D-01, -1.5000000000000000D+00, 488 : $ -1.5000000000000000D+00, 1.0000000000000000D+00, 489 : $ 2.0000000000000000D+00, 1.0000000000000000D+00, 490 : $ 1.4999999999999967D+00 491 : $ / 492 : DATA ((PNM21(I,J),I=1, 21),J= 2, 2)/ 493 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 494 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 495 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 496 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 497 : $ -2.0000000000000000D+00, 0.0000000000000000D+00, 498 : $ 2.0000000000000000D+00, 2.0000000000000000D+00, 499 : $ 0.0000000000000000D+00, 5.9999999999999991D+00, 500 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 501 : $ 0.0000000000000000D+00, -2.0000000000000000D+00, 502 : $ -2.0000000000000000D+00, 0.0000000000000000D+00, 503 : $ -5.9999999999999991D+00 504 : $ / 505 : DATA ((PNM21(I,J),I=1, 21),J= 3, 3)/ 506 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 507 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 508 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 509 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 510 : $ 5.0000000000000000D-01, 0.0000000000000000D+00, 511 : $ -1.0000000000000000D+00, 0.0000000000000000D+00, 512 : $ 0.0000000000000000D+00, -1.4999999999999998D+00, 513 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 514 : $ 0.0000000000000000D+00, 1.0000000000000000D+00, 515 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 516 : $ 1.4999999999999998D+00 517 : $ / 518 : DATA ((PNM21(I,J),I=1, 21),J= 4, 4)/ 519 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 520 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 521 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 522 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 523 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 524 : $ 0.0000000000000000D+00, -2.0000000000000000D+00, 525 : $ 0.0000000000000000D+00, 5.9999999999999991D+00, 526 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 527 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 528 : $ 2.0000000000000000D+00, 0.0000000000000000D+00, 529 : $ -5.9999999999999991D+00 530 : $ / 531 : DATA ((PNM21(I,J),I=1, 21),J= 5, 5)/ 532 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 533 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 534 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 535 : $ 1.5543122344752192D-15, 0.0000000000000000D+00, 536 : $ 0.0000000000000000D+00, 5.0000000000000000D-01, 537 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 538 : $ -1.0000000000000000D+00, -1.4999999999999998D+00, 539 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 540 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 541 : $ 0.0000000000000000D+00, 1.0000000000000000D+00, 542 : $ 1.4999999999999982D+00 543 : $ / 544 : DATA ((PNM21(I,J),I=1, 21),J= 6, 6)/ 545 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 546 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 547 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 548 : $ -2.6645352591003757D-15, 0.0000000000000000D+00, 549 : $ 0.0000000000000000D+00, -2.0000000000000000D+00, 550 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 551 : $ 2.0000000000000000D+00, 5.9999999999999991D+00, 552 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 553 : $ 2.0000000000000000D+00, 0.0000000000000000D+00, 554 : $ -2.0000000000000000D+00, -2.0000000000000000D+00, 555 : $ -5.9999999999999964D+00 556 : $ / 557 : DATA ((PNM21(I,J),I=1, 21),J= 7, 7)/ 558 : $ 1.0000000000000000D+00, -3.0000000000000000D+00, 559 : $ -3.0000000000000000D+00, 2.0000000000000000D+00, 560 : $ 4.0000000000000000D+00, 2.0000000000000000D+00, 561 : $ 2.9999999999999933D+00, 0.0000000000000000D+00, 562 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 563 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 564 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 565 : $ -1.0000000000000000D+00, 3.0000000000000000D+00, 566 : $ 3.0000000000000000D+00, -2.0000000000000000D+00, 567 : $ -4.0000000000000000D+00, -2.0000000000000000D+00, 568 : $ -2.9999999999999933D+00 569 : $ / 570 : DATA ((PNM21(I,J),I=1, 21),J= 8, 8)/ 571 : $ 0.0000000000000000D+00, -1.0000000000000000D+00, 572 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 573 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 574 : $ 2.9999999999999996D+00, 0.0000000000000000D+00, 575 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 576 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 577 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 578 : $ 0.0000000000000000D+00, 1.0000000000000000D+00, 579 : $ 0.0000000000000000D+00, -2.0000000000000000D+00, 580 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 581 : $ -2.9999999999999996D+00 582 : $ / 583 : DATA ((PNM21(I,J),I=1, 21),J= 9, 9)/ 584 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 585 : $ -1.0000000000000000D+00, 0.0000000000000000D+00, 586 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 587 : $ 2.9999999999999964D+00, 0.0000000000000000D+00, 588 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 589 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 590 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 591 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 592 : $ 1.0000000000000000D+00, 0.0000000000000000D+00, 593 : $ 0.0000000000000000D+00, -2.0000000000000000D+00, 594 : $ -2.9999999999999964D+00 595 : $ / 596 : DATA ((PNM21(I,J),I=1, 21),J= 10, 10)/ 597 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 598 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 599 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 600 : $ 1.5543122344752192D-15, 5.0000000000000000D-01, 601 : $ -1.5000000000000000D+00, -1.5000000000000000D+00, 602 : $ 1.0000000000000000D+00, 2.0000000000000000D+00, 603 : $ 1.0000000000000000D+00, 1.4999999999999982D+00, 604 : $ 5.0000000000000000D-01, -1.5000000000000000D+00, 605 : $ -1.5000000000000000D+00, 1.0000000000000000D+00, 606 : $ 2.0000000000000000D+00, 1.0000000000000000D+00, 607 : $ 1.4999999999999967D+00 608 : $ / 609 : DATA ((PNM21(I,J),I=1, 21),J= 11, 11)/ 610 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 611 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 612 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 613 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 614 : $ 2.0000000000000000D+00, 0.0000000000000000D+00, 615 : $ -2.0000000000000000D+00, -2.0000000000000000D+00, 616 : $ 0.0000000000000000D+00, -5.9999999999999991D+00, 617 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 618 : $ 0.0000000000000000D+00, -2.0000000000000000D+00, 619 : $ -2.0000000000000000D+00, 0.0000000000000000D+00, 620 : $ -5.9999999999999991D+00 621 : $ / 622 : DATA ((PNM21(I,J),I=1, 21),J= 12, 12)/ 623 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 624 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 625 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 626 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 627 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 628 : $ 1.0000000000000000D+00, 0.0000000000000000D+00, 629 : $ 0.0000000000000000D+00, 1.4999999999999998D+00, 630 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 631 : $ 0.0000000000000000D+00, 1.0000000000000000D+00, 632 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 633 : $ 1.4999999999999998D+00 634 : $ / 635 : DATA ((PNM21(I,J),I=1, 21),J= 13, 13)/ 636 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 637 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 638 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 639 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 640 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 641 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 642 : $ 0.0000000000000000D+00, -5.9999999999999991D+00, 643 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 644 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 645 : $ 2.0000000000000000D+00, 0.0000000000000000D+00, 646 : $ -5.9999999999999991D+00 647 : $ / 648 : DATA ((PNM21(I,J),I=1, 21),J= 14, 14)/ 649 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 650 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 651 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 652 : $ 1.5543122344752192D-15, 0.0000000000000000D+00, 653 : $ 0.0000000000000000D+00, -5.0000000000000000D-01, 654 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 655 : $ 1.0000000000000000D+00, 1.4999999999999998D+00, 656 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 657 : $ -5.0000000000000000D-01, 0.0000000000000000D+00, 658 : $ 0.0000000000000000D+00, 1.0000000000000000D+00, 659 : $ 1.4999999999999982D+00 660 : $ / 661 : DATA ((PNM21(I,J),I=1, 21),J= 15, 15)/ 662 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 663 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 664 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 665 : $ -3.5527136788005009D-15, 0.0000000000000000D+00, 666 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 667 : $ 0.0000000000000000D+00, -2.0000000000000000D+00, 668 : $ -2.0000000000000000D+00, -5.9999999999999991D+00, 669 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 670 : $ 2.0000000000000000D+00, 0.0000000000000000D+00, 671 : $ -2.0000000000000000D+00, -2.0000000000000000D+00, 672 : $ -5.9999999999999956D+00 673 : $ / 674 : DATA ((PNM21(I,J),I=1, 21),J= 16, 16)/ 675 : $ 0.0000000000000000D+00, 4.0000000000000000D+00, 676 : $ 0.0000000000000000D+00, -4.0000000000000000D+00, 677 : $ -4.0000000000000000D+00, 0.0000000000000000D+00, 678 : $ -1.1999999999999998D+01, 0.0000000000000000D+00, 679 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 680 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 681 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 682 : $ 0.0000000000000000D+00, -4.0000000000000000D+00, 683 : $ 0.0000000000000000D+00, 4.0000000000000000D+00, 684 : $ 4.0000000000000000D+00, 0.0000000000000000D+00, 685 : $ 1.1999999999999998D+01 686 : $ / 687 : DATA ((PNM21(I,J),I=1, 21),J= 17, 17)/ 688 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 689 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 690 : $ 4.0000000000000000D+00, 0.0000000000000000D+00, 691 : $ -1.1999999999999998D+01, 0.0000000000000000D+00, 692 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 693 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 694 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 695 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 696 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 697 : $ -4.0000000000000000D+00, 0.0000000000000000D+00, 698 : $ 1.1999999999999998D+01 699 : $ / 700 : DATA ((PNM21(I,J),I=1, 21),J= 18, 18)/ 701 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 702 : $ 4.0000000000000000D+00, 0.0000000000000000D+00, 703 : $ -4.0000000000000000D+00, -4.0000000000000000D+00, 704 : $ -1.1999999999999991D+01, 0.0000000000000000D+00, 705 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 706 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 707 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 708 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 709 : $ -4.0000000000000000D+00, 0.0000000000000000D+00, 710 : $ 4.0000000000000000D+00, 4.0000000000000000D+00, 711 : $ 1.1999999999999991D+01 712 : $ / 713 : DATA ((PNM21(I,J),I=1, 21),J= 19, 19)/ 714 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 715 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 716 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 717 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 718 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 719 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 720 : $ 0.0000000000000000D+00, -1.3499999999999998D+01, 721 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 722 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 723 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 724 : $ 1.3499999999999998D+01 725 : $ / 726 : DATA ((PNM21(I,J),I=1, 21),J= 20, 20)/ 727 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 728 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 729 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 730 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 731 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 732 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 733 : $ 0.0000000000000000D+00, 1.3499999999999998D+01, 734 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 735 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 736 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 737 : $ 1.3499999999999998D+01 738 : $ / 739 : DATA ((PNM21(I,J),I=1, 21),J= 21, 21)/ 740 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 741 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 742 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 743 : $ 2.6999999999999996D+01, 0.0000000000000000D+00, 744 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 745 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 746 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 747 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 748 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 749 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 750 : $ -2.6999999999999996D+01 751 : $ / 752 : DATA ((PNM15(I,J),I=1, 15),J= 1, 1)/ 753 : $ 1.0000000000000000D+00, -3.0000000000000000D+00, 754 : $ -3.0000000000000000D+00, -3.0000000000000000D+00, 755 : $ 2.0000000000000000D+00, 4.0000000000000000D+00, 756 : $ 4.0000000000000000D+00, 2.0000000000000000D+00, 757 : $ 4.0000000000000000D+00, 2.0000000000000000D+00, 758 : $ -3.9999999999999334D+00, -5.6477649705392858D-15, 759 : $ 2.9999999999999964D+00, 2.9999999999999964D+00, 760 : $ 2.9999999999999964D+00 761 : $ / 762 : DATA ((PNM15(I,J),I=1, 15),J= 2, 2)/ 763 : $ 0.0000000000000000D+00, 4.0000000000000000D+00, 764 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 765 : $ -4.0000000000000000D+00, -4.0000000000000000D+00, 766 : $ -4.0000000000000000D+00, 0.0000000000000000D+00, 767 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 768 : $ 3.1999999999999964D+01, 6.8821901302862365D-15, 769 : $ 0.0000000000000000D+00, -1.1999999999999998D+01, 770 : $ -1.1999999999999998D+01 771 : $ / 772 : DATA ((PNM15(I,J),I=1, 15),J= 3, 3)/ 773 : $ 0.0000000000000000D+00, -1.0000000000000000D+00, 774 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 775 : $ 2.0000000000000000D+00, 0.0000000000000000D+00, 776 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 777 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 778 : $ -3.9999999999999920D+00, 2.9999999999999982D+00, 779 : $ 0.0000000000000000D+00, 2.9999999999999996D+00, 780 : $ 2.9999999999999996D+00 781 : $ / 782 : DATA ((PNM15(I,J),I=1, 15),J= 4, 4)/ 783 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 784 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 785 : $ 0.0000000000000000D+00, 4.0000000000000000D+00, 786 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 787 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 788 : $ 3.1999999999999989D+01, -1.1999999999999998D+01, 789 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 790 : $ -1.1999999999999998D+01 791 : $ / 792 : DATA ((PNM15(I,J),I=1, 15),J= 5, 5)/ 793 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 794 : $ -1.0000000000000000D+00, 0.0000000000000000D+00, 795 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 796 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 797 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 798 : $ -3.9999999999999920D+00, 2.9999999999999982D+00, 799 : $ 2.9999999999999996D+00, 0.0000000000000000D+00, 800 : $ 2.9999999999999996D+00 801 : $ / 802 : DATA ((PNM15(I,J),I=1, 15),J= 6, 6)/ 803 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 804 : $ 4.0000000000000000D+00, 0.0000000000000000D+00, 805 : $ 0.0000000000000000D+00, -4.0000000000000000D+00, 806 : $ 0.0000000000000000D+00, -4.0000000000000000D+00, 807 : $ -4.0000000000000000D+00, 0.0000000000000000D+00, 808 : $ 3.1999999999999964D+01, 6.8821901302862365D-15, 809 : $ -1.1999999999999998D+01, 0.0000000000000000D+00, 810 : $ -1.1999999999999998D+01 811 : $ / 812 : DATA ((PNM15(I,J),I=1, 15),J= 7, 7)/ 813 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 814 : $ 0.0000000000000000D+00, 4.0000000000000000D+00, 815 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 816 : $ -4.0000000000000000D+00, 0.0000000000000000D+00, 817 : $ -4.0000000000000000D+00, -4.0000000000000000D+00, 818 : $ 3.1999999999999964D+01, 6.8821901302862365D-15, 819 : $ -1.1999999999999998D+01, -1.1999999999999998D+01, 820 : $ 0.0000000000000000D+00 821 : $ / 822 : DATA ((PNM15(I,J),I=1, 15),J= 8, 8)/ 823 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 824 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 825 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 826 : $ 4.0000000000000000D+00, 0.0000000000000000D+00, 827 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 828 : $ 3.1999999999999989D+01, -1.1999999999999998D+01, 829 : $ 0.0000000000000000D+00, -1.1999999999999998D+01, 830 : $ 0.0000000000000000D+00 831 : $ / 832 : DATA ((PNM15(I,J),I=1, 15),J= 9, 9)/ 833 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 834 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 835 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 836 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 837 : $ 4.0000000000000000D+00, 0.0000000000000000D+00, 838 : $ 3.1999999999999989D+01, -1.1999999999999998D+01, 839 : $ -1.1999999999999998D+01, 0.0000000000000000D+00, 840 : $ 0.0000000000000000D+00 841 : $ / 842 : DATA ((PNM15(I,J),I=1, 15),J= 10, 10)/ 843 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 844 : $ 0.0000000000000000D+00, -1.0000000000000000D+00, 845 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 846 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 847 : $ 0.0000000000000000D+00, 2.0000000000000000D+00, 848 : $ -3.9999999999999920D+00, 2.9999999999999982D+00, 849 : $ 2.9999999999999996D+00, 2.9999999999999996D+00, 850 : $ 0.0000000000000000D+00 851 : $ / 852 : DATA ((PNM15(I,J),I=1, 15),J= 11, 11)/ 853 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 854 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 855 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 856 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 857 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 858 : $ -1.0799999999999996D+02, -7.1183903498257817D-15, 859 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 860 : $ 2.6999999999999996D+01 861 : $ / 862 : DATA ((PNM15(I,J),I=1, 15),J= 12, 12)/ 863 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 864 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 865 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 866 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 867 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 868 : $ -1.0799999999999996D+02, -7.1183903498257817D-15, 869 : $ 0.0000000000000000D+00, 2.6999999999999996D+01, 870 : $ 0.0000000000000000D+00 871 : $ / 872 : DATA ((PNM15(I,J),I=1, 15),J= 13, 13)/ 873 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 874 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 875 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 876 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 877 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 878 : $ -1.0800000000000003D+02, 2.7000000000000007D+01, 879 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 880 : $ 0.0000000000000000D+00 881 : $ / 882 : DATA ((PNM15(I,J),I=1, 15),J= 14, 14)/ 883 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 884 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 885 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 886 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 887 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 888 : $ -1.0799999999999996D+02, -7.1183903498257817D-15, 889 : $ 2.6999999999999996D+01, 0.0000000000000000D+00, 890 : $ 0.0000000000000000D+00 891 : $ / 892 : DATA ((PNM15(I,J),I=1, 15),J= 15, 15)/ 893 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 894 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 895 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 896 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 897 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 898 : $ 2.5600000000000000D+02, -2.2189281662088338D-15, 899 : $ 0.0000000000000000D+00, 0.0000000000000000D+00, 900 : $ 0.0000000000000000D+00 901 : $ / 902 : * 903 : * 904 : C======================================================================= 905 : C QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS 906 : C 907 : C SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE FORME DU ND INO 908 : C SHP(2,INO) DERIVEES PR RPPRT A QSI 909 : C SHP(3,INO) DERIVEES PR RPPRT A ETA 910 : C SHP(4,INO) DERIVEES PR RPPRT A DZE 911 : C SHP(5,INO) DERIVEES PR RPPRT A 4 912 : C DERIVEES BARYCENTRIQUE 913 : C CAS DE DIMENSION 3 914 : C 15 NOEUDS TETRAEDRE 915 : C 19 NOEUDS PYRAMIDE 916 : C 21 NOEUDS PRISME 917 : C 27 NOEUDS CUBE 918 : C======================================================================= 919 : C 920 : C CUB27 921 : C 922 : C======================================================================= 923 : IF (IELE.EQ.33) THEN 924 : NBNOD=27 925 : * Init exposant des polynomes de la base polynomiale 926 : INOD=0 927 : DO IDZE=0,2 928 : DO IETA=0,2 929 : DO IKSI=0,2 930 : INOD=INOD+1 931 : NMON27(INOD)=1 932 : CMON27(1,INOD)=1.D0 933 : XMON27(1,1,INOD)=IKSI 934 : XMON27(2,1,INOD)=IETA 935 : XMON27(3,1,INOD)=IDZE 936 : ENDDO 937 : ENDDO 938 : ENDDO 939 : ORDDER(1)=0 940 : ORDDER(2)=0 941 : ORDDER(3)=0 942 : * Calcul des fonctions de forme et des dérivées premières 943 : * au point QSI,ETA,DZE 944 : * WRITE(6,*) 'QSI=',QSI,' ETA=',ETA,' DZE=',DZE 944 : CALL NI2(NDIM,NBNOD,QSI,ETA,DZE, 945 : $ NMMO27,NMON27,CMON27,XMON27,PNM27,XCO1,ORDDER, 946 : $ VECTP, 947 : $ SHP) 948 : C======================================================================= 949 : C 950 : C PR21 951 : C 952 : C======================================================================= 953 : ELSEIF (IELE.EQ.34) THEN 954 : NBNOD=21 955 : * Init exposant des polynomes de la base polynomiale 956 : INOD=0 957 : DO IDZE=0,2 958 : * D'abord la base polynomiale complète 959 : DO IDEG=0,2 960 : DO IEXP=IDEG,0,-1 961 : INOD=INOD+1 962 : NMON21(INOD)=1 963 : CMON21(1,INOD)=1.D0 964 : XMON21(1,1,INOD)=IEXP 965 : XMON21(2,1,INOD)=IDEG-IEXP 966 : XMON21(3,1,INOD)=IDZE 967 : ENDDO 968 : ENDDO 969 : * Puis la bulle 970 : * On génère un polynôme du type : 971 : * x_1.x_2 (1-x_1-x_2) 972 : INOD=INOD+1 973 : NMON21(INOD)=3 974 : * On initialise les monômes à x_1.x_2 975 : DO IMON=1,NMON21(INOD) 976 : XMON21(1,IMON,INOD)=1 977 : XMON21(2,IMON,INOD)=1 978 : XMON21(3,IMON,INOD)=IDZE 979 : ENDDO 980 : * On ajuste les coeff. et les monômes 981 : CMON21(1,INOD)=1.D0 982 : DO IMON=2,NMON21(INOD) 983 : CMON21(IMON,INOD)=-1.D0 984 : XMON21(IMON-1,IMON,INOD)=XMON21(IMON-1,IMON,INOD)+1 985 : ENDDO 986 : ENDDO 987 : ORDDER(1)=0 988 : ORDDER(2)=0 989 : ORDDER(3)=0 990 : * Calcul des fonctions de forme et des dérivées premières 991 : * au point QSI,ETA,DZE 992 : * WRITE(6,*) 'QSI=',QSI,' ETA=',ETA,' DZE=',DZE 992 : CALL NI2(NDIM,NBNOD,QSI,ETA,DZE, 993 : $ NMMO21,NMON21,CMON21,XMON21,PNM21,XCO1,ORDDER, 994 : $ VECTP, 995 : $ SHP) 996 : C======================================================================= 997 : C 998 : C TE15 999 : C 1000 : C======================================================================= 1001 : ELSEIF (IELE.EQ.35) THEN 1002 : NBNOD=15 1003 : * Init exposant des polynomes de la base polynomiale 1004 : INOD=0 1005 : * D'abord la base polynomiale complète 1006 : DO IDEG=0,2 1007 : DO IEXP=IDEG,0,-1 1008 : DO IEXP2=IDEG-IEXP,0,-1 1009 : INOD=INOD+1 1010 : NMON15(INOD)=1 1011 : CMON15(1,INOD)=1.D0 1012 : XMON15(1,1,INOD)=IEXP 1013 : XMON15(2,1,INOD)=IEXP2 1014 : XMON15(3,1,INOD)=IDEG-(IEXP+IEXP2) 1015 : ENDDO 1016 : ENDDO 1017 : ENDDO 1018 : * Puis la bulle du centre 1019 : * On génère un polynôme du type : 1020 : * x_1.x_2.x_3 (1-x_1-x_2-x_3) 1021 : INOD=INOD+1 1022 : NMON15(INOD)=4 1023 : * On initialise les monômes à x_1.x_2.x_3 1024 : DO IMON=1,NMON15(INOD) 1025 : XMON15(1,IMON,INOD)=1 1026 : XMON15(2,IMON,INOD)=1 1027 : XMON15(3,IMON,INOD)=1 1028 : ENDDO 1029 : * On ajuste les coeff. et les monômes 1030 : CMON15(1,INOD)=1.D0 1031 : DO IMON=2,NMON15(INOD) 1032 : CMON15(IMON,INOD)=-1.D0 1033 : XMON15(IMON-1,IMON,INOD)=XMON15(IMON-1,IMON,INOD)+1 1034 : ENDDO 1035 : * Puis les bulles des faces 1036 : * On génère le polynôme : x_1.x_2...x_n 1037 : INOD=INOD+1 1038 : NMON15(INOD)=1 1039 : CMON15(1,INOD)=1.D0 1040 : XMON15(1,1,INOD)=1 1041 : XMON15(2,1,INOD)=1 1042 : XMON15(3,1,INOD)=1 1043 : * On génère les polynômes du type : 1044 : * ((x_1.x_2,x_3)/x_i) (1-x_1...-x_3) 1045 : DO JDIM=1,NDIM 1046 : INOD=INOD+1 1047 : NMON15(INOD)=4 1048 : * On initialise les monômes à x_1.x_2.x_3 1049 : DO IMON=1,NMON15(INOD) 1050 : XMON15(1,IMON,INOD)=1 1051 : XMON15(2,IMON,INOD)=1 1052 : XMON15(3,IMON,INOD)=1 1053 : ENDDO 1054 : * On divise par x_i 1055 : DO IMON=1,NMON15(INOD) 1056 : XMON15(JDIM,IMON,INOD)=XMON15(JDIM,IMON,INOD)-1 1057 : ENDDO 1058 : * On ajuste les coeff. et les monômes 1059 : CMON15(1,INOD)=UN 1060 : DO IMON=2,NMON15(INOD) 1061 : CMON15(IMON,INOD)=-UN 1062 : XMON15(IMON-1,IMON,INOD)=XMON15(IMON-1,IMON,INOD)+1 1063 : ENDDO 1064 : ENDDO 1065 : ORDDER(1)=0 1066 : ORDDER(2)=0 1067 : ORDDER(3)=0 1068 : * Calcul des fonctions de forme et des dérivées premières 1069 : * au point QSI,ETA,DZE 1070 : * WRITE(6,*) 'QSI=',QSI,' ETA=',ETA,' DZE=',DZE 1070 : CALL NI2(NDIM,NBNOD,QSI,ETA,DZE, 1071 : $ NMMO15,NMON15,CMON15,XMON15,PNM15,XCO1,ORDDER, 1072 : $ VECTP, 1073 : $ SHP) 1074 : ELSE 1075 : * WRITE(6,*) 'IELE=',IELE,' non implem' 1075 : * CALL ERREUR(5) 1076 : IRET=0 1077 : RETURN 1078 : ENDIF 1079 : IRET=1 1080 : RETURN 1081 : END 1082 : 1083 : 1084 : 1085 : 1086 :
© Cast3M 2003 - All rights reserved.
Disclaimer