1 : C PRE61 SOURCE BECC 11/01/05 21:15:19 6836 2 : SUBROUTINE PRE61() 3 : C************************************************************************ 4 : C 5 : C PROJET : CASTEM 2000 6 : C 7 : C NOM : PRE61 8 : C 9 : C DESCRIPTION : Voir PRE6 10 : C 11 : C Discrete Equations Method 12 : C 13 : C 2me ordre en espace 14 : C 15 : C Creation of the MCHAMLs IALPF, IROF, IVITF, IPF, 16 : C for the two phases. 17 : C 18 : C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI) 19 : C 20 : C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF 21 : C 22 : C************************************************************************ 23 : C 24 : C 25 : C************************************************************************ 26 : C 27 : C HISTORIQUE (Anomalies et modifications éventuelles) 28 : C 29 : C HISTORIQUE : Crée le 03.12.2009 30 : C Estension au 3D le 21.12.2010 31 : C 32 : C************************************************************************ 33 : C 34 : IMPLICIT INTEGER(I-N) 35 : C 36 : C**** Les variables 37 : C 38 : INTEGER ICOND, IRETOU, MMODEL, ICELL 39 : & , IDOMA, ICEN, IFACE, IFACEL, INORM 40 : & , IAL1, IGRAL1, ILIAL1 41 : & , IAL2, IGRAL2, ILIAL2 42 : & , IRN1, IGRRN1, ILIRN1 43 : & , IRN2, IGRRN2, ILIRN2 44 : & , IVN1, IGRVN1, ILIVN1 45 : & , IVN2, IGRVN2, ILIVN2 46 : & , IPN1, IGRPN1, ILIPN1 47 : & , IPN2, IGRPN2, ILIPN2 48 : & , IAL1F, IAL2F, IRN1F, IRN2F, IVN1F, IVN2F, IPN1F, IPN2F 49 : C 50 : CHARACTER*(4) NOMGRA(27),NOMLIM(9) 51 : CHARACTER*(8) MTYPR 52 : CHARACTER*(40) MESERR 53 : C 54 : C**** Les Includes 55 : C 56 : -INC CCOPTIO 57 : INTEGER JGM, JGN 58 : -INC SMLMOTS 59 : POINTEUR MLMCOM.MLMOTS, MLMVIT.MLMOTS, MLMTEN.MLMOTS 60 : C 61 : C**** Nom de composantes de gradients (HP. <= 9 composantes) 62 : C 63 : DATA NOMGRA /'P1DX','P1DY','P1DZ', 64 : & 'P2DX','P2DY','P2DZ', 65 : & 'P3DX','P3DY','P3DZ', 66 : & 'P4DX','P4DY','P4DZ', 67 : & 'P5DX','P5DY','P5DZ', 68 : & 'P6DX','P6DY','P6DZ', 69 : & 'P7DX','P7DY','P7DZ', 70 : & 'P8DX','P8DY','P8DZ', 71 : & 'P9DX','P9DY','P9DZ'/ 72 : C 73 : DATA NOMLIM /'P1 ', 74 : & 'P2 ', 75 : & 'P3 ', 76 : & 'P4 ', 77 : & 'P5 ', 78 : & 'P6 ', 79 : & 'P7 ', 80 : & 'P8 ', 81 : & 'P9 '/ 82 : C 83 : C**** Initialisation of some segment 84 : C 85 : JGN=4 86 : JGM=1 87 : SEGINI MLMCOM 88 : JGN=4 89 : JGM=IDIM 90 : SEGINI MLMVIT 91 : JGN=4 92 : JGM=IDIM*IDIM 93 : SEGINI MLMTEN 94 : C 95 : C**** Lecture de l'objet MODELE 96 : C 97 : ICOND = 1 98 : MTYPR = 'MMODEL ' 98 : CALL LIROBJ(MTYPR, MMODEL, ICOND, IRETOU) 99 : IF(IERR.NE.0)GOTO 9999 99 : CALL LEKMOD(MMODEL, IDOMA, ICELL) 100 : IF(IERR.NE.0)GOTO 9999 101 : C 102 : C**** Lecture du MELEME SPG des points CENTRE. 103 : C 104 : C 104 : C CALL LEKTAB(IDOMA,'CENTRE',IP) 105 : C 106 : C**** Probleme du LEKTAB: si IDOMA.'CENTRE' n'existe pas, 107 : C il crèe IDOMA.'CENTRE' sans recrèer 'FACEL' 108 : C -> la correspondance global des noeuds saut! 109 : C 110 : C On peut utilizer ACCTAB ou ACMO 111 : C 112 : MTYPR = 'MAILLAGE' 112 : CALL ACMO(IDOMA,'CENTRE',MTYPR,ICEN) 113 : IF(IERR.NE.0)GOTO 9999 114 : C 115 : C**** Lecture du MELEME 'FACE' 116 : C 117 : MTYPR = 'MAILLAGE' 117 : CALL ACMO(IDOMA,'FACE',MTYPR,IFACE) 118 : IF(IERR.NE.0)GOTO 9999 119 : C 120 : C**** Lecture du MELEME 'FACEL' 121 : C 122 : MTYPR = 'MAILLAGE' 122 : CALL ACMO(IDOMA,'FACEL',MTYPR,IFACEL) 123 : IF(IERR.NE.0)GOTO 9999 124 : C 125 : C**** Lecture du CHPOINT contenant les normales (tangentes) aux faces 126 : C 127 : IF(IDIM .EQ. 2)THEN 128 : C Que les normales 128 : CALL LEKTAB(IDOMA,'XXNORMAF',INORM) 129 : IF(IERR .NE. 0) GOTO 9999 130 : MLMVIT.MOTS(1) = 'UX ' 131 : MLMVIT.MOTS(2) = 'UY ' 131 : CALL QUEPO1(INORM, IFACE, MLMVIT) 132 : IF(IERR.NE.0)GOTO 9999 133 : ELSE 134 : C Les normales et les tangentes 135 : MTYPR = ' ' 135 : CALL ACMO(IDOMA,'MATROT',MTYPR,INORM) 136 : IF (MTYPR .NE. 'CHPOINT ') THEN 136 : CALL MATRAN(IDOMA,INORM) 137 : IF(IERR .NE. 0) GOTO 9999 138 : ENDIF 139 : MLMTEN.MOTS(1) = 'UX ' 140 : MLMTEN.MOTS(2) = 'UY ' 141 : MLMTEN.MOTS(3) = 'UZ ' 142 : MLMTEN.MOTS(4) = 'RX ' 143 : MLMTEN.MOTS(5) = 'RY ' 144 : MLMTEN.MOTS(6) = 'RZ ' 145 : MLMTEN.MOTS(7) = 'MX ' 146 : MLMTEN.MOTS(8) = 'MY ' 147 : MLMTEN.MOTS(9) = 'MZ ' 147 : CALL QUEPO1(INORM, IFACE, MLMTEN) 148 : IF(IERR.NE.0)GOTO 9999 149 : ENDIF 150 : C 151 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 152 : C**** Lecture des CHPOINTs alpha ****C 153 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 154 : C 155 : C**** IAL1 156 : C 157 : ICOND = 1 158 : MTYPR = 'CHPOINT ' 158 : CALL LIROBJ(MTYPR,IAL1,ICOND,IRETOU) 159 : IF(IERR .NE. 0)GOTO 9999 160 : C Control du CHPOINT: QUEPO1 161 : SEGACT MLMCOM*MOD 162 : MLMCOM.MOTS(1)='SCAL' 162 : CALL QUEPO1(IAL1, ICEN, MLMCOM) 163 : SEGDES MLMCOM 164 : IF(IERR .NE. 0)THEN 165 : GOTO 9999 166 : ENDIF 167 : C 168 : C**** Lecture du CHPOINT IGRAL1 169 : C 170 : ICOND = 1 171 : MTYPR = 'CHPOINT ' 171 : CALL LIROBJ(MTYPR,IGRAL1,ICOND,IRETOU) 172 : IF (IERR.NE.0) GOTO 9999 173 : C Control du CHPOINT: QUEPO1 174 : SEGACT MLMVIT*MOD 175 : MLMVIT.MOTS(1)=NOMGRA(1) 176 : MLMVIT.MOTS(2)=NOMGRA(2) 177 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMGRA(3) 177 : CALL QUEPO1(IGRAL1, ICEN, MLMVIT) 178 : SEGDES MLMVIT 179 : IF(IERR .NE. 0)THEN 180 : GOTO 9999 181 : ENDIF 182 : C 183 : C**** Lecture du CHPOINT ILIAL1 184 : C 185 : ICOND = 1 186 : MTYPR = 'CHPOINT ' 186 : CALL LIROBJ(MTYPR,ILIAL1,ICOND,IRETOU) 187 : IF (IERR.NE.0) GOTO 9999 188 : C Control du CHPOINT: QUEPO1 189 : SEGACT MLMCOM*MOD 190 : MLMCOM.MOTS(1)= NOMLIM(1) 190 : CALL QUEPO1(ILIAL1, ICEN, MLMCOM) 191 : SEGDES MLMCOM 192 : IF(IERR .NE. 0)THEN 193 : GOTO 9999 194 : ENDIF 195 : C 196 : C**** IAL2 197 : C 198 : ICOND = 1 199 : MTYPR = 'CHPOINT ' 199 : CALL LIROBJ(MTYPR,IAL2,ICOND,IRETOU) 200 : IF(IERR .NE. 0)GOTO 9999 201 : C Control du CHPOINT: QUEPO1 202 : SEGACT MLMCOM*MOD 203 : MLMCOM.MOTS(1)='SCAL' 203 : CALL QUEPO1(IAL2, ICEN, MLMCOM) 204 : SEGDES MLMCOM 205 : IF(IERR .NE. 0)THEN 206 : GOTO 9999 207 : ENDIF 208 : C 209 : C**** Lecture du CHPOINT IGRAL2 210 : C 211 : ICOND = 1 212 : MTYPR = 'CHPOINT ' 212 : CALL LIROBJ(MTYPR,IGRAL2,ICOND,IRETOU) 213 : IF (IERR.NE.0) GOTO 9999 214 : C Control du CHPOINT: QUEPO1 215 : SEGACT MLMVIT*MOD 216 : MLMVIT.MOTS(1)=NOMGRA(1) 217 : MLMVIT.MOTS(2)=NOMGRA(2) 218 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMGRA(3) 218 : CALL QUEPO1(IGRAL2, ICEN, MLMVIT) 219 : SEGDES MLMVIT 220 : IF(IERR .NE. 0)THEN 221 : GOTO 9999 222 : ENDIF 223 : C 224 : C**** Lecture du CHPOINT ILIAL2 225 : C 226 : ICOND = 1 227 : MTYPR = 'CHPOINT ' 227 : CALL LIROBJ(MTYPR,ILIAL2,ICOND,IRETOU) 228 : IF (IERR.NE.0) GOTO 9999 229 : C Control du CHPOINT: QUEPO1 230 : SEGACT MLMCOM*MOD 231 : MLMCOM.MOTS(1)= NOMLIM(1) 231 : CALL QUEPO1(ILIAL2, ICEN, MLMCOM) 232 : SEGDES MLMCOM 233 : IF(IERR .NE. 0)THEN 234 : GOTO 9999 235 : ENDIF 236 : C 237 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 238 : C**** Lecture des CHPOINTs rho ****C 239 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 240 : C 241 : C**** RN1 242 : C 243 : ICOND = 1 244 : MTYPR = 'CHPOINT ' 244 : CALL LIROBJ(MTYPR,IRN1,ICOND,IRETOU) 245 : IF(IERR .NE. 0)GOTO 9999 246 : C Control du CHPOINT: QUEPO1 247 : SEGACT MLMCOM*MOD 248 : MLMCOM.MOTS(1)='SCAL' 248 : CALL QUEPO1(IRN1, ICEN, MLMCOM) 249 : SEGDES MLMCOM 250 : IF(IERR .NE. 0)THEN 251 : GOTO 9999 252 : ENDIF 253 : C 254 : C**** Lecture du CHPOINT IGRRN1 255 : C 256 : ICOND = 1 257 : MTYPR = 'CHPOINT ' 257 : CALL LIROBJ(MTYPR,IGRRN1,ICOND,IRETOU) 258 : IF (IERR.NE.0) GOTO 9999 259 : C Control du CHPOINT: QUEPO1 260 : SEGACT MLMVIT*MOD 261 : MLMVIT.MOTS(1)=NOMGRA(1) 262 : MLMVIT.MOTS(2)=NOMGRA(2) 263 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMGRA(3) 263 : CALL QUEPO1(IGRRN1, ICEN, MLMVIT) 264 : SEGDES MLMVIT 265 : IF(IERR .NE. 0)THEN 266 : GOTO 9999 267 : ENDIF 268 : C 269 : C**** Lecture du CHPOINT ILIRN1 270 : C 271 : ICOND = 1 272 : MTYPR = 'CHPOINT ' 272 : CALL LIROBJ(MTYPR,ILIRN1,ICOND,IRETOU) 273 : IF (IERR.NE.0) GOTO 9999 274 : C Control du CHPOINT: QUEPO1 275 : SEGACT MLMCOM*MOD 276 : MLMCOM.MOTS(1)= NOMLIM(1) 276 : CALL QUEPO1(ILIRN1, ICEN, MLMCOM) 277 : SEGDES MLMCOM 278 : IF(IERR .NE. 0)THEN 279 : GOTO 9999 280 : ENDIF 281 : C 282 : C**** IRN2 283 : C 284 : ICOND = 1 285 : MTYPR = 'CHPOINT ' 285 : CALL LIROBJ(MTYPR,IRN2,ICOND,IRETOU) 286 : IF(IERR .NE. 0)GOTO 9999 287 : C Control du CHPOINT: QUEPO1 288 : SEGACT MLMCOM*MOD 289 : MLMCOM.MOTS(1)='SCAL' 289 : CALL QUEPO1(IRN2, ICEN, MLMCOM) 290 : SEGDES MLMCOM 291 : IF(IERR .NE. 0)THEN 292 : GOTO 9999 293 : ENDIF 294 : C 295 : C**** Lecture du CHPOINT IGRRN2 296 : C 297 : ICOND = 1 298 : MTYPR = 'CHPOINT ' 298 : CALL LIROBJ(MTYPR,IGRRN2,ICOND,IRETOU) 299 : IF (IERR.NE.0) GOTO 9999 300 : C Control du CHPOINT: QUEPO1 301 : SEGACT MLMVIT*MOD 302 : MLMVIT.MOTS(1)=NOMGRA(1) 303 : MLMVIT.MOTS(2)=NOMGRA(2) 304 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMGRA(3) 304 : CALL QUEPO1(IGRRN2, ICEN, MLMVIT) 305 : SEGDES MLMVIT 306 : IF(IERR .NE. 0)THEN 307 : GOTO 9999 308 : ENDIF 309 : C 310 : C**** Lecture du CHPOINT ILIRN2 311 : C 312 : ICOND = 1 313 : MTYPR = 'CHPOINT ' 313 : CALL LIROBJ(MTYPR,ILIRN2,ICOND,IRETOU) 314 : IF (IERR.NE.0) GOTO 9999 315 : C Control du CHPOINT: QUEPO1 316 : SEGACT MLMCOM*MOD 317 : MLMCOM.MOTS(1)= NOMLIM(1) 317 : CALL QUEPO1(ILIRN2, ICEN, MLMCOM) 318 : SEGDES MLMCOM 319 : IF(IERR .NE. 0)THEN 320 : GOTO 9999 321 : ENDIF 322 : C 323 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 324 : C**** Lecture des CHPOINTs vitesse ****C 325 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 326 : C 327 : C**** VN1 328 : C 329 : ICOND = 1 330 : MTYPR = 'CHPOINT ' 330 : CALL LIROBJ(MTYPR,IVN1,ICOND,IRETOU) 331 : IF(IERR .NE. 0)GOTO 9999 332 : C Control du CHPOINT: QUEPO1 333 : SEGACT MLMVIT*MOD 334 : MLMVIT.MOTS(1)='UX ' 335 : MLMVIT.MOTS(2)='UY ' 336 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = 'UZ ' 336 : CALL QUEPO1(IVN1, ICEN, MLMVIT) 337 : SEGDES MLMVIT 338 : IF(IERR .NE. 0)THEN 339 : GOTO 9999 340 : ENDIF 341 : C 342 : C**** Lecture du CHPOINT IGRVN1 343 : C 344 : ICOND = 1 345 : MTYPR = 'CHPOINT ' 345 : CALL LIROBJ(MTYPR,IGRVN1,ICOND,IRETOU) 346 : IF (IERR.NE.0) GOTO 9999 347 : C Control du CHPOINT: QUEPO1 348 : SEGACT MLMTEN*MOD 349 : IF (IDIM .EQ. 2)THEN 350 : MLMTEN.MOTS(1)=NOMGRA(1) 351 : MLMTEN.MOTS(2)=NOMGRA(2) 352 : MLMTEN.MOTS(3)=NOMGRA(4) 353 : MLMTEN.MOTS(4)=NOMGRA(5) 354 : ELSEIF(IDIM .EQ. 3) THEN 355 : MLMTEN.MOTS(1)=NOMGRA(1) 356 : MLMTEN.MOTS(2)=NOMGRA(2) 357 : MLMTEN.MOTS(3)=NOMGRA(3) 358 : MLMTEN.MOTS(4)=NOMGRA(4) 359 : MLMTEN.MOTS(5)=NOMGRA(5) 360 : MLMTEN.MOTS(6)=NOMGRA(6) 361 : MLMTEN.MOTS(7)=NOMGRA(7) 362 : MLMTEN.MOTS(8)=NOMGRA(8) 363 : MLMTEN.MOTS(9)=NOMGRA(9) 364 : ENDIF 364 : CALL QUEPO1(IGRVN1, ICEN, MLMTEN) 365 : SEGDES MLMTEN 366 : IF(IERR .NE. 0)THEN 367 : GOTO 9999 368 : ENDIF 369 : C 370 : C 371 : C**** Lecture du CHPOINT ILIVN1 372 : C 373 : ICOND = 1 374 : MTYPR = 'CHPOINT ' 374 : CALL LIROBJ(MTYPR,ILIVN1,ICOND,IRETOU) 375 : IF (IERR.NE.0) GOTO 9999 376 : C Control du CHPOINT: QUEPO1 377 : SEGACT MLMVIT*MOD 378 : MLMVIT.MOTS(1)=NOMLIM(1) 379 : MLMVIT.MOTS(2)=NOMLIM(2) 380 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMLIM(3) 380 : CALL QUEPO1(ILIVN1, ICEN, MLMVIT) 381 : SEGDES MLMVIT 382 : IF(IERR .NE. 0)THEN 383 : GOTO 9999 384 : ENDIF 385 : C 386 : C**** VN2 387 : C 388 : ICOND = 1 389 : MTYPR = 'CHPOINT ' 389 : CALL LIROBJ(MTYPR,IVN2,ICOND,IRETOU) 390 : IF(IERR .NE. 0)GOTO 9999 391 : C Control du CHPOINT: QUEPO1 392 : SEGACT MLMVIT*MOD 393 : MLMVIT.MOTS(1)='UX ' 394 : MLMVIT.MOTS(2)='UY ' 395 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = 'UZ ' 395 : CALL QUEPO1(IVN2, ICEN, MLMVIT) 396 : SEGDES MLMVIT 397 : IF(IERR .NE. 0)THEN 398 : GOTO 9999 399 : ENDIF 400 : C 401 : C**** Lecture du CHPOINT IGRVN2 402 : C 403 : ICOND = 1 404 : MTYPR = 'CHPOINT ' 404 : CALL LIROBJ(MTYPR,IGRVN2,ICOND,IRETOU) 405 : IF (IERR.NE.0) GOTO 9999 406 : C Control du CHPOINT: QUEPO1 407 : SEGACT MLMTEN*MOD 408 : IF (IDIM .EQ. 2)THEN 409 : MLMTEN.MOTS(1)=NOMGRA(1) 410 : MLMTEN.MOTS(2)=NOMGRA(2) 411 : MLMTEN.MOTS(3)=NOMGRA(4) 412 : MLMTEN.MOTS(4)=NOMGRA(5) 413 : ELSEIF(IDIM .EQ. 3) THEN 414 : MLMTEN.MOTS(1)=NOMGRA(1) 415 : MLMTEN.MOTS(2)=NOMGRA(2) 416 : MLMTEN.MOTS(3)=NOMGRA(3) 417 : MLMTEN.MOTS(4)=NOMGRA(4) 418 : MLMTEN.MOTS(5)=NOMGRA(5) 419 : MLMTEN.MOTS(6)=NOMGRA(6) 420 : MLMTEN.MOTS(7)=NOMGRA(7) 421 : MLMTEN.MOTS(8)=NOMGRA(8) 422 : MLMTEN.MOTS(9)=NOMGRA(9) 423 : ENDIF 423 : CALL QUEPO1(IGRVN2, ICEN, MLMTEN) 424 : SEGDES MLMTEN 425 : IF(IERR .NE. 0)THEN 426 : GOTO 9999 427 : ENDIF 428 : C 429 : C 430 : C**** Lecture du CHPOINT ILIVN2 431 : C 432 : ICOND = 1 433 : MTYPR = 'CHPOINT ' 433 : CALL LIROBJ(MTYPR,ILIVN2,ICOND,IRETOU) 434 : IF (IERR.NE.0) GOTO 9999 435 : C Control du CHPOINT: QUEPO1 436 : SEGACT MLMVIT*MOD 437 : MLMVIT.MOTS(1)=NOMLIM(1) 438 : MLMVIT.MOTS(2)=NOMLIM(2) 439 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMLIM(3) 439 : CALL QUEPO1(ILIVN2, ICEN, MLMVIT) 440 : SEGDES MLMVIT 441 : IF(IERR .NE. 0)THEN 442 : GOTO 9999 443 : ENDIF 444 : C 445 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 446 : C**** Lecture des CHPOINTs rho ****C 447 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 448 : C 449 : C**** PN1 450 : C 451 : ICOND = 1 452 : MTYPR = 'CHPOINT ' 452 : CALL LIROBJ(MTYPR,IPN1,ICOND,IRETOU) 453 : IF(IERR .NE. 0)GOTO 9999 454 : C Control du CHPOINT: QUEPO1 455 : SEGACT MLMCOM*MOD 456 : MLMCOM.MOTS(1)='SCAL' 456 : CALL QUEPO1(IPN1, ICEN, MLMCOM) 457 : SEGDES MLMCOM 458 : IF(IERR .NE. 0)THEN 459 : GOTO 9999 460 : ENDIF 461 : C 462 : C**** Lecture du CHPOINT IGRPN1 463 : C 464 : ICOND = 1 465 : MTYPR = 'CHPOINT ' 465 : CALL LIROBJ(MTYPR,IGRPN1,ICOND,IRETOU) 466 : IF (IERR.NE.0) GOTO 9999 467 : C Control du CHPOINT: QUEPO1 468 : SEGACT MLMVIT*MOD 469 : MLMVIT.MOTS(1)=NOMGRA(1) 470 : MLMVIT.MOTS(2)=NOMGRA(2) 471 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMGRA(3) 471 : CALL QUEPO1(IGRPN1, ICEN, MLMVIT) 472 : SEGDES MLMVIT 473 : IF(IERR .NE. 0)THEN 474 : GOTO 9999 475 : ENDIF 476 : C 477 : C**** Lecture du CHPOINT ILIPN1 478 : C 479 : ICOND = 1 480 : MTYPR = 'CHPOINT ' 480 : CALL LIROBJ(MTYPR,ILIPN1,ICOND,IRETOU) 481 : IF (IERR.NE.0) GOTO 9999 482 : C Control du CHPOINT: QUEPO1 483 : SEGACT MLMCOM*MOD 484 : MLMCOM.MOTS(1)= NOMLIM(1) 484 : CALL QUEPO1(ILIPN1, ICEN, MLMCOM) 485 : SEGDES MLMCOM 486 : IF(IERR .NE. 0)THEN 487 : GOTO 9999 488 : ENDIF 489 : C 490 : C**** IPN2 491 : C 492 : ICOND = 1 493 : MTYPR = 'CHPOINT ' 493 : CALL LIROBJ(MTYPR,IPN2,ICOND,IRETOU) 494 : IF(IERR .NE. 0)GOTO 9999 495 : C Control du CHPOINT: QUEPO1 496 : SEGACT MLMCOM*MOD 497 : MLMCOM.MOTS(1)='SCAL' 497 : CALL QUEPO1(IPN2, ICEN, MLMCOM) 498 : SEGDES MLMCOM 499 : IF(IERR .NE. 0)THEN 500 : GOTO 9999 501 : ENDIF 502 : C 503 : C**** Lecture du CHPOINT IGRPN2 504 : C 505 : ICOND = 1 506 : MTYPR = 'CHPOINT ' 506 : CALL LIROBJ(MTYPR,IGRPN2,ICOND,IRETOU) 507 : IF (IERR.NE.0) GOTO 9999 508 : C Control du CHPOINT: QUEPO1 509 : SEGACT MLMVIT*MOD 510 : MLMVIT.MOTS(1)=NOMGRA(1) 511 : MLMVIT.MOTS(2)=NOMGRA(2) 512 : IF(IDIM .EQ. 3) MLMVIT.MOTS(3) = NOMGRA(3) 512 : CALL QUEPO1(IGRPN2, ICEN, MLMVIT) 513 : SEGDES MLMVIT 514 : IF(IERR .NE. 0)THEN 515 : GOTO 9999 516 : ENDIF 517 : C 518 : C**** Lecture du CHPOINT ILIPN2 519 : C 520 : ICOND = 1 521 : MTYPR = 'CHPOINT ' 521 : CALL LIROBJ(MTYPR,ILIPN2,ICOND,IRETOU) 522 : IF (IERR.NE.0) GOTO 9999 523 : C Control du CHPOINT: QUEPO1 524 : SEGACT MLMCOM*MOD 525 : MLMCOM.MOTS(1)= NOMLIM(1) 525 : CALL QUEPO1(ILIPN2, ICEN, MLMCOM) 526 : SEGDES MLMCOM 527 : IF(IERR .NE. 0)THEN 528 : GOTO 9999 529 : ENDIF 530 : C 531 : C write(*,*) 'Fin qui' 532 : C write(*,*) IAL1, IGRAL1, ILIAL1 533 : C write(*,*) IAL2, IGRAL2, ILIAL2 534 : C write(*,*) IRN1, IGRRN1, ILIRN1 535 : C write(*,*) IRN2, IGRRN2, ILIRN2 536 : C write(*,*) IVN1, IGRVN1, ILIVN1 537 : C write(*,*) IVN2, IGRVN2, ILIVN2 538 : C write(*,*) IPN1, IGRPN1, ILIPN1 539 : C write(*,*) IPN2, IGRPN2, ILIPN2 540 : C goto 9999 541 : C 541 : CALL PRE611( 542 : & ICEN,IFACE,IFACEL,INORM, 543 : & IAL1, IGRAL1, ILIAL1, 544 : & IAL2, IGRAL2, ILIAL2, 545 : & IRN1, IGRRN1, ILIRN1, 546 : & IRN2, IGRRN2, ILIRN2, 547 : & IVN1, IGRVN1, ILIVN1, 548 : & IVN2, IGRVN2, ILIVN2, 549 : & IPN1, IGRPN1, ILIPN1, 550 : & IPN2, IGRPN2, ILIPN2, 551 : & IAL1F, IAL2F, IRN1F, IRN2F, IVN1F, IVN2F, IPN1F, IPN2F) 552 : C 553 : SEGSUP MLMCOM 554 : SEGSUP MLMVIT 555 : SEGSUP MLMTEN 556 : C 557 : C**** Ecriture de 558 : C IAL1F, IAL2F, IRN1F, IRN2F, IVN1F, IVN2F, IPN1F, IPN2F 559 : C 560 : MTYPR = 'MCHAML ' 560 : CALL ECROBJ(MTYPR,IPN2F) 560 : CALL ECROBJ(MTYPR,IPN1F) 560 : CALL ECROBJ(MTYPR,IVN2F) 560 : CALL ECROBJ(MTYPR,IVN1F) 560 : CALL ECROBJ(MTYPR,IRN2F) 560 : CALL ECROBJ(MTYPR,IRN1F) 560 : CALL ECROBJ(MTYPR,IAL2F) 560 : CALL ECROBJ(MTYPR,IAL1F) 561 : C 562 : 9999 CONTINUE 563 : C 564 : RETURN 565 : END 566 : 567 :
© Cast3M 2003 - All rights reserved.
Disclaimer