Télécharger prigfm.eso

Retour à la liste

Numérotation des lignes :

  1. C PRIGFM SOURCE BECC 11/05/18 21:15:35 6973
  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**** Variables de COOPTIO
  94. C
  95. C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  96. C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
  97. C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
  98. C & ,IECHO, IIMPI, IOSPI
  99. C & ,IDIM
  100. C & ,MCOORD
  101. C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  102. C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  103. C & ,NORINC,NORVAL,NORIND,NORVAD
  104. C & ,NUCROU, IPSAUV
  105. C
  106. C**** Les variables
  107. C
  108. IMPLICIT INTEGER(I-N)
  109. INTEGER NESP, NESP1
  110. INTEGER ICOND, IRETOU, INDIC, NBCOMP
  111. & , ICEN, IRO, IROVIT, IROET
  112. & , IPRES, IVIT, JGN, JGM, IPHI
  113. & , IPGAS, IROY, IALP, IY
  114. REAL*8 VALER(2)
  115. CHARACTER*(40) MESERR(2),MESCEL
  116. CHARACTER*(4) MOT(1)
  117. CHARACTER*(6) NOMTRI
  118. CHARACTER*(8) MTYPR, MTYPI
  119. LOGICAL LOGNEG, LOGTRI
  120. C
  121. C**** Variables en ACCTAB
  122. C
  123. INTEGER IVALI, IRETI,IVALR, IRETR
  124. REAL*8 XVALI,XVALR
  125. LOGICAL LOGII, LOGIR
  126. CHARACTER*(8) CHARR
  127. C
  128. C**** Les Includes
  129. C
  130. -INC CCOPTIO
  131. -INC SMCHPOI
  132. -INC SMLMOTS
  133. -INC SMLREEL
  134. POINTEUR MLRMGA.MLREEL, MLRPGA.MLREEL,
  135. & MLRMPI.MLREEL, MLRPPI.MLREEL
  136. C
  137. C
  138. C**** Initialisation des parametres d'erreur
  139. C
  140. LOGNEG = .FALSE.
  141. MESCEL = ' '
  142. MESERR(1) = MESCEL
  143. MESERR(2) = MESCEL
  144. MOTERR(1:40) = MESCEL
  145. VALER(1) = 0.0D0
  146. VALER(2) = 0.0D0
  147. C
  148. C**** Initialisation des MOT(1)
  149. C
  150. MOT(1) = ' '
  151. C
  152. C**** Initialisation des variables en ACCTAB
  153. C
  154. IVALI = 0
  155. IVALR = 0
  156. XVALI = 0.0D0
  157. XVALR = 0.0D0
  158. LOGII = .FALSE.
  159. LOGIR = .FALSE.
  160. IRETI = 0
  161. IRETR = 0
  162. CHARR = ' '
  163. C
  164. C**************************************************
  165. C**** Lecture de la table des proprietes du gaz ***
  166. C**************************************************
  167. C
  168. ICOND = 1
  169. CALL QUETYP(MTYPR,ICOND,IRETOU)
  170. IF(IERR .NE. 0)GOTO 9999
  171. IF(MTYPR .NE. 'TABLE ')THEN
  172. C
  173. C******* Message d'erreur standard
  174. C 37 2
  175. C On ne trouve pas d'objet de type %m1:8
  176. C
  177. MOTERR(1:8) = 'TABLE '
  178. CALL ERREUR(37)
  179. GOTO 9999
  180. ELSE
  181. ICOND = 1
  182. CALL LIROBJ(MTYPR,IPGAS,ICOND,IRETOU)
  183. IF(IERR .NE. 0)GOTO 9999
  184. ENDIF
  185. C
  186. C**** Nom de l'espece qui n'est pas dans les equations d'Euler
  187. C
  188. MTYPI = 'MOT '
  189. MTYPR = ' '
  190. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'ESPNEULE', LOGII,IRETI,
  191. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  192. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'MOT ')) THEN
  193. C
  194. C******* Message d'erreur standard
  195. C -301 0 %m1:40
  196. C
  197. MOTERR = 'TAB1 . ESPNEULE = ??? '
  198. WRITE(IOIMP,*) MOTERR(1:40)
  199. C
  200. C********** Message d'erreur standard
  201. C 21 2
  202. C Données incompatibles
  203. C
  204. CALL ERREUR(21)
  205. GOTO 9999
  206. ENDIF
  207. C
  208. C**** Les especes qui sont dans les Equations d'Euler
  209. C
  210. MTYPR = ' '
  211. CALL ACMO(IPGAS,'ESPEULE',MTYPR,MLMOT1)
  212. IF(MTYPR .NE. 'LISTMOTS')THEN
  213. C
  214. C******* Message d'erreur standard
  215. C -301 0 %m1:40
  216. C
  217. MOTERR(1:40) = 'TAB1 . ESPEULE = ??? '
  218. WRITE(IOIMP,*) MOTERR(1:40)
  219. C
  220. C******* Message d'erreur standard
  221. C 21 2
  222. C Données incompatibles
  223. C
  224. CALL ERREUR(21)
  225. GOTO 9999
  226. ENDIF
  227. SEGACT MLMOT1
  228. NESP = MLMOT1.MOTS(/2)
  229. SEGDES MLMOT1
  230. C
  231. C JGN = 4
  232. C JGM = NESP + 1
  233. C SEGINI MLMOT2
  234. C DO I1 = 1, NESP, 1
  235. C MLMOT2.MOTS(I1) = MLMOT1.MOTS(I1)
  236. C ENDDO
  237. C MLMOT2.MOTS(NESP+1)=CHARR(1:4)
  238. C DO I1 = 1, (NESP + 1)
  239. C write(*,*) MLMOT2.MOTS(I1)
  240. C ENDDO
  241. C SEGSUP MLMOT2
  242. C stop
  243. C
  244. C
  245. C**** List de gamma, PHI < 0 (MLRMGA, MLRrel Minus GAMMA)
  246. C
  247. MTYPI = 'MOT '
  248. MTYPR = ' '
  249. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'MGAM', LOGII,IRETI,
  250. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  251. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  252. C
  253. C********** Message d'erreur standard
  254. C -301 0 %m1:40
  255. C
  256. MOTERR = ' TAB1 . MGAM = ??? '
  257. WRITE(IOIMP,*) MOTERR
  258. C
  259. C********** Message d'erreur standard
  260. C 21 2
  261. C Données incompatibles
  262. C
  263. CALL ERREUR(21)
  264. GOTO 9999
  265. ENDIF
  266. MLRMGA = IRETR
  267. SEGACT MLRMGA
  268. NESP1 = MLRMGA.PROG(/1)
  269. IF(NESP1 .NE. (NESP + 1))THEN
  270. C
  271. C********** Message d'erreur standard
  272. C -301 0 %m1:40
  273. C
  274. MOTERR = 'DIME(TAB1.MGAM) != NESP '
  275. WRITE(IOIMP,*) MOTERR(1:40)
  276. C
  277. C********** Message d'erreur standard
  278. C 21 2
  279. C Données incompatibles
  280. C
  281. CALL ERREUR(21)
  282. GOTO 9999
  283. ENDIF
  284. C
  285. C**** List de gamma, PHI > 0 (MLRPGA, MLRrel Plus GAMMA)
  286. C
  287. MTYPI = 'MOT '
  288. MTYPR = ' '
  289. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'PGAM', LOGII,IRETI,
  290. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  291. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  292. C
  293. C********** Message d'erreur standard
  294. C -301 0 %m1:40
  295. C
  296. MOTERR = ' TAB1 . PGAM = ??? '
  297. WRITE(IOIMP,*) MOTERR
  298. C
  299. C********** Message d'erreur standard
  300. C 21 2
  301. C Données incompatibles
  302. C
  303. CALL ERREUR(21)
  304. GOTO 9999
  305. ENDIF
  306. MLRPGA = IRETR
  307. SEGACT MLRPGA
  308. NESP1 = MLRPGA.PROG(/1)
  309. IF(NESP1 .NE. (NESP + 1))THEN
  310. C
  311. C********** Message d'erreur standard
  312. C -301 0 %m1:40
  313. C
  314. MOTERR = 'DIME(TAB1.PGAM) != NESP '
  315. WRITE(IOIMP,*) MOTERR(1:40)
  316. C
  317. C********** Message d'erreur standard
  318. C 21 2
  319. C Données incompatibles
  320. C
  321. CALL ERREUR(21)
  322. GOTO 9999
  323. ENDIF
  324. C
  325. C**** List de PINF, PHI < 0 (MLRMPI)
  326. C
  327. MTYPI = 'MOT '
  328. MTYPR = ' '
  329. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'MPIN', LOGII,IRETI,
  330. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  331. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  332. C
  333. C********** Message d'erreur standard
  334. C -301 0 %m1:40
  335. C
  336. MOTERR = ' TAB1 . MPIN = ??? '
  337. WRITE(IOIMP,*) MOTERR
  338. C
  339. C********** Message d'erreur standard
  340. C 21 2
  341. C Données incompatibles
  342. C
  343. CALL ERREUR(21)
  344. GOTO 9999
  345. ENDIF
  346. MLRMPI = IRETR
  347. SEGACT MLRMPI
  348. NESP1 = MLRMPI.PROG(/1)
  349. IF(NESP1 .NE. (NESP + 1))THEN
  350. C
  351. C********** Message d'erreur standard
  352. C -301 0 %m1:40
  353. C
  354. MOTERR = 'DIME(TAB1.MPIN) != NESP '
  355. WRITE(IOIMP,*) MOTERR(1:40)
  356. C
  357. C********** Message d'erreur standard
  358. C 21 2
  359. C Données incompatibles
  360. C
  361. CALL ERREUR(21)
  362. GOTO 9999
  363. ENDIF
  364. C
  365. C**** List de PINF, PHI > 0 (MLRPPI)
  366. C
  367. MTYPI = 'MOT '
  368. MTYPR = ' '
  369. CALL ACCTAB(IPGAS,MTYPI,IVALI,XVALI,'PPIN', LOGII,IRETI,
  370. & MTYPR,IVALR,XVALR ,CHARR,LOGIR,IRETR)
  371. IF((IERR .NE. 0) .OR. (MTYPR .NE. 'LISTREEL')) THEN
  372. C
  373. C********** Message d'erreur standard
  374. C -301 0 %m1:40
  375. C
  376. MOTERR = ' TAB1 . PPIN = ??? '
  377. WRITE(IOIMP,*) MOTERR
  378. C
  379. C********** Message d'erreur standard
  380. C 21 2
  381. C Données incompatibles
  382. C
  383. CALL ERREUR(21)
  384. GOTO 9999
  385. ENDIF
  386. MLRPPI = IRETR
  387. SEGACT MLRPPI
  388. NESP1 = MLRPPI.PROG(/1)
  389. IF(NESP1 .NE. (NESP + 1))THEN
  390. C
  391. C********** Message d'erreur standard
  392. C -301 0 %m1:40
  393. C
  394. MOTERR = 'DIME(TAB1.PPIN) != NESP '
  395. WRITE(IOIMP,*) MOTERR(1:40)
  396. C
  397. C********** Message d'erreur standard
  398. C 21 2
  399. C Données incompatibles
  400. C
  401. CALL ERREUR(21)
  402. GOTO 9999
  403. ENDIF
  404. C
  405. C***********************************************************************
  406. C******************** Lecture du reste *********************************
  407. C***********************************************************************
  408. C
  409. C**** Lecture du CHPOINT ROPHI
  410. C
  411. ICOND = 1
  412. CALL LIROBJ('CHPOINT',IPHI,ICOND,IRETOU)
  413. IF (IERR.NE.0) GOTO 9999
  414. C
  415. C**** On cherche le pointeur de son maillage et on l'impose sur les
  416. C autres CHPOINT
  417. C
  418. MCHPOI = IPHI
  419. SEGACT MCHPOI
  420. MSOUPO = MCHPOI.IPCHP(1)
  421. SEGACT MSOUPO
  422. ICEN = MSOUPO.IGEOC
  423. SEGDES MSOUPO
  424. SEGDES MCHPOI
  425. C
  426. C**** Control du CHPOINT: QUEPOI
  427. C
  428. C On controlle que le chpoint est non-partitione
  429. C
  430. C INDIC = 1 -> on impose le pointeur du support geometrique (ICEN)
  431. C INDIC = 0 -> on ne fait que verifier le support geometrique (ICEN)
  432. C
  433. C NBCOMP > 0 -> numero des composantes
  434. C
  435. C MOT(1) = ' ' obligatoire s'on connais pas les noms des composantes
  436. C
  437. INDIC = 1
  438. NBCOMP = 1
  439. MOT(1) = 'SCAL'
  440. CALL QUEPOI(IPHI, ICEN, INDIC, NBCOMP, MOT)
  441. IF(IERR .NE. 0)THEN
  442. C
  443. C******** Message d'erreur standard
  444. C -301 0 %m1:40
  445. C
  446. MOTERR = 'CHPO0 = ??? '
  447. WRITE(IOIMP,*) MOTERR(1:40)
  448.  
  449. GOTO 9999
  450. ENDIF
  451. C
  452. C**** Lecture du CHPOINT RO
  453. C
  454. ICOND = 1
  455. CALL LIROBJ('CHPOINT',IRO,ICOND,IRETOU)
  456. IF (IERR.NE.0) GOTO 9999
  457. C
  458. C**** Control du CHPOINT
  459. C
  460. INDIC = 1
  461. NBCOMP = 1
  462. MOT(1) = 'SCAL'
  463. CALL QUEPOI(IRO, ICEN, INDIC, NBCOMP, MOT)
  464. IF(IERR .NE. 0)THEN
  465. C
  466. C******** Message d'erreur standard
  467. C -301 0 %m1:40
  468. C
  469. MOTERR = 'CHPO1 = ??? '
  470. WRITE(IOIMP,*) MOTERR(1:40)
  471.  
  472. GOTO 9999
  473. ENDIF
  474. C
  475. C**** Lecture du CHPOINT DEBIT.
  476. C
  477. ICOND = 1
  478. CALL LIROBJ('CHPOINT',IROVIT,ICOND,IRETOU)
  479. IF (IERR.NE.0) GOTO 9999
  480. C
  481. C**** Control du CHPOINT
  482. C
  483. INDIC = 1
  484. NBCOMP = IDIM
  485. JGN = 4
  486. JGM = IDIM
  487. SEGINI MLMOTS
  488. MLMOTS.MOTS(1) = 'UX '
  489. MLMOTS.MOTS(2) = 'UY '
  490. IF(IDIM .EQ. 3) MLMOTS.MOTS(3) = 'UZ '
  491. CALL QUEPO1(IROVIT, ICEN, MLMOTS)
  492. IF(IERR .NE. 0)THEN
  493. C
  494. C******** Message d'erreur standard
  495. C -301 0 %m1:40
  496. C
  497. MOTERR = 'CHPO2 = ??? '
  498. WRITE(IOIMP,*) MOTERR(1:40)
  499.  
  500. GOTO 9999
  501. ENDIF
  502. SEGSUP MLMOTS
  503. C
  504. C**** Lecture du CHPOINT ROET
  505. C
  506. ICOND = 1
  507. CALL LIROBJ('CHPOINT',IROET,ICOND,IRETOU)
  508. IF (IERR.NE.0) GOTO 9999
  509. C
  510. C**** Control du CHPOINT
  511. C
  512. INDIC = 1
  513. NBCOMP = 1
  514. MOT(1) = 'SCAL'
  515. CALL QUEPOI(IROET, ICEN, INDIC, NBCOMP, MOT)
  516. IF(IERR .NE. 0)THEN
  517. C
  518. C******** Message d'erreur standard
  519. C -301 0 %m1:40
  520. C
  521. MOTERR = 'CHPO3 = ??? '
  522. WRITE(IOIMP,*) MOTERR(1:40)
  523.  
  524. GOTO 9999
  525. ENDIF
  526. C
  527. C**** Lecture du CHPOINT IROY and IALP
  528. C
  529. IF (NESP .GE. 1) THEN
  530. ICOND = 1
  531. CALL LIROBJ('CHPOINT',IROY,ICOND,IRETOU)
  532. IF (IERR.NE.0) GOTO 9999
  533. C
  534. C
  535. C**** Control du CHPOINT
  536. C
  537. CALL QUEPO1(IROY , ICEN , MLMOT1)
  538. IF(IERR .NE. 0)THEN
  539.  
  540. C
  541. C******* Message d'erreur standard
  542. C -301 0 %m1:40
  543. C
  544. MOTERR = 'CHPO4 = ??? '
  545. WRITE(IOIMP,*) MOTERR(1:40)
  546. C
  547. C******* Message d'erreur standard
  548. C 21 2
  549. C Données incompatibles
  550. C
  551. CALL ERREUR(21)
  552. GOTO 9999
  553. ENDIF
  554. C
  555. ICOND = 1
  556. CALL LIROBJ('CHPOINT',IALP,ICOND,IRETOU)
  557. IF (IERR.NE.0) GOTO 9999
  558. C
  559. C
  560. C**** Control du CHPOINT
  561. C
  562. CALL QUEPO1(IALP , ICEN , MLMOT1)
  563. IF(IERR .NE. 0)THEN
  564.  
  565. C
  566. C******* Message d'erreur standard
  567. C -301 0 %m1:40
  568. C
  569. MOTERR = 'CHPO5 = ??? '
  570. WRITE(IOIMP,*) MOTERR(1:40)
  571. C
  572. C******* Message d'erreur standard
  573. C 21 2
  574. C Données incompatibles
  575. C
  576. CALL ERREUR(21)
  577. GOTO 9999
  578. ENDIF
  579. ELSE
  580. IROY = 0
  581. IALP = 0
  582. ENDIF
  583. C
  584. C
  585. C******* Option TRICHE
  586. C
  587. ICOND = 0
  588. CALL LIRCHA(NOMTRI,ICOND,IRETOU)
  589. IF(IERR .NE. 0)GOTO 9999
  590. IF(IRETOU .EQ. 0)THEN
  591. LOGTRI = .FALSE.
  592. ELSEIF(NOMTRI .EQ. 'TRICHE')THEN
  593. LOGTRI = .TRUE.
  594. ELSE
  595. LOGTRI = .FALSE.
  596. CALL ECRCHA(NOMTRI)
  597. ENDIF
  598. C
  599. C**** Calcul des sorties.
  600. C
  601. CALL PRGFM1(NESP,
  602. & ICEN,IPHI,IRO,IROVIT,IROET,IROY,IALP,
  603. & MLRMGA,MLRPGA,MLRMPI,MLRPPI,
  604. & IVIT,IPRES,IY,
  605. & LOGNEG,MESERR,
  606. & VALER)
  607. C
  608. IF(LOGNEG)THEN
  609. C
  610. C******* Pression (energie thermique) ou densité negative
  611. C
  612. C
  613. C******* Message d'erreur standard
  614. C 41 2
  615. C %m1:8 = %r1 inférieur à %r2
  616. C
  617. MESCEL = MESERR(1)
  618. MOTERR(1:8) = MESCEL(1:8)
  619. REAERR(1) = REAL(VALER(1))
  620. REAERR(2) = 0.0
  621. CALL ERREUR(41)
  622. IF(LOGTRI)THEN
  623. IERR = 0
  624. ELSE
  625. GOTO 9999
  626. ENDIF
  627. ENDIF
  628. IF (NESP .GE. 1) THEN
  629. CALL ECROBJ('CHPOINT',IY)
  630. ENDIF
  631. C
  632. C**** Ecriture du CHPOINT contenant la pression.
  633. C
  634. CALL ECROBJ('CHPOINT',IPRES)
  635. C
  636. C**** Ecriture du CHPOINT contenant la vitesse.
  637. C
  638. CALL ECROBJ('CHPOINT',IVIT)
  639. C
  640. SEGDES MLRMPI
  641. SEGDES MLRMGA
  642. SEGDES MLRPPI
  643. SEGDES MLRPGA
  644. C
  645. 9999 CONTINUE
  646. C
  647. RETURN
  648. END
  649.  
  650.  

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