Télécharger feqpr.eso

Retour à la liste

Numérotation des lignes :

  1. C FEQPR SOURCE MB234859 16/12/14 21:15:07 9253
  2. SUBROUTINE FEQPR(IPMODL,IPCHM1,IPCHM2,IPCHP4,IRET)
  3. C_______________________________________________________________________
  4. C
  5. C ENTREES:
  6. C ________
  7. C
  8. C IPMODL Pointeur sur un MMODEL
  9. C IPCHM1 Pointeur sur un MCHAML de CONTRAINTES
  10. C IPCHM2 Pointeur sur un MCHAML de CARACTERISTIQUES
  11. C
  12. C SORTIES:
  13. C ________
  14. C
  15. C IPCHP4 Pointeur sur un CHPOINT de forces aux noeuds
  16. C IRET = 1 OU 0 suivant succes ou pas (Message d'erreur
  17. C imprime dans ce cas)
  18. C_______________________________________________________________________
  19. C
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22. C
  23. -INC CCOPTIO
  24. -INC CCHAMP
  25. -INC SMMODEL
  26. -INC SMCHAML
  27. -INC SMCHPOI
  28. -INC SMELEME
  29. -INC SMINTE
  30. -INC SMLENTI
  31. C
  32. SEGMENT NOTYPE
  33. CHARACTER*16 TYPE(NBTYPE)
  34. ENDSEGMENT
  35. C
  36. SEGMENT MPTVAL
  37. INTEGER IPOS(NS) ,NSOF(NS)
  38. INTEGER IVAL(NCOSOU)
  39. CHARACTER*16 TYVAL(NCOSOU)
  40. ENDSEGMENT
  41. C
  42. SEGMENT LIMODL(0)
  43. C
  44. PARAMETER ( NINF=3 )
  45. INTEGER INFOS(NINF)
  46. CHARACTER*8 CMATE
  47. CHARACTER*(NCONCH) CONM
  48. C_______________________________________________________________________
  49. C
  50. C A T-ON BIEN UN MMODEL CHARGEMENT PRESSION
  51. C_______________________________________________________________________
  52. C
  53. MMODEL=IPMODL
  54. SEGACT,MMODEL
  55. NSOUS = MMODEL.KMODEL(/1)
  56. SEGINI, LIMODL
  57. DO ISOUS = 1,NSOUS
  58. IMODEL = MMODEL.KMODEL(ISOUS)
  59. SEGACT, IMODEL
  60. IF (FORMOD(1).EQ.'CHARGEMENT') THEN
  61. LIMODL(**) = IMODEL
  62. ELSE
  63. SEGDES, IMODEL
  64. ENDIF
  65. ENDDO
  66. C
  67. NSOUS = LIMODL(/1)
  68. IF (NSOUS.EQ.0) THEN
  69. SEGDES, MMODEL
  70. MOTERR(1:10)='un MMODEL '
  71. MOTERR(11:20)='CHARGEMENT'
  72. MOTERR(21:30)=' PRESSION '
  73. CALL ERREUR(881)
  74. RETURN
  75. ENDIF
  76. C
  77. C TEST DE NON REDONDANCES DES SOUS-MODELES
  78. C
  79. N1 = 1
  80. DO I = NSOUS,2,-1
  81. IMODE1 = LIMODL(I)
  82. DO J = (I-1),1,-1
  83. IMODE2 = LIMODL(J)
  84. IF (IMODE1.EQ.IMODE2) THEN
  85. LIMODL(I) = 0
  86. GOTO 10
  87. ELSE IF (IMODE1.IMAMOD.EQ.IMODE2.IMAMOD .AND.
  88. & IMODE1.CONMOD.EQ.IMODE2.CONMOD) THEN
  89. LIMODL(I) = 0
  90. GOTO 10
  91. ENDIF
  92. SEGDES,IMODE1
  93. ENDDO
  94. N1 = N1 + 1
  95. 10 CONTINUE
  96. ENDDO
  97. C
  98. C CREATION DU MMODEL
  99. C
  100. J = 0
  101. SEGINI,MMODE1
  102. DO i = 1,NSOUS
  103. IF (LIMODL(I).GT.0) THEN
  104. j = j + 1
  105. MMODE1.KMODEL(J) = LIMODL(I)
  106. ENDIF
  107. ENDDO
  108. IPMOD0 = MMODE1
  109. SEGDES, MMODE1
  110. SEGSUP, LIMODL
  111. C_______________________________________________________________________
  112. C
  113. C QUELQUES INITIALISATIONS
  114. C_______________________________________________________________________
  115. C
  116. ISUP1 = 0
  117. ISUP2 = 0
  118. IRET = 0
  119. IPCHP4 = 0
  120. MCHELM = 0
  121. MCHAML = 0
  122. IPCHE1 = 0
  123. IPCHE2 = 0
  124. C_______________________________________________________________________
  125. C
  126. C REDUCTION DES MCHAML EN ENTREE SUR LE MODELE
  127. C_______________________________________________________________________
  128. C
  129. IF (IPCHM1.NE.0) THEN
  130. CALL REDUAF(IPCHM1,IPMOD0,IPCHE1,0,IR,KER)
  131. IF (IR.NE.1) CALL ERREUR(KER)
  132. IF (IERR.NE.0) RETURN
  133. ENDIF
  134. C
  135. IF (IPCHM2.NE.0) THEN
  136. CALL REDUAF(IPCHM2,IPMOD0,IPCHE2,0,IR,KER)
  137. IF (IR.NE.1) CALL ERREUR(KER)
  138. IF (IERR.NE.0) RETURN
  139. ENDIF
  140. C_______________________________________________________________________
  141. C
  142. C VERIFICATION DES LIEUX SUPPORT DES MCHAML
  143. C_______________________________________________________________________
  144. C
  145. IF (IPCHE1.NE.0) THEN
  146. CALL QUESUP(IPMOD0,IPCHE1,3,0,ISUP1,IRET1C)
  147. IF (ISUP1.GT.1) RETURN
  148. ENDIF
  149. C
  150. IF (IPCHE2.NE.0) THEN
  151. CALL QUESUP(IPMOD0,IPCHE2,3,0,ISUP2,IR)
  152. IF (ISUP2.GT.1) RETURN
  153. ENDIF
  154. C_______________________________________________________________________
  155. C
  156. C ACTIVATION DU MODELE
  157. C_______________________________________________________________________
  158. C
  159. MMODEL=IPMOD0
  160. SEGACT, MMODEL
  161. NSOUS = MMODEL.KMODEL(/1)
  162. DO IM = 1, NSOUS
  163. IMODEL = MMODEL.KMODEL(IM)
  164. SEGACT, IMODEL
  165. ENDDO
  166. C
  167. C ACTIVATION DU MCHELM CONTENANT 'PRES'
  168. C
  169. IF (IPCHE1.NE.0) THEN
  170. MCHEL1 = IPCHE1
  171. ELSE
  172. MCHEL1 = IPCHE2
  173. ENDIF
  174. SEGACT, MCHEL1
  175. C
  176. C INITIALISATION DU MCHELM DE FORCES
  177. C
  178. N1 = NSOUS
  179. L1 = 6
  180. N3 = 5
  181. SEGINI, MCHELM
  182. IPCHE5 = MCHELM
  183. MCHELM.IFOCHE = IFOUR
  184. MCHELM.TITCHE = 'FORCES'
  185. C
  186. C=======================================================================
  187. C
  188. C BOUCLE SUR LES MODELES ELEMENTAIRES
  189. C
  190. C=======================================================================
  191. C
  192. ISOUS = 0
  193. C
  194. DO 200 KISOUS = 1, NSOUS
  195. C
  196. C INITIALISATION
  197. C
  198. IVAMAT=0
  199. IVACAR=0
  200. IVASTR=0
  201. IVAFOR=0
  202. MOMATR=0
  203. MOCARA=0
  204. MOSTRS=0
  205. MOFORC=0
  206. IPMINT=0
  207. IPMIN1=0
  208. C
  209. C TRAITEMENT DU MODELE
  210. C
  211. IMODEL = MMODEL.KMODEL(KISOUS)
  212. ISOUS = ISOUS+1
  213. MELE = IMODEL.NEFMOD
  214. C
  215. C PETITE VERIFICATION SUR LE TYPE D'ELEMENT
  216. C
  217. IF (MELE.EQ.0) THEN
  218. C
  219. C ERREUR : IMPOSSIBLE D UTILISER L OPERATEUR PRESSI POUR
  220. C LES ELEMENTS DE FORMULATION MELM
  221. C
  222. MOTERR(1:8)=NOMTP(MELM)
  223. CALL ERREUR(193)
  224. GOTO 9992
  225. ENDIF
  226. IPMAIL = IMODEL.IMAMOD
  227. CONM = IMODEL.CONMOD
  228. C
  229. C CREATION DU TABLEAU INFOS
  230. C
  231. CALL IDENT(IPMAIL,CONM,IPCHE1,IPCHE2,INFOS,IRTD)
  232. IF (IRTD.EQ.0) GOTO 9992
  233. C
  234. C NATURE DU MATERIAU
  235. C
  236. CMATE = CMATEE
  237. MATE = IMATEE
  238. INAT = INATUU
  239. C_______________________________________________________________________
  240. C
  241. C ACTIVATION DU MELEME
  242. C_______________________________________________________________________
  243. C
  244. MELEME = IPMAIL
  245. SEGACT, MELEME
  246. NBNN = MELEME.NUM(/1)
  247. NBELEM = MELEME.NUM(/2)
  248. C_______________________________________________________________________
  249. C
  250. C INFORMATIONS SUR L'ELEMENT FINI
  251. C_______________________________________________________________________
  252. C
  253. NBPGAU= INFELE(4)
  254. MINTE = INFMOD(5)
  255. MINTE1= INFMOD(8)
  256. MFR = INFELE(13)
  257. MFR2 = NUMMFR(MELE)
  258. NSTRS = INFELE(16)
  259. C
  260. IPMINT= MINTE
  261. IPMIN1= MINTE1
  262. SEGACT, MINTE
  263. IPPORE= 0
  264. IF (MFR2.EQ.33.OR.MFR2.EQ.57.OR.MFR2.EQ.59) IPPORE = NBNN
  265. C
  266. IMACHE(ISOUS) = IPMAIL
  267. INFCHE(ISOUS,1)=0
  268. INFCHE(ISOUS,2)=0
  269. INFCHE(ISOUS,3)=NIFOUR
  270. INFCHE(ISOUS,4)=0
  271. INFCHE(ISOUS,5)=0
  272. C_______________________________________________________________________
  273. C
  274. C NOMS DE COMPOSANTES DE CONTRAINTES
  275. C_______________________________________________________________________
  276. C
  277. IF (IPCHE1.NE.0) THEN
  278. MOSTRS=LNOMID(4)
  279. NOMID=MOSTRS
  280. SEGACT, NOMID
  281. NSTR=LESOBL(/2)
  282. NFAC=LESFAC(/2)
  283. C
  284. NBTYPE=1
  285. SEGINI, NOTYPE
  286. TYPE(1)='REAL*8'
  287. MOTYPE=NOTYPE
  288. C
  289. C VERIFICATION DE LEUR PRESENCE
  290. C
  291. CALL KOMCHA(IPCHE1,IPMAIL,CONM,MOSTRS,MOTYPE,1,INFOS,3,IVASTR)
  292. SEGSUP, NOTYPE
  293. SEGDES, NOMID
  294. IF (IERR.NE.0) GOTO 9991
  295. C
  296. IF (ISUP1.EQ.1) THEN
  297. NSTRS = 1
  298. CALL VALCHE(IVASTR,NSTRS,IPMINT,IPPORE,MOSTRS,MELE)
  299. ENDIF
  300. ELSE
  301. MOMATR=LNOMID(6)
  302. NOMID=MOMATR
  303. SEGACT, NOMID
  304. NMATR=LESOBL(/2)
  305. NMATF=LESFAC(/2)
  306. C
  307. NBTYPE=1
  308. SEGINI, NOTYPE
  309. TYPE(1)='REAL*8'
  310. MOTYPE=NOTYPE
  311. C
  312. C VERIFICATION DE LEUR PRESENCE
  313. C
  314. CALL KOMCHA(IPCHE2,IPMAIL,CONM,MOMATR,MOTYPE,1,INFOS,3,IVAMAT)
  315. SEGSUP, NOTYPE
  316. SEGDES, NOMID
  317. IF (IERR.NE.0) GOTO 9991
  318. C
  319. IF (ISUP2.EQ.1) THEN
  320. CALL VALCHE(IVAMAT,NMATR,IPMINT,IPPORE,MOMATR,MELE)
  321. ENDIF
  322. ENDIF
  323. C_______________________________________________________________________
  324. C
  325. C NOMS DE COMPOSANTES DE FORCES ET CREATION DU MCHAML DE FORCE
  326. C (CE MCHAML SERA TRANSFORME EN FIN DE SUBROUTINE EN CHPOINT)
  327. C_______________________________________________________________________
  328. C
  329. MOFORC = LNOMID(2)
  330. NOMID=MOFORC
  331. SEGACT, NOMID
  332. NFORC=LESOBL(/2)
  333. NFACF=LESFAC(/2)
  334. C
  335. N2=NFORC
  336. SEGINI, MCHAML
  337. ICHAML(ISOUS)=MCHAML
  338. C
  339. DO 110 ICOMP=1,NFORC
  340. NOMCHE(ICOMP)=LESOBL(ICOMP)
  341. TYPCHE(ICOMP)='REAL*8'
  342. 110 CONTINUE
  343. C
  344. IF (NFACF.NE.0) THEN
  345. IFAC = 0
  346. DO 111 ICOMP=(NFORC+1),N2
  347. IFAC = IFAC + 1
  348. NOMCHE(ICOMP)=LESFAC(IFAC)
  349. TYPCHE(ICOMP)='REAL*8'
  350. 111 CONTINUE
  351. ENDIF
  352. C
  353. SEGDES, NOMID
  354. C
  355. C TAILLES DE MELVAL
  356. C
  357. N1EL=NBELEM
  358. N1PTEL=NBNN
  359. NBPTEL=NBPGAU
  360. NEL =N1EL
  361. C
  362. C CREATION DU MELVAL DE FORCES
  363. C
  364. NS=1
  365. NCOSOU=NFORC+NFACF
  366. SEGINI, MPTVAL
  367. IVAFOR=MPTVAL
  368. DO 100 ICOMP=1,NCOSOU
  369. N2PTEL=0
  370. N2EL=0
  371. SEGINI, MELVAL
  372. IELVAL(ICOMP)=MELVAL
  373. IVAL(ICOMP)=MELVAL
  374. 100 CONTINUE
  375. C_______________________________________________________________________
  376. C
  377. C NOMS DE COMPOSANTES DE CARACTERISTIQUES (FACULTATIF)
  378. C_______________________________________________________________________
  379. C
  380. NBROBL=0
  381. NBRFAC=0
  382. NOMID=0
  383. IVECT=0
  384. C
  385. C EPAISSEUR DANS LE CAS MASSIF EN CONTRAINTES PLANES
  386. C
  387. IF(MFR.EQ.72.AND.IFOUR.EQ.-2)THEN
  388. C
  389. NBRFAC=1
  390. SEGINI, NOMID
  391. LESFAC(1)='DIM3'
  392. C
  393. NBTYPE=1
  394. SEGINI, NOTYPE
  395. TYPE(1)='REAL*8'
  396. C
  397. C EPAISSEUR ET EXCENTREMENT DANS LE CAS DES COQUES
  398. C
  399. ELSEIF (MFR.EQ.74.AND.MFR2.EQ.3) THEN
  400. NBROBL=0
  401. NBRFAC=1
  402. IF (IFOUR.EQ.-2) NBRFAC=2
  403. SEGINI, NOMID
  404. LESFAC(1)='EXCE'
  405. IF(IFOUR.EQ.-2) LESFAC(2)='DIM3'
  406. C
  407. NBTYPE=1
  408. SEGINI, NOTYPE
  409. TYPE(1)='REAL*8'
  410. C
  411. C EPAISSEUR ET RAYON EXTERNE DANS LE CAS DES TUYAUX
  412. C
  413. ELSEIF (MFR.EQ.74.AND.MFR2.EQ.13) THEN
  414. NBROBL=2
  415. NBRFAC=2
  416. SEGINI, NOMID
  417. LESOBL(1)='EPAI'
  418. LESOBL(2)='RAYO'
  419. LESFAC(1)='RACO'
  420. LESFAC(2)='VECT'
  421. C
  422. NBTYPE=4
  423. SEGINI, NOTYPE
  424. TYPE(1)='REAL*8'
  425. TYPE(2)='REAL*8'
  426. TYPE(3)='REAL*8'
  427. TYPE(4)='POINTEURPOINT '
  428. C
  429. ENDIF
  430. MOCARA=NOMID
  431. MOTYPE=NOTYPE
  432. C
  433. NCARA=NBROBL
  434. NCARF=NBRFAC
  435. NCARR=NCARA+NCARF
  436. C
  437. IF (IPCHE2.NE.0) THEN
  438. IF (MOCARA.NE.0) THEN
  439. CALL KOMCHA(IPCHE2,IPMAIL,CONM,MOCARA,MOTYPE,1,INFOS,3,
  440. + IVACAR)
  441. SEGSUP, NOTYPE
  442. SEGSUP, NOMID
  443. IF (IERR.NE.0) GOTO 9990
  444. IF (ISUP2.EQ.1) THEN
  445. CALL VALCHE(IVACAR,NCARR,IPMINT,IPPORE,MOCARA,MELE)
  446. IF (IERR.NE.0)THEN
  447. ISUP2=0
  448. GOTO 9990
  449. ENDIF
  450. ENDIF
  451. ELSE
  452. SEGSUP, NOTYPE
  453. SEGSUP, NOMID
  454. ENDIF
  455. ELSE IF (NCARA.GT.0) THEN
  456. SEGSUP, NOTYPE
  457. SEGSUP, NOMID
  458. MOTERR(1:8)='CARACTER'
  459. MOTERR(9:12)=NOMTP(MELE)
  460. MOTERR(13:20)='FEQPR'
  461. CALL ERREUR(145)
  462. GOTO 9990
  463. ENDIF
  464. C_______________________________________________________________________
  465. C
  466. C CALCUL DES FORCES EQUIVALENTES
  467. C_______________________________________________________________________
  468. C
  469. IF (IPCHE1.NE.0) THEN
  470. MPTVAL=IVASTR
  471. ELSE
  472. MPTVAL=IVAMAT
  473. ENDIF
  474. IVAPRE = IVAL(1)
  475. C
  476. C - ELEMENTS DE FORMULATION MASSIF
  477. C
  478. IF (MFR2.EQ.1) THEN
  479. C
  480. IF (MELE.EQ.2.OR.MELE.EQ.3.OR.MELE.EQ.79.OR.MELE.EQ.80) THEN
  481. CALL FPMA2D(IVAPRE,IPMAIL,IPMINT,IVAFOR,IVACAR,0)
  482. ELSE IF(MELE.EQ.4.OR.MELE.EQ.6.OR.MELE.EQ.8.OR.
  483. + MELE.EQ.10.OR.MELE.EQ.81.OR.MELE.EQ.82.OR.
  484. + MELE.EQ.83) THEN
  485. IF (IDIM.EQ.3) THEN
  486. CALL FPMA3D(IVAPRE,IPMAIL,IPMINT,IVAFOR,0)
  487. ELSE
  488. C ERREUR, APPLICATION PRESSION SUR SURFACE DANS UN PB 2D
  489. CALL ERREUR(820)
  490. GOTO 9990
  491. ENDIF
  492. ELSE IF (MELE.EQ.45) THEN
  493. CALL FPMA1D(IVAPRE,IPMAIL,IPMINT,IVAFOR,0)
  494. ELSE
  495. C ERREUR, ELEMENT NON IMPLEMENTE
  496. MOTERR(1:4)=NOMTP(MELE)
  497. MOTERR(5:12)='FEQPR '
  498. CALL ERREUR(86)
  499. GOTO 9990
  500. ENDIF
  501. C
  502. C - ELEMENTS DE FORMULATION COQUES
  503. C
  504. ELSE IF (MFR2.EQ.3.OR.MFR2.EQ.5.OR.MFR2.EQ.9) THEN
  505. C
  506. IF (MELE.EQ.27.OR.MELE.EQ.28.OR.MELE.EQ.45.OR.
  507. + MELE.EQ.93.OR.MELE.EQ.44.OR.MELE.EQ.49.OR.
  508. + MELE.EQ.41.OR.MELE.EQ.56) THEN
  509. C VERIFICATION DE L ORIENTATION DU IPMAIL
  510. CALL ECROBJ('MAILLAGE',IPMAIL)
  511. CALL VERSEN
  512. CALL LIROBJ('MAILLAGE',IPMAIL,1,IRETOU)
  513. IF (IERR.NE.0) GOTO 9990
  514. IF (MELE.EQ.27.OR.MELE.EQ.28.OR.MELE.EQ.45
  515. + .OR.MELE.EQ.93) THEN
  516. CALL FPCO3D(IVAPRE,IPMAIL,IVAFOR)
  517. ELSE IF (MELE.EQ.44) THEN
  518. CALL FPCO2D(IVAPRE,IPMAIL,IVAFOR,IVACAR)
  519. ELSE IF (MELE.EQ.49) THEN
  520. IPT1 = IPMAIL
  521. SEGACT, IPT1
  522. CALL FPCOQ4(IVAPRE,IPMAIL,IPMINT,IVAFOR)
  523. SEGDES, IPT1
  524. ELSE IF (MELE.EQ.41.OR.MELE.EQ.56) THEN
  525. IPT1 = IPMAIL
  526. SEGACT, IPT1
  527. CALL FPCOQ8(IVAPRE,IPMAIL,IPMINT,IVACAR,IPMIN1,IVAFOR)
  528. SEGDES, IPT1
  529. ENDIF
  530. ELSE
  531. C ERREUR, ELEMENT NON IMPLEMENTE
  532. MOTERR(1:4)=NOMTP(MELE)
  533. MOTERR(5:12)='FEQPR '
  534. CALL ERREUR(86)
  535. GOTO 9990
  536. ENDIF
  537. C
  538. ELSE IF (MFR2.EQ.13) THEN
  539. C
  540. C - ELEMENTS TUYAU
  541. C
  542. CALL FPELTU(IVAPRE,IVACAR,IPMAIL,I,IVAFOR)
  543. C
  544. C - SINON TENTATIVE D'UTILISATION D'UNE OPTION NON IMPLEMENTEE
  545. C
  546. ELSE
  547. CALL ERREUR(251)
  548. GOTO 9990
  549. ENDIF
  550. C
  551. C DESACTIVATION AVANT DE PASSER A LA SOUS ZONE SUIVANTE
  552. C
  553. SEGDES, MELEME
  554. SEGDES, MINTE
  555. C
  556. IF(ISUP1.EQ.1)THEN
  557. CALL DTMVAL(IVASTR,3)
  558. ELSE
  559. CALL DTMVAL(IVASTR,1)
  560. ENDIF
  561. C
  562. CALL DTMVAL(IVAFOR,1)
  563. C
  564. IF(ISUP2.EQ.1)THEN
  565. CALL DTMVAL(IVAMAT,3)
  566. CALL DTMVAL(IVACAR,3)
  567. ELSE
  568. CALL DTMVAL(IVAMAT,1)
  569. CALL DTMVAL(IVACAR,1)
  570. ENDIF
  571. C
  572. IF (IERR.NE.0) GO TO 9990
  573. C
  574. C=======================================================================
  575. C
  576. C FIN DE BOUCLE SUR LES MODELES ELEMENTAIRES
  577. C
  578. C=======================================================================
  579. 200 CONTINUE
  580. C_______________________________________________________________________
  581. C
  582. C TRANSFORMATION DU CHAMELEM EN CHPOINT
  583. C_______________________________________________________________________
  584. C
  585. CALL CHAMPO(IPCHE5,0,IPCHP4,IRETOU)
  586. CALL DTCHAM(IPCHE5)
  587. IF (IRETOU.EQ.0) GOTO 9000
  588. C
  589. C FIN NORMALE
  590. C
  591. IRET = 1
  592. GOTO 9000
  593. C_______________________________________________________________________
  594. C
  595. C ERREUR DANS UNE SOUS ZONE, DESACTIVATION ET RETOUR
  596. C_______________________________________________________________________
  597. C
  598. 9990 CONTINUE
  599. IF (MCHAML.NE.0) SEGSUP, MCHAML
  600. CALL DTMVAL(IVAFOR,3)
  601. C
  602. 9991 CONTINUE
  603. SEGDES, MINTE
  604. C
  605. 9992 CONTINUE
  606. IF (MCHELM.NE.0) SEGSUP, MCHELM
  607. SEGDES, MELEME
  608. C
  609. IRET = 0
  610. SEGSUP, MCHELM
  611. C_______________________________________________________________________
  612. C
  613. C DERNIERES DESACTIVATION AVANT DE QUITTER
  614. C_______________________________________________________________________
  615. C
  616. 9000 CONTINUE
  617. MMODEL = IPMOD0
  618. DO IM = 1, NSOUS
  619. IMODEL = MMODEL.KMODEL(IM)
  620. SEGDES, IMODEL
  621. ENDDO
  622. SEGDES, MCHEL1
  623. C
  624. RETURN
  625. END
  626.  
  627.  
  628.  
  629.  

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