Télécharger fpfiss.eso

Retour à la liste

Numérotation des lignes :

  1. C FPFISS SOURCE AM 16/04/12 21:15:42 8903
  2. SUBROUTINE FPFISS(P,IPCHE1,IPMODL,IPVECT,IPPOIN,IPCHE2,
  3. 1 IPTFP,IRET)
  4. C_____________________________________________________________________
  5. C
  6. C CALCULE LES FORCES DE PRESSIONS APPLIQUEES SUR DES LEVRES D UNE
  7. C FISSURE (ELT LINESPRING)
  8. C
  9. C ENTREES :
  10. C ---------
  11. C
  12. C P VALEUR DE LA PRESSION SI ELLE EST CONSTANTE
  13. C IPCHE1 CHPOINT CONTENANT LES VALEURS DES PRESSIONS AUX NOEUDS
  14. C IPMODL OBJET MMODEL SUR LEQUEL S APPLIQUE LA PRESSION
  15. C IPVECT VECTEUR INDIQUANT LA DIRECTION DANS LAQUELLE
  16. C S APPLIQUE LA PRESSION
  17. C IPPOIN POINT OU SE RAPPORTE LE VECTEUR
  18. C IPCHE2 MCHAML CONTENANT LES CARACTERISTIQUES
  19. C
  20. C SORTIE :
  21. C --------
  22. C
  23. C IPTFP CHPOINT DES FORCES NODALES EQUIVALENTES
  24. C IRET 1 OU 0 SUIVANT SUCCES OU NON
  25. C
  26. C REVISION JACQUELINE BROCHARD SEPTEMBRE 86
  27. C PASSAGE AUX NOUVEAUX MCHAMLS PAR JM CAMPENON LE 05 09 90
  28. C
  29. C_____________________________________________________________________
  30. IMPLICIT INTEGER(I-N)
  31. IMPLICIT REAL*8(A-H,O-Z)
  32. C
  33. -INC CCOPTIO
  34. -INC CCHAMP
  35. C
  36. -INC SMCOORD
  37. -INC SMELEME
  38. -INC SMMODEL
  39. -INC SMCHAML
  40. -INC SMCHPOI
  41. -INC SMINTE
  42. C
  43. C SEGMENT DONNANT LE POINTEUR DE MAILLAGE CORRECTE AU MCHAML DE
  44. C CARACTERISTIQUE APRES CREATION D'UN MMODEL
  45. C
  46. logical ltelq
  47. SEGMENT JPMAIL
  48. INTEGER MAIL1 (NSOUS1)
  49. INTEGER MAIL2 (NSOUS1)
  50. ENDSEGMENT
  51. *
  52. SEGMENT NOTYPE
  53. CHARACTER*16 TYPE(NBTYPE)
  54. ENDSEGMENT
  55. C
  56. SEGMENT MPTVAL
  57. INTEGER IPOS(NS) ,NSOF(NS)
  58. INTEGER IVAL(NCOSOU)
  59. CHARACTER*16 TYVAL(NCOSOU)
  60. ENDSEGMENT
  61. C
  62. DIMENSION V(3),XP(3)
  63. DIMENSION BPSS(3,3),XE(3,4),XEL(3,3),V1(3),V2(3),H1(3),H2(3)
  64. CHARACTER*8 MOT
  65. CHARACTER*(NCONCH) CONM
  66. PARAMETER ( NINF=3 )
  67. INTEGER INFOS(NINF)
  68. LOGICAL lsupfo
  69. C
  70. DATA X774/.774596669241483D0/
  71. DATA UN,UNDEMI,ZERO/1.D0,.5D0,0.D0/
  72. DATA MOT/'NOEUD '/
  73.  
  74. IRET=0
  75. C
  76. C VERIFICATION DU LIEU SUPPORT DU MCHAML DE CARACTERISTIQUES
  77. C
  78. CALL QUESUP(IPMODL,IPCHE2,3,0,ISUP,IRETOU)
  79. IF (ISUP.GT.1) RETURN
  80. C
  81. IFLAG=0
  82. NHRM=NIFOUR
  83. C
  84. C ON RECUPERE LES COORDONNEES DU VECTEUR
  85. C
  86. IREF=(IPVECT-1)*(IDIM+1)
  87. V(1)=XCOOR(IREF+1)
  88. V(2)=XCOOR(IREF+2)
  89. IF (IDIM.EQ.2) THEN
  90. VN=SQRT(V(1)**2+V(2)**2)
  91. IF (VN.EQ.0.) THEN
  92. CALL ERREUR(277)
  93. RETURN
  94. ENDIF
  95. V(1)=V(1)/VN
  96. V(2)=V(2)/VN
  97. ELSE
  98. V(3)=XCOOR(IREF+3)
  99. VN=SQRT(V(1)**2+V(2)**2+V(3)**2)
  100. IF (VN.EQ.0.) THEN
  101. CALL ERREUR(277)
  102. RETURN
  103. ENDIF
  104. V(1)=V(1)/VN
  105. V(2)=V(2)/VN
  106. V(3)=V(3)/VN
  107. ENDIF
  108. C
  109. C LE FLAG SERT A INDIQUER SI L'ON DOIT OU NON DETRUIRE LE MODELE
  110. C EN CAS DE CREATION ( 0 : DESTRUCTION D'UN MMODEL CREE )
  111. C
  112. JPMAIL=0
  113. IF (IPCHE1.NE.0) THEN
  114. C
  115. C ON CREE LE MMODEL S'ACCROCHANT AU CHPOINT
  116. C
  117. CALL NOMCOM(IPCHE1,'SCAL',IPCHE,IRETOU)
  118. IF (IERR.NE.0) RETURN
  119. C
  120. C ON CREE L OBJET MAILLAGE CONTENANT TOUS LES POINT DU CHPOINT
  121. C
  122. MCHPOI=IPCHE
  123. SEGACT MCHPOI
  124. NSOUPO=IPCHP(/1)
  125. IPGEOM = 0
  126. DO 1140 I=1,NSOUPO
  127. MSOUPO=IPCHP(I)
  128. SEGACT MSOUPO
  129. IF (IPGEOM.EQ.0) THEN
  130. IPGEOM = IGEOC
  131. ELSE
  132. IPP2 = IGEOC
  133. ltelq=.false.
  134. CALL FUSE (IPGEOM,IPP2,IPPT,ltelq)
  135. IPGEOM = IPPT
  136. ENDIF
  137. SEGDES MSOUPO
  138. 1140 CONTINUE
  139. SEGDES MCHPOI
  140. C
  141. N1=0
  142. SEGINI MMODEL
  143. IPMOD=MMODEL
  144. C
  145. MMODE1=IPMODL
  146. SEGACT MMODE1
  147. NSOUS1=MMODE1.KMODEL(/1)
  148. C
  149. C BOUCLE SUR LES SOUS ZONE GEOMETRIQUE ELEMENTAIRE
  150. C
  151. IRRT=0
  152. DO 50 ISOUS=1,NSOUS1
  153. IMODE1=MMODE1.KMODEL(ISOUS)
  154. SEGACT IMODE1
  155. ITGEOM=IMODE1.IMAMOD
  156. CALL ECROBJ('MAILLAGE',IPGEOM)
  157. CALL ECRCHA('STRI')
  158. CALL ECRCHA('APPU')
  159. CALL ECROBJ('MAILLAGE',ITGEOM)
  160. CALL EXTREL(IRR,0,IBNOR)
  161. IF (IRR.EQ.0) THEN
  162. C
  163. C ON A VERIFIER L ADHERENCE DU CHPOINT A CE MAILLAGE
  164. C
  165. CALL LIROBJ('MAILLAGE',IPOGEO,1,IRETOU)
  166. IF (IERR.NE.0) THEN
  167. SEGDES MMODE1
  168. SEGDES IMODE1
  169. SEGSUP MMODEL
  170. RETURN
  171. ENDIF
  172. N1=N1+1
  173. SEGADJ MMODEL
  174. C
  175. C CREATION DE L'OBJET IMODEL DE CETTE SOUS ZONE
  176. C
  177. NFOR=IMODE1.FORMOD(/2)
  178. NMAT=IMODE1.MATMOD(/2)
  179. MN3 =IMODE1.INFMOD(/1)
  180. NPARMO=0
  181. nobmod=0
  182. C
  183. SEGINI IMODEL
  184. conmod(17:24)=' '
  185. IMAMOD=IPOGEO
  186. NEFMOD=IMODE1.NEFMOD
  187. CONMOD=IMODE1.CONMOD
  188. IPDPGE=IMODE1.IPDPGE
  189. C
  190. C CREATION D'UN TABLEAU DE CORRESPONDANCE LE IMAMOD DU
  191. C MMODEL (IPMODL) ET DU IMAMOD DU NVX MMODEL QUE L'ON CREE
  192. C
  193. IF (JPMAIL.EQ.0) SEGINI JPMAIL
  194. MAIL1(ISOUS)=ITGEOM
  195. MAIL2(ISOUS)=IPOGEO
  196. DO 47 I=1,MN3
  197. INFMOD(I)=IMODE1.INFMOD(I)
  198. 47 CONTINUE
  199. CONMOD=IMODE1.CONMOD
  200. DO 48 I=1,NFOR
  201. FORMOD(I)=IMODE1.FORMOD(I)
  202. 48 CONTINUE
  203. DO 49 I=1,NMAT
  204. MATMOD(I)=IMODE1.MATMOD(I)
  205. 49 CONTINUE
  206. KMODEL(N1)=IMODEL
  207. SEGDES IMODEL
  208. ELSE
  209. C
  210. C LE CHPOINT N'ADHERE PAS A CETTE ZONE
  211. C
  212. IRRT=IRRT+1
  213. ENDIF
  214. SEGDES IMODE1
  215. 50 CONTINUE
  216. SEGDES MMODE1
  217. SEGDES MMODEL
  218. C
  219. IF (NSOUPO.GT.1) THEN
  220. MELEME=IPGEOM
  221. SEGSUP MELEME
  222. ENDIF
  223. C
  224. IF (IRRT.EQ.NSOUS1) THEN
  225. C
  226. C L'OBJET MAILLAGE ET LE CHPOINT SONT INCOMPATIBLES
  227. C
  228. MOTERR(1:8)='MAILLAGE'
  229. MOTERR(9:16)='CHPOINT'
  230. CALL ERREUR(135)
  231. MMODEL=IPMOD
  232. SEGSUP MMODEL
  233. RETURN
  234. ENDIF
  235. C
  236. CALL CHAME1(0,IPMOD,IPCHE,' ',IPCH1,3)
  237. IF (IERR.NE.0) THEN
  238. CALL DTMODL(IPMOD)
  239. SEGSUP JPMAIL
  240. RETURN
  241. ENDIF
  242. ELSE
  243. IFLAG=1
  244. IPMOD=IPMODL
  245. CALL ZEROP(IPMOD,MOT,IPCH1)
  246. IF (IERR.NE.0) RETURN
  247. MCHEL1=IPCH1
  248. SEGACT MCHEL1
  249. NSOUS=MCHEL1.ICHAML(/1)
  250. DO 11 ISOUS=1,NSOUS
  251. MCHAM1=MCHEL1.ICHAML(ISOUS)
  252. SEGACT MCHAM1
  253. MELVA1=MCHAM1.IELVAL(1)
  254. SEGACT MELVA1
  255. N1PTEL=MELVA1.VELCHE(/1)
  256. N1EL =MELVA1.VELCHE(/2)
  257. DO 9 IGAU=1,N1PTEL
  258. DO 9 IB=1,N1EL
  259. MELVA1.VELCHE(IGAU,IB)=P
  260. 9 CONTINUE
  261. SEGDES MELVA1
  262. SEGDES MCHAM1
  263. 11 CONTINUE
  264. SEGDES MCHEL1
  265. ENDIF
  266. C
  267. C ACTIVATION DU MODEL
  268. C
  269. MMODEL=IPMOD
  270. SEGACT MMODEL
  271. NSOUS=KMODEL(/1)
  272. C
  273. C CREATION DU MCHELM DES FORCES NODALES
  274. C
  275. N1=NSOUS
  276. L1=5
  277. N3=5
  278. SEGINI MCHELM
  279. IPCHEL=MCHELM
  280. TITCHE='FORCE'
  281. IFOCHE=IFOUR
  282. C_______________________________________________________________________
  283. C
  284. C BOUCLE SUR LES SOUS ZONES DU MAILLAGE
  285. C_______________________________________________________________________
  286. C
  287. DO 500 ISOUS=1,NSOUS
  288. C
  289. C ON RECUPERE L INFORMATION GENERALE
  290. C
  291. IMODEL=KMODEL(ISOUS)
  292. SEGACT IMODEL
  293. IPMAIL=IMAMOD
  294. CONM =CONMOD
  295. IMACHE(ISOUS)=IPMAIL
  296. C
  297. C TRAITEMENT DU MODEL
  298. C
  299. MELE=NEFMOD
  300. C
  301. C ERREUR L ELEMENT N EST PAS ENCORE IMPLEMENTE
  302. IF (MELE.NE.30) THEN
  303. MOTERR(1:4)=NOMTP(MELE)
  304. MOTERR(5:12)='FPFISS'
  305. CALL ERREUR(86)
  306. SEGDES IMODEL,MMODEL
  307. SEGSUP MCHELM
  308. IF (IFLAG.EQ.0) CALL DTMODL (IPMOD)
  309. IF (JPMAIL.NE.0) SEGSUP JPMAIL
  310. CALL DTCHAM(IPCH1)
  311. RETURN
  312. ENDIF
  313. C
  314. MELEME=IMAMOD
  315. IPTGEO=MELEME
  316. C
  317. C INFORMATION SUR L'ELEMENT FINI
  318. C
  319. * CALL ELQUOI(MELE,0,3,IPINF,IMODEL)
  320. * IF (IERR.NE.0) THEN
  321. * SEGDES IMODEL,MMODEL
  322. * SEGSUP MCHELM
  323. * IF (IFLAG.EQ.0) CALL DTMODL (IPMOD)
  324. * IF (JPMAIL.NE.0) SEGSUP JPMAIL
  325. * CALL DTCHAM(IPCH1)
  326. * RETURN
  327. * ENDIF
  328. * INFO=IPINF
  329. MFR =INFELE(13)
  330. * IPTINT=INFELE(11)
  331. IPTINT=infmod(5)
  332. MINTE=IPTINT
  333. SEGACT,MINTE
  334. * SEGSUP INFO
  335. C
  336. C CREATION DU TABLEAU INFOS
  337. C
  338. CALL IDENT(IPMAIL,CONM,IPCH1,IPCHE2,INFOS,IRTD)
  339. IF (IRTD.EQ.0) THEN
  340. SEGDES,MINTE
  341. SEGDES IMODEL,MMODEL
  342. SEGSUP MCHELM
  343. IF (IFLAG.EQ.0) CALL DTMODL (IPMOD)
  344. IF (JPMAIL.NE.0) SEGSUP JPMAIL
  345. CALL DTCHAM(IPCH1)
  346. RETURN
  347. ENDIF
  348. C
  349. INFCHE(ISOUS,1)=0
  350. INFCHE(ISOUS,2)=0
  351. INFCHE(ISOUS,3)=NHRM
  352. INFCHE(ISOUS,4)=IPTINT
  353. INFCHE(ISOUS,5)=0
  354. C
  355. C RECHERCHE DU MELVAL DU CHAMELEM DE PRESSION
  356. C
  357. NCARA=0
  358. NCARF=0
  359. MOCARA=0
  360. NFOR=0
  361. MOFORC=0
  362. C
  363. NBROBL=1
  364. NBRFAC=0
  365. SEGINI NOMID
  366. MOSCAL=NOMID
  367. LESOBL(1)='SCAL'
  368. *
  369. NBTYPE=1
  370. SEGINI NOTYPE
  371. MOTYPE=NOTYPE
  372. TYPE(1)='REAL*8'
  373. *
  374. CALL KOMCHA(IPCH1,IPMAIL,CONM,MOSCAL,MOTYPE,1,INFOS,3,IVASCA)
  375. SEGSUP NOTYPE
  376. IF (IERR.NE.0) GOTO 9990
  377. SEGDES NOMID
  378. MPTVAL=IVASCA
  379. IPTVPR=IVAL(1)
  380. C
  381. C CALCUL DES FORCES NODALES EQUIVALENTES
  382. C BRANCHEMENT SUIVANT LE TYPE DES ELEMENTS
  383. C
  384. C RECHERCHE DES NOM DE COMPOSANTES
  385. C
  386. if(lnomid(2).ne.0) then
  387. nomid=lnomid(2)
  388. segact nomid
  389. moforc=nomid
  390. nfor=lesobl(/2)
  391. nfac=0
  392. lsupfo=.false.
  393. else
  394. lsupfo=.true.
  395. CALL IDFORC(MFR,IFOUR,MOFORC,NFOR,NFAC)
  396. endif
  397. C
  398. C ELEMENT LINESPRING
  399. C
  400. SEGACT MELEME
  401. NBNN =NUM(/1)
  402. NBELEM=NUM(/2)
  403. IPPORE=0
  404. IF(MFR.EQ.33) IPPORE=NBNN
  405.  
  406. C CREATION DU MCHAML DE LA SOUS ZONE
  407. C
  408. C INIT DU MELVAL DEVANT CONTENIR LES FORCES DE PRESSION
  409. C
  410. N1PTEL=4
  411. N1EL=NBELEM
  412. N2PTEL=0
  413. N2EL=0
  414. C
  415. N2=NFOR
  416. SEGINI MCHAML
  417. ICHAML(ISOUS)=MCHAML
  418. NS=1
  419. NCOSOU=NFOR
  420. SEGINI MPTVAL
  421. IVAFOR=MPTVAL
  422. NOMID=MOFORC
  423. SEGACT NOMID
  424. DO 1100 ICOMP=1,NFOR
  425. NOMCHE(ICOMP)=LESOBL(ICOMP)
  426. TYPCHE(ICOMP)='REAL*8'
  427. SEGINI MELVAL
  428. IELVAL(ICOMP)=MELVAL
  429. IVAL(ICOMP)=MELVAL
  430. 1100 CONTINUE
  431. SEGDES NOMID
  432. C
  433. C TRAITEMENT DES CHAMPS DE CARACTERISTIQUES POUR LES LINESPRING
  434. C
  435. NBROBL=5
  436. NBRFAC=0
  437. SEGINI NOMID
  438. MOCARA=NOMID
  439. LESOBL(1)='EPAI'
  440. LESOBL(2)='FISS'
  441. LESOBL(3)='VX '
  442. LESOBL(4)='VY '
  443. LESOBL(5)='VZ '
  444. IF (JPMAIL.NE.0) THEN
  445. C
  446. C ON RECUPERE LE IMAMOD DU MMODEL D'ORIGINE POUR QUE LE
  447. C DONNE CORRESPONDE A CELUI DE IPCHE21
  448. C
  449. DO 60 KISOUS=1,NSOUS1
  450. IF (IPMAIL.EQ.MAIL2(KISOUS)) THEN
  451. IPMAI1=MAIL1(KISOUS)
  452. GOTO 61
  453. ENDIF
  454. 60 CONTINUE
  455. C
  456. C NE DOIT NORMALEMENT JAMAIS SE PRODUIRE
  457. C
  458. CALL ERREUR (472)
  459. GOTO 9990
  460. ELSE
  461. IPMAI1=IPMAIL
  462. ENDIF
  463. 61 CONTINUE
  464. C
  465. NBTYPE=1
  466. SEGINI NOTYPE
  467. MOTYPE=NOTYPE
  468. TYPE(1)='REAL*8'
  469. *
  470. CALL KOMCHA(IPCHE2,IPMAI1,CONM,MOCARA,MOTYPE,
  471. 1 1,INFOS,3,IVACAR)
  472. SEGSUP NOTYPE
  473. IF (IERR.NE.0) GOTO 9990
  474. C
  475. NCARA=NBROBL
  476. NCARF=NBRFAC
  477. NCARR=NCARA+NCARF
  478. C
  479. IF (ISUP.EQ.1) THEN
  480. CALL VALCHE(IVACAR,NCARR,IPTINT,IPPORE,MOCARA,MELE)
  481. ENDIF
  482. C
  483. SEGDES NOMID
  484. C
  485. C ELEMENT LINESPRING
  486. C
  487. CALL FPLISP(IPTVPR,IPTGEO,IPTINT,IVACAR,IVAFOR)
  488. C
  489. C DESACTIVATION DES SEGMENT PROPRE A LA GEOMETRIE ISOUS
  490. C
  491. SEGDES,MINTE
  492. SEGDES IMODEL
  493. SEGDES MCHAML
  494. C
  495. IF (ISUP.EQ.1) THEN
  496. CALL DTMVAL(IVACAR,3)
  497. ELSE
  498. CALL DTMVAL(IVACAR,1)
  499. ENDIF
  500. C
  501. CALL DTMVAL(IVAFOR,1)
  502. C
  503. CALL DTMVAL(IVASCA,1)
  504. C
  505. NOMID=MOFORC
  506. if(lsupfo)SEGSUP NOMID
  507. NOMID=MOCARA
  508. SEGSUP NOMID
  509. NOMID=MOSCAL
  510. SEGSUP NOMID
  511. C
  512. SEGDES MELEME
  513. C
  514. 500 CONTINUE
  515. SEGDES MMODEL
  516. CALL DTCHAM(IPCH1)
  517. IF (IFLAG.EQ.0) CALL DTMODL(IPMOD)
  518. IF (JPMAIL.NE.0) SEGSUP JPMAIL
  519. C
  520. C ON TRANSFORME LE CHAM/ELEM EN CHAM/POIN
  521. C
  522. C* SEGDES MCHELM
  523. CALL CHAMPO(IPCHEL,0,IPTFP,IRETOU)
  524. CALL DTCHAM(IPCHEL)
  525. IF (IRETOU.EQ.0) RETURN
  526. C
  527. C ON COMPARE LE SENS DE LA FORCE AU SENS DU VECTEUR AU POINT INDIQUE
  528. C
  529. MCHPOI=IPTFP
  530. SEGACT MCHPOI
  531. DO 201 I=1,IPCHP(/1)
  532. MSOUPO=IPCHP(I)
  533. SEGACT MSOUPO
  534. MELEME=IGEOC
  535. SEGACT MELEME
  536. DO 202 K=1,NUM(/2)
  537. IF (NUM(1,K).EQ.IPPOIN) GO TO 205
  538. 202 CONTINUE
  539. SEGDES MSOUPO,MELEME
  540. 201 CONTINUE
  541. C
  542. C LE POINT DONNE N APPARTIENT PAS A LA STRUCTURE
  543. C
  544. INTERR(1)=IPPOIN
  545. MOTERR(1:8)=' '
  546. CALL ERREUR(64)
  547. SEGDES MCHPOI
  548. RETURN
  549. C
  550. 205 CONTINUE
  551. SEGDES MELEME
  552. MPOVAL=IPOVAL
  553. SEGACT MPOVAL
  554. FN2=ZERO
  555. DO 210 J=1,IDIM
  556. r_z = VPOCHA(K,J)
  557. FN2=FN2 + r_z*r_z
  558. TEST=TEST+ V(J)*r_z
  559. 210 CONTINUE
  560. FN=SQRT(FN2)
  561. SEGDES MPOVAL,MSOUPO,MCHPOI
  562. C
  563. C ERREUR IMPOSSIBLE D ORIENTER LES FORCES DE PRESSION
  564. C
  565. IF (ABS(TEST).LE.0.025*FN) THEN
  566. CALL ERREUR(192)
  567. RETURN
  568. ENDIF
  569. IF (TEST.LE.0.) THEN
  570. XFLOT=-UN
  571. CALL MUCHPO(IPTFP,XFLOT,IPTFP0,1)
  572. CALL DTCHPO(IPTFP)
  573. IPTFP=IPTFP0
  574. ENDIF
  575. IRET = 1
  576. RETURN
  577. C
  578. C ERREUR DANS UNE SOUS ZONE / DESACTIVATION ET RETOUR
  579. C
  580. 9990 CONTINUE
  581. IRET=0
  582. IF (IFLAG.EQ.0) CALL DTMODL(IPMOD)
  583. IF (JPMAIL.NE.0) SEGSUP JPMAIL
  584. C
  585. SEGSUP MCHELM
  586. C
  587. IF (ISUP.EQ.1) THEN
  588. CALL DTMVAL(IVACAR,3)
  589. ELSE
  590. CALL DTMVAL(IVACAR,1)
  591. ENDIF
  592. C
  593. CALL DTMVAL(IVAFOR,3)
  594. C
  595. CALL DTMVAL(IVASCA,1)
  596. C
  597. NOMID=MOCARA
  598. IF (MOCARA.NE.0) SEGSUP NOMID
  599. NOMID=MOFORC
  600. IF (lsupfo.and.MOFORC.NE.0) SEGSUP NOMID
  601. NOMID=MOSCAL
  602. SEGSUP NOMID
  603. C
  604. SEGDES,MINTE
  605. SEGDES IMODEL
  606. SEGDES MMODEL
  607. C
  608. CALL DTCHAM(IPCH1)
  609. C
  610. RETURN
  611. END
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  

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