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