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.  
  141. -INC PPARAM
  142. -INC CCOPTIO
  143. -INC SMCHAML
  144. POINTEUR MELALP.MELVAL, MELP.MELVAL,
  145. & MELTV.MELVAL, MELTL.MELVAL,
  146. & MELRV.MELVAL, MELRL.MELVAL
  147. POINTEUR MELUVN.MELVAL, MELUVT.MELVAL,
  148. & MELULN.MELVAL, MELULT.MELVAL,
  149. & MELVNX.MELVAL, MELVNY.MELVAL, MELVTX.MELVAL,
  150. & MELVTY.MELVAL, MELLNX.MELVAL, MELLNY.MELVAL,
  151. & MELLTX.MELVAL, MELLTY.MELVAL
  152. -INC SMCHPOI
  153. POINTEUR MPOVSU.MPOVAL, MPOVDI.MPOVAL
  154. & , MPOFLU.MPOVAL
  155. POINTEUR MCHAMY.MCHAML
  156. -INC SMELEME
  157. -INC SMLMOTS
  158. -INC SMLENTI
  159. C
  160. C**** Flux at each interface and each reference framework
  161. C
  162. SEGMENT IFLUX
  163. REAL*8 FLUX(8)
  164. ENDSEGMENT
  165. POINTEUR IFLU1.IFLUX
  166. C
  167. C**** Initialisation des MCHAMLs
  168. C
  169. C**** Void fraction
  170. C
  171. MCHEL1 = IALPF
  172. SEGACT MCHEL1
  173. MCHAM1 = MCHEL1.ICHAML(1)
  174. SEGACT MCHAM1
  175. MELALP = MCHAM1.IELVAL(1)
  176. SEGDES MCHEL1
  177. SEGDES MCHAM1
  178. C
  179. C**** Vapour velocity and director cosines (cosenos directores vamos??)
  180. C
  181. MCHEL1 = IUVF
  182. SEGACT MCHEL1
  183. C
  184. C**** La vitesse a comme SPG MELEFE
  185. C Le cosinus directeurs ont comme SPG MELEMF
  186. C
  187. C MCHAM1 -> Cosinus directeurs
  188. C MCHAM2 -> Vitesse
  189. C
  190. SPG1 = MCHEL1.IMACHE(1)
  191. SPG2 = MCHEL1.IMACHE(2)
  192. IF((SPG1 .EQ. MELEMF) .AND. (SPG2 .EQ. MELEFE))THEN
  193. MCHAM1 = MCHEL1.ICHAML(1)
  194. MCHAM2 = MCHEL1.ICHAML(2)
  195. ELSEIF((SPG1 .EQ. MELEFE) .AND. (SPG2 .EQ. MELEMF))THEN
  196. MCHAM1 = MCHEL1.ICHAML(2)
  197. MCHAM2 = MCHEL1.ICHAML(1)
  198. ELSE
  199. LOGAN = .TRUE.
  200. GOTO 9999
  201. ENDIF
  202. SEGACT MCHAM1
  203. MELVNX = MCHAM1.IELVAL(1)
  204. MELVNY = MCHAM1.IELVAL(2)
  205. MELVTX = MCHAM1.IELVAL(3)
  206. MELVTY = MCHAM1.IELVAL(4)
  207. SEGDES MCHAM1
  208. SEGACT MCHAM2
  209. MELUVN = MCHAM2.IELVAL(1)
  210. MELUVT = MCHAM2.IELVAL(2)
  211. SEGDES MCHAM2
  212. SEGDES MCHEL1
  213. C
  214. C**** Liquid velocity and director cosines (cosenos directores vamos??)
  215. C Claramente tengo que eliminar los cosenos directores de la
  216. C velocidad del liquido!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  217. C
  218. MCHEL1 = IULF
  219. SEGACT MCHEL1
  220. C
  221. C**** La vitesse a comme SPG MELEFE
  222. C Le cosinus directeurs ont comme SPG MELEMF
  223. C
  224. C MCHAM1 -> Cosinus directeurs
  225. C MCHAM2 -> Vitesse
  226. C
  227. SPG1 = MCHEL1.IMACHE(1)
  228. SPG2 = MCHEL1.IMACHE(2)
  229. IF((SPG1 .EQ. MELEMF) .AND. (SPG2 .EQ. MELEFE))THEN
  230. MCHAM1 = MCHEL1.ICHAML(1)
  231. MCHAM2 = MCHEL1.ICHAML(2)
  232. ELSEIF((SPG1 .EQ. MELEFE) .AND. (SPG2 .EQ. MELEMF))THEN
  233. MCHAM1 = MCHEL1.ICHAML(2)
  234. MCHAM2 = MCHEL1.ICHAML(1)
  235. ELSE
  236. LOGAN = .TRUE.
  237. GOTO 9999
  238. ENDIF
  239. SEGACT MCHAM1
  240. MELLNX = MCHAM1.IELVAL(1)
  241. MELLNY = MCHAM1.IELVAL(2)
  242. MELLTX = MCHAM1.IELVAL(3)
  243. MELLTY = MCHAM1.IELVAL(4)
  244. SEGDES MCHAM1
  245. SEGACT MCHAM2
  246. MELULN = MCHAM2.IELVAL(1)
  247. MELULT = MCHAM2.IELVAL(2)
  248. SEGDES MCHAM2
  249. SEGDES MCHEL1
  250. C
  251. C**** Pressure
  252. C
  253. MCHEL1 = IPF
  254. SEGACT MCHEL1
  255. MCHAM1 = MCHEL1.ICHAML(1)
  256. SEGACT MCHAM1
  257. MELP = MCHAM1.IELVAL(1)
  258. SEGDES MCHEL1
  259. SEGDES MCHAM1
  260. C
  261. C**** Vapour temperature
  262. C
  263. MCHEL1 = ITVF
  264. SEGACT MCHEL1
  265. MCHAM1 = MCHEL1.ICHAML(1)
  266. SEGACT MCHAM1
  267. MELTV = MCHAM1.IELVAL(1)
  268. SEGDES MCHEL1
  269. SEGDES MCHAM1
  270. C
  271. C**** Liquid temperature
  272. C
  273. MCHEL1 = ITLF
  274. SEGACT MCHEL1
  275. MCHAM1 = MCHEL1.ICHAML(1)
  276. SEGACT MCHAM1
  277. MELTL = MCHAM1.IELVAL(1)
  278. SEGDES MCHEL1
  279. SEGDES MCHAM1
  280. C
  281. C**** Vapour density
  282. C
  283. MCHEL1 = IRVF
  284. SEGACT MCHEL1
  285. MCHAM1 = MCHEL1.ICHAML(1)
  286. SEGACT MCHAM1
  287. MELRV = MCHAM1.IELVAL(1)
  288. SEGDES MCHEL1
  289. SEGDES MCHAM1
  290. C
  291. C**** Liquid density
  292. C
  293. MCHEL1 = IRLF
  294. SEGACT MCHEL1
  295. MCHAM1 = MCHEL1.ICHAML(1)
  296. SEGACT MCHAM1
  297. MELRL = MCHAM1.IELVAL(1)
  298. SEGDES MCHEL1
  299. SEGDES MCHAM1
  300. C
  301. C**** Initialisation des MELEMEs
  302. C
  303. C 'CENTRE', 'FACEL'
  304. C
  305. IPT2 = MELEFE
  306. SEGACT IPT2
  307. NFAC = IPT2.NUM(/2)
  308. C
  309. C**** KRIPAD pour la correspondance global/local de centre
  310. C
  311. CALL KRIPAD(MELEMC,MLENT1)
  312. C
  313. C**** MLENTI1 a MCORD.XCOORD(/1)/(IDIM+1) elements
  314. C
  315. C Si i est le numero global d'un noeud de ICEN,
  316. C MLENT1.LECT(i) contient sa position, i.e.
  317. C
  318. C I = numero global du noeud centre
  319. C MLENT1.LECT(i) = numero local du noeud centre
  320. C
  321. C MLENT1 déjà activé, i.e.
  322. C
  323. C SEGACT MLENT1
  324. C
  325. C
  326. C**** KRIPAD pour la correspondance global/local de 'FACE'
  327. C
  328. CALL KRIPAD(MELEMF,MLENT2)
  329. C SEGACT MELNT2
  330. C
  331. C**** Initialisation de flux
  332. C
  333. SEGINI IFLU1
  334. C
  335. C**** CHPOINTs de la table DOMAINE
  336. C
  337. CALL LICHT(ICHPSU,MPOVSU,TYPE,IGEOMF)
  338. CALL LICHT(ICHPDI,MPOVDI,TYPE,IGEOMC)
  339. C
  340. C**** LICHT active les MPOVALs en *MOD
  341. C
  342. C i.e.
  343. C
  344. C SEGACT MPOVSU*MOD
  345. C SEGACT MPOVDI*MOD
  346. C
  347. IF(IGEOMF .NE. MELEMF)THEN
  348. MESERR = 'Il ne faut pas jouer avec la console. '
  349. LOGAN = .TRUE.
  350. GOTO 9999
  351. ENDIF
  352. IF(IGEOMC .NE. MELEMC)THEN
  353. MESERR = 'Il ne faut pas jouer avec la console. '
  354. LOGAN = .TRUE.
  355. GOTO 9999
  356. ENDIF
  357. C
  358. C**** Les FLUX aux face
  359. C
  360. C La densité
  361. C
  362. CALL LICHT(ICHFLU,MPOFLU,TYPE,IGEOMF)
  363. C
  364. C SEGACT MPOFLU*MOD
  365. C
  366. IF(IGEOMC .NE. MELEMC)THEN
  367. MESERR = 'Il ne faut pas jouer avec la console. '
  368. LOGAN = .TRUE.
  369. GOTO 9999
  370. ENDIF
  371. C
  372. C**** Activation des MCHAMLs
  373. C
  374. SEGACT MELALP
  375. SEGACT MELP
  376. SEGACT MELTV
  377. SEGACT MELTL
  378. SEGACT MELRV
  379. SEGACT MELRL
  380. SEGACT MELUVN
  381. SEGACT MELUVT
  382. SEGACT MELULN
  383. SEGACT MELULT
  384. SEGACT MELVNX
  385. SEGACT MELVNY
  386. SEGACT MELVTX
  387. SEGACT MELVTY
  388. SEGACT MELLNX
  389. SEGACT MELLNY
  390. SEGACT MELLTX
  391. SEGACT MELLTY
  392. C
  393. C**** Initialisation de 1/DT
  394. C
  395. UNSDT = 0.0D0
  396. C
  397. C**** BOUCLE SUR FACEL pour le calcul du FLUX
  398. C
  399. DO NLCF = 1, NFAC
  400. C
  401. C******* NLCF = numero local du centre de facel
  402. C NGCF = numero global du centre de facel
  403. C NLCF1 = numero local du centre de face
  404. C NGCEG = numero global du centre ELT "gauche"
  405. C NLCEG = numero local du centre ELT "gauche"
  406. C NGCED = numero global du centre ELT "droite"
  407. C NLCED = numero local du centre ELT "droite"
  408. C
  409. NGCEG = IPT2.NUM(1,NLCF)
  410. NGCED = IPT2.NUM(3,NLCF)
  411. NGCF = IPT2.NUM(2,NLCF)
  412. NLCF1 = MLENT2.LECT(NGCF)
  413. NLCEG = MLENT1.LECT(NGCEG)
  414. NLCED = MLENT1.LECT(NGCED)
  415. C
  416. C******* NLCF != NLCF1 -> l'auteur (MOI) n'a rien compris.
  417. C
  418. IF(NLCF .NE. NLCF1)THEN
  419. MESERR = 'Il ne faut pas jouer avec la console. '
  420. LOGAN = .TRUE.
  421. GOTO 9999
  422. ENDIF
  423. C
  424. C******* Recuperation des Etats "gauche" et "droite"
  425. C
  426. ALPG = MELALP.VELCHE(1,NLCF)
  427. UVNG = MELUVN.VELCHE(1,NLCF)
  428. UVTG = MELUVT.VELCHE(1,NLCF)
  429. ULNG = MELULN.VELCHE(1,NLCF)
  430. ULTG = MELULT.VELCHE(1,NLCF)
  431. PG = MELP.VELCHE(1,NLCF)
  432. TVG = MELTV.VELCHE(1,NLCF)
  433. TLG = MELTL.VELCHE(1,NLCF)
  434. RVG = MELRV.VELCHE(1,NLCF)
  435. RLG = MELRL.VELCHE(1,NLCF)
  436. C
  437. ALPD = MELALP.VELCHE(3,NLCF)
  438. UVND = MELUVN.VELCHE(3,NLCF)
  439. UVTD = MELUVT.VELCHE(3,NLCF)
  440. ULND = MELULN.VELCHE(3,NLCF)
  441. ULTD = MELULT.VELCHE(3,NLCF)
  442. PD = MELP.VELCHE(3,NLCF)
  443. TVD = MELTV.VELCHE(3,NLCF)
  444. TLD = MELTL.VELCHE(3,NLCF)
  445. RVD = MELRV.VELCHE(3,NLCF)
  446. RLD = MELRL.VELCHE(3,NLCF)
  447. C
  448.  
  449. CNX = MELVNX.VELCHE(1,NLCF)
  450. CNY = MELVNY.VELCHE(1,NLCF)
  451. CTX = MELVTX.VELCHE(1,NLCF)
  452. CTY = MELVTY.VELCHE(1,NLCF)
  453. C
  454. C******* Calcul du flux aux interfaces
  455. C
  456. C AUSM+ (1) and "preconditioned" AUSM+ (2)
  457. IF((INDMET .EQ. 1) .OR. (INDMET .EQ. 2))THEN
  458. CALL fAP22F(INDMET,
  459. & ALPG, UVNG, UVTG, ULNG, ULTG, PG, TVG, TLG,
  460. & RVG, RLG,
  461. & ALPD, UVND, UVTD, ULND, ULTD, PD, TVD, TLD,
  462. & RVD, RLD,
  463. & IFLU1.FLUX)
  464. C AUSMDV (3) and "preconditioned" AUSMDV (3)
  465. ELSE IF((INDMET .EQ. 3) .OR. (INDMET .EQ. 4))THEN
  466. CALL fDV22F(INDMET,
  467. & ALPG, UVNG, UVTG, ULNG, ULTG, PG, TVG, TLG,
  468. & RVG, RLG,
  469. & ALPD, UVND, UVTD, ULND, ULTD, PD, TVD, TLD,
  470. & RVD, RLD,
  471. & IFLU1.FLUX)
  472. ELSE
  473. C
  474. C******** Message d'erreur standard
  475. C 251 2
  476. C Tentative d'utilisation d'une option non implémentée
  477. C
  478. CALL ERREUR(251)
  479. ENDIF
  480. IF(LOGAN) GOTO 9999
  481. IF(LOGNC) GOTO 9999
  482. C
  483. C******* Ecriture des flux
  484. C
  485. C FLUX(1) = alpha rv uvn
  486. C FLUX(2) = (1-alpha) rl uln
  487. C FLUX(3) = alpha rv uvx uvn + alpha p
  488. C FLUX(4) = alpha rv uvy uvn + alpha p
  489. C FLUX(5) = (1-alpha) rl ulx uln + (1-alpha) p
  490. C FLUX(6) = (1-alpha) rl uly uln + (1-alpha) p
  491. C FLUX(7) = alpha rv Hv uvn
  492. C FLUX(8) = (1-alpha) rl Hl uln
  493. C
  494. SURF = MPOVSU.VPOCHA(NLCF,1)
  495.  
  496. MPOFLU.VPOCHA(NLCF,1) = (IFLU1.FLUX(1) * SURF)
  497.  
  498. MPOFLU.VPOCHA(NLCF,2) = (IFLU1.FLUX(2) * SURF)
  499.  
  500. MPOFLU.VPOCHA(NLCF,3) = ((IFLU1.FLUX(3)*CNX+IFLU1.FLUX(5)*CTX)
  501. & * SURF)
  502.  
  503. MPOFLU.VPOCHA(NLCF,4) = ((IFLU1.FLUX(3)*CNY+IFLU1.FLUX(5)*CTY)
  504. & * SURF)
  505.  
  506. MPOFLU.VPOCHA(NLCF,5) = ((IFLU1.FLUX(4)*CNX+IFLU1.FLUX(6)*CTX)
  507. & * SURF)
  508.  
  509. MPOFLU.VPOCHA(NLCF,6) = ((IFLU1.FLUX(4)*CNY+IFLU1.FLUX(6)*CTY)
  510. & * SURF)
  511.  
  512. MPOFLU.VPOCHA(NLCF,7) = (IFLU1.FLUX(7) * SURF)
  513.  
  514. MPOFLU.VPOCHA(NLCF,8) = (IFLU1.FLUX(8) * SURF)
  515. C
  516. C******* CALCUL DU PAS DU TEMPS (CFL)
  517. C
  518. C****** A) ETAT A L'ITNERFACE
  519. C
  520. DIAMG = MPOVDI.VPOCHA(NLCEG,1)
  521. DIAMD = MPOVDI.VPOCHA(NLCED,1)
  522. DIAM = MIN(DIAMG,DIAMD)
  523. C CELL = 1.0D0/DIAM/CELLT
  524. C IF(CELL .GT. UNSDT)THEN
  525. C UNSDT = CELL
  526. C ENDIF
  527. C
  528. C****** B) ETAT GAUCHE
  529. C gammav = 1.4D0
  530. C Cpv = 1008.D0 J/(kg·K)
  531. C gammal = 2.8D0
  532. C Cpl = 4186.D0 J/(kg·K)
  533. C pil = 8.5D8 Pa
  534.  
  535. AVSON = SQRT(1008.d0*(1.4d0 - 1.d0)*TVG)
  536. LAMDAV = ABS(UVNG) + AVSON
  537. ALSON = SQRT(4186.D0*(2.8d0 - 1.d0)*TLG)
  538. LAMDAL = ABS(ULNG) + ALSON
  539. CELL = 1.D0/(((1.d0 - ALPG)/LAMDAL + ALPG/LAMDAV)* DIAM)
  540. IF(CELL .GT. UNSDT)THEN
  541. UNSDT = CELL
  542. ENDIF
  543. C
  544. C****** C) ETAT DROITE
  545. C
  546. AVSON = SQRT(1008.d0*(1.4d0 - 1.d0)*TVD)
  547. LAMDAV = ABS(UVND) + AVSON
  548. ALSON = SQRT(4186.D0*(2.8d0 - 1.d0)*TLD)
  549. LAMDAL = ABS(ULND) + ALSON
  550. CELL = 1.D0/(((1.d0 - ALPD)/LAMDAL + ALPD/LAMDAV)* DIAM)
  551. IF(CELL .GT. UNSDT)THEN
  552. UNSDT = CELL
  553. ENDIF
  554. C
  555. C
  556. C**** FIN BOUCLE SUR FACEL
  557. C
  558. ENDDO
  559. C
  560. C**** PAS DU TEMPS (CONDITION DE NON INTERACTION EN 1D)
  561. C
  562. DT = 1.D0 / UNSDT
  563. C
  564. C**** DESACTIVATION DES SEGMENTS ET
  565. C ON DETRUIT LES MCHAMLS
  566. C
  567. C
  568. C**** SEGSUP FRAMAG
  569. C SEGSUP FRAMAD
  570. C
  571. C MEME SI LOGME = .FALSE.
  572. C
  573. C
  574. SEGSUP MLENT1
  575. SEGSUP MLENT2
  576. SEGDES IPT2
  577. C
  578. SEGSUP IFLU1
  579. C
  580. SEGDES MPOVSU
  581. SEGDES MPOVDI
  582. C
  583. SEGDES MPOFLU
  584. C
  585.  
  586. SEGDES MELALP
  587. SEGDES MELP
  588. SEGDES MELTV
  589. SEGDES MELTL
  590. SEGDES MELRV
  591. SEGDES MELRL
  592. SEGDES MELUVN
  593. SEGDES MELUVT
  594. SEGDES MELULN
  595. SEGDES MELULT
  596. SEGDES MELVNX
  597. SEGDES MELVNY
  598. SEGDES MELVTX
  599. SEGDES MELVTY
  600. SEGDES MELLNX
  601. SEGDES MELLNY
  602. SEGDES MELLTX
  603. SEGDES MELLTY
  604. CC
  605. 9999 CONTINUE
  606. C
  607. RETURN
  608. END
  609. C
  610.  
  611.  
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  
  619.  
  620.  
  621.  
  622.  

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