Télécharger konf21.eso

Retour à la liste

Numérotation des lignes :

  1. C KONF21 SOURCE CB215821 16/04/21 21:17:32 8920
  2. SUBROUTINE KONF21(INDMET,
  3. & IALPF,IUVF,IULF,IPF,
  4. & ITVF,ITLF,IRVF,IRLF,
  5. & ICHPSU,ICHPDI,
  6. & MELEMC,MELEMF,MELEFE,
  7. & ICHFLU,DT,
  8. & LOGNC,LOGAN,MESERR)
  9. C************************************************************************
  10. C
  11. C PROJET : CASTEM 2000
  12. C
  13. C NOM : KONFL8
  14. C
  15. C DESCRIPTION : Voir KONV2F
  16. C
  17. C 2D two fluid flow
  18. C
  19. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  20. C
  21. C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
  22. C Modified for two phase flow by
  23. C José R.
  24. C
  25. C************************************************************************
  26. C
  27. C
  28. C APPELES (Outils
  29. C CASTEM) : KRIPAD, LICHT
  30. C
  31. C APPELES (Calcul) : FLURIE, FLUXVL, FLUVLH, FLHUS1, FLHUS2
  32. C
  33. C
  34. C************************************************************************
  35. C
  36. C ENTREES
  37. C
  38. C
  39. C 1) PARAMETRES
  40. C
  41. C INDMET : 6 AUSM+
  42. C
  43. C
  44. C 2) Pointeurs des MCHAMLs
  45. C
  46. C IALPF : MCHAML sur "FACEL" containing void fraction
  47. C
  48. C IUVF : MCHAML sur "FACEL" containing vapour velocity
  49. C
  50. C IULF : MCHAML sur "FACEL" containing liquid velocity
  51. C
  52. C IPF : MCHAML sur "FACEL" containing pressure
  53. C
  54. C ITVF : MCHAML sur "FACEL" containing vapour temperature
  55. C
  56. C ITLF : MCHAML sur "FACEL" containing liquid temperature
  57. C
  58. C IRVF : MCHAML sur "FACEL" containing vapour density
  59. C
  60. C IRLF : MCHAML sur "FACEL" containing liquid density
  61. C
  62. C
  63. C 3) Pointeurs de CHPOINTs de la table DOMAINE
  64. C
  65. C ICHPSU : CHPOINT "FACE" contenant la surface des faces
  66. C
  67. C ICHPDI : CHPOINT "CENTRE" contenant le diametre minimum
  68. C de chaque element
  69. C
  70. C
  71. C 4) Pointeurs de MELEME de la table DOMAINE
  72. C
  73. C MELEMC : MELEME 'CENTRE' du SPG des CENTRES
  74. C
  75. C MELEMF : MELEME 'FACE' du SPG des FACES
  76. C
  77. C MELEFE : MELEME 'FACEL' du connectivité FACES -> ELEM
  78. C
  79. C SORTIES
  80. C
  81. C ICHFLU : pointeurs de CHPOINTs "FACE" des flux aux interfaces:
  82. C
  83. C DT : pas de temps pour le respect de la CFL-like condition
  84. C DT < DIAMMIN /2 /max(Lambda_i)
  85. C En maillage regulier cette condition garantie la
  86. C non-interaction des ondes
  87. C
  88. C
  89. C LOGNC : (LOGICAL): si .TRUE. la methode de Newton-Rapson, utilisée
  90. C dans pour la solution du probleme Riemann exacte ou dans
  91. C l'algorithm HUS, n'a pas bien marchéee; MESERR = 'Goudunov'
  92. C ou 'HUS'.
  93. C
  94. C LOGAN : (LOGICAL): si .TRUE. une anomalie à été detectée
  95. C
  96. C MESERR : pour l'ecriture des messages d'erreurs
  97. C
  98. C************************************************************************
  99. C
  100. C HISTORIQUE (Anomalies et modifications éventuelles)
  101. C
  102. C HISTORIQUE :
  103. C
  104. C************************************************************************
  105. C
  106. C
  107. C N.B.: On suppose qu'on a déjà controllé RO, P > 0
  108. C GAMMA \in (1,3)
  109. C Y \in (0,1)
  110. C Si non il faut le faire!!!
  111. C
  112. C************************************************************************
  113. C
  114. IMPLICIT INTEGER(I-N)
  115. IMPLICIT REAL*8(A-H,O-Z)
  116.  
  117. INTEGER I1
  118. & ,INDMET
  119. & ,IALPF,IUVF,IULF,IPF,ITVF,ITLF,IRVF,IRLF
  120. & ,ICHPSU,ICHPDI,MELEMC,MELEMF,MELEFE
  121. & ,IGEOMC,IGEOMF
  122. & ,ICHFLU
  123. & ,NESP, NFAC
  124. & ,NLCF, NGCEG, NGCED, NLCEG, NLCED
  125. & ,NGCF, NLCF1, SPG1, SPG2
  126. REAL*8 DT, UNSDT
  127. & , ALPG,UVNG, UVTG, ULNG, ULTG, PG, TVG, TLG
  128. & , RVG, RLG
  129. & , ALPD, UVND, UVTD, ULND, ULTD, PD, TVD, TLD
  130. & , RVD, RLD
  131. & , SURF,CNX, CNY, CTX , CTY
  132. & , CELL, DIAMG, DIAMD, DIAM
  133. & , AVSON, LAMDAV, ALSON, LAMDAL
  134. LOGICAL LOGNC, LOGAN
  135. CHARACTER*(40) MESERR
  136. CHARACTER*(8) TYPE
  137. C
  138. C**** LES INCLUDES
  139. C
  140. -INC CCOPTIO
  141. -INC SMCHAML
  142. POINTEUR MELALP.MELVAL, MELP.MELVAL,
  143. & MELTV.MELVAL, MELTL.MELVAL,
  144. & MELRV.MELVAL, MELRL.MELVAL
  145. POINTEUR MELUVN.MELVAL, MELUVT.MELVAL,
  146. & MELULN.MELVAL, MELULT.MELVAL,
  147. & MELVNX.MELVAL, MELVNY.MELVAL, MELVTX.MELVAL,
  148. & MELVTY.MELVAL, MELLNX.MELVAL, MELLNY.MELVAL,
  149. & MELLTX.MELVAL, MELLTY.MELVAL
  150. -INC SMCHPOI
  151. POINTEUR MPOVSU.MPOVAL, MPOVDI.MPOVAL
  152. & , MPOFLU.MPOVAL
  153. POINTEUR MCHAMY.MCHAML
  154. -INC SMELEME
  155. -INC SMLMOTS
  156. -INC SMLENTI
  157. C
  158. C**** Flux at each interface and each reference framework
  159. C
  160. SEGMENT IFLUX
  161. REAL*8 FLUX(8)
  162. ENDSEGMENT
  163. POINTEUR IFLU1.IFLUX
  164. C
  165. C**** Initialisation des MCHAMLs
  166. C
  167. C**** Void fraction
  168. C
  169. MCHEL1 = IALPF
  170. SEGACT MCHEL1
  171. MCHAM1 = MCHEL1.ICHAML(1)
  172. SEGACT MCHAM1
  173. MELALP = MCHAM1.IELVAL(1)
  174. SEGDES MCHEL1
  175. SEGDES MCHAM1
  176. C
  177. C**** Vapour velocity and director cosines (cosenos directores vamos??)
  178. C
  179. MCHEL1 = IUVF
  180. SEGACT MCHEL1
  181. C
  182. C**** La vitesse a comme SPG MELEFE
  183. C Le cosinus directeurs ont comme SPG MELEMF
  184. C
  185. C MCHAM1 -> Cosinus directeurs
  186. C MCHAM2 -> Vitesse
  187. C
  188. SPG1 = MCHEL1.IMACHE(1)
  189. SPG2 = MCHEL1.IMACHE(2)
  190. IF((SPG1 .EQ. MELEMF) .AND. (SPG2 .EQ. MELEFE))THEN
  191. MCHAM1 = MCHEL1.ICHAML(1)
  192. MCHAM2 = MCHEL1.ICHAML(2)
  193. ELSEIF((SPG1 .EQ. MELEFE) .AND. (SPG2 .EQ. MELEMF))THEN
  194. MCHAM1 = MCHEL1.ICHAML(2)
  195. MCHAM2 = MCHEL1.ICHAML(1)
  196. ELSE
  197. LOGAN = .TRUE.
  198. GOTO 9999
  199. ENDIF
  200. SEGACT MCHAM1
  201. MELVNX = MCHAM1.IELVAL(1)
  202. MELVNY = MCHAM1.IELVAL(2)
  203. MELVTX = MCHAM1.IELVAL(3)
  204. MELVTY = MCHAM1.IELVAL(4)
  205. SEGDES MCHAM1
  206. SEGACT MCHAM2
  207. MELUVN = MCHAM2.IELVAL(1)
  208. MELUVT = MCHAM2.IELVAL(2)
  209. SEGDES MCHAM2
  210. SEGDES MCHEL1
  211. C
  212. C**** Liquid velocity and director cosines (cosenos directores vamos??)
  213. C Claramente tengo que eliminar los cosenos directores de la
  214. C velocidad del liquido!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  215. C
  216. MCHEL1 = IULF
  217. SEGACT MCHEL1
  218. C
  219. C**** La vitesse a comme SPG MELEFE
  220. C Le cosinus directeurs ont comme SPG MELEMF
  221. C
  222. C MCHAM1 -> Cosinus directeurs
  223. C MCHAM2 -> Vitesse
  224. C
  225. SPG1 = MCHEL1.IMACHE(1)
  226. SPG2 = MCHEL1.IMACHE(2)
  227. IF((SPG1 .EQ. MELEMF) .AND. (SPG2 .EQ. MELEFE))THEN
  228. MCHAM1 = MCHEL1.ICHAML(1)
  229. MCHAM2 = MCHEL1.ICHAML(2)
  230. ELSEIF((SPG1 .EQ. MELEFE) .AND. (SPG2 .EQ. MELEMF))THEN
  231. MCHAM1 = MCHEL1.ICHAML(2)
  232. MCHAM2 = MCHEL1.ICHAML(1)
  233. ELSE
  234. LOGAN = .TRUE.
  235. GOTO 9999
  236. ENDIF
  237. SEGACT MCHAM1
  238. MELLNX = MCHAM1.IELVAL(1)
  239. MELLNY = MCHAM1.IELVAL(2)
  240. MELLTX = MCHAM1.IELVAL(3)
  241. MELLTY = MCHAM1.IELVAL(4)
  242. SEGDES MCHAM1
  243. SEGACT MCHAM2
  244. MELULN = MCHAM2.IELVAL(1)
  245. MELULT = MCHAM2.IELVAL(2)
  246. SEGDES MCHAM2
  247. SEGDES MCHEL1
  248. C
  249. C**** Pressure
  250. C
  251. MCHEL1 = IPF
  252. SEGACT MCHEL1
  253. MCHAM1 = MCHEL1.ICHAML(1)
  254. SEGACT MCHAM1
  255. MELP = MCHAM1.IELVAL(1)
  256. SEGDES MCHEL1
  257. SEGDES MCHAM1
  258. C
  259. C**** Vapour temperature
  260. C
  261. MCHEL1 = ITVF
  262. SEGACT MCHEL1
  263. MCHAM1 = MCHEL1.ICHAML(1)
  264. SEGACT MCHAM1
  265. MELTV = MCHAM1.IELVAL(1)
  266. SEGDES MCHEL1
  267. SEGDES MCHAM1
  268. C
  269. C**** Liquid temperature
  270. C
  271. MCHEL1 = ITLF
  272. SEGACT MCHEL1
  273. MCHAM1 = MCHEL1.ICHAML(1)
  274. SEGACT MCHAM1
  275. MELTL = MCHAM1.IELVAL(1)
  276. SEGDES MCHEL1
  277. SEGDES MCHAM1
  278. C
  279. C**** Vapour density
  280. C
  281. MCHEL1 = IRVF
  282. SEGACT MCHEL1
  283. MCHAM1 = MCHEL1.ICHAML(1)
  284. SEGACT MCHAM1
  285. MELRV = MCHAM1.IELVAL(1)
  286. SEGDES MCHEL1
  287. SEGDES MCHAM1
  288. C
  289. C**** Liquid density
  290. C
  291. MCHEL1 = IRLF
  292. SEGACT MCHEL1
  293. MCHAM1 = MCHEL1.ICHAML(1)
  294. SEGACT MCHAM1
  295. MELRL = MCHAM1.IELVAL(1)
  296. SEGDES MCHEL1
  297. SEGDES MCHAM1
  298. C
  299. C**** Initialisation des MELEMEs
  300. C
  301. C 'CENTRE', 'FACEL'
  302. C
  303. IPT2 = MELEFE
  304. SEGACT IPT2
  305. NFAC = IPT2.NUM(/2)
  306. C
  307. C**** KRIPAD pour la correspondance global/local de centre
  308. C
  309. CALL KRIPAD(MELEMC,MLENT1)
  310. C
  311. C**** MLENTI1 a MCORD.XCOORD(/1)/(IDIM+1) elements
  312. C
  313. C Si i est le numero global d'un noeud de ICEN,
  314. C MLENT1.LECT(i) contient sa position, i.e.
  315. C
  316. C I = numero global du noeud centre
  317. C MLENT1.LECT(i) = numero local du noeud centre
  318. C
  319. C MLENT1 déjà activé, i.e.
  320. C
  321. C SEGACT MLENT1
  322. C
  323. C
  324. C**** KRIPAD pour la correspondance global/local de 'FACE'
  325. C
  326. CALL KRIPAD(MELEMF,MLENT2)
  327. C SEGACT MELNT2
  328. C
  329. C**** Initialisation de flux
  330. C
  331. SEGINI IFLU1
  332. C
  333. C**** CHPOINTs de la table DOMAINE
  334. C
  335. CALL LICHT(ICHPSU,MPOVSU,TYPE,IGEOMF)
  336. CALL LICHT(ICHPDI,MPOVDI,TYPE,IGEOMC)
  337. C
  338. C**** LICHT active les MPOVALs en *MOD
  339. C
  340. C i.e.
  341. C
  342. C SEGACT MPOVSU*MOD
  343. C SEGACT MPOVDI*MOD
  344. C
  345. IF(IGEOMF .NE. MELEMF)THEN
  346. MESERR = 'Il ne faut pas jouer avec la console. '
  347. LOGAN = .TRUE.
  348. GOTO 9999
  349. ENDIF
  350. IF(IGEOMC .NE. MELEMC)THEN
  351. MESERR = 'Il ne faut pas jouer avec la console. '
  352. LOGAN = .TRUE.
  353. GOTO 9999
  354. ENDIF
  355. C
  356. C**** Les FLUX aux face
  357. C
  358. C La densité
  359. C
  360. CALL LICHT(ICHFLU,MPOFLU,TYPE,IGEOMF)
  361. C
  362. C SEGACT MPOFLU*MOD
  363. C
  364. IF(IGEOMC .NE. MELEMC)THEN
  365. MESERR = 'Il ne faut pas jouer avec la console. '
  366. LOGAN = .TRUE.
  367. GOTO 9999
  368. ENDIF
  369. C
  370. C**** Activation des MCHAMLs
  371. C
  372. SEGACT MELALP
  373. SEGACT MELP
  374. SEGACT MELTV
  375. SEGACT MELTL
  376. SEGACT MELRV
  377. SEGACT MELRL
  378. SEGACT MELUVN
  379. SEGACT MELUVT
  380. SEGACT MELULN
  381. SEGACT MELULT
  382. SEGACT MELVNX
  383. SEGACT MELVNY
  384. SEGACT MELVTX
  385. SEGACT MELVTY
  386. SEGACT MELLNX
  387. SEGACT MELLNY
  388. SEGACT MELLTX
  389. SEGACT MELLTY
  390. C
  391. C**** Initialisation de 1/DT
  392. C
  393. UNSDT = 0.0D0
  394. C
  395. C**** BOUCLE SUR FACEL pour le calcul du FLUX
  396. C
  397. DO NLCF = 1, NFAC
  398. C
  399. C******* NLCF = numero local du centre de facel
  400. C NGCF = numero global du centre de facel
  401. C NLCF1 = numero local du centre de face
  402. C NGCEG = numero global du centre ELT "gauche"
  403. C NLCEG = numero local du centre ELT "gauche"
  404. C NGCED = numero global du centre ELT "droite"
  405. C NLCED = numero local du centre ELT "droite"
  406. C
  407. NGCEG = IPT2.NUM(1,NLCF)
  408. NGCED = IPT2.NUM(3,NLCF)
  409. NGCF = IPT2.NUM(2,NLCF)
  410. NLCF1 = MLENT2.LECT(NGCF)
  411. NLCEG = MLENT1.LECT(NGCEG)
  412. NLCED = MLENT1.LECT(NGCED)
  413. C
  414. C******* NLCF != NLCF1 -> l'auteur (MOI) n'a rien compris.
  415. C
  416. IF(NLCF .NE. NLCF1)THEN
  417. MESERR = 'Il ne faut pas jouer avec la console. '
  418. LOGAN = .TRUE.
  419. GOTO 9999
  420. ENDIF
  421. C
  422. C******* Recuperation des Etats "gauche" et "droite"
  423. C
  424. ALPG = MELALP.VELCHE(1,NLCF)
  425. UVNG = MELUVN.VELCHE(1,NLCF)
  426. UVTG = MELUVT.VELCHE(1,NLCF)
  427. ULNG = MELULN.VELCHE(1,NLCF)
  428. ULTG = MELULT.VELCHE(1,NLCF)
  429. PG = MELP.VELCHE(1,NLCF)
  430. TVG = MELTV.VELCHE(1,NLCF)
  431. TLG = MELTL.VELCHE(1,NLCF)
  432. RVG = MELRV.VELCHE(1,NLCF)
  433. RLG = MELRL.VELCHE(1,NLCF)
  434. C
  435. ALPD = MELALP.VELCHE(3,NLCF)
  436. UVND = MELUVN.VELCHE(3,NLCF)
  437. UVTD = MELUVT.VELCHE(3,NLCF)
  438. ULND = MELULN.VELCHE(3,NLCF)
  439. ULTD = MELULT.VELCHE(3,NLCF)
  440. PD = MELP.VELCHE(3,NLCF)
  441. TVD = MELTV.VELCHE(3,NLCF)
  442. TLD = MELTL.VELCHE(3,NLCF)
  443. RVD = MELRV.VELCHE(3,NLCF)
  444. RLD = MELRL.VELCHE(3,NLCF)
  445. C
  446.  
  447. CNX = MELVNX.VELCHE(1,NLCF)
  448. CNY = MELVNY.VELCHE(1,NLCF)
  449. CTX = MELVTX.VELCHE(1,NLCF)
  450. CTY = MELVTY.VELCHE(1,NLCF)
  451. C
  452. C******* Calcul du flux aux interfaces
  453. C
  454. C AUSM+ (1) and "preconditioned" AUSM+ (2)
  455. IF((INDMET .EQ. 1) .OR. (INDMET .EQ. 2))THEN
  456. CALL fAP22F(INDMET,
  457. & ALPG, UVNG, UVTG, ULNG, ULTG, PG, TVG, TLG,
  458. & RVG, RLG,
  459. & ALPD, UVND, UVTD, ULND, ULTD, PD, TVD, TLD,
  460. & RVD, RLD,
  461. & IFLU1.FLUX)
  462. C AUSMDV (3) and "preconditioned" AUSMDV (3)
  463. ELSE IF((INDMET .EQ. 3) .OR. (INDMET .EQ. 4))THEN
  464. CALL fDV22F(INDMET,
  465. & ALPG, UVNG, UVTG, ULNG, ULTG, PG, TVG, TLG,
  466. & RVG, RLG,
  467. & ALPD, UVND, UVTD, ULND, ULTD, PD, TVD, TLD,
  468. & RVD, RLD,
  469. & IFLU1.FLUX)
  470. ELSE
  471. C
  472. C******** Message d'erreur standard
  473. C 251 2
  474. C Tentative d'utilisation d'une option non implémentée
  475. C
  476. CALL ERREUR(251)
  477. ENDIF
  478. IF(LOGAN) GOTO 9999
  479. IF(LOGNC) GOTO 9999
  480. C
  481. C******* Ecriture des flux
  482. C
  483. C FLUX(1) = alpha rv uvn
  484. C FLUX(2) = (1-alpha) rl uln
  485. C FLUX(3) = alpha rv uvx uvn + alpha p
  486. C FLUX(4) = alpha rv uvy uvn + alpha p
  487. C FLUX(5) = (1-alpha) rl ulx uln + (1-alpha) p
  488. C FLUX(6) = (1-alpha) rl uly uln + (1-alpha) p
  489. C FLUX(7) = alpha rv Hv uvn
  490. C FLUX(8) = (1-alpha) rl Hl uln
  491. C
  492. SURF = MPOVSU.VPOCHA(NLCF,1)
  493.  
  494. MPOFLU.VPOCHA(NLCF,1) = (IFLU1.FLUX(1) * SURF)
  495.  
  496. MPOFLU.VPOCHA(NLCF,2) = (IFLU1.FLUX(2) * SURF)
  497.  
  498. MPOFLU.VPOCHA(NLCF,3) = ((IFLU1.FLUX(3)*CNX+IFLU1.FLUX(5)*CTX)
  499. & * SURF)
  500.  
  501. MPOFLU.VPOCHA(NLCF,4) = ((IFLU1.FLUX(3)*CNY+IFLU1.FLUX(5)*CTY)
  502. & * SURF)
  503.  
  504. MPOFLU.VPOCHA(NLCF,5) = ((IFLU1.FLUX(4)*CNX+IFLU1.FLUX(6)*CTX)
  505. & * SURF)
  506.  
  507. MPOFLU.VPOCHA(NLCF,6) = ((IFLU1.FLUX(4)*CNY+IFLU1.FLUX(6)*CTY)
  508. & * SURF)
  509.  
  510. MPOFLU.VPOCHA(NLCF,7) = (IFLU1.FLUX(7) * SURF)
  511.  
  512. MPOFLU.VPOCHA(NLCF,8) = (IFLU1.FLUX(8) * SURF)
  513. C
  514. C******* CALCUL DU PAS DU TEMPS (CFL)
  515. C
  516. C****** A) ETAT A L'ITNERFACE
  517. C
  518. DIAMG = MPOVDI.VPOCHA(NLCEG,1)
  519. DIAMD = MPOVDI.VPOCHA(NLCED,1)
  520. DIAM = MIN(DIAMG,DIAMD)
  521. C CELL = 1.0D0/DIAM/CELLT
  522. C IF(CELL .GT. UNSDT)THEN
  523. C UNSDT = CELL
  524. C ENDIF
  525. C
  526. C****** B) ETAT GAUCHE
  527. C gammav = 1.4D0
  528. C Cpv = 1008.D0 J/(kg·K)
  529. C gammal = 2.8D0
  530. C Cpl = 4186.D0 J/(kg·K)
  531. C pil = 8.5D8 Pa
  532.  
  533. AVSON = SQRT(1008.d0*(1.4d0 - 1.d0)*TVG)
  534. LAMDAV = ABS(UVNG) + AVSON
  535. ALSON = SQRT(4186.D0*(2.8d0 - 1.d0)*TLG)
  536. LAMDAL = ABS(ULNG) + ALSON
  537. CELL = 1.D0/(((1.d0 - ALPG)/LAMDAL + ALPG/LAMDAV)* DIAM)
  538. IF(CELL .GT. UNSDT)THEN
  539. UNSDT = CELL
  540. ENDIF
  541. C
  542. C****** C) ETAT DROITE
  543. C
  544. AVSON = SQRT(1008.d0*(1.4d0 - 1.d0)*TVD)
  545. LAMDAV = ABS(UVND) + AVSON
  546. ALSON = SQRT(4186.D0*(2.8d0 - 1.d0)*TLD)
  547. LAMDAL = ABS(ULND) + ALSON
  548. CELL = 1.D0/(((1.d0 - ALPD)/LAMDAL + ALPD/LAMDAV)* DIAM)
  549. IF(CELL .GT. UNSDT)THEN
  550. UNSDT = CELL
  551. ENDIF
  552. C
  553. C
  554. C**** FIN BOUCLE SUR FACEL
  555. C
  556. ENDDO
  557. C
  558. C**** PAS DU TEMPS (CONDITION DE NON INTERACTION EN 1D)
  559. C
  560. DT = 1.D0 / UNSDT
  561. C
  562. C**** DESACTIVATION DES SEGMENTS ET
  563. C ON DETRUIT LES MCHAMLS
  564. C
  565. C
  566. C**** SEGSUP FRAMAG
  567. C SEGSUP FRAMAD
  568. C
  569. C MEME SI LOGME = .FALSE.
  570. C
  571. C
  572. SEGSUP MLENT1
  573. SEGSUP MLENT2
  574. SEGDES IPT2
  575. C
  576. SEGSUP IFLU1
  577. C
  578. SEGDES MPOVSU
  579. SEGDES MPOVDI
  580. C
  581. SEGDES MPOFLU
  582. C
  583.  
  584. SEGDES MELALP
  585. SEGDES MELP
  586. SEGDES MELTV
  587. SEGDES MELTL
  588. SEGDES MELRV
  589. SEGDES MELRL
  590. SEGDES MELUVN
  591. SEGDES MELUVT
  592. SEGDES MELULN
  593. SEGDES MELULT
  594. SEGDES MELVNX
  595. SEGDES MELVNY
  596. SEGDES MELVTX
  597. SEGDES MELVTY
  598. SEGDES MELLNX
  599. SEGDES MELLNY
  600. SEGDES MELLTX
  601. SEGDES MELLTY
  602. CC
  603. 9999 CONTINUE
  604. C
  605. RETURN
  606. END
  607. C
  608.  
  609.  
  610.  
  611.  
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  
  619.  
  620.  

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