Télécharger prigfm.eso

Retour à la liste

Numérotation des lignes :

prigfm
  1. C PRIGFM SOURCE OF166741 24/12/13 21:17:11 12097
  2. SUBROUTINE PRIGFM()
  3. C************************************************************************
  4. C
  5. C PROJET : CASTEM 2000
  6. C
  7. C NOM : PRIGFM (OPERATEUR GIBIANE)
  8. C
  9. C DESCRIPTION : Voir PRIMIT
  10. C
  11. C GFMP. Calcul des variables primitives.
  12. C
  13. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  14. C
  15. C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
  16. C
  17. C************************************************************************
  18. C
  19. C
  20. C APPELES (E/S) : LIROBJ, QUEPOI, QUEPO1, ERREUR, ECROBJ, LIRCHA,
  21. C ECRCHA
  22. C
  23. C APPELES (Calcul) : PRGFM1
  24. C
  25. C
  26. C************************************************************************
  27. C
  28. C PHRASE D'APPEL (GIBIANE) :
  29. C
  30. C RCHPO0 [=vn] RCHPO1[=pn] (RCHPO2 [=yn]) = 'PRIM' 'GFMP' TAB1
  31. C CHPO0 [=phi] CHPO1 [=rn] CHPO2 [=gn]
  32. C CHPO3 [=ren] (CHPO4 [=ryn]) (CHPO5 [=aln]) ;
  33. C
  34. C ENTREES :
  35. C
  36. C TAB1 : TABLE qui contient :
  37. C * les noms des espèces qui apparaissent explicitement
  38. C dans les équations d'Euler en TAB1 . 'ESPEULE'
  39. C (LISTMOTS) ;
  40. C * le nom de l'espèce qui n'y est pas
  41. C (TAB1 . 'ESPNEULE' (MOT)) ;
  42. C * les gamma et leS pinf dans la zone phi < 0
  43. C TAB1 . 'MGAM' (LISTREEL) ;
  44. C TAB1 . 'MPIN' (LISTREEL) ;
  45. C * les gamma et les pinf dans la zone phi > 0
  46. C TAB1 . 'PGAM' (LISTREEL) ;
  47. C TAB1 . 'PPIN' (LISTREEL) ;
  48. C NB
  49. C La première valeur dans les objets LISTEEL
  50. C TAB1 . 'MGAM', ... est celle de l'espèce
  51. C TAB1 . 'ESPNEULE'; les autres sont celles des
  52. C espèces TAB1 . 'ESPEULE'.
  53. C
  54. C CHPO0 : CHPOINT contenant la fonction phi
  55. C (une composante, 'SCAL').
  56. C
  57. C CHPO1 : CHPOINT contenant la masse volumique
  58. C (une composante, 'SCAL').
  59. C
  60. C CHPO2 : CHPOINT contenant les débits
  61. C (2 composantes en 2D, 'UX ','UY ');
  62. C
  63. C CHPO3 : CHPOINT contenant l'énergie totale per
  64. C unité de volume (RHO Et),
  65. C (une composante, 'SCAL').
  66. C
  67. C CHPO4 : CHPOINT contenant les masses volumiques
  68. C des espèces en (TAB1. 'ESPEULE')
  69. C (composante dans TAB1. 'ESPEULE').
  70. C
  71. C CHPO5 : CHPOINT contenant les fractions volumiques
  72. C des espèces en (TAB1. 'ESPEULE')
  73. C (composante dans TAB1. 'ESPEULE').
  74. C
  75. C
  76. C MCLE2 : Option personnelle: pas dans la notice
  77. C officielle!!!
  78. C Mot clé, 'TRICHE' (s'il y a un erreur,
  79. C les objets RCHPO1 et RCHPO2 ne sont pas
  80. C des type ANNULE et le programme
  81. C ne s'arrête pas!!!)
  82. C
  83. C SORTIES :
  84. C
  85. C RCHPO0 : CHPOINT contenant la vitesse
  86. C
  87. C RCHPO1 : CHPOINT contenant la pression du gaz;
  88. C
  89. C RCHPO2 : CHPOINT contenant les fractions massiques
  90. C des espèces en (TAB1. 'ESPEULE')
  91. C (composante dans TAB1. 'ESPEULE').
  92. C
  93. C**** Les variables
  94. C
  95. IMPLICIT INTEGER(I-N)
  96. INTEGER NESP, NESP1
  97. INTEGER ICOND, IRETOU, INDIC, NBCOMP
  98. & , ICEN, IRO, IROVIT, IROET
  99. & , IPRES, IVIT, JGN, JGM, IPHI
  100. & , IPGAS, IROY, IALP, IY
  101. REAL*8 VALER(2)
  102. CHARACTER*(40) MESERR(2),MESCEL
  103. CHARACTER*(4) MOT(1)
  104. CHARACTER*(6) NOMTRI
  105. CHARACTER*(8) MTYPR, MTYPI
  106. LOGICAL LOGNEG, LOGTRI
  107. C
  108. C**** Variables en ACCTAB
  109. C
  110. INTEGER IVALI, IRETI,IVALR, IRETR
  111. REAL*8 XVALI,XVALR
  112. LOGICAL LOGII, LOGIR
  113. CHARACTER*(8) CHARR
  114. C
  115. C**** Les Includes
  116. C
  117.  
  118. -INC PPARAM
  119. -INC CCOPTIO
  120. -INC SMCHPOI
  121. -INC SMLMOTS
  122. -INC SMLREEL
  123. POINTEUR MLRMGA.MLREEL, MLRPGA.MLREEL,
  124. & MLRMPI.MLREEL, MLRPPI.MLREEL
  125. C
  126. C
  127. C**** Initialisation des parametres d'erreur
  128. C
  129. LOGNEG = .FALSE.
  130. MESCEL = ' '
  131. MESERR(1) = MESCEL
  132. MESERR(2) = MESCEL
  133. MOTERR(1:40) = MESCEL
  134. VALER(1) = 0.0D0
  135. VALER(2) = 0.0D0
  136. C
  137. C**** Initialisation des MOT(1)
  138. C
  139. MOT(1) = ' '
  140. C
  141. C**** Initialisation des variables en ACCTAB
  142. C
  143. IVALI = 0
  144. IVALR = 0
  145. XVALI = 0.0D0
  146. XVALR = 0.0D0
  147. LOGII = .FALSE.
  148. LOGIR = .FALSE.
  149. IRETI = 0
  150. IRETR = 0
  151. CHARR = ' '
  152. C
  153. C**************************************************
  154. C**** Lecture de la table des proprietes du gaz ***
  155. C**************************************************
  156. C
  157. ICOND = 1
  158. CALL QUETYP(MTYPR,ICOND,IRETOU)
  159. IF(IERR .NE. 0)GOTO 9999
  160. IF(MTYPR .NE. 'TABLE ')THEN
  161. C
  162. C******* Message d'erreur standard
  163. C 37 2
  164. C On ne trouve pas d'objet de type %m1:8
  165. C
  166. MOTERR(1:8) = 'TABLE '
  167. CALL ERREUR(37)
  168. GOTO 9999
  169. ELSE
  170. ICOND = 1
  171. CALL LIROBJ(MTYPR,IPGAS,ICOND,IRETOU)
  172. CALL ACTOBJ(MTYPR,IPGAS,1)
  173. IF(IERR .NE. 0)GOTO 9999
  174. ENDIF
  175. C
  176. C**** Nom de l'espece qui n'est pas dans les equations d'Euler
  177. C
  178. MTYPI = 'MOT '
  179. MTYPR = ' '
  180. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'ESPNEULE', LOGII,IRETI,
  181. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  182. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'MOT ')) THEN
  183. C
  184. C******* Message d'erreur standard
  185. C -301 0 %m1:40
  186. C
  187. MOTERR = 'TAB1 . ESPNEULE = ??? '
  188. WRITE(IOIMP,*) MOTERR(1:40)
  189. C
  190. C********** Message d'erreur standard
  191. C 21 2
  192. C Données incompatibles
  193. C
  194. CALL ERREUR(21)
  195. GOTO 9999
  196. ENDIF
  197. C
  198. C**** Les especes qui sont dans les Equations d'Euler
  199. C
  200. MTYPR = ' '
  201. CALL ACMO(IPGAS,'ESPEULE',MTYPR,MLMOT1)
  202. IF(MTYPR .NE. 'LISTMOTS')THEN
  203. C
  204. C******* Message d'erreur standard
  205. C -301 0 %m1:40
  206. C
  207. MOTERR(1:40) = 'TAB1 . ESPEULE = ??? '
  208. WRITE(IOIMP,*) MOTERR(1:40)
  209. C
  210. C******* Message d'erreur standard
  211. C 21 2
  212. C Données incompatibles
  213. C
  214. CALL ERREUR(21)
  215. GOTO 9999
  216. ENDIF
  217. SEGACT MLMOT1
  218. NESP = MLMOT1.MOTS(/2)
  219. SEGDES MLMOT1
  220. C
  221. C JGN = 4
  222. C JGM = NESP + 1
  223. C SEGINI MLMOT2
  224. C DO I1 = 1, NESP, 1
  225. C MLMOT2.MOTS(I1) = MLMOT1.MOTS(I1)
  226. C ENDDO
  227. C MLMOT2.MOTS(NESP+1)=CHARR(1:4)
  228. C DO I1 = 1, (NESP + 1)
  229. C write(*,*) MLMOT2.MOTS(I1)
  230. C ENDDO
  231. C SEGSUP MLMOT2
  232. C stop
  233. C
  234. C
  235. C**** List de gamma, PHI < 0 (MLRMGA, MLRrel Minus GAMMA)
  236. C
  237. MTYPI = 'MOT '
  238. MTYPR = ' '
  239. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'MGAM', LOGII,IRETI,
  240. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  241. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  242. C
  243. C********** Message d'erreur standard
  244. C -301 0 %m1:40
  245. C
  246. MOTERR = ' TAB1 . MGAM = ??? '
  247. WRITE(IOIMP,*) MOTERR
  248. C
  249. C********** Message d'erreur standard
  250. C 21 2
  251. C Données incompatibles
  252. C
  253. CALL ERREUR(21)
  254. GOTO 9999
  255. ENDIF
  256. MLRMGA = IRETR
  257. SEGACT MLRMGA
  258. NESP1 = MLRMGA.PROG(/1)
  259. IF(NESP1 .NE. (NESP + 1))THEN
  260. C
  261. C********** Message d'erreur standard
  262. C -301 0 %m1:40
  263. C
  264. MOTERR = 'DIME(TAB1.MGAM) != NESP '
  265. WRITE(IOIMP,*) MOTERR(1:40)
  266. C
  267. C********** Message d'erreur standard
  268. C 21 2
  269. C Données incompatibles
  270. C
  271. CALL ERREUR(21)
  272. GOTO 9999
  273. ENDIF
  274. C
  275. C**** List de gamma, PHI > 0 (MLRPGA, MLRrel Plus GAMMA)
  276. C
  277. MTYPI = 'MOT '
  278. MTYPR = ' '
  279. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'PGAM', LOGII,IRETI,
  280. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  281. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  282. C
  283. C********** Message d'erreur standard
  284. C -301 0 %m1:40
  285. C
  286. MOTERR = ' TAB1 . PGAM = ??? '
  287. WRITE(IOIMP,*) MOTERR
  288. C
  289. C********** Message d'erreur standard
  290. C 21 2
  291. C Données incompatibles
  292. C
  293. CALL ERREUR(21)
  294. GOTO 9999
  295. ENDIF
  296. MLRPGA = IRETR
  297. SEGACT MLRPGA
  298. NESP1 = MLRPGA.PROG(/1)
  299. IF(NESP1 .NE. (NESP + 1))THEN
  300. C
  301. C********** Message d'erreur standard
  302. C -301 0 %m1:40
  303. C
  304. MOTERR = 'DIME(TAB1.PGAM) != NESP '
  305. WRITE(IOIMP,*) MOTERR(1:40)
  306. C
  307. C********** Message d'erreur standard
  308. C 21 2
  309. C Données incompatibles
  310. C
  311. CALL ERREUR(21)
  312. GOTO 9999
  313. ENDIF
  314. C
  315. C**** List de PINF, PHI < 0 (MLRMPI)
  316. C
  317. MTYPI = 'MOT '
  318. MTYPR = ' '
  319. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'MPIN', LOGII,IRETI,
  320. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  321. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  322. C
  323. C********** Message d'erreur standard
  324. C -301 0 %m1:40
  325. C
  326. MOTERR = ' TAB1 . MPIN = ??? '
  327. WRITE(IOIMP,*) MOTERR
  328. C
  329. C********** Message d'erreur standard
  330. C 21 2
  331. C Données incompatibles
  332. C
  333. CALL ERREUR(21)
  334. GOTO 9999
  335. ENDIF
  336. MLRMPI = IRETR
  337. SEGACT MLRMPI
  338. NESP1 = MLRMPI.PROG(/1)
  339. IF(NESP1 .NE. (NESP + 1))THEN
  340. C
  341. C********** Message d'erreur standard
  342. C -301 0 %m1:40
  343. C
  344. MOTERR = 'DIME(TAB1.MPIN) != NESP '
  345. WRITE(IOIMP,*) MOTERR(1:40)
  346. C
  347. C********** Message d'erreur standard
  348. C 21 2
  349. C Données incompatibles
  350. C
  351. CALL ERREUR(21)
  352. GOTO 9999
  353. ENDIF
  354. C
  355. C**** List de PINF, PHI > 0 (MLRPPI)
  356. C
  357. MTYPI = 'MOT '
  358. MTYPR = ' '
  359. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'PPIN', LOGII,IRETI,
  360. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  361. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  362. C
  363. C********** Message d'erreur standard
  364. C -301 0 %m1:40
  365. C
  366. MOTERR = ' TAB1 . PPIN = ??? '
  367. WRITE(IOIMP,*) MOTERR
  368. C
  369. C********** Message d'erreur standard
  370. C 21 2
  371. C Données incompatibles
  372. C
  373. CALL ERREUR(21)
  374. GOTO 9999
  375. ENDIF
  376. MLRPPI = IRETR
  377. SEGACT MLRPPI
  378. NESP1 = MLRPPI.PROG(/1)
  379. IF(NESP1 .NE. (NESP + 1))THEN
  380. C
  381. C********** Message d'erreur standard
  382. C -301 0 %m1:40
  383. C
  384. MOTERR = 'DIME(TAB1.PPIN) != NESP '
  385. WRITE(IOIMP,*) MOTERR(1:40)
  386. C
  387. C********** Message d'erreur standard
  388. C 21 2
  389. C Données incompatibles
  390. C
  391. CALL ERREUR(21)
  392. GOTO 9999
  393. ENDIF
  394. C
  395. C***********************************************************************
  396. C******************** Lecture du reste *********************************
  397. C***********************************************************************
  398. C
  399. C**** Lecture du CHPOINT ROPHI
  400. C
  401. ICOND = 1
  402. CALL LIROBJ('CHPOINT',IPHI,ICOND,IRETOU)
  403. CALL ACTOBJ('CHPOINT',IPHI,1)
  404. IF (IERR.NE.0) GOTO 9999
  405. C
  406. C**** On cherche le pointeur de son maillage et on l'impose sur les
  407. C autres CHPOINT
  408. C
  409. MCHPOI = IPHI
  410. SEGACT MCHPOI
  411. MSOUPO = MCHPOI.IPCHP(1)
  412. SEGACT MSOUPO
  413. ICEN = MSOUPO.IGEOC
  414. SEGDES MSOUPO
  415. SEGDES MCHPOI
  416. C
  417. C**** Control du CHPOINT: QUEPOI
  418. C
  419. C On controlle que le chpoint est non-partitione
  420. C
  421. C INDIC = 1 -> on impose le pointeur du support geometrique (ICEN)
  422. C INDIC = 0 -> on ne fait que verifier le support geometrique (ICEN)
  423. C
  424. C NBCOMP > 0 -> numero des composantes
  425. C
  426. C MOT(1) = ' ' obligatoire s'on connais pas les noms des composantes
  427. C
  428. INDIC = 1
  429. NBCOMP = 1
  430. MOT(1) = 'SCAL'
  431. CALL QUEPOI(IPHI, ICEN, INDIC, NBCOMP, MOT)
  432. IF(IERR .NE. 0)THEN
  433. C
  434. C******** Message d'erreur standard
  435. C -301 0 %m1:40
  436. C
  437. MOTERR = 'CHPO0 = ??? '
  438. WRITE(IOIMP,*) MOTERR(1:40)
  439.  
  440. GOTO 9999
  441. ENDIF
  442. C
  443. C**** Lecture du CHPOINT RO
  444. C
  445. ICOND = 1
  446. CALL LIROBJ('CHPOINT ',IRO,ICOND,IRETOU)
  447. CALL ACTOBJ('CHPOINT ',IRO,1)
  448. IF (IERR.NE.0) GOTO 9999
  449. C
  450. C**** Control du CHPOINT
  451. C
  452. INDIC = 1
  453. NBCOMP = 1
  454. MOT(1) = 'SCAL'
  455. CALL QUEPOI(IRO, ICEN, INDIC, NBCOMP, MOT)
  456. IF(IERR .NE. 0)THEN
  457. C
  458. C******** Message d'erreur standard
  459. C -301 0 %m1:40
  460. C
  461. MOTERR = 'CHPO1 = ??? '
  462. WRITE(IOIMP,*) MOTERR(1:40)
  463.  
  464. GOTO 9999
  465. ENDIF
  466. C
  467. C**** Lecture du CHPOINT DEBIT.
  468. C
  469. ICOND = 1
  470. CALL LIROBJ('CHPOINT ',IROVIT,ICOND,IRETOU)
  471. CALL ACTOBJ('CHPOINT ',IROVIT,1)
  472. IF (IERR.NE.0) GOTO 9999
  473. C
  474. C**** Control du CHPOINT
  475. C
  476. INDIC = 1
  477. NBCOMP = IDIM
  478. JGN = 4
  479. JGM = IDIM
  480. SEGINI MLMOTS
  481. MLMOTS.MOTS(1) = 'UX '
  482. MLMOTS.MOTS(2) = 'UY '
  483. IF(IDIM .EQ. 3) MLMOTS.MOTS(3) = 'UZ '
  484. CALL QUEPO1(IROVIT, ICEN, MLMOTS)
  485. IF(IERR .NE. 0)THEN
  486. C
  487. C******** Message d'erreur standard
  488. C -301 0 %m1:40
  489. C
  490. MOTERR = 'CHPO2 = ??? '
  491. WRITE(IOIMP,*) MOTERR(1:40)
  492.  
  493. GOTO 9999
  494. ENDIF
  495. SEGSUP MLMOTS
  496. C
  497. C**** Lecture du CHPOINT ROET
  498. C
  499. ICOND = 1
  500. CALL LIROBJ('CHPOINT ',IROET,ICOND,IRETOU)
  501. CALL ACTOBJ('CHPOINT ',IROET,1)
  502. IF (IERR.NE.0) GOTO 9999
  503. C
  504. C**** Control du CHPOINT
  505. C
  506. INDIC = 1
  507. NBCOMP = 1
  508. MOT(1) = 'SCAL'
  509. CALL QUEPOI(IROET, ICEN, INDIC, NBCOMP, MOT)
  510. IF(IERR .NE. 0)THEN
  511. C
  512. C******** Message d'erreur standard
  513. C -301 0 %m1:40
  514. C
  515. MOTERR = 'CHPO3 = ??? '
  516. WRITE(IOIMP,*) MOTERR(1:40)
  517.  
  518. GOTO 9999
  519. ENDIF
  520. C
  521. C**** Lecture du CHPOINT IROY and IALP
  522. C
  523. IF (NESP .GE. 1) THEN
  524. ICOND = 1
  525. CALL LIROBJ('CHPOINT ',IROY,ICOND,IRETOU)
  526. CALL ACTOBJ('CHPOINT ',IROY,1)
  527. IF (IERR.NE.0) GOTO 9999
  528. C
  529. C
  530. C**** Control du CHPOINT
  531. C
  532. CALL QUEPO1(IROY , ICEN , MLMOT1)
  533. IF(IERR .NE. 0)THEN
  534.  
  535. C
  536. C******* Message d'erreur standard
  537. C -301 0 %m1:40
  538. C
  539. MOTERR = 'CHPO4 = ??? '
  540. WRITE(IOIMP,*) MOTERR(1:40)
  541. C
  542. C******* Message d'erreur standard
  543. C 21 2
  544. C Données incompatibles
  545. C
  546. CALL ERREUR(21)
  547. GOTO 9999
  548. ENDIF
  549. C
  550. ICOND = 1
  551. CALL LIROBJ('CHPOINT ',IALP,ICOND,IRETOU)
  552. CALL ACTOBJ('CHPOINT ',IALP,1)
  553. IF (IERR.NE.0) GOTO 9999
  554. C
  555. C
  556. C**** Control du CHPOINT
  557. C
  558. CALL QUEPO1(IALP , ICEN , MLMOT1)
  559. IF(IERR .NE. 0)THEN
  560.  
  561. C
  562. C******* Message d'erreur standard
  563. C -301 0 %m1:40
  564. C
  565. MOTERR = 'CHPO5 = ??? '
  566. WRITE(IOIMP,*) MOTERR(1:40)
  567. C
  568. C******* Message d'erreur standard
  569. C 21 2
  570. C Données incompatibles
  571. C
  572. CALL ERREUR(21)
  573. GOTO 9999
  574. ENDIF
  575. ELSE
  576. IROY = 0
  577. IALP = 0
  578. ENDIF
  579. C
  580. C
  581. C******* Option TRICHE
  582. C
  583. ICOND = 0
  584. CALL LIRCHA(NOMTRI,ICOND,IRETOU)
  585. IF(IERR .NE. 0)GOTO 9999
  586. IF(IRETOU .EQ. 0)THEN
  587. LOGTRI = .FALSE.
  588. ELSEIF(NOMTRI .EQ. 'TRICHE')THEN
  589. LOGTRI = .TRUE.
  590. ELSE
  591. LOGTRI = .FALSE.
  592. CALL ECRCHA(NOMTRI)
  593. ENDIF
  594. C
  595. C**** Calcul des sorties.
  596. C
  597. CALL PRGFM1(NESP,
  598. & ICEN,IPHI,IRO,IROVIT,IROET,IROY,IALP,
  599. & MLRMGA,MLRPGA,MLRMPI,MLRPPI,
  600. & IVIT,IPRES,IY,
  601. & LOGNEG,MESERR,
  602. & VALER)
  603. C
  604. IF(LOGNEG)THEN
  605. C
  606. C******* Pression (energie thermique) ou densité negative
  607. C
  608. C
  609. C******* Message d'erreur standard
  610. C 41 2
  611. C %m1:8 = %r1 inférieur à %r2
  612. C
  613. MESCEL = MESERR(1)
  614. MOTERR(1:8) = MESCEL(1:8)
  615. REAERR(1) = REAL(VALER(1))
  616. REAERR(2) = 0.0
  617. CALL ERREUR(41)
  618. IF(LOGTRI)THEN
  619. IERR = 0
  620. ELSE
  621. GOTO 9999
  622. ENDIF
  623. ENDIF
  624. IF (NESP .GE. 1) THEN
  625. CALL ACTOBJ('CHPOINT ',IY,1)
  626. CALL ECROBJ('CHPOINT ',IY)
  627. ENDIF
  628. C
  629. C**** Ecriture du CHPOINT contenant la pression.
  630. C
  631. CALL ACTOBJ('CHPOINT ',IPRES,1)
  632. CALL ECROBJ('CHPOINT ',IPRES)
  633. C
  634. C**** Ecriture du CHPOINT contenant la vitesse.
  635. C
  636. CALL ACTOBJ('CHPOINT ',IVIT,1)
  637. CALL ECROBJ('CHPOINT ',IVIT)
  638. C
  639. C
  640. 9999 CONTINUE
  641. END
  642.  
  643.  
  644.  
  645.  
  646.  

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